Параметры региональных МСК для MapInfo

[dab]

Андрей, спасибо.

[Донецков]

Пытаюсь составить описание МСК на одно из муниципальных образований Сахалинской области, земельный кадастр в нем ведется в "своей" местной СК. Хотел сделать строку описания для пересчета данных из СК-95 в эту МСК, естественно, без аффинных преобразований на этой локальной территории не обойдешься...
Вычислили параметры трансформации:

Спойлер

┌─────────────────────────────────┐
│ Коэффициенты преобразования: │
╞═════════════════════════════════╡
│ X'=Mx*cos(U)*X-My*sin(U)*Y+Dx │
│ Y'=Mx*sin(U)*X+My*cos(U)*Y+Dy │
├─────────────────────────────────┤
│ U = 2 21'27.2" │
│ Mx= 0.99975167 │
│ My= 0.99975167 │
│ Dx=-ххххххх.ххх │
│ Dy=-yyyyyy.yyyy │
└─────────────────────────────────┘

после соответствующих вычислений и с учетом того, что МСК, по косвенным данным, разрабатывалась когда-то на базе СК-42, получил

Спойлер

"0000000, Сахалинская обл. Affine (Пулково 1995 ГОСТ 51794-2008)", 1008, 9999, 3, 24.47, -130.89, -81.56, 0, 0, -0.13, -0.22, 0, 7, 141, 0, 1, 24500000, 0, 0.998905457427707, -0.041125281627371, -24yyyyyy.yyyy, 0.041125281627371, 0.998905457427707, -xxxxxxx.xxx

и получается совсем не то, что ожидалось...
В чем может быть ошибка, т.к. ранее никогда не писал параметры с аффинными преобразованиями...

[Boris]

поскольку спойлеры пропали (да и что бы они дали), по опыту в mapinfo, могу сказать, что дело в направлении поворота, со всем остальным трудно ошибиться.
самое простое лечение - заставить мапинфо вычислить эти параметры самой. я делаю это так:
1) совмещение точки 0 - выполняется за пределами мапинфо, на пример в автокаде. в нем получаем смещение и угол поворота.
2) создаем точку в 0 локальной СК и выводим ее на экран, центруем экран по этой (точке Ctrl+A - для этого хорошо подходит, если изначально точки не было на экране, можно и руками набить цетр в "показать по другому" окна)
3) находим в "программах" программу "поворот окна" и в ней говорим повернуть окно и указываем вычисленные градусы. в окне мапбасика мапинфо напишет параметры трансформации, вернее новую проекцию с аффинным преобразованием, которую она рассчитала для того, что бы окно повернуть. вот параметры поворота без смещения надо взять именно оттуда. по моему это 1,2,4,5 параметры.

[Донецков]

Спасибо, но спойлеры на месте, проверил сейчас, все открывается...

дело в направлении поворота

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

[dab]

На первый взгляд, читая документацию:

Спецификатор аффинных единиц и коэффициенты следуют после основных параметров для системы...

попробуйте добавить 7 перед шестью параметрами аффинного

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

"0000000, Сахалинская обл. Affine (Пулково 1995 ГОСТ 51794-2008)", 1008, 9999, 3, 24.47, -130.89, -81.56, 0, 0, -0.13, -0.22, 0, 7, 141, 0, 1, 24500000, 0, 7, 0.998905457427707, -0.041125281627371, -24yyyyyy.yyyy, 0.041125281627371, 0.998905457427707, -xxxxxxx.xxx

[Донецков]

Спасибо, результат на несколько порядков лучше, но все равно далеко от идеала..., остальное в ЛС...

[juffin_h]

Позволю себе высказать сомнения.

но все равно далеко от идеала.

А на сколько далеко? И откуда взяты данные для вычисления параметров?
Если из параметров МСК, то аффинное преобразование не годится. Должно быть выполнено перепроецирование.

[Донецков]

Данные каталожные...
Непосредственно параметров МСК не знаю...

[juffin_h]

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

[Донецков]

