Карты участков в план-схеме на карте Росреестра в MapInfo

MapInfo, MapBasic
Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2227
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1500
Откуда: Казань

Re: Карты Yandex и Google в MapInfo

Сообщение Игорь Белов » 24 дек 2015, 09:21

fvgbh писал(а):а поделитесь секретом как это вы вычисляете? :)
Непременно. Если дадите знать, устраивает ли результат. :wink:
The purpose of computing is insight, not numbers

fvgbh
Новоприбывший
Сообщения: 7
Зарегистрирован: 20 дек 2015, 15:35
Репутация: 0

Re: Карты Yandex и Google в MapInfo

Сообщение fvgbh » 26 дек 2015, 21:33

не совсем то что надо, вот результаты перевода:
исходный участок в МСК
29787.45 25219.59
29807.85 25220.29
29806.39 25260.62
29786.21 25259.79
29786.88 25238.98
29787.45 25219.59

он же переведенный в СК63 (эталон)
4935962.00000 4382893.96000
4935982.41000 4382894.54000
4935981.19000 4382934.88000
4935961.01000 4382934.17000
4935961.55000 4382913.36000
4935962.00000 4382893.96000

а при использовании параметров перевода "MCK-91", 8, 1001, 7, 33, 0, 1, -18000, -4915000 в "1963_X3 зона 4", 8, 1001, 7, 32.5, 0, 1, 4300000, -9414.70
получается вот так:
4382892.15 4935759.4
4382892.73 4935779.8
4382933.07 4935778.59
4382932.37 4935758.4
4382911.56 4935758.95
4382892.15 4935759.4

