Форумы GIS-Lab.info

Выложите координаты эталонных точек.

Evgeniy30 писал(а): ↑

24 янв 2020, 11:36

просто помогут правильно применить 7 параметров ортогонального преобразования

Раз вам нужна геодезическая точность, то и надо нанимать геодезистов. Для каждой точки свои 7-параметров, полученные на основе ближайших геодезических пунктов. С дивана это не делается.

Если такая точность преобразования устраивает - обращайтесь dr_gis @ mail.ru

Upd. Сократил информативность сообщения. Для оценки и одной точки достаточно.

freeExec писал(а): ↑

24 янв 2020, 12:49

Для каждой точки свои 7-параметров, полученные на основе ближайших геодезических пунктов. С дивана это не делается.

Даже если, как советовали выше, используется NTv2 — где координаты ближайших пунктов известны в обеих системах?

drrrrr писал(а): ↑

25 янв 2020, 11:52

Если такая точность преобразования устраивает

Каким образом получен результат?

Evgeniy30 писал(а): ↑

30 янв 2020, 22:47

Даже если, как советовали выше, используется NTv2 — где координаты ближайших пунктов известны в обеих системах?

Нет там координат, это все лишь поле поправок, которые были рассчитаны имея координаты в двух системах. Но вы его днём с огнём не сыщице.

Evgeniy30 писал(а): ↑

30 янв 2020, 22:47

Каким образом получен результат?

Результат получен преобразованием координат приведенных точек из местной системы координат (параметры подобрал сам, понятно, они отличаются от приведенных в http://sngo.mggt.ru/documents/sngo/PMSK_Moscow.pdf или используемых в ГИС вроде московской isogd.mos.ru и т.п.) в WGS-84. Расхождение при использовании этой СК для перепроецирования границ земельных участков, построенных из координат, приведенных в кадастровых выписках в местных координатах в прямоугольную WGS-84-Pseudo Mercator/Google Mercator (EPSG:3857) в которой отдает границы ЗУ WMS сервис ПКК (http://pkk5.rosreestr.ru/arcgis/service ... ervice=WMS) не превышает 2 см.
Для сопоставления выложите координаты нескольких своих эталонных точек (желательно пространственно разнесенных) в местной СК.

Evgeniy30 писал(а): ↑

24 янв 2020, 11:36

Правда, в начале топика надеялся, что просто помогут правильно применить 7 параметров ортогонального преобразования

Не раз эту тему обсуждал (хотя я тоже совсем далек от точного понимания работы с СК), это засекречено. Передо мной сидели специалисты МГГТ и просто улыбались на мои наивные вопросы "дескать, да, засекречено".

Непонятно одно. Как работу с такими допусками по точности, могли предложить выполнить ТС без оформления более других допусков?

evgeny.shirinyan писал(а): ↑

01 фев 2020, 19:09

Не раз эту тему обсуждал (хотя я тоже совсем далек от точного понимания работы с СК), это засекречено.

tikhpetr писал(а): ↑

01 фев 2020, 20:25

Как работу с такими допусками по точности, могли предложить выполнить ТС без оформления более других допусков?

Все данные для решения задачи, как ни странно, есть. Но что с ними делать — отдельный вопрос. Я просмотрел ГОСТы, описывающие алгоритмы преобразований (сложение матриц, вот это все), но не то чтобы все понял. И уж точно не понял, как реализовать программу для перевода понятными мне средствами (какой-нибудь PROJ4 или PostGIS).

Насколько я понимаю, приведенные данные для ортогонального преобразования — эмпирические (усредненные), а для действительных геодезических работ нужен расчет с привязкой к геодезическим пунктам. Тут советуют NTv2, но как реализовать этот алгоритм существующими средствами — куда больший вопрос.

Evgeniy30 писал(а): ↑

02 фев 2020, 14:04

Все данные для решения задачи, как ни странно, есть. Но что с ними делать — отдельный вопрос. Я просмотрел ГОСТы, описывающие алгоритмы преобразований (сложение матриц, вот это все), но не то чтобы все понял. И уж точно не понял, как реализовать программу для перевода понятными мне средствами (какой-нибудь PROJ4 или PostGIS).

Я тут посмотрел несколько роликов от Esri CIS, они делали адаптацию Аркгиса под работу с Росреестром. Ответ по точности: около 4 м. Все набегает из-за этих трансформаций из одной СК в другую через третью. https://youtu.be/-tUWPhG8F-M , https://blogs.esri-cis.ru/2018/02/14/ge ... ormations/ , https://blogs.esri-cis.ru/2018/08/06/ca ... X4I0CBThq0 Может, уже это видели.

В кугисе я вставил еще раз 1)параметры с mapbasic, 2) параметры с форума геодезистов из файла от МКА. И там, и там сдвиг векторных участков в МГГТ относительно WMS Росреестра около 3-4 м, в разные стороны.

