Відповіді:
Завжди існує метод "грубої сили":
Візьміть шар із відомою системою координат, який повинен накладатися на ваш невідомий шар.
Тепер зробіть кілька освічених здогадок про те, якою може бути проекція невідомого шару. (UTM, Plate Carree тощо). Проектуйте свій відомий рівень системи координат у кожну проекцію, поки ви не знайдете той, який максимально відповідає невідомому шару.
Я просуваю коментарі Mapperz та Brandon Copeland, додаючи відповідь, яка використовує їх техніку.
Ця методика використовує той факт, що ArcMap може перетворити дані в нову систему координат в пам'ять. Дані без відомих-coordsys неможливо перетворити, тому вони просто відображаються. Проектуючи відомі дані в пам'ять, ми можемо швидко побачити, які координати вони мають у різних можливих системах координат.
Оскільки це питання ніколи не старіє, я створив сайт, який використовує метод Brute Force . Якщо ви перетягніть на карту блискавку shp + shx, вона відобразить її у кожній системі координат, доступній у PostGIS. Припускаючи, що ви знаєте, як виглядає "правильний", ви можете збільшити масштаб до цієї області та натиснути на полігон, щоб отримати файл .prj з epsg.io.
Є два чудові посилання від Есрі, які детально розглядають це питання:
На моїй останній роботі у мене був шар даних про геологію (багатокутник) під назвою "FSU_Geol.shp". Мій начальник дав мені це і попросив дізнатися низку речей. По-перше, клієнтові він передав цей файл файлу форми, а файлу .prj не було, тому він хотів, щоб я це зрозумів. Він також хотів, щоб я зрозумів, що таке категорії геології. Я міг би продовжувати і продовжувати, але давайте вирізати погоню .... Я погуглив її, і я закінчився ТУТ .
Скорочення в конвенції про іменування означало "геологію поверхні колишнього Радянського Союзу", і пошук Google привів мене "безпосередньо" до джерела (USGS). Все, про що я міг коли-небудь знати про це формату, було на верхньому посиланні, яке я потрапив. Я не кажу, що Google може знайти що завгодно і все , але я був недалеко від університету і просто робив «постріл у темряві», і подивіться на відгуки, які я отримав!
В іншому випадку хтось в офісі завантажив декілька файлів форм за допомогою пакетної функції. Наразі не пам’ятаю назву програмного пакету, але файли .prj відсутні. Я просто зайшов до метаданих в ArcCatalog, і я там дійсно знайшов свою відповідь. Знову ж таки, не цитуйте мене, але я думаю, що там була URL-адреса до джерела.
З цього дня я починаю з пошуку в Google, якщо в метаданих немає підказки!
Я дуже люблю цей веб-інструмент: http://projfinder.com/ . Перевірте свій файл і знайдіть деякі координати. Наведіть масштаб приблизно на те саме місце на земній кулі і дайте йому здогадуватися.
Не інструмент (я не знаю жодного з існуючих, який дозволив би вам це зробити), але перевірити відповідь @ mkennedy на те, як би я перетворив цю точку на WKID 4326? . Вона пояснює, як дійшла до правильної просторової орієнтації. SpatialReference.org і терпіння будуть вашими друзями.
Крім того, Esri пропонує посібник про те, як відгадати систему координат (хоча я віддаю перевагу методу mkennedy, якщо ви знаєте трохи більше про дані).
Найкраща відповідь, яку я знайшов на це питання, - нетехнічна: дізнайтеся, звідки беруться ваші дані. Агенції та організації, як правило, узгоджуються з використанням їх прогнозів. Знаєте, це походить від вашого штату ДОТ? Подивіться на решту їх даних і подивіться, що вони вам повідомляють. Не знаєте, звідки це? Освічена здогадка настільки ж ймовірна, що направить вас в потрібний шлях.
Принаймні, це дозволяє вирішити проблему з грубою силою трохи більше можливостей!
На цей момент серйозно застаріли, але Вернер Флакке і Біргіт Клаус розмістили Find Projection на ArcScripts в 2007 році. Я не думаю, що вихідний код є, на жаль. Це на основі VBA, тому його можна використовувати лише в ArcGIS Desktop v9.2 і, можливо, 9.3. Він включає два файли форми з цікавими областями з набору даних геодезичних параметрів EPSG, які можуть бути використані для звуження можливостей.
Географічний калькулятор блакитного мармуру та географічний трансформатор мають інструменти відновлення системи координат.
Спробуйте програму ogrinfo, що постачається як частина GDAL.
Див. Як отримати доступ до метаданих Shapefile за допомогою OGR?
Таким чином, ogrinfo може не надати вам інформацію про проекцію без файлу .prj, але це все ще корисний інструмент, який допоможе вам дослідити список потенційних прогнозів.
Наприклад:
Geometry: Polygon
Feature Count: 269
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000)
Ця відповідь підказує, що у файлі форм використовується система відліку, заснована на метрах, а не на градусах.
Якщо припустити, що ви знаєте, звідки приблизно дані, тепер у вас є коротший список можливих прогнозів.
Іншими найкращими практиками можуть бути:
Напевно, краще починати з якоїсь позиції знань, ніж намагатися з грубої сили з нуля.
Щоб додати до інших відповідей у цій публікації, я додам наступне:
Як: Ідентифікувати невідому проектовану систему координат за допомогою ArcMap
Надайте спеціальний чек на NAD1927
Якщо координати знаходяться в десяткових градусах, наприклад, між довжиною -180 та +180 та широтою -90 та +90, визначте Географічну систему координат (дату), що використовується для даних. До версії 9.2 ArcMap за замовчуванням призначає дані GCS_Assumed_Geographic_1. Це розміщує дані на даті NAD_1927
В додаток
Якщо дані знаходяться у Сполучених Штатах і показують ступінь, в якому координати зліва від десяткової значень складають 6, 7 або 8 цифр, дані, ймовірно, прогнозуються на системи координат Штату Плану або UTM.
FAQ: Основи прогнозування: Що повинен знати професіонал ГІС
На цій сторінці представлено кілька перевірок щодо визначення проекцій, що важливо допомагати в спробі виявити невідому систему координат.
Наступні концепції є основоположними для розуміння використання картографічних проекцій в ArcGIS. Зауважте, що тема прогнозів надзвичайно широка, і ця стаття може не лише торкнутися кількох важливих тем.
Системи координат, також відомі як картні проекції, є довільними позначеннями просторових даних. Їх мета - створити загальну основу для спілкування про певне місце чи місце на земній поверхні. Найважливішим питанням, пов’язаним з проекцією карт, є знання, що таке проекція, і наявність правильної інформації системи координат, пов'язаної з набором даних.
Коли були складені перші прогнози на карті, помилково було припущено, що земля рівна. Пізніше припущення було переглянуто, і земля вважалася ідеальною сферою. У 18 столітті люди почали розуміти, що земля не ідеально кругла. Це було початком концепції картографічного сфероїда.
Для більш точного представлення розташувань на земній поверхні, картографісти вивчили форму землі (геодезія) та створили концепцію сфероїда. Потім були розроблені географічні системи координат (GCS), які включають в себе дату, одиниці вимірювання та основний меридіан. Дана пов'язує сфероїд із певною частиною земної поверхні. Останні дані розроблені так, щоб добре вмістити всю земну поверхню.
Найпоширенішими даними в Північній Америці є:
• NAD 1927 (північноамериканська дата 1927) за допомогою сфероїда Clarke 1866
• NAD 1983 (North American Date 1983) з використанням сфероїда GRS 1980
• WGS 1984 (World Geodetic Survey 1984) з використанням сфероїда WGS 1984Новіші сфероїди розроблені за допомогою супутникових вимірювань і є більш точними, ніж ті, що були розроблені Кларком у 1866 р. Терміни "географічна система координат" та "дата" використовуються взаємозамінно, але, як зазначено вище, GCS включає дату, сфероїд, одиниці вимірювання. і основний меридіан.
- Координати змінюються в залежності від дати та сфероїда, на яких базуються ці координати, навіть якщо вони використовують однакові проекції карти та параметри.
Наприклад, наведені нижче географічні координати призначені для однієї точки, розташованої в межах міста Беллінгем, штат Вашингтон, з використанням 3 різних даних:
Код: X-координат DATUM Y-координата NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.7438798534299
- Принцип хорошого управління даними полягає в отриманні параметрів проекції з джерела даних, що надає дані. Не робіть обґрунтованих здогадок щодо прогнозування даних, оскільки результатом стане неточна база даних ГІС. Необхідні параметри:
• Проектування
• Одиниці вимірювання
• ЗОНА (для UTM)
• зона FIPS (для площини штату)
• ДатаІнші параметри можуть знадобитися, залежно від проекції. Наприклад, для проекції Альберса та Ламберта потрібні такі параметри:
• 1-а стандартна паралель, у градусах, хвилинах і секундах (DMS)
• 2-а стандартна паралель (DMS)
• Центральний меридіан (DMS)
• Широта походження проекцій (DMS)
• Неправдивий схід і одиниці вимірювання
• Неправдиві норми та одиниці виміру
• Зсув X та одиниці вимірювання
• Зсув Y та одиниці виміру
Подивіться на цей веб-сайт:
http://www.egger-gis.at/shapefile-projectionfinder/
Я розробник цього безкоштовного інструменту. Можливо, цей інструмент може допомогти вам знайти та визначити правильну проекцію вашої форми форми. Він ґрунтується на ідеї http://projfinder.com/ Аарона Рачико.
Ви можете також спробувати ці спеціальні рішення для:
Австрія (німецька): https://www.data.gv.at/anwendungen/checkaustrianprojection/
Австралія: https://maegger.github.io/map_australia.html
UTM - Зони: https://maegger.github.io/map_utm.html
Ви можете це зробити з Фіоною.
import fiona as f
a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
print(a.crs)
Ви отримаєте відповідь як
{'init': 'epsg:4326'}