georeferencja z plikiem points

[ksnich]

Witam, z quantum GISem mam już jakiś czas do czynienia, ale przed takim problemem stanąłem po raz pierwszy.

Mam fragmenty map w png, bmp oraz niektóre w zwykłych tiff. Do wszystkich mam dołączony plik points. Z tego co wiem, mapy jakimi dysponuję są bezpośrednimi obrazami z WGS 84 (epsg 4326).
W QGis wykorzystuję georeferencję, wczytuję mapę, punkty współrzędnych dodają się przeważnie same (gdy mają te same nazwy), wybieram transformację, układ współrzędnych (oczywiście WGS 84) i ok.
Powstała mapa (geotiff) jest zawsze zniekształcona - węższa i dłuższa (np. z pliku 5000px x 3750px robi 5590px x 2795px). Wynik jest taki sam przy wszystkich rodzajach transformacji, chociaż przy Helmercie zniekształcenie jest odrobinę inne.

Czy to co robi QGis jest prawidłowe skoro mapa była w WGS 84, a ja układ pozostawiam bez zmian?

Jeżeli jakimś cudem okazałoby się to działaniem prawidłowym to jak ręcznie zamienić tiff'a na geotiff'a. QGis zamieniając plik 5000x3750px na 5590x2795px bardzo robi to niedokładnie, jakość takiej zamiany jest mierna. Więc wolałbym zrobić ją sobie ręcznie w jakimkolwiek programie graficznym, tylko jak później stworzyć z tego geotiff'a?

Pozdrawiam.

[Kristoforos]

Może spróbować sobie to zarejestrować w programie Map Window.

Generalnie w WGS 84 mapa powinna wyglądać nieco spłaszczona w stosunku do naszych, lokalnych współrzędnych. A może coś z tymi punktami orientacyjnymi nie gra? Może pomoże ręczne wprowadzenie tych punktów na nowo?

[ksnich]

Racja, WGS 84 jest spłaszczony, o czym zapomniałem.
Poświęciłem dzisiaj kolejny cały dzień na mój problem, no i prześledziłem jeszcze raz po kolei co się dzieje z tymi mapami, a więc:

Mapy mają dołączony plik "points", niektóre plik "map". Z tego drugiego idzie odczytać odwzorowanie i ewidentnie jest WGS 84
1 ,Map Code,
WGS 84,WGS 84, 0.0000, 0.0000,WGS 84
Reserved 1
Reserved 2
Magnetic Variation,,,E

Tworzę nowy projekt, definiuje odwzorowanie, oczywiście Poland CS92, żeby nie było wspomnianego zniekształcenia. Wczytuje sobie kontur kraju (wektor w CS92) i jest wszystko ok.
Następnie zabieram się za którykolwiek kawałek mapy (te z plikiem points) i dokonuję ich georeferencji.
Przy transformacji liniowej, wielomianu 1, sklejania i projective (niezależnie czy wybiorę WGS 84, czy CS92) mapa zawsze jest spłaszczona i w innym miejscu przestrzeni niż obrys Polski.
Przy transformacji Helmerta mapa wczytana jest bez żadnych zniekształceń (wizualnie 100% ok), ale umieszczona jest w dokładnie tym samym miejscu przestrzeni, co przy pozostałych transformacjach. Miejsce jest to poza obrysem kraju, a podziałka wskazuje 0,1m, a powinna pomiędzy 20-40km...

Jakieś pomysły?

Poniżej zawartość jednego z plików points (mapa ma 5000x3750px)
mapX,mapY,pixelX,pixelY,enable
18.527999999999999,49.903999999999996,0.000000000000000,0.000000000000000,1
19.007999999999999,49.903999999999996,5000.000000000000000,0.000000000000000,1
19.007999999999999,49.664000000000001,5000.000000000000000,-3750.000000000000000,1
18.527999999999999,49.664000000000001,0.000000000000000,-3750.000000000000000,1

W sumie jak przeczytałem to wszystko jeszcze raz, to przyszło mi banalne rozwiązanie, ale nie wiem czy się ono sprawdzi.
Wczytanie mapy transformacją Helmerta i zdefiniowanie jej punktów na nowo, wtedy powinna chyba umieścić się prawidłowo.
Jednak będzie to proces troszkę żmudny, bo kawałki map są różnych wielkości (przeważnie bardzo małe). Może jest jakiś matematyczny wzór pod który mógłbym podstawić wielkości z pliku points i je po prostu podmienić?