Přesná PostGIS transformace S-JTSK (5514) <-> WGS-84 (4326)

[motyc]

Cílem je najít řešení nezávislé na službě ČÚZK a její dostupnosti.

Např. zde lze najít transformace s centimetrovou přesností: http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid (podle autorů to lze implementovat v PostGIS, je tam i návod).

Nejdříve je třeba provést analýzu, zda je tato nebo jiná obdobná implementace realistická a na základě toho navrhnout update stávajícího řešení.

Problém mimo jiné je, že nyní se nám souřadnice transformují s odchylkou. Testoval jsem to na http://192.168.254.30:8080/pas/zapsat a srovnával s https://geoportal.cuzk.cz/(S(go2kkco2unz1m1ccb0vssijc))/Default.aspx?mode=TextMeta&text=wcts&menu=191. Přitom údajně tuto službu využíváme (což tím pádem asi není pravda).

13.8098558 -> -785478.58 (CUZK) 13,8098558 -> -785477.88 (AMCR)

50.1278954 -> -1032422.39 (CUZK) 50,1278954 -> -1032422.01 (AMCR)

-785478.58 -> 13.8098558 (CUZK) -785478.58 -> 13,8098553 (AMCR)

-1032422.38 -> 50.1278955 (CUZK) -1032422.38 -> 50,1278952 (AMCR)

Pokud to nejde na naší straně v PostGIS, alespoň správně napojit službu ČÚZK.

SELECT ident_cely, geom_system, ST_asText(geom) as wgs, ST_asText(geom_sjtsk) as sjtsk FROM public.samostatny_nalez
ORDER BY id DESC LIMIT 30;
SELECT ident_cely, geom_system, ST_asText(geom) as wgs, ST_asText(geom_sjtsk) as sjtsk FROM pian ORDER BY id DESC LIMIT 30;

[motyc]

Test v0.4.36

Transformace probíhá perfektně, je třeba pouze dořešit související úpravy:

• [ ] Na některých místech v kódu stále pracujeme s 5514* (označovalo situaci, kdy transformace proběhla nepřesným způsobem kvůli nedostupnosti služby - to už nebudeme potřebovat). Bylo vhodné kód pročistit a související řešení odstranit. Jde minimálně o: https://github.com/ARUP-CAS/aiscr-webamcr/blob/92607ec1379163c15443b1c1137a972ed9e39709/webclient/pas/models.py#L271 https://github.com/ARUP-CAS/aiscr-webamcr/blob/92607ec1379163c15443b1c1137a972ed9e39709/webclient/pian/models.py#L122 https://github.com/ARUP-CAS/aiscr-webamcr/blob/92607ec1379163c15443b1c1137a972ed9e39709/webclient/cron/tasks.py#L75 https://github.com/ARUP-CAS/aiscr-webamcr/blob/92607ec1379163c15443b1c1137a972ed9e39709/webclient/cron/tasks.py#L238
• [ ] Je třeba migraci (https://github.com/ARUP-CAS/aiscr-webamcr/blob/dev/migrace/copy_data.py) doplnit o krok, kdy se vygenerují data SJTSK pro staré samostatné nálezy a PIAN.
• [ ] pian.geom_sjtsk se systematicky nevyplňují data při vytvoření nového PIAN katastru (viz např. P-2334-900077). Je tak třeba upravit: https://github.com/ARUP-CAS/aiscr-webamcr/blob/92607ec1379163c15443b1c1137a972ed9e39709/webclient/pian/models.py#L272

[jhavrlant]

@motyc Pro migraci jsem upravil funkce pian_to_sjtsk() a nalez_to_sjtsk() v task.py. Po migraci by se spustily v Djangu shell. Stačí to tak? Nechce se mi to integrovat do copy_data.py, neznám jeho parametry spuštění a taky bych to těžko testoval.

[motyc]

@jhavrlant Zítra bych se s Tebou spojil a nějak se domluvíme.

[jhavrlant]

@motyc Migrace transformace souřadnic s spustí: sudo docker exec -it $(sudo docker ps -q -f name=swarm_webamcr_web) python3 manage.py shell --settings=webclient.settings.production

import cron.tasks; cron.tasks.nalez_to_sjtsk() cron.tasks.pian_to_sjtsk() Pokud jsi jako administrátor můžeš smazat sudo

[motyc]

@jhavrlant Tak na pokus na .30 skončil chybou:

EDIT: Teď koukám, že tam je ještě jeden související pull, tak to asi bude oprava a testoval jsem příliš brzy...