Сравнение 2 шейп шайлов по "прочим" признакам

[Soitanen]

Добрый день!
У меня нестандартная проблема. Есть два шейп файла:
- нарезка побережья из OSM - http://yadi.sk/d/BhKGJMby6xSmG
- тот же район из v0_map - http://yadi.sk/d/PqqOm2js74JD8
Я вставляю один из файлов для дальнейшей обработки в одну утилиту (для симулятора FlightGear). Заданный квадрат для обработки - между 29 и 30 г в.д. и 59 и 60 г с.ш. Почему-то на файле v0_map все работает нормально, а вот на файле из ОСМ в итоге пустой результат. Дело вроде бы не в геометрии, потому что я могу задать квадрат со степью Ставропольской (где шейп из ОСМ будет представлять собой просто квадрат), но он также без всякой ошибки выдаст пустой результат.
Соответственно, у меня появилось подозрение, что может влиять какой-либо тип геометрии, или может проекции какие-то... Я пользуюсь QGIS, он прекрасно открывает оба файла, но разницу между ними я понять не могу. Помогите, пожалуйста, найти принципиальную разницу между двумя этими файлами.
Спасибо!

[Александр Мурый]

Посмотрим вывод утилиты <ogrinfo> для этих двух шейпов.

Спойлер

INFO: Open of `land_polygons-cut.shp'
using driver `ESRI Shapefile' successful.

Layer name: land_polygons-cut
Geometry: Polygon
Feature Count: 41
Extent: (29.000000, 59.000000) - (30.000000, 60.000000)
Layer SRS WKT:
PROJCS["WGS_84_Pseudo_Mercator",
GEOGCS["GCS_WGS_1984",
DATUM["WGS_1984",
SPHEROID["WGS_84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["Meter",1],
PARAMETER["latitude_of_origin",0.0]]

Спойлер

INFO: Open of `v0_landmass.shp'
using driver `ESRI Shapefile' successful.

Layer name: v0_landmass
Geometry: Polygon
Feature Count: 7
Extent: (27.322111, 55.701248) - (36.240273, 61.323326)
Layer SRS WKT:
GEOGCS["GCS_WGS_1984",
DATUM["WGS_1984",
SPHEROID["WGS_84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]]
OGC_FID: Real (11.0)
CAT: Real (11.0)

В шейпе "land_polygons-cut.shp" прописана проекция Меркатора (Mercator_1SP), а координаты явно географические в WGS84.

Так можно переназначить проекцию:

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

ogr2ogr -a_srs "EPSG:4326" land_polygons-cut_WGS84.shp  land_polygons-cut.shp

Вот архив с полученным шейпом.

land_polygons-cut_WGS84.zip

(67.96 КБ) 313 скачиваний

[Soitanen]

Огромноейшее спасибо! Все заработало!