GoldenLine
Zaloguj się

konto usunięte

Temat: Kalibracja dawnych planów

Piotr Tracz:
Ta metoda nazywa się TPS (dostępna zarówno w ArcGIS, darmowym QGIS jak i w innych programach), wg mnie nie ma sensu kombinować jeśli są już dostępne odpowiednie metody i narzędzia.
Oczywiście, że nie ma, tylko ja jej nie stosowałem, więc nie wiedziałem :).
Czy transformacja następuje tylko w zadanym obszarze rastra, a jego pozostała część jest nieruszana? Z definicji metody tak to wygląda (pobieżnie poczytałem), ale...
Czy to nie jest metoda "globalna", czyli obliczająca parametry transformacji dla całego zestawu punktów? A przecież chodzi o to, aby skalibrować każdą strefę oddzielnie w oparciu o własne, zindywidualizowane parametry.Kamil Nieścioruk edytował(a) ten post dnia 24.02.10 o godzinie 11:51
Link do wypowiedziLink
Piotr T.

Piotr T. programista GIS,
analityk GIS,
(www.cgis.pl)

Temat: Kalibracja dawnych planów

Kamil Nieścioruk:
Piotr Tracz:
Ta metoda nazywa się TPS (dostępna zarówno w ArcGIS, darmowym QGIS jak i w innych programach), wg mnie nie ma sensu kombinować jeśli są już dostępne odpowiednie metody i narzędzia.
Oczywiście, że nie ma, tylko ja jej nie stosowałem, więc nie wiedziałem :).
Czy transformacja następuje tylko w zadanym obszarze rastra, a jego pozostała część jest nieruszana? Z definicji metody tak to wygląda (pobieżnie poczytałem), ale...
Czy to nie jest metoda "globalna", czyli obliczająca parametry transformacji dla całego zestawu punktów? A przecież chodzi o to, aby skalibrować każdą strefę oddzielnie w oparciu o własne, zindywidualizowane parametry.Kamil Nieścioruk edytował(a) ten post dnia 24.02.10 o godzinie 11:51

Metoda TPS gwarantuje, że współrzędna (GCP), którą przyporządkowałeś do piksela rastra przed transformacją znajdzie się również w tym pixelu i po transformacji - nie ma wyliczanej wiec żadnej globalnej transformacji (tak jak w wielomianie gdzie wiązanie współrzędna (GCP)->pixel przed i po transformacji może być różna). Można wiec uogólnić, że TPS dzieli raster na takie jakby sekcje - jeśli wiesz jak działa tzw "Rubber Sheeting" dla rastrów to jest to identyczne. Oczywiście, że TPS przetwarza również obszary rastra znajdującego się poza zasięgiem punktów kontrolnych (GCP) ale ponieważ tam dane są ekstrapolowane to nie ma gwarancji co do jakości wyników (zresztą nie mając danych to i człowiek ręcznie nie da rady no chyba, że na oko i przypuszczenie).

Jakby co to ESRI jak zwykle wszystko ładnie tłumaczy:
http://webhelp.esri.com/arcgisdesktop/9.3/body.cfm?toc...
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

a co wy na to:
znajac polozenie punktow stabilnych XiY i dodajac do nich kolumne z wartoscia bledu RMS (odczytanego z transformacji np. 2 wielomianu) otrzymamy .shp z trzema kolumnami X,Y,RMS_Error. Nastepnie mozemy przeprowadzic interpolacje np. IDW i otrzymamy w efekcie rastra z rozmieszczeniem bledu RMS i podzialem na klasy bledu RMS (np. klasyfikacja manualna).
Link do wypowiedziLink

konto usunięte

Temat: Kalibracja dawnych planów

