Пересчет коорд. из местной 3-градусн.сист.коорд. в СК42

[127.0.0.1]

Подскажите пожалуйста,
возможно ли пересчитать координаты местной системы координат с параметрами в "MAPINFOW.PRJ":

"МСК-А", 8, 1001, 7, 66, 0.999999999, 1, 850000, -8900000,
Здесь:
8 - проекция Гаусса-Крюгера
1001 - система координат Пулково 1942, эллипсоид Крассовского
7 - единица измерения - метры
66 – нулевая долгота
0.999999999 – ???
1 – коэффициент сжатия
850000, -8900000 – восточное/северное смещение
и других подобных МСК:
"МСК-Б", 8, 1001, 7, 20.55, 0.99999999999, 1, 850000, 8900000,
"МСК-В", 8, 1001, 7, 23.55, 0.9999999999, 1, 950000, 8900000,

в СК42 (зона 9)
"GK Зона 9 (Пулково 1942)", 8, 1001, 7, 101, 0, 1, 9500000, 0
(значение параметров аналогично)

MapInfo это делает при выполнении "Файл/ Сохранить копию...", с выбором соответствующей проекции. Но мне нужен независимый от MapInfo инструмент.

Достаточно ли указанных параметров для пересчета координат из, например, "МСК-А" в "СК42 (зона 9)" ?

Если можно пример, или ссылку на методический материал. Спасибо.

[a.temushkin]

0.133333333333333 - это начальная широта проекции.

то есть в сторонних конверторах параметры МСК-А будут такими:
Название=МСК-А
Тип проекции=8=Универсальная Поперечная Меркатора
Эллипсоид=1001=Крассовский(в зависимости от софта для пересчёта может потребоваться указание его параметров - a; b;1/f а также датума)
Единицы измерения проекции=7=метры
Центральный меридиан(нулевая широта проекции)=50.55
Начальная параллель(Нулевая долгота проекции)=0.133333333333333
Масштабный коэффициент=1
Восточное смещение=2250000
Северное смещение=-5900000

указав их можно будет пересчитать координаты точек из мскА в СК42 (зона 9), аналогичным образом указав её параметры.
как-то так.
ps - мск знакомая)

[127.0.0.1]

спасибо.
да, разумеется, геодезические калькуляторы такой пересчет делают...
но мне самому! нужно написать хранимую процедуру в БД для пересчета коорд. из МСК в СК42!
мне бы алгоритм, методику или пример пересчета...

[a.temushkin]

а что за СУБД?

формулы:
по буржуйски с их обозначениями.

по русски:
Конформное преобразование координат из местной системы в государственную:
X = (X0+x') + Q1',
Y = (Y0+y') + Q2' - Q3',
где x' = cos(a)/m*dx - sin(a)/m*dy,
y' = cos(a)/m*dy + sin(a)/m*dx,
dx = x - x0,
dy = y - y0,
Q1'=x'*Y0*(2y'+Y0)*f,
Q2'=y'*Y02*f,
Q3'=Y0*(x'+y') *(x'-y').
Здесь:
x0, y0, X0,Y0- координаты начального пункта соответственно в местной и государственной системах координат;
m - масштабный коэффициент,(отношение длин линий в местной системе к линиям в государственной системе);
(a) - угол разворота местной системы относительно государственной;
X,Y и x, y - преобразуемые координаты соответственно в государственной и местной системах;
f=1/2R0, где R0 - радиус кривизны эллипсоида в точке X0, Y0.
Поправка за высоту поверхности гск относительно местной системе координат должна быть учтена в масштабном коэффициенте m.
За начальный пункт обычно принимают начало координат местной системы
---------------------------------------------------------
Афинное преобразование (иногда достачтно его.)
производится общим по формулам афинного преобразования:
x' = x2 + a1*dX + b1*dY,
y' = y2 + a2*dX + b2*dY,
где
dX = x - x1,
dY = y - y1,

[127.0.0.1]

[quote a.temushkin]а что за СУБД?[/quote]
Oracle10gR2

спасибо,
а откуда переменная "Y02" здесь: Q2'=y'*Y02*f, - наверно опечатка?

[a.temushkin]

Y02 - читать как "Y0 в квадрате".

[Boris]

Может вы знаете и ответ на вопрос как по известным точкам в обоих системах координат узнать все переменные параметры в МСК? А то на мой такой вопрос, что-то никто раньше не дал ответа...

[gstof]

В руководстве написано, чтобы работать в координатах плана небходимо выполнить оператор CoordSys Nonearth.
Но при компиляции дал давать ошибку на слудующую за этим выражением строку. Хотя при операторе CoordSys earth компилировалось без ошибок.

[Boris]

В руководстве написано, чтобы работать в координатах плана небходимо выполнить оператор CoordSys Nonearth.
Но при компиляции дал давать ошибку на слудующую за этим выражением строку. Хотя при операторе CoordSys earth компилировалось без ошибок.

Давайте начнем новую тему и в разделе про mapinfo.

[a.temushkin]

Может вы знаете и ответ на вопрос как по известным точкам в обоих системах координат узнать все переменные параметры в МСК? А то на мой такой вопрос, что-то никто раньше не дал ответа...

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

[Boris]

если не трудно дайте здесь список этих программ. как не выйдет за пару часов подобрать - обращусь к специализированному софту. а то в прок его изучать - дело нудное

[gstof]

http://autojp1.narod.ru