Identyfikowanie układu współrzędnych pliku kształtu, gdy jest nieznany?

Odpowiedzi:

49

Zawsze istnieje metoda „brutalnej siły”:

  1. Weź warstwę ze znanym układem współrzędnych, który powinien nakładać się na twoją nieznaną warstwę.

  2. Teraz dokonaj pewnych świadomych domysłów na temat tego, jaką projekcją może być nieznana warstwa. (UTM, Plate Carree itp.). Projektuj znaną warstwę układu współrzędnych na każdą projekcję, aż znajdziesz taką, która pasuje do nieznanej warstwy w jak największym stopniu.

brenth
źródło
1
to także moja metoda. Używam dwóch różnych programów, aby to „udowodnić”. Możesz również zapoznać się z tym pytaniem w celu uzyskania podstawowych wskazówek dotyczących projekcji. gis.stackexchange.com/questions/2769/…
Brad Nesom
6
Możesz trochę uprościć krok 2, po prostu zmieniając rzut ramki danych (ArcMap) i pozwalając znanej zmianie warstwy na bieżąco.
Brandon Copeland
13
Pamiętaj, że nie powinieneś wyświetlać danych. Powinieneś definiować projekcję. Rzutowanie danych zmienia współrzędne. Definiowanie projekcji po prostu przypisuje projekcję do tych danych i informuje oprogramowanie, w jakiej projekcji są już dane. Są to dwa różne narzędzia w ArcMap, „Zdefiniuj projekcję” jest tym, którego chcesz użyć.
Sean
Wolę czytać wartości współrzędnych, zwykle mogę stwierdzić, czy są to stopnie dziesiętne, metry czy stopy (trudniej rozróżnić), a następnie przejść do najlepszego układu. Zacznij od WGS 84, jeśli wydaje się dziesiętny, UTM, jeśli liczniki, stan płaszczyzny, jeśli stopy. Kiedyś miałem taki, którego nikt nie mógł naprawić (w przypadku anomalii grawitacyjnych). Był to amerykański samolot państwowy o jednostce liniowej ...... kilometrów. Z faktu, że było blisko, mogłem stwierdzić, że system miał rację, ale stan był niewielki, ahh km to jednostka liniowa.
Jeśli nie wiesz - po prostu GIS
33

Promuję komentarze Mapperza i Brandona Copelanda, dodając odpowiedź, która wykorzystuje ich technikę.

  1. Jeśli masz ArcGIS Desktop, dodaj do ArcMap dane referencyjne obejmujące ten sam obszar. Te dane odniesienia muszą mieć poprawną definicję układu współrzędnych (rzutowania).
  2. Dodaj dane do nieznanego układu współrzędnych.
  3. Zbadaj, jakie możliwe układy współrzędnych są stosowane w tym obszarze. Jednym miejscem do sprawdzenia jest rejestr EPSG .
  4. Ustaw układ współrzędnych ramki danych zgodnie z możliwościami. Jeśli dane ze znanymi współrzędnymi pokrywają się z danymi bez znanego współrzędnych, prawdopodobnie zidentyfikowano układ współrzędnych.

Ta technika wykorzystuje fakt, że ArcMap może konwertować dane do nowego układu współrzędnych w pamięci. Dane bez znanego współrzędnego nie mogą być konwertowane, więc są po prostu wyświetlane. Projektując znane dane w pamięci, możemy szybko zobaczyć, jakie współrzędne ma w różnych możliwych układach współrzędnych.

Mkennedy
źródło
4
Uwaga dodatkowa: W QGIS przerzutowanie w pamięci nazywa się „w locie”.
podmroku
3
W ArcGIS nazywa się to również „w locie”! Ciągle znajduję wielu ludzi, którzy nie rozumieją tego wyrażenia, więc zamiast tego próbowałem użyć słowa „w pamięci”.
mkennedy
Czy ten link jest nadal aktualny? Wydaje się nieco bardziej skomplikowany niż to, jak go pamiętam, jednak istnieje epsg-registry.org
GISKid
1
@GISKid Masz rację. Wydaje mi się, że strona epsg-registry.org jest późniejsza, kiedy napisałem tę odpowiedź. Zaktualizuję.
mkennedy
1
Kolejny wielki zasób: epsg.io
GISKid
29

Ponieważ to pytanie nigdy się nie zestarzeje, stworzyłem witrynę, która wykonuje metodę Brute Force . Przeciągnięcie spakowanego shp + shx na mapę spowoduje odwzorowanie go w każdym układzie współrzędnych dostępnym w PostGIS. Zakładając, że wiesz, jak wygląda „poprawny”, możesz powiększyć ten obszar i kliknąć wielokąt, aby pobrać plik .prj z epsg.io.