Też ciekawe. Takie mapy wykorzystuje się w analizach - wektory przesunięć lub właśnie izolinie wektorów przesunięć (zobacz np. http://jenny.cartography.ch/pdf/2007_Jenny_etal_MapAna... lub http://www.unibw.de/ipk/karto-en/publications/pubbeine... - to drugie to ciekawy doktorat, zwłaszcza gdy zna się niemiecki :), ale nawet bez tego sobie jakoś skorzystałem...). Przeważnie wartością są jednostki metryczne, ale co to szkodzi spróbować dla RMS? Jeśli zrobiłeś takie coś, to możesz się pochwalić, chętnie zobaczę :).Kamil Nieścioruk edytował(a) ten post dnia 01.03.10 o godzinie 09:47
Link do wypowiedziLink
Piotr T.

Piotr T. programista GIS,
analityk GIS,
(www.cgis.pl)

Temat: Kalibracja dawnych planów

Piotr Noszczyński:
a co wy na to:
znajac polozenie punktow stabilnych XiY i dodajac do nich kolumne z wartoscia bledu RMS (odczytanego z transformacji np. 2 wielomianu) otrzymamy .shp z trzema kolumnami X,Y,RMS_Error. Nastepnie mozemy przeprowadzic interpolacje np. IDW i otrzymamy w efekcie rastra z rozmieszczeniem bledu RMS i podzialem na klasy bledu RMS (np. klasyfikacja manualna).

Wg mnie lepsza będzie proponowana wcześniej przez Kamila siatka zniekształceń. Zrobienie interpolacji z punktów RMS jest banalnie proste i pewnie będzie wyglądało ok (szczególnie jak się ładnie dobierze przedziały klasowe) ale obawiam się, ze merytorycznie łatwo się będzie do tego przyczepić. Jeśli RMS'y będą duże to kalibracja będzie delikatnie mówiąc przybliżona, a jeśli małe to odchyłki pokazane na takiej wyinterpolowanej mapie nie będą miały znaczenia dużego (nie wiem więc co tak naprawdę będzie przedstawiane na takiej mapie). Kolejny problem wg mnie to różnice w działaniu interpolacji metodą IDW oraz wielomianu - każdy działa inaczej (ma inne podstawy matematyczne) i wyniki raczej nie będą spójne...
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

IDW jako raster z X,Y,RMS_Error dla mapy z 1700 roku.


Obrazek


Obrazek
Link do wypowiedziLink

konto usunięte

Temat: Kalibracja dawnych planów

Ciekawe, ale...
:)
Sam wskaźnik - jak mówił Piotr - jest nieco podejrzany. Co przedstawia ta mapa? Wartość RMS jest trudna w odbiorze, zwłaszcza w formie izolinii i niewiele mówi odbiorcy o charakterystyce geometrycznej mapy. Nie da się z niej wyczytać np. kierunku zniekształcenia, a jedynie wartość w punkcie.
Dziwnie też wyglądają wartości przedziałów dobrane z automatu, które dodatkowo się zazębiają (brak informacji o charakterze przedziału - otwarty/zamknięty - w związku z czym jedna wartość jest w dwu klasach).
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

OK, ale warto było spróbować :)
Teraz spróbuje się zabrać za siatkę zniekształceń, tylko jeszcze nie wiem jak. Najpierw poczytam artykuły, które zamieścił Kamil w formie linków.
Link do wypowiedziLink

konto usunięte

Temat: Kalibracja dawnych planów

Zerknij na mój poprzedni post w temacie. Najprościej nanieść siatkę na mapę współczesną i przekalibrować (naciągnąć) do mapy dawnej. Wtedy siatka na współczesnej ("kwadraty") sama się zniekształci :) i możesz ją przerysować na mapę dawną. Sposób dość banalny, ale wystarcza, bo siatka zniekształceń jest pewnym "obrazem koncepcyjnym", a nie dokładnym obrazem matematycznym (teoretycznie).
Link do wypowiedziLink
Piotr T.

Piotr T. programista GIS,
analityk GIS,
(www.cgis.pl)

Temat: Kalibracja dawnych planów

