Разный результат при перепроецировании в ogr2ogr и ArcGIS

Вопросы по нескольким пакетам сразу, или вопросы, которые непонятно к какой ГИС отнести
Ответить
denis27
Интересующийся
Сообщения: 39
Зарегистрирован: 18 фев 2008, 10:08
Репутация: 1

Разный результат при перепроецировании в ogr2ogr и ArcGIS

Сообщение denis27 » 27 фев 2008, 12:23

Исходный файл MIF 100000 лист K52-011

Вот что про него говорит ArcCatalog:

Код: Выделить всё

Horizontal coordinate system
Projected coordinate system name: _MI_0
Geographic coordinate system name: PULKOVO-1942
Details 
Map Projection Name: Transverse Mercator
Scale Factor at Central Meridian: 1.000000
Longitude of Central Meridian: 129.000000
Latitude of Projection Origin: 0.000000
False Easting: 22500000.000000
False Northing: 0.000000

Planar Coordinate Information
Planar Distance Units: meters
Coordinate Encoding Method: coordinate pair

Geodetic Model
Horizontal Datum Name: D_PULKOVO-1942
Ellipsoid Name: unnamed
Semi-major Axis: 6378137.000000
Denominator of Flattening Ratio: 298.257224
А вот что написано в нем самом:

Код: Выделить всё

CoordSys Earth Projection 8, 1001, "m", 129, 0, 1, 22500000, 0
Перепроецирую его в Pulkovo 1942.prj

Код: Выделить всё

GEOGCS["GCS_Pulkovo_1942",
DATUM["D_Pulkovo_1942",
SPHEROID["Krasovsky_1940",6378245,298.3]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]]
Сначала с помощью ArcGIS потом с помощью FWTools2.1.0 вот такой командой

Код: Выделить всё

ogr2ogr.exe -nlt POLYGON -f "ESRI Shapefile" -skipfailures -t_srs "Pulkovo 1942.prj" "admin.shp" "admin.mif"
Результат получается разный, если открыть в ArcMap все три слоя исходный mif и shap-ы видно что AcrGIS-овский уезжает на северо восток примерно на 85 метров. А тот что спроецирован ogr2ogr.exe совпадает с исходным мифом. Пробовал также открывать в GlobalMapper - показывает ту же картинку.

Что я делаю не так?

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 28 фев 2008, 01:52

а зачем проецировать mif сначала в шейп, а потом обратно в mif?

ArcGIS-овский уезжает относительно чего?
пристегивайтесь, турбулентность прямо по курсу

denis27
Интересующийся
Сообщения: 39
Зарегистрирован: 18 фев 2008, 10:08
Репутация: 1

Сообщение denis27 » 28 фев 2008, 04:42

Нет обратно в миф не надо.

В общем ситуация такая - есть много мифов в разных зонах Гаусса крюгера разбитых по листам. Хочу собрать их в один большой шейп в географической проекции.

Для примера я отдельно нарисовал полигон на границе листа.
В архиве исходник и результат
http://dvride.ru/arcgis_ogr2ogr_difference.rar

При попытке открыть все три слоя в ArcMap получаем такую картинку
http://dvride.ru/arcmap_screenshot.gif

Сетка 1:100000 (1км) отсюда http://gis-lab.info/qa/topogrids.html

Может быть причина в том что аркинфо и ogr2ogr по разному интерпретируют содержимое исходного mif файла?
Последний раз редактировалось denis27 28 фев 2008, 05:09, всего редактировалось 1 раз.

denis27
Интересующийся
Сообщения: 39
Зарегистрирован: 18 фев 2008, 10:08
Репутация: 1

Сообщение denis27 » 28 фев 2008, 04:50

Похоже действительно проблема с интерпретаций MIF в Data Interoperability.

Сконвертировал MIF в TAB через MapInfo.
TAB - открылся в ArcMap там же где ogr2ogr-овский.
И напрямую и перепроицированный.

Перегнал TAB обратно в MIF через мапинфо - открываю в ArcMap уезжает.

Ответить

Вернуться в «Общий - ПО»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 2 гостя