Перевести карту для мапинфо из плансхемы в систему координат

[Boris]

Пришлите побольше точек, но только в одной зоне. Я на досуге запишу видео или набор картинок.

[donmisha]

Boris, подобрал 22 пары точек. Все они находятся в пределах 8 зоны проекции Гаусса-Крюгера
1166890 385880 8576310 6200520
1177170 384580 8586600 6199570
1164160 364590 8574850 6179190
1178340 362270 8588650 6177290
1165300 376620 8575090 6191170
1172070 387310 8581410 6202150
1169160 398080 8578250 6212690
1180130 386290 8589460 6201400
1180250 381150 8589790 6196290
1165510 393860 8574550 6208300
1166310 356380 8576860 6170990
1179710 337110 8590940 6152260
1182050 321290 8593920 6136480
1218610 405800 8627210 6222330
1221700 399910 8630540 6216490
1229530 425600 8637490 6242450
1242290 409130 8650740 6226530
1214530 401950 8623330 6218280
1200790 398710 8609880 6214440
1194230 397960 8601630 6212870
1201900 406280 8610540 6222120
1218430 397100 8627360 6213640

[gamm]

1) уточните, это план-схема (откуда она взялась - территория большая) или МСК?
2) некоторые точки явно ошибочные (номер 20 по крайней мере), ошибки для аффинного преобразования:

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

[,1][,2]
1     39.920263   86.9033635
2      7.926676   48.4286570
3    464.309701  100.5764278
4      5.522866  -16.4819236
5     41.753809   35.3456865
6     21.329086   95.9241333
7    195.139244   -0.3135888
8    -21.625124   61.3401333
9    -15.060119   75.0820607
10   -25.342568  -43.0917837
11     1.175044   11.1703703
12   -69.557337    4.5879305
13   -53.875223  -79.3514807
14    56.488565   82.9671591
15    66.094167    3.6181624
16   210.652729   38.7599302
17    60.210673   72.5827151
18   100.262054   27.1603984
19   249.277062  -65.0001256
20 -1476.477058 -640.7152362
21   100.075978   20.6916147
22    41.799513   79.8153958

после удаления 3 и 20 точек:

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

         [,1]        [,2]
1   -7.365171   56.341277
2  -37.937176   20.486890
4    9.520364  -23.994285
5   15.465452   12.972403
6  -29.915879   64.788413
7  120.164822  -41.639848
8  -71.753844   32.259537
9  -53.679554   50.687253
10 -90.323053  -81.095839
11  20.129200    7.323777
12  -9.336274   20.121813
13  41.483012  -49.119640
14 -42.983707   41.552530
15 -20.618461  -32.013535
16  65.191175  -19.111222
17 -50.170127   31.459008
18  10.016456  -11.328332
19 168.293828 -102.469691
21   1.955330  -23.498831
22 -38.136393   46.278320

[Boris]

Я не смотрел данные, но если есть большие погрешности в одной точке, то это не правильно, т.к. весь алгоритм состоит в постепенном приближении начала координат глобальной СК к началу координат локальной СК путем расчета параметров переноса через аффинное преобразование. По тексту статьи, должно быть ортогональное смещение-поворот-масштабирование, но т.к. в реализации этого алгоритма меня посетили сомнения в правильности выписанной подпрограммы, то он был заменен на аффинное преобразование с усреднением масштабирования по осям.

[Boris]

Да, меня тут посетили сомнения - вопросы. А как записаны координаты? В моем XLS они пишутся как в мапинфо - X-восток, Y-север.
И первый столбец, если не путаю координаты в глобальной СК.

[Донецков]

Если не ошибаюсь:
Умск 1166890 Хмск 385880
Уск42 8576310 Хск42 6200520

и зачем Вам дополнительно подбирать параметры, если МСК-21 региона для Чувашии построена на основе СК-63 блок D и имеет вид для MapInfo:

MSK-21_D3 зона 1", 8, 1001, 7, ЦМ, смещение от экватора, 1, 1(значение У по ЦМ), -5800000

ЦМ, смещение от экватора, значение У по ЦМ - аналогичны параметрам СК-63 (в открытую не пишу т.к. "С"), если у кого из знакомых есть Кадастровый офис, там эти параметры были приведены....

[donmisha]

1) уточните, это план-схема (откуда она взялась - территория большая) или МСК?[/code]

Данные были взяты из Республиканского отдела картографии. У них территория всей республики оцифрована в Мапинфо в одной системе координат (размер территории Чувашии 150 на 200 км). Когда я смотрю проекцию таблиц Мапинфо - вывод "план-схема, метры".

[donmisha]

ххх

2) некоторые точки явно ошибочные (номер 20 по крайней мере)[/code]

Да, проверил 20 точку - она не верная.
Должно быть так
1192610 397330 8601630 6212870

[donmisha]

Да, меня тут посетили сомнения - вопросы. А как записаны координаты? В моем XLS они пишутся как в мапинфо - X-восток, Y-север.

Да, кординаты, записаны таким образом: 1 столбец - долгота (Х) в план-схеме, 2 столбец - широта (У) в план-схеме, 3 столбец - долгота (Х) в Гаусса-Крюгера, 4 столбец - широта (У) в Гауса-Крюгера.

И первый столбец, если не путаю координаты в глобальной СК.

Когда я сохраняю таблицу в mif файле, то при открытии в блокноте выводится:
Version 450
Charset "WindowsCyrillic"
Delimiter ","
CoordSys NonEarth Units "m" Bounds (1100000, 200000) (1400000, 500000).

Я правильно понимаю, что параметр "Bounds (1100000, 200000) (1400000, 500000)" описывает границы карты?

[donmisha]

Если не ошибаюсь:
Умск 1166890 Хмск 385880
Уск42 8576310 Хск42 6200520

Нет, наоборот:
Хмск 1166890 Умск 385880
Хск42 8576310 Уск42 6200520,
при этом Х - долгота, У-широта

и зачем Вам дополнительно подбирать параметры, если МСК-21 региона для Чувашии построена на основе СК-63 блок D и имеет вид для MapInfo:

MSK-21_D3 зона 1", 8, 1001, 7, ЦМ, смещение от экватора, 1, 1(значение У по ЦМ), -5800000

ЦМ, смещение от экватора, значение У по ЦМ - аналогичны параметрам СК-63 (в открытую не пишу т.к. "С"), если у кого из знакомых есть Кадастровый офис, там эти параметры были приведены....

Кадастрового офиса пока в шаговой доступности нет.

[gamm]

видимо, тов.Донецков прав. Найдите это "С" полишинеля, добавьте эту "С" проекцию в Мапинфо, загрузите свои точки в UTM и переведите в нее. Потом посчитайте свои dX, dY - поиск в гугле и на данном форме вам поможет.

[Донецков]

Естественно, я приводил X и Y геодезические, аналогичные получаемым выпискам из местной кадастровой палаты...

[Boris]

Я попробовал посчитать - ошибка схождения порядка 3-х километров. Либо
а) точки неправильные или взятые с большой долей фантазии
б) Чувашия не 8-ой зоне СК-42/95
в) все таки попутан порядок точек в столбцах - точки должны быть в стиле мапинфо : X-восток, Y-север. Расчет предназначен для мапинфо - не стоит впутывать сюда российский ГОСТ.
г) локальная СК с масштабным коэффициентом отличным от 1-цы?

[Донецков]

По указанным параметрам результаты следующие:
Уск42 8576310 Хск42 6200520
вычисленные:
Умск 1166885.24 Хмск 385901.94
т.е. расхождение несколько десятков метров, но не 3 км...

[donmisha]

По указанным параметрам результаты следующие:
Уск42 8576310 Хск42 6200520
вычисленные:
Умск 1166885.24 Хмск 385901.94
т.е. расхождение несколько десятков метров, но не 3 км...

Что показывают вычисленные значения? Куда их подставлять?