1. Tworzysz pusty TIFF o rozmiarach (w pixelach) jak skan mapy
2. Rysujesz na nim oczka siatki
3. Wczytujesz swój pusty TIFF
4. Wczytujesz punkty kalibracji
5. Uruchamiasz rektyfikacje
6. Bawisz się usuwaniem koloru tła i w ten sposób możesz nakładać swoją siatkę na zrektyfikowany raster lub zamieniasz siatkę raster na siatkę wektor
Link do wypowiedziLink

konto usunięte

Temat: Kalibracja dawnych planów

Nawet jeszcze prościej!
Chyba, że masz już gotową mapę współczesną z siatką (o czym ja pisałem) - a często właśnie na takową się kalibruje.
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

rozumiem ze punkty 1 i 2 tworze przy pomocy programu graficznego (np.Photoshop) a cala resztę (pkt.3-6) np. w QGIS lub ArcGIS (?).
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

spróbowałem zrobić siatkę zniekształceń w darmowym MapAnalyst, i taki jest tego efekt:

Obrazek
Link do wypowiedziLink

konto usunięte

Temat: Kalibracja dawnych planów

Podejrzanie regularna :). Czy te "wąsy" przy punktach to wektory przesunięć?
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

tak, to są wektory przesunięć wraz z okręgami.
Link do wypowiedziLink

konto usunięte

Temat: Kalibracja dawnych planów

No właśnie.
W siatce zniekształceń, którą mam na myśli nie może być takiej sytuacji, że między punktem, a jego zniekształconym obrazem są dodatkowe linie siatki. Co to znaczy dodatkowe? To znaczy, że zniekształcona jest przestrzeń oczka siatki, w którym znajduje się badany punkt, a zatem obszar między punktem, a jego zniekształconym obrazem musi nadal pozostać jednym oczkiem siatki. Zobacz na ilustracje - górna to sytuacja wyjściowa (A', B', C' - punkty na planie dawnym), dolna - przekształcona. Zakrzywieniu ulegają linie siatki, ale zachowana zostaje topologia pomiędzy punktami. Tymczasem na tym, co załączyłeś, widać że przesunięty punkt przekracza swoje oczko siatki (a właściwie odwrotnie - oczko jest za małe i za słabo zniekształcone).

Obrazek
Kamil Nieścioruk edytował(a) ten post dnia 17.03.10 o godzinie 21:46
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

witam,
tez bylem zdziwiony tym ze siatka jest stosunkowo regularna.
Korzystałem ze darmowego MapAnalyst (http://mapanalyst.cartography.ch/) i raczej caly proces zostal przeze mnie przeprowadzony poprawnie - wzorowalem sie na dosc dobrym tutorialu, ktory jest tam umieszczony.
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

witam,
a tak wygląda siatka zniekształceń oparta o metodę opisana przez Piotrka:

Obrazek
Piotr Noszczyński edytował(a) ten post dnia 19.03.10 o godzinie 19:52
Link do wypowiedziLink

konto usunięte

Temat: Kalibracja dawnych planów

Tendencja taka sama, ale druga wersja jakoś lepiej wygląda. Oczka siatki powyginane sensownie... Widać też z tej siatki, że punktów stabilnych nie ma za dużo (ale cóż - takie czasem życie przy dawnych planach).
Kiedyś czytałem o MapAnalyst, ale na sobie nie sprawdzałem. Może warto by było. Strona programu niestety podaje mało informacji o metodologii.
Link do wypowiedziLink
Piotr Noszczyński

Piotr Noszczyński Archeolog,
Forschungsstelle
Kaiserpfalz
Ingelheim

Temat: Kalibracja dawnych planów

a tak okiem specjalisty patrząc, to co właściwie można odczytać (wnioski) z powyższej siatki zniekształceń?
Link do wypowiedziLink
« Wróć do tematów
Odpowiedz

Podobne tematy


Następna dyskusja:

kalibracja WIG 1936




Wyślij zaproszenie do
Anuluj