Używasz niestandardowego układu współrzędnych w ArcGIS Desktop?

16

Nie wiem zbyt wiele o układach współrzędnych ...

W moim biurze zajmujemy się danymi przestrzennymi pochodzącymi ze stanowisk archeologicznych. Każda strona ma swój własny układ współrzędnych xyz (GCS). Trzy proste ortogonalne osie kartezjańskie . W ostatnich latach zarządzaliśmy tymi danymi przestrzennymi za pomocą oprogramowania GIS (ArcGIS), bez użycia określonego układu współrzędnych (po prostu pozostaw to jako „niezdefiniowane”)

Chciałbym wiedzieć, czy istnieje jakiś GCS zaprojektowany do obsługi takich zestawów danych przy użyciu prostej kartezjańskiej osi ortogonalnej, bez zniekształceń siatki typowego GCS. Ponadto chciałbym wiedzieć, czy ten system jest odpowiedni do używania go w aplikacji do mapowania online.

Nawiasem mówiąc, zarządzamy środowiskami 2D (ArcMap) i 3D (ArcScene) i pracujemy z „mm” jako jednostką bazową długości.

Jeśli coś takiego nie istnieje, być może ktoś wie, jak go stworzyć.

Maik
źródło

Odpowiedzi:

17

Zakładając, że twoje obszary zainteresowania są raczej małe w porównaniu do globu, możesz ustawić niestandardową projekcję poprzecznego mercatora.

Musisz znać współrzędne geograficzne lat_0 i lon_0 początku twojego CRS, a także kierunek osi xiy:

+proj=tmerc +lat_0=51.4 +lon_0=7 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

Do użytku z Arcgis plik .prj powinien zawierać:

PROJCS["Transverse_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",51.4],PARAMETER["central_meridian",7],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

Jeśli nie są one równoległe do kierunków północnego i wschodniego, konieczny byłby obrót za pomocą skośnego rzutu mercatora:

+proj=omerc +lat_0=51.4 +lonc=7 +alpha=-10 +k=1 +x_0=0 +y_0=0 +gamma=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

Plik .prj dla tej projekcji to:

PROJCS["Hotine_Oblique_Mercator_Azimuth_Center",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],PARAMETER["latitude_of_center",51.4],PARAMETER["longitude_of_center",7],PARAMETER["azimuth",-10],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

Oto oba razem (tmerc na niebiesko i omerc na czerwono):

wprowadź opis zdjęcia tutaj

Lat_0 i lon_0 (lub lonc) są współrzędnymi lokalnego pochodzenia CRS. Musisz je eksplorować za pomocą GPS lub z innych map. Kąt, który musisz obliczyć na podstawie współrzędnych dwóch znanych punktów na lokalnej osi X.

Jedyne, co nie działa zgodnie z oczekiwaniami, to milimetr jako jednostki. Powinien działać, ale nie działa w obecnym QGIS :-(

http://hub.qgis.org/issues/9414


Aktualizacja

Obecne wersje QGIS są teraz w stanie pracować z lokalnymi rzutami przy użyciu mm jako jednostek.

AndreJ
źródło