Dzień dobry wszystkim
Nie radzę sobie z transformacją Helmerta w QGIS - znając parametry transformacji nie potrafię napisać poprawnej definicji układu odniesienia w PROJ (lub WKT).
Jest to niejako kontynuacja mojego poprzedniego
tematu dotyczącego georeferencji map w układzie lokalnym ale kilka lat sprawiło, że mam już więcej wiedzy i tamte sprawy mam ogarnięte, a ten wątek dotyczy już samej transformacji i jej użycia w praktyce.
Ale do brzegu:
Posiadam mapy, ze współrzędnymi w metrach. Sekcje map są w lokalnym układzie odniesienia, obrócone o 77,87st.
Potrafię wyliczyć parametry transformacji Helmerta [LOK -> ukł. 2000] na podstawie punktów ze znanymi współrzędnymi w obu układach.
Potrafię przeliczać współrzędne układu lokalnego na układ 2000 (i odwrotnie)
ale robię to w excelu...
Potrafię je poprawnie wyświetlić na warstwie po imporcie CSV do QGIS
ale:
- chciałbym móc pominąć excela, czyli importować pliki CSV ze współrzędnymi w układzie lokalnym bezpośrednio do QGIS
Posiadam współrzędne środków ciężkości dla obydwu układów i jeszcze dwa parametry transformacji u=0,9776 i v=0,21
Stworzyłem w QGIS nowy układ odniesienia użytkownika (
na podstawie dokumentacji PROJ), ale nie jestem pewien które wartości dokładnie powinienem wpisać.
Kod:
+proj=helmert +convention=coordinate_frame +x=5586823.219 +y=6577924.687 +s=0.999958259058672 +theta=280338.3
gdzie x i y to współrzędne (w układzie 2000) początku układu lokalnego (0, 0), s to skala a theta to obrót (77,87st * 3600)
ale...
to nie działa...
Próbowałem też po omacku:
Kod:
+proj=pipeline
+step +proj=tmerc +lat_0=0 +lon_0=18 +k=0.999923 +x_0=6500000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
+step +proj=helmert +convention=coordinate_frame +x=5586823.219 +y=6577924.687 +s=0.999958259058672 +theta=-280338.3 +units=m
Jak ustawię jeden lub drugi układ współrzędnych dla warstwy, to QGIS wyrzuca błąd że nie posiada parametrów transformacji. Jak wejdę w ustawienia QGISa i w Transformacje, to mogę tam tylko dodać nową transformację MOJ_LOK -> 2000 ale nie mogę w żaden sposób 'edytować' jej parametrów... I ostatecznie nic to nie zmienia...
Szukam tego w internetach od dłuższego czasu, przejrzałem wiele zagranicznych wątków i wiele prób za mną, ale... może ja to w ogóle źle robię?
Cytując "georobot" z 2018 roku:
Cytat:A potem to już "bułka z masłem"
Niestety nie dla mnie... Proszę o pomoc.
Jeśli to co napisałem powyżej jest nie do końca zrozumiałe, to moje pytanie brzmi:
Jak w praktyce wykorzystać posiadane parametry transformacji Helmerta do poprawnego wyświetlania na mapie współrzędnych w układzie lokalnym bez przeliczania ich poza QGISem (w excelu).
Z góry dziękuję za zainteresowanie i pomoc!
RE: [QGIS][PROJ] Transformacja Helmerta - lokalny układ odniesienia