Если делать, как говорилось выше, с дивана. А нам и надо с дивана, если это делается для визуализации (в InfraWorks можно было бы задавать этой план-схеме пользовательскую СК). По идее мы должны были бы работать с план-схемой, и тем же shp не задавать вообще никакой СК (это было возможно в старом ArcView, есть в Map3d, Mapinfo, но не получается нормально никак даже в новых версиях QGIS)

Тут советуют NTv2, но как реализовать этот алгоритм существующими средствами — куда больший вопрос.

Это вопрос к геодезистам - надо идти в поле и коорденировать все пункты. Я думал у вас есть уже такие материалы

Передо мной сидели специалисты МГГТ и просто улыбались на мои наивные вопросы "дескать, да, засекречено"

Просто денег хотят - http://sngo.mggt.ru/coord-systems.html

evgeny.shirinyan писал(а): ↑

02 фев 2020, 18:01

Ответ по точности: около 4 м. Все набегает из-за этих трансформаций из одной СК в другую через третью.

Да, сплошная путаница. Я специально выше спрашивал — один ли это переход; меня заверили, что один:

trir писал(а): ↑

24 янв 2020, 10:44

в каком порядке выполняется перевод? Формула, приведенная выше — это сумма параметров нескольких переходов (вроде ПМСК—СК-95—ПЗ90—WGS84), или переход всего один?

переход всего один

Хотя ГОСТы как раз описывают преобразования поэтапно.

trir писал(а): ↑

03 фев 2020, 06:45

Это вопрос к геодезистам - надо идти в поле и коорденировать все пункты. Я думал у вас есть уже такие материалы

Есть данные, которые формируют грид. По ним, насколько я понимаю, нужно будет пересчитывать координаты объектов — искать ближайшие точки, выполнять триангуляцию, вот это все.

Наверняка уже существуют какие-то программные инструменты, которым можно скормить грид в определенном формате, и сразу получить работающий инструмент.

evgeny.shirinyan писал(а): ↑

02 фев 2020, 18:01

параметры с mapbasic

А почему свет прямо-таки сошелся клином на параметрах Mapbasic? (То есть, я понимаю, что они позволяют получать очень и очень хорошие результаты, но подобная монополия выглядит странно.) Откуда у Mapbasic эти параметры? Как были получены?

В исходниках сайта можно найти строку PROJ4:

"+proj=tmerc +lat_0=55.66666666667 +lon_0=37.5 +k=1 +x_0=11.86143 +y_0=12.1761150 +ellps=bessel +towgs84=316.151,78.924,589.65,-1.57273,2.69209,2.34693,8.4507 +units=m +no_defs"

Люди с глазами говорят, что указанные параметры мало похожи на оригинальные

Evgeniy30,
Ну да, отличается от официальной МГГТ на +x_0=11.86143 +y_0=12.1761150, что уменьшает разницу с wms росреестра до ~3 метров, а если еще чуть подвинуть +x_0=14.86143 +y_0=12.7761150 - разница в пределах МКАД и вовсе уложится в полметра.
Росреестр, кстати, вчера свой WMS сервис прикрыл. Странные люди, этим сервисом многие пользовались, в частности, федеральный ФГИС ТП и популярный модуль для Qgis QMS.