geopaparazzi

geopaparazziを土地被覆地上検証に使うには。

geopaparazzi 6.1.0 https://www.geopaparazzi.org/geopaparazzi/#_handling_datasets_from_and_for_geopaparazzi

まずQGISにHRLULCを読み込み, geopaparazziがインポートできるmbtilesという形式に変換する!

プロセシング → ツールボックス → ラスターツール − XYZタイルの生成（MBTiles形式）

次にこれをスマホにコピーする（私はUbuntu Linux 20.04のパソコンにUSBケーブルでAndroidスマホをつなぎ, 勝手にマウントしてくれたスマホのファイルシステムにコピーしました）:

$ cp ~/LC_N36E140.mbtiles /run/user/1000/gvfs/(...)/geopaparazzi/projects/
$ cp ~/LC_N36E140.mbtiles /run/user/1000/gvfs/mtp:host=OPPO_OPPO_Reno5_A_93c88c27/内部共有ストレージ/geopaparazzi/projects/

MBTiles は SQLite を利用したラスタ形式のひとつ

$ gdal_translate -co "RESAMPLING=NEAREST" -co "ZOOM_LEVEL_STRATEGY=UPPER" -of MBTiles -ot Byte -strict -scale 0 255 bogota.tif bogota.mbtiles
$ gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -r bilinear JAXA_HRLULC_Japan_v21.11/LC_N36E140.tif dummy.tif

RGBA画像にする (結果はヴァーチャルラスターファイル):

$ gdal_translate -of vrt -expand rgba dummy.tif temp.vrt

タイルにする:

gdal2tiles.py temp.vrt dummy2
vi dummy2/tilemapresource.xml
gdal2tiles.py dummy.tif
gdal_translate -co "RESAMPLING=NEAREST" -co "ZOOM_LEVEL_STRATEGY=UPPER" -of MBTiles dummy.tif LC_N36E140.mbtiles 
gdal_merge.py L*.tif -o test.vrt -ul_lr 122 46 154 20 -n 0 -a_nodata 0