jmk
źródło
15
świetny pomysł, ale myślę, że strona tak naprawdę nie działa
Fábio Dias,
22

Podczas mojej ostatniej pracy miałem warstwę danych geologicznych (wielokąt) o nazwie „FSU_Geol.shp”. Mój szef dał mi go i poprosił, abym dowiedział się wielu rzeczy. Po pierwsze, klient przekazał mu ten plik kształtu i nie było pliku .prj, więc chciał, żebym to zrozumiał. Chciał również, abym zorientował się, jakie są kategorie geologii. Mógłbym iść naprzód, ale przejdźmy do sedna ... Przejrzałem go i skończyłem TUTAJ .

Akronim w konwencji nazewnictwa oznaczał „geologię powierzchniową byłego Związku Radzieckiego”, a wyszukiwanie Google doprowadziło mnie „bezpośrednio” do źródła (USGS). Wszystko, co mogłem kiedykolwiek wiedzieć o tym pliku kształtu, znajdowało się na górnym linku, który trafiłem. Nie twierdzę, że Google może znaleźć wszystko i wszystko , ale byłem świeżo po studiach i po prostu zrobiłem „zastrzyk w ciemności” i popatrzyłem na opinie, które otrzymałem!

W innym przypadku ktoś w biurze pobrał wiele plików kształtów za pomocą funkcji wsadowej. W tej chwili nie pamiętam nazwy pakietu oprogramowania, ale brakowało plików .prj. Po prostu przeszedłem do metadanych w ArcCatalog i znalazłem tam swoją odpowiedź. Ponownie, nie cytuj mnie, ale myślę, że był tam adres URL źródła.

Od tego dnia „zaczynam” wyszukiwanie w Google, jeśli w metadanych nie ma żadnych wskazówek!

Dano
źródło
18

Zapytaj producenta.

Jeśli znasz wspólne CRS dla swojego regionu geograficznego, możesz wypróbować niektóre z nich. Ale pytanie jest lepsze.

podmrok
źródło
Jedynym problemem jest to, że nie mogę go zapytać ... To była również pierwsza myśl, którą miałem.
com
Potem będziesz musiał użyć brutalnej siły, jak powiedział brenth.
podmrok
18

Naprawdę uwielbiam to narzędzie internetowe: http://projfinder.com/ . Sprawdź plik i wyszukaj współrzędne. Przybliż do mniej więcej tego samego miejsca na kuli ziemskiej i niech zgadnie.

johanvdw
źródło
To właściwie najlepsze narzędzie, jakie znalazłem. Zgadywanie w ArcMap nie jest naprawdę przydatne na dłuższą metę lub jeśli nie wiesz, co powinno być. Nie może pomóc w układach współrzędnych homebrew i braku jednego, ale kto może!
Martin
Szkoda, że ​​to już nie działa - ktoś wie, z kim się skontaktować?
sventechie
16

Nie narzędzie (nie wiem od jednego w istnieniu, który pozwoli Ci zrobić), ale sprawdzić odpowiedź @ mkennedy za do Jak przekonwertować ten punkt do WKID 4326? . Wyjaśnia, jak doszła do właściwego odniesienia przestrzennego. SpatialReference.org i cierpliwość będą twoimi przyjaciółmi.

