Sk:Generovanie dlaždíc Ortofotomozaika SR
Jump to navigation
Jump to search
Postup:
- Stiahnúť časti Slovenska z https://www.geoportal.sk/sk/udaje/ortofotomozaika/ a rozbaliť
- Vygenerovať vrt - poradie od najstaršieho po najnovšie:
gdalbuildvrt -r cubic -b 1 -b 2 -b 3 -resolution highest sk.vrt Ortofoto_2020_zapad_jtsk_rgb/*.tif Ortofoto_2021_stred_jtsk_rgb/*.tif Ortofoto_2022_vychod_jtsk_rgb/*.tif
- Pripraviť si systém na presnú reprojekciu z S-JTSK: https://www.geoportal.sk/files/gz/s-jtsk_jtsk03_v_qgis.pdf
- Vygenerovať reprojektnutý GeoTIFF optimalizovaný pre Z19 (pixel-aligned):
gdalwarp -s_srs 'EPSG:5514' -t_srs 'EPSG:3857' -ct "+proj=pipeline +step +inv +proj=krovak +lat_0=49.5 +lon_0=24.8333333333333 +alpha=30.2881397527778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +step +inv +proj=hgridshift +grids=Slovakia_JTSK03_to_JTSK.gsb +step +proj=krovak +lat_0=49.5 +lon_0=24.8333333333333 +alpha=30.2881397527778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +step +inv +proj=krovak +lat_0=49.5 +lon_0=24.8333333333333 +alpha=30.2881397527778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +step +proj=push +v_3 +step +proj=cart +ellps=bessel +step +proj=helmert +x=485.021 +y=169.465 +z=483.839 +rx=-7.786342 +ry=-4.397554 +rz=-4.102655 +s=0 +convention=coordinate_frame +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=webmerc +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84" -co JPEG_QUALITY=95 -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co NUM_THREADS=24 -co TILED=YES -co BIGTIFF=YES -wo NUM_THREADS=ALL_CPUS -multi -r cubic -te 1873395.1262413848 6061762.091127617 2512485.1197337434 6380199.049545804 -ts 2140416 1066497 sk.vrt ../sk-warped.tif
Výpočet parametrov-te
a-ts
pre zvolenú oblasť (ak chceme oblasť inú ako celé Slovensko):- Na https://tools.geofabrik.de/calc/#type=geofabrik_standard&bbox=16.828995,47.72824,22.570038,49.616939&tab=1&proj=EPSG:3857 si zvoľte želanú oblasť a skopírujte údaje z poľa Simple copy. Príklad:
1873395 6061762 2512486 6380199
- Spustite skript
computeTeTs.js [zoom] [minX] [minY] [maxX] [maxY]
. Príklad:node computeTeTs.js 19 1873395 6061762 2512485 6380199
- Na https://tools.geofabrik.de/calc/#type=geofabrik_standard&bbox=16.828995,47.72824,22.570038,49.616939&tab=1&proj=EPSG:3857 si zvoľte želanú oblasť a skopírujte údaje z poľa Simple copy. Príklad:
- Voliteľne orezať: https://groups.google.com/g/osm_sk/c/po6X5eMkxUI (TODO pridať do repozitára vygenerovaný shapefile na orezanie a príkazy)
- Vygenerovať tiles (pre JPEG viď. https://github.com/OSGeo/gdal/pull/7655):
gdal2tiles.py --tiledriver=JPEG -r near -z 19 --xyz -x --processes=$(getconf _NPROCESSORS_ONLN) -n -w none sk-warped.tif ofmozaika
- Skontrolovať či nie sú niektoré súbory vadné:
find -name '*.jpg' -print0 | parallel -0 -j24 -m jpeginfo --check 2>&1 | grep -v '\[OK\]$'
- Odstrániť okrajové prázdne taily a taily s artefaktami
- Odstrániť prázdne taily (upraviť veľkosť):
find . -type f -size -800c -delete
- Odstrániť malé čiernobiele taily:
find . -type f -size -10000c | parallel 'convert {} -colorspace HSB -channel G -separate +channel -format "%[filename] avg=%[mean] peak=%[max]\n" info:' | grep 'avg=0 peak=0' | awk '{print $1}' | xargs rm
- Odstrániť prázdne taily (upraviť veľkosť):
- Zoomy 18 a nižšie vygenerovať pomocou Pyramid Builder.
Nejaké poznámky sú spomenuté aj na https://groups.google.com/g/osm_sk/c/W3Ayv1OflF4/m/zOJk9kWxAAAJ.