если из "MCK-91", 8, 1001, 7, 33, 0, 1, -18001.82, -4915002.58 в "1963_X3 зона 4", 8, 1001, 7, 32.5, 0, 1, 4300000, -9414.70
то вот так:
4382893.96 4935761.99
4382894.54 4935782.39
4382934.88 4935781.18
4382934.17 4935760.99
4382913.36 4935761.54
4382893.96 4935761.99
т.е по Y все ок совпадает, а по Х нет( (Х и Y тут переставлены местами 49... это Х 43.... это Y)

а если использовать параметры такие: "MCK-91", 8, 1001, 7, 33, 0, 1, -18001.82, -4915002.58 в "1963_X3 зона 4 вариант", 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, 7, 32.5, 0.083333333333333333333333333333333, 1, 4300000, 0
то получается:
4382895.23 4935960.45
4382895.81 4935980.84
4382936.16 4935979.63
4382935.45 4935959.46
4382914.63 4935960
4382895.23 4935960.45
здесь по иксам и игрекам чуть-чуть не дотягивает до эталона

т.е параметры "MCK-91", 8, 1001, 7, 33, 0, 1, -18000, -4915000 нуждаются в корректировке

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2227
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1500
Откуда: Казань

Re: Карты Yandex и Google в MapInfo

Сообщение Игорь Белов » 27 дек 2015, 00:10

В топографии всё и всегда нуждается в корректировке. Границы земельных участков между соседними районами не стыкуются зачастую не на пару метров, а много больше. Вот и в Вашей задачке тоже встретились СК города и республики.

Кроме того, возможно использование особого датума, созданного специально для геодезических и прочих работ на объекте. В таком случае вместо номера 1001 следует использовать хитрые цифры, одни для СК-63, другие для МСК.

В любом случае что-то корректировать придётся. Вряд ли от этого всё станет точнее. Просто наборы координат будут как бы согласованы.

Между прочим, у меня из СК-63 ("CS63 zone X4", 8, 1001, 7, 32.5, 0.0833333333333333, 1, 4300000, 0) в местную с подправленными параметрами сдвига ("MCK-91", 8, 1001, 7, 33, 0, 1, -18001.82, -4915002.58) все точки "эталона" пересчитываются сантиметр в сантиметр.

Кстати об эталоне. Вы сказали, что координаты получены переводом. Было бы здорово выйти с вопросами на этого переводчика.
The purpose of computing is insight, not numbers

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2227
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1500
Откуда: Казань

Ответ на вопрос «как»

Сообщение Игорь Белов » 27 дек 2015, 13:07

Просьба к модераторам перенести последние посты в отдельную тему.
====
Ответ на вопрос «как».

Предоставленные данные: слой кадастровых кварталов в местной системе координат.
Задача: определить параметры проекции.

Имеются программы, которые умеют делать:
  • конвертирование координат в картографических проекциях
  • вычисление параметров конформного преобразования на плоскости
Открываем слой кварталов в карте. Нужны те же данные на глобусе. В другой карте открываем кадастровый слой сервиса веб-карт Росреестра. Находим эти кварталы. В каждом окне создаём новый векторный слой. Проекция первого — МСК, второго — EPSG:3857. Равномерно по всей площади исходного слоя создаём точки на углах кварталов. В исходной карте они должны прилипать к узлам полигонов. В карте Росрееста тщательно прицеливаемся при зуме до масштаба 1:141.

Экспортируем результаты в текстовый файл с координатами. Первый слой в координатах МСК:
44728.174996844267298 3320.422987725022267
48620.485990035311261 15983.358001437964049
14098.864999825218547 22277.095007283980522
31131.646000856319006 45996.24199719344324
30919.928509850578848 26201.10451288792683
29525.749993967954651 17983.320004405981308
26314.190010227906896 30790.110006726576103

Второй слой выводим в двух системах. "СК-42 зона 6":
6562729.793782969 4918323.8591784192
6566622.108865004 4930986.8215941722
6532100.4826228367 4937280.5503298976
6549133.2689905353 4960999.6918238997
6548921.5379769998 4941204.5575788934
6547527.3665311718 4932986.7749165995
6544315.803866528 4945793.5579636665

"СК-63 зона X4":
4402565.3337153988 4909613.8265660787
4406380.4877342982 4922301.3944166964
4371818.6427516174 4928384.0862910813
4388706.5577096445 4952208.700989048
4388616.3314933265 4932411.3121750057
4387272.4270656202 4924184.5980842886
4383982.2493979409 4936972.2852538209

В обеих системах для перехода от Пулково 1942 к WGS 84 использован набор параметров ГОСТ Р 51794-2001:
23.92, -141.27, -80.9, 0, -0.35, -0.82, -0.12

Вычисляем параметры конформного преобразования из СК-42 в МСК и получаем невязки:
0.000 -0.012
0.003 0.014
0.003 0.003
0.005 -0.005
-0.008 0.001
-0.000 0.004
-0.003 -0.005

Невязки преобразования из СК-63 в МСК:
-0.106 -0.169
0.109 -0.028
0.204 0.075
-0.235 -0.006
0.007 0.050
0.012 0.081
0.009 -0.003

Выводы:
  1. городская система координат построена на основе СК-42;
  2. предоставленный слой кварталов однороден с данными на сайте Росреестра.
Смотрим на параметры преобразования из СК-42 в МСК:

Разворот,градусы: -1.9971332226009363e-06
Масштаб: 0.99999985703266503
Сдвиг по X: -6518000.5088546816
Сдвиг по Y: -4915002.9738647649

Практически разворот ноль, масштаб единица, величины сдвига целые километры. Это просто подарок! МСК — это СК-42 с простым смещением. Взглянем на координаты в МСК и в СК-42. Так и есть, с точностью до сдвига в первые метры они отличаются на 6518000 и 4915000 метров. Средние значения:
-6518001.61744778 -4915003.45155256

Оставшиеся метры можно объяснить тем, что использованный Росреестром датум отличается от ГОСТ Р 51794-2001. Имеются и более мрачные объяснения. В общем, можно смело округлить параметры сдвига. Параметры проекции:

Latitude of natural origin: 0
Longitude of natural origin: 33
Scale factor at natural origin: 1
False Easting: -18000
False Northing: -4915000

На всякий случай тестировал СК-95 с датумом по ГОСТ 32453-2013:
24.47, -130.89, -81.56, 0, 0, -0.13, -0.22

Получилось не очень:

Разворот,градусы: 0.00017057098897460825
Масштаб: 0.99999922761820415
Сдвиг по X: -6518013.991140632
Сдвиг по Y: -4914975.986209875
====
The purpose of computing is insight, not numbers

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: Ответ на вопрос «как»

Сообщение Александр Мурый » 27 дек 2015, 13:19

ErnieBoyd писал(а):Просьба к модераторам перенести последние посты в отдельную тему.
Модераторы тоже нуждаются в помощи. Скажите конкретнее, с какого именно поста переносить и как лучше назвать новую тему. Вам виднее — вы же здесь пишете.
Редактор материалов, модератор форума

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2227
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1500
Откуда: Казань

Re: Ответ на вопрос «как»

Сообщение Игорь Белов » 27 дек 2015, 13:59

Александр Мурый писал(а):Скажите конкретнее, с какого именно поста переносить и как лучше назвать новую тему.
Предлагаю создать тему "Карты участков в план-схеме на карте Росреестра в MapInfo" и перенести туда все сообщения, начиная с viewtopic.php?f=17&t=8038&start=60#p138245.
The purpose of computing is insight, not numbers

Ответить

Вернуться в «MapInfo»

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

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