Ponadto Esri zawiera przewodnik na temat odgadywania układu współrzędnych (chociaż wolę metodę Mkennedy'ego, jeśli wiesz trochę więcej na temat danych).

Sasa Ivetic
źródło
13

Najlepsza odpowiedź, jaką znalazłem na to pytanie, jest nietechniczna: dowiedz się, skąd pochodzą Twoje dane. Agencje i organizacje zwykle zachowują spójność ze swoimi prognozami. Wiesz, że pochodzi z twojego stanu DOT? Spójrz na resztę ich danych i zobacz, co ci mówi. Nie wiesz skąd się wziął? Wykształcone przypuszczenie równie dobrze może skierować cię właściwą drogą.

Przynajmniej sprawia, że ​​radzenie sobie z problemem brutalnej siły jest trochę bardziej wykonalne!

Jason Baker
źródło
11

Poważnie nieaktualne w tym momencie, ale Werner Flacke i Birgit Klaus opublikowali Find Projection na ArcScripts w 2007 roku. Niestety nie sądzę, że kod źródłowy istnieje. Jest oparty na VBA, więc można go używać tylko w ArcGIS Desktop v9.2 i ewentualnie 9.3. Obejmuje on dwa pliki kształtów z obszarami zainteresowania z zestawu danych parametrów geodezyjnych EPSG, które można wykorzystać do zawężenia możliwości.

Blue Marble Geographic Calculator i Geographic Transformer mają narzędzia do odzyskiwania systemu współrzędnych.

Mkennedy
źródło
Dzięki za wskazówkę Melita. Kod jest tam (* .frm). Licencja jest własnością publiczną. PrjFinder jest towarzyszem książki „Working with Projections and Datum Transformations in ArcGIS” (2006).
matt wilkie 30.08.11
4

Wypróbuj program ogrinfo dostarczony jako część GDAL.

Zobacz Jak uzyskać dostęp do metadanych Shapefile przy użyciu OGR?

Więc ogrinfo może nie dać ci informacji o projekcji bez pliku .prj, ale nadal jest użytecznym narzędziem, które pomoże ci zbadać listę potencjalnych projekcji.

Na przykład:

Geometry: Polygon
Feature Count: 269 
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000) 

Ta odpowiedź sugeruje mi, że plik kształtu używa systemu odniesienia opartego na metrach, a nie na stopniach.

Zakładając, że wiesz, skąd pochodzą dane, masz teraz krótszą listę możliwych prognoz.

Inne najlepsze praktyki mogą obejmować:

  • Zapytaj przede wszystkim osoby, które przekazały Ci dane, nawet jeśli nie wiedzą dokładnie, że powinny mieć listę możliwych projekcji, które mogłyby zostać wyświetlone w pierwszej kolejności.
  • Sprawdź, czy są jakieś powiązane metadane.

Prawdopodobnie lepiej zacząć od pewnej pozycji znanych, niż próbować brutalnej siły od zera.

nmtoken
źródło
1
To nie odpowiada na zadane pytanie.
Dan C
1
Jeśli nie ma pliku prj, hen ogrinfo nie zwróci żadnych wyników dla informacji o układzie współrzędnych.
Devdatta Tengshe,
4

Aby dodać do innych odpowiedzi w tym poście, dodałbym następujące:

Instrukcje: zidentyfikowanie nieznanego rzutowanego układu współrzędnych za pomocą ArcMap

Podaj szczególną kontrolę na NAD1927

Jeśli współrzędne są w stopniach dziesiętnych, na przykład między długością -180 a +180 oraz szerokością geograficzną -90 i +90, zidentyfikuj układ współrzędnych geograficznych (układ odniesienia) użyty dla danych. Przed wersją 9.2 ArcMap domyślnie przypisuje GCS_Assumed_Geographic_1 do danych. To umieszcza dane w układzie odniesienia NAD_1927

Dodatkowo

Jeśli dane znajdują się w Stanach Zjednoczonych i pokazują zakres, w którym współrzędne po lewej stronie dziesiętnej wynoszą 6, 7 lub 8 cyfr, dane prawdopodobnie są rzutowane na układ współrzędnych płaszczyzny stanu lub UTM.

FAQ: Podstawy projekcji: Co powinien wiedzieć specjalista GIS

Ta strona zawiera kilka sprawdzeń dotyczących określania rzutów, co jest ważne podczas próby zidentyfikowania nieznanego układu współrzędnych.

Poniższe pojęcia są fundamentalne dla zrozumienia zastosowania rzutów map w ArcGIS. Należy jednak pamiętać, że temat prognoz jest bardzo szeroki, a ten artykuł może jedynie poruszyć kilka ważnych tematów.

  1. Układy współrzędnych, znane również jako rzuty mapy, są dowolnymi oznaczeniami danych przestrzennych. Ich celem jest zapewnienie wspólnej podstawy komunikacji na temat określonego miejsca lub obszaru na powierzchni ziemi. Najważniejszym problemem związanym z rzutowaniem mapy jest wiedza o tym, czym jest rzut i posiadanie prawidłowych informacji o układzie współrzędnych powiązanych z zestawem danych.

  2. Kiedy opracowano pierwsze rzuty mapy, błędnie założono, że ziemia była płaska. Później założenie to zostało zmienione i przyjęto, że Ziemia jest idealną kulą. W XVIII wieku ludzie zaczęli zdawać sobie sprawę, że ziemia nie była idealnie okrągła. To był początek koncepcji sferoidu kartograficznego.

  3. Aby dokładniej przedstawić miejsca na powierzchni Ziemi, twórcy map zbadali jej kształt (geodezję) i stworzyli koncepcję sferoidy. Następnie opracowano układy współrzędnych geograficznych (GCS), które obejmują układ odniesienia, jednostki miary i główny południk. Układ odniesienia łączy sferoidę z określoną częścią powierzchni ziemi. Najnowsze dane zostały zaprojektowane tak, aby dobrze pasowały do ​​całej powierzchni ziemi.

  4. Najczęściej używane dane w Ameryce Północnej to:

• NAD 1927 (North American Datum 1927) przy użyciu sferoidy Clarke 1866
• NAD 1983 (North American Datum 1983) przy użyciu sferoidy GRS 1980
• WGS 1984 (World Geodetic Survey 1984) przy użyciu sferoidy WGS 1984

Nowsze sferoidy są opracowywane na podstawie pomiarów satelitarnych i są bardziej dokładne niż te opracowane przez Clarke'a w 1866 r. Terminy „układ współrzędnych geograficznych” i „układ odniesienia” są używane zamiennie, ale jak wspomniano powyżej, GCS obejmuje układ odniesienia, sferoidę, jednostki miary i główny południk.

  1. Współrzędne danych zmieniają się w zależności od układu odniesienia i sferoidy, na których oparte są te współrzędne, nawet jeśli używają tych samych rzutów mapy i parametrów.

Na przykład poniższe współrzędne geograficzne dotyczą jednego punktu położonego w mieście Bellingham w stanie Waszyngton przy użyciu 3 różnych punktów odniesienia:

Kod: DATUM Współrzędna X Współrzędna Y NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.7438798534299

  1. Zasadą dobrego zarządzania danymi jest uzyskanie parametrów projekcji ze źródła danych dostarczającego dane. Nie zgaduj w sposób przemyślany na temat projekcji danych, ponieważ wynikiem tego będzie niedokładna baza danych GIS. Niezbędne parametry są następujące:

• Projekcja
• Jednostki miary
• STREFA (dla UTM)
• Strefa FIPS (dla Płaszczyzny stanu)
• Punkt odniesienia

W zależności od projekcji mogą być wymagane inne parametry. Na przykład projekcje Albersa i Lamberta wymagają następujących parametrów:

• 1. standardowa równoległa, w stopniach, minutach i sekundach (DMS)
• 2. standardowa równoległa (DMS)
• Południk środkowy (DMS)
• Początek szerokości geograficznej rzutów (DMS)
• Fałszywy wschód i jednostki miary
• Fałszywy północ i jednostki miary
• Przesunięcie X i jednostki miary
• Przesunięcie Y i jednostki miary

whyzar
źródło
Co mam zrobić, jeśli mam warstwę siatki, która pokazuje nieznany układ współrzędnych i nie mogę mu nic przypisać? Dzięki.
AndrewLebron,
2

Miałem ten sam problem podczas pracy nad moim projektem (DEM). Jedna z DEM była nieznanymi współrzędnymi i nie można jej rzutować, więc zrobiłem „georeferencję” za pomocą narzędzia do georeferencji w ArcMap . Po georeferencji otrzymasz plik prj.file.

ndr
źródło
0

Spójrz na tę stronę:

http://www.egger-gis.at/shapefile-projectionfinder/

Jestem twórcą tego bezpłatnego narzędzia. Być może to narzędzie może pomóc w znalezieniu i zdefiniowaniu prawidłowej projekcji pliku kształtu. Opiera się na pomyśle http://projfinder.com/ autorstwa Aarona Racicota.

Możesz również wypróbować te specjalne rozwiązania dla:

Austria (niemiecki): https://www.data.gv.at/anwendungen/checkaustrianprojection/

Australia: https://maegger.github.io/map_australia.html

UTM - Strefy: https://maegger.github.io/map_utm.html

Manfred Egger
źródło
0

Możesz rozróżnić współrzędne geograficzne i kartograficzne, pierwsze łaty / lon są w stopniach, a drugie północy / wschody są w metrach. Jeśli plik ma kształt, możesz zrzucić współrzędne geometrii za pomocą shapelib-1.3.0 i zobaczyć wartości. ..

natan
źródło
0

Możesz to zrobić z Fioną.

    import fiona as f
    a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
    print(a.crs)

Otrzymasz odpowiedź jako

    {'init': 'epsg:4326'}
Sandip Nepal
źródło
Czy na pewno zadziała bez pliku .prj? skąd Fiona się dowie?
JR
Nie bez pliku .prj, nie będzie działać
Sandip Nepal