Система координат МСК в квантуме - кто имел дело?

[Slinger]

В QGIS вроде преобразование Гельмерта используется. Может значения поворотов нужно с обратным знаком указывать ?

[lam]

При такой форме записи все абсолютно не в ту степь улетает. О координатах начальной точки проекции я помню.

выложите координаты одной и той же точки в WGS84 и СК-42

[Игорь Белов]

В QGIS вроде преобразование Гельмерта используется. Может значения поворотов нужно с обратным знаком указывать ?

Верно. Если параметры из Pinnacle в Бурша-Вольфа, так и следует поступить.

При такой форме записи все абсолютно не в ту степь улетает. О координатах начальной точки проекции я помню.

выложите координаты одной и той же точки в WGS84 и СК-42

И полную строку Ваших параметров, пожалуйста.

[Charaunica]

При такой форме записи все абсолютно не в ту степь улетает. О координатах начальной точки проекции я помню.

выложите координаты одной и той же точки в WGS84 и СК-42

При написании параметров использовала данную статью: http://gis-lab.info/qa/local-cs.html#.D ... 1.80.D0.B0
Соответственно, здесь k - масштабный коэффициент;x_0 и y_0 - координаты начальной точки проекции, полученные при нажатии кнопки "Рассчитать" в Установки-Ввод системы координат

+proj=tmerc +lat_0=53.90311667 +lon_0=27.55481667 +k=0.7990658317 +x_0=-93.1395 +y_0=110.1021 +ellps=krass +towgs84=-9.160106617,-136.8964304,-60.6945809,-0,-3.359180803,-3.248368344,-5.900805023 +units=m +no_defs
Используемые точки:
1 5976614,588 546617,874 232,041 53°54'47.71309"N 27°42'27.30552"E 255,441
2 5976290,593 546584,215 229,59 53°54'37.24560"N 27°42'25.28473"E 252,993
3 5977005,786 546806,919 234,465 53°55'00.30386"N 27°42'37.87595"E 257,865
4 5977372,793 546921,291 236,301 53°55'12.13610"N 27°42'44.34383"E 259,701
5 5976919,989 544984,713 232,226 53°54'58.10558"N 27°40'58.00085"E 255,659

моя МСК жестко привязана к СК-42, поэтому и использую параметры перехода туда. У точек должны получиться такие координаты:
1 1134,322 10036,046 232,038
2 810,606 9999,78 229,59
3 1523,989 10228,235 234,465
4 1890,067 10345,558 236,301
5 1452,861 8405,381 232,226

Есть подозрение, что слишком близкие друг к другу точки взяты для пересчета, хоть это и МСК(

[Донецков]

Близость не должна влиять, если есть жесткая связи между СК-42 и МСК...
а вот взаимосвязь с WGS84 все от "лукавого" и имеет локальный характер для каждого участка...

[Игорь Белов]

Вы путаетесь в различных сущностях. «Scale: 0.7990658317 ppm» — это вовсе не масштабный коэффициент проекции k_0, а записанный в специфической форме «хвостик» преобразования датума: ∆S = (S − 1) × 10⁺⁶. Он должен быть седьмым параметром в +towgs84=…

Итак, перед Вами стоят две задачи:

• определить параметры датума;
• определить параметры проекции.

Первую задачу Вы решили с помощью Pinnacle, вторую руководствуясь статьёй. Осталось свести результаты в одном наборе параметров, но, к сожалению, где-то Ваше решение оказалось ошибочным.
Выполним задачи с использованием программ Ракурса.

Вычисление 7 параметров

Подготовим два файла.
В первом файле с координатами в СК-42 к значениям координат Y добавим по 5000000 (пятая зона):

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

1,5976614.588,5546617.874,232.041
2,5976290.593,5546584.215,229.59
3,5977005.786,5546806.919,234.465
4,5977372.793,5546921.291,236.301
5,5976919.989,5544984.713,232.226

Второй файл содержит координаты в WGS 84:

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

1,53 54 47.71309,27 42 27.30552,255.441
2,53 54 37.24560,27 42 25.28473,252.993
3,53 55 00.30386,27 42 37.87595,257.865
4,53 55 12.13610,27 42 44.34383,259.701
5,53 54 58.10558,27 40 58.00085,255.659

Запустим программу ComputeSevenParams. Нажимая каждую из четырёх кнопок в верхнем ряду, укажем для первой и второй координатных систем тип данных и эллипсоид. Для первой выберем, конечно, эллипсоид Красовского.
Загрузим данные в соответствующие окошечки. Нажмём кнопку [Вычислить].

Вычислить 7 параметров

ComputeSevenParams.png (35.38 КБ) 7233 просмотра

Соберём строку параметров датума:

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

+towgs84=-8.6433,-136.8101,-59.7188,3.353,3.251,5.904,0.629

Здесь компоненты разворота меняют знак, а масштаб вычисляется по вышеприведённой формуле.

Вычисление параметров проекции

Придётся пересчитать координаты из пятой зоны СК-42 в широту/долготу СК-42:

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

1,53 54 48.221832,27 42 33.922494
2,53 54 37.754395,27 42 31.901214
3,53 55 00.812251,27 42 44.493620
4,53 55 12.644388,27 42 50.961935
5,53 54 58.618126,27 41 04.619818

Из файла координат МСК удалим высоты:

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

1,1134.322,10036.046
2,810.606,9999.78
3,1523.989,10228.235
4,1890.067,10345.558
5,1452.861,8405.381

Запустим программу ComputeDatumParams. Загрузим данные, вычислим параметры:

Вычислить параметры проекции

ComputeDatumParams.png (31.43 КБ) 7233 просмотра

Сформируем окончательный набор параметров:

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

+proj=tmerc +lat_0=0 +lon_0=27.5705643917 +k_0=1 +x_0=911.2154 +y_0=-5975255.9586 +towgs84=-8.6433,-136.8101,-59.7188,3.353,3.251,5.904,0.629 +ellps=krass +no_defs

Здесь «недостающие» параметры +lat_0=0 и +k_0=1 принимаются программой ComputeDatumParams в качестве постоянных.

Удачи!

[Charaunica]

Итак, перед Вами стоят две задачи:

• определить параметры датума;
• определить параметры проекции.

Первую задачу Вы решили с помощью Pinnacle, вторую руководствуясь статьёй. Осталось свести результаты в одном наборе параметров, но, к сожалению, где-то Ваше решение оказалось ошибочным.

.........

Удачи!

ErnieBoyd, спасибо Вам огромное за подсказку, достойную отдельной статьи! Все получилось, МСК садится нормально, вылетают только максимум 30 см по сравнению с данными съемки