Odpowiedzi:
Zawsze istnieje metoda „brutalnej siły”:
Weź warstwę ze znanym układem współrzędnych, który powinien nakładać się na twoją nieznaną warstwę.
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.
Promuję komentarze Mapperza i Brandona Copelanda, dodając odpowiedź, która wykorzystuje ich technikę.
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.
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.
Istnieją dwa świetne linki od Esri, które szczegółowo omawiają:
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!
Zapytaj producenta.
Jeśli znasz wspólne CRS dla swojego regionu geograficznego, możesz wypróbować niektóre z nich. Ale pytanie jest lepsze.
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.
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).
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!
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.
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ć:
Prawdopodobnie lepiej zacząć od pewnej pozycji znanych, niż próbować brutalnej siły od zera.
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.
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.
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.
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.
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 1984Nowsze 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.
- 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
- 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 odniesieniaW 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
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
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'}