Sk:Generovanie dlaždíc Ortofotomozaika SR

From OpenStreetMap Wiki
Jump to navigation Jump to search

Postup:

  1. Stiahnúť časti Slovenska z https://www.geoportal.sk/sk/udaje/ortofotomozaika/ a rozbaliť
  2. 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
  3. Pripraviť si systém na presnú reprojekciu z S-JTSK: https://www.geoportal.sk/files/gz/s-jtsk_jtsk03_v_qgis.pdf
  4. 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):
    1. 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
    2. Spustite skript computeTeTs.js [zoom] [minX] [minY] [maxX] [maxY]. Príklad: node computeTeTs.js 19 1873395 6061762 2512485 6380199
  5. Voliteľne orezať: https://groups.google.com/g/osm_sk/c/po6X5eMkxUI (TODO pridať do repozitára vygenerovaný shapefile na orezanie a príkazy)
  6. 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
  7. Skontrolovať či nie sú niektoré súbory vadné: find -name '*.jpg' -print0 | parallel -0 -j24 -m jpeginfo --check 2>&1 | grep -v '\[OK\]$'
  8. Odstrániť okrajové prázdne taily a taily s artefaktami
    1. Odstrániť prázdne taily (upraviť veľkosť): find . -type f -size -800c -delete
    2. 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
  9. 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.