Геометрич. размеры shp в QGIS, Google Map vs Matlab

[ДимонС]

Всем привет.

Я хотел задать вопрос профессионалам в этой области. Меня интересует конкретный вопрос, который я опишу по шагам.
1. Я открыл Google Map карту и нашел рандомное здание. Измерил ее длину прям там же, получил 292 м.
2. Открыл Openstreet Map плагин через QGIS, измерил длину там же, получил те же 292 м.
3. Сгенерировал shp файл как написано во многих источниках литературы и сохранил его в EPSG:3857 (OTF) стандарт. Везде написано, что проекция на плоскость с Google Map осуществляется с помощью этого стандарта. Ниже приведена расшифровка:

PROJCS["WGS_84_Pseudo_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Mercator"],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],PARAMETER["standard_parallel_1",6.64456744726493],UNIT["Meter",1]]


4. Я открываю в матлабе сгенерированный ранее shp файл, используя функцию shaperead.
5. Затем я мерю габариты здания, визуализированного через mapshow. Получаю число ~434м.

В течении 4 дней я возился с этой проблемой и до сих пор ее не решил. Откуда берется 434м, не пойму, когда все остальные браузерные карты показывают ~292м. Объясните, пожалуйста.

Я спрашивал этот вопрос на зарубежных сайтах. Мне сказали, следующее:

cos(47) = 0.68. 434*cos(47) = 295 m. My guess is that one software (Matlab) is calculating a 2D Cartesian length/distance. The more GIS-aware software is calculating the great circle/geodesic distance which is more correct.

Так как мне все-таки конвертировать данный в Матлабе, чтобы получить правильные габариты здания.

Благодарю за помощь!

[Denis Rykov]

Вполне ожидаемое поведение. Проекция 3857 не подходит для измерений. Либо измеряйте на эллипсоиде (это умеет делать почти весь ГИС софт), либо, если хотите мерить в единицах измерения слоя, то сохраните в более подходящую для этих задач проекцию.

[ДимонС]

Спасибо за ответ. А не могли бы вы ответить на два вопроса:
1. Как измерять дистанцию на проекции 3857? Harvestine formua для этого подходит?
2. На самом деле мне действительно нужна проекция на плоскости так, чтобы можно было мерить расстояние между объектами как dist = sqrt(x.^2 + y.^2). Не подскажите какую лучше выбрать?

[trir]

2. UTM, СК42

[Denis Rykov]

Сетку зон UTM можно скачать тут. Если ваши объекты разбросаны на несколько зон, то лучше проводить измерения на эллипсоиде или на сфере (подробнее тут).