Исходными было 5 точек с известными координатами в одной и другой системе, которые располагаются вокруг интересующего участка, параметры вычислялись в нескольких программах, можно в том же R, значения получаются схожие, но не оптимизированы, по моему для использования в MI, т.к. в вычислениях использованы данные по Y без номера зоны (24), то добавление 24 значительно ухудшает результат.
Возможно, стоит из Y убрать 500000 и потом выполнить расчет параметров...

[Institor]

МСК, по косвенным данным, разрабатывалась когда-то на базе СК-42, получил

Может, косвенные данные врут.
Если не врут - проверьте датум

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

3, 24.47, -130.89, -81.56, 0, 0, -0.13, -0.22, 0,

Это точно СК-42?
Данные в известной системе - это в СК-42?

[Донецков]

24.47, -130.89, -81.56, 0, 0, -0.13, -0.22

- это параметры СК95 согласно ГОСТ от 2008...
Данные в известной СК - это данные в СК95...

косвенные данные врут

- это офф. полученные выписки, по ним работали, все в пределах допусков...
Суть в том, что даже если найти параметры МСК созданной на базе СК42, это не решит проблемы точного пересчета с СК95 в МСК, т.к. нет точных параметров перехода СК95 - СК42, а если их вычислять, то и получим опять масштабный К и угол поворота, т.к. простым смещением это описать не возможно...
Вот и задался целью написать параметры с аффинными коэфф. описывающими МСК, но с основой на СК95, чтобы минимизировать погрешности перехода СК95-WGS84-СК42 (или как производная - МСК, что равнозначно), т.к. MI как и большинство др. программ производит все вычисления через WGS84, а это не есть хорошо по точности из-за неоднозначности пересчета с Красовского на WGS84 и обратно...
Пока можно судить только о том, что для получения параметров МСК предварительно нужно данные в СК95 как-то еще подготовить, возможно вычесть из Y 500 000, а может и еще что-то....

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

...для получения параметров МСК предварительно нужно данные в СК95 как-то еще подготовить, возможно вычесть из Y 500 000, а может и еще что-то...

Афинное преобразование сохраняет формы и размеры при переносе координат, так что это бесполезно.

...MI ... производит все вычисления через WGS84, а это не есть хорошо по точности из-за неоднозначности пересчета с Красовского на WGS84 и обратно...

Пересчёт вполне однозначен. Неоднозначность в параметрах, которые Вам неизвестны.

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

Если МСК построена на СК-42/63/95 без разворота и масштабирования, афинное преобразование - собаке пятая нога.
Считайте координаты СК-95 точными, а для МСК подберите параметры проекции, и будет Вам счастье. Проекцию назовите по имени участка, поскольку для другого участка километрах в тридцати, обеспеченного другим набором исходных пунктов, параметры будут другие. Так учтём деформации старой сети. Сим победиши.
Ежели захочется терять время и силы на параметры конформного преобразования, вспомните о тех же деформациях сети и бросьте.
Другое дело, если имеешь дело с МСК населённого пункта с имманентными разворотом и масштабом. Но это совсем другая тема.

[juffin_h]

Если Вы хотите увязать 4 системы координат, то, как мне кажется надо разделить задачи.
Согласен с ErnieBoyd. Если МСК разрабатывалась на базе СК42, то это значит, что датум у них один. Здесь важно не то, что он (датум) СК42, а то, что он один для обеих систем. Следовательно, датум оставляем стандартный (для СК42), а определяем параметры проекции, как это делает ComputeDatumParams (РАКУРС). И в последнюю очередь (если есть не преодолимое желание ) используем аффинные преобразования.
Связь между СК42 и СК95 математически не описывается. Можно получить усредненные параметры. Но если упоминается WGS84, то логично все увязать через нее.
Итого:
МСК - определяются параметры проекции на датуме СК42.
СК42 –> WGS84 и СК95 –> WGS84 - определяются 7 параметров на основе тех же опорных точек с известными координатами для каждой СК. Параметры по ГОСТ не годятся.
После получения приемлемых результатов по параметрам СК42 – WGS84, их можно использовать в качестве датума для МСК.

[gimran]

Из опыта — ряд предприятий нарушают принятые ими же МСК, например смещение вместо 100000 бывает 99999, иногда забывают угол смещения произвести и тп.