Страница 1 из 1
GEOGCS - простой вопрос
Добавлено: 01 фев 2009, 08:47
Boris
У меня есть данные в SHP-файлах, и у большинства из них имеется файл PRJ вот такого содержания:
GEOGCS["Lat Long for MAPINFO type 0 Datum",DATUM["D_MAPINFO",SPHEROID["World_Geodetic_System_of_1984",
6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]]
Проясните два вопроса:
1. где глянуть описание формата записи проекции? а то как-то сразу не нашел
2. такое мудреное название соответствует просто Lat/Long на WGS84? Смутило обилие названий MAPINFO - это просто названия системы координат и датума?
Re: GEOGCS - простой вопрос
Добавлено: 01 фев 2009, 09:27
Максим Дубинин
1. у
ESRI
2. да, там где Mapinfo можно хоть неприлично писать.
Re: GEOGCS - простой вопрос
Добавлено: 01 фев 2009, 11:33
Boris
Спасибо. Прочитаю.
Re: GEOGCS - простой вопрос
Добавлено: 02 фев 2009, 11:31
geologic
Mapinfo всего лишь название, но оно означает, что ваши prj родом оттуда. Бывают случаи, когда такие файлы не всегда читаются правильно в ArcMap. Проблемы совместимости при загрузке в Arc таких файлов могут быть, лучше заменить prj на родные Arc-овские, или хотя бы сопоставить название датума.
Если захотите наоборот, конвертить данные в Mapinfo, то учтите, что ogr, например, умеет в Tab-Map вписывать проекцию из shp-prj, после чего Arc-проекция не всегда разпознается Мапинфой. В этот момент подставляете "мапинфовский" prj и дело в шляпе.
Re: GEOGCS - простой вопрос
Добавлено: 02 фев 2009, 18:13
Boris
2 geologic
geologic писал(а):Mapinfo всего лишь название...
...
и дело в шляпе.
Я бы сказал, что вы меня сильно-сильно запутали.
Теперь ответьте - судя по тексту PRJ файла - это система координат LL на WGS84? Желательно, да или нет.
Дело в том, что у Mapinfo есть такая система координат:
"Longitude / Latitude", 1, 0
а есть такая
"Longitude / Latitude (WGS 84)\p4326", 1, 104
Мне надо было быть уверенным, что их не перепутали. А то по данной работе много разных данных сходится от разных организаций, и все как правило слабо документированные. Большей частью в SHP-ах, многие без PRJ, и практически по всем видно, что SHP - это не их родной формат.
Re: GEOGCS - простой вопрос
Добавлено: 02 фев 2009, 18:32
Boris
Во истину вышло - спроси умных людей, будешь знать ответ на свой вопрос. Причем именно на тот, который нужно узнать, а не в тот, что сумел спросить.
Последний совет надоумил меня проверить наконец-то, что такое в Mapinfo проекция
"Longitude / Latitude", 1, 0
Файлы настроек FME Universal Translator описывают его вот так:
Coordinate System Parameters
============================
CS_NAME: LL-MAPINFO
DESC_NM: MapInfo Generic Lat/Long
DT_NAME: MAPINFO
GROUP: MAPINFO
MAP_SCL: 1
PARM1: -180
PARM2: 180
PROJ: LL
QUAD: 1
SCL_RED: 1
SOURCE: MapInfo
UNIT: DEGREE
Datum Parameters
================
DELTA_Z: 1e-098
DESC_NM: MapInfo Default Datum (type 0)
ELLIPSOID: WGS84
SOURCE: MapInfo
USE: WGS84
Ellipsoid Parameters
====================
DESC_NM: World Geodetic System of 1984
E_RAD: 6378137
P_RAD: 6356752.3142
SOURCE: US Defense Mapping Agency, TR-8350.2-B, December 1987
OGC WKT Description
===================
GEOGCS["MapInfo Generic Lat/Long",
DATUM["MAPINFO",
SPHEROID["World Geodetic System of 1984",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]]
ESRI WKT Description
===================
GEOGCS["MapInfo Generic Lat/Long",DATUM["D_MAPINFO",SPHEROID["World_Geodetic_System_of_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
В какой-то литературе эта проекция описывалась как "геоцентрический шар", правда параметры шара не приводились.
Пришлось построить в Mapinfo две таблицы в проекциях
"Longitude / Latitude", 1, 0
"Longitude / Latitude (WGS 84)\p4326", 1, 104
соответственно, состоящие из точек от 0 до 80 градуса широты и от 35 до 41 долготы с шагом в 1 градус.
Потом высчитать для каждой из них координаты X;Y в UTM37/WGS84.
И вычислить квадратичное отклонение между точками в разных таблицах.
Пересчет координат получился с точностью в дециметрах. Разность координат равна 0.
Что и требовалось узнать. Обе системы координат - идентичны. Удивительно, что для первой из них не указан номер ESPG. И как я понимаю, проекция № 0, в MapInfo появилась раньше, чем была окончательно принята WGS84.
Re: GEOGCS - простой вопрос
Добавлено: 02 фев 2009, 19:01
Каминский Вадим
Во истину не понятно в чем была проблема! Boris вы все говорите о проекции, но из указанного prj файла делается вывод о том, что о проекции нет нислова, данные неспроецированы и хранятся в геодезических координатах (широта, долгота) на эллипсоиде WGS-84. Соответственно самое главное в названии ваших файлов привязки - это параметры WGS-84, остальное ни на что не влияет. Если пришли данные с разными названиями WGS-84 проверьте значения большой полуоси и эксцентриситета, вдруг кто-то ошибся

Re: GEOGCS - простой вопрос
Добавлено: 02 фев 2009, 19:42
Boris
Проекция, да, извините не точное название, - но это не мое изобретение. Mapinfo их все называет Projection. А если честно, и не знаю как их правильно обобщенно называть. Система координат, то же на мой взгляд не очень полное название, и проекция короче в написании.
А вот на счет полуосей - это очень ценный совет. Сейчас проверю, теперь я понимаю, что в том что я спрашивал - не два ничего не значащих названия, три.
Я уже имел по этой работе карту Минрегиона в формате MapInfo, в которую банда государственных векторизаторов запихнула первую попавшуюся географическую систему координат. Естественно "Longitude / Latitude", 1, 0
К счастью было с чем сравнивать, что бы понять, что это был эллипсоид Красовского.