Как преобразовать координату 4175503.20 в 37.510882?

Пожалуйста:
4035419.20 7233396.30 36.252500 54.514300
4035424.70 7233396.30 36.252500 54.514400
4036036.00 7233401.90 36.258004 54.514445
4036076.20 7233395.60 36.258365 54.514412
4037358.50 7238289.80 36.269885 54.539983

Дома посмотрю.. на работе автокад не стоит.

 

Да, это похоже на Меркатора. В первом приближении выполняется формула
Y=A(Y1)*180/pi*ln(tan((Y1+90)*pi/360), но коэффициент A(Y1) для разных широт разный (для Москвы - 110796.5, для Калуги - 110786.8). Наверное, бралась хорошая модель геоида.

Правда, разного смещения по долготе это не объясняет.

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

 

Оффтоп

мужики, а вам не кажется, что вас нагло эксплуатируют в темную? задачка-то денег стОит, судя по количеству и дислокации точек - хоть и не шибко сложная и, к тому же, довольно интересная. а рекламщик вполне способен заплатить - на коньяк к очному заседанию заочного клуба хватило бы... а то и на билеты к месту заседания...

 

А хрен его знает, всю голову сломал... Тык эти координаты откудато ведь взялись? А спросить\прочитать, что это за система координат нельзя? А то чую я мы тут бьемся над названием, а ответ уже известен...

 

"Дело было вечером, делать было нечего". Решил тоже покопаться.
Анализирую то, что было дано.
1 4175503.20 7484666.00 37.510882 55.806533
2 4175497.30 7484682.50 37.510800 55.806600
3 4176233.60 7484330.30 37.517400 55.804800
4 4177683.80 7483672.70 37.530469 55.801508
5 4180080.00 7481841.70 37.552000 55.792200
6 4178257.70 7483006.50 37.535600 55.798100
В 3 и 4 колонках геодезические координаты, не известно в какой системе. При этом, в колонке 3 дана долгота, а в колонке 4 широта (иначе это координаты не в Москве).
Ищу максимальное значение B°. Это строка 2 (точка2). Минимальное значение B° в строке 5 (точка 5). То же самое проделываю в колонке с L°.
Далее в колонках 1 и 2 аналогично нахожу максимальные и минимальные значения прямоугольных координат. Сопоставляю результаты в колонках 1, 2 и 3, 4. Максимальное B° во второй строке и ему должно соответствовать максимальное X. Такое максимальное X находится тоже во второй строке, во второй колонке. Следовательно, в колонке 2 значения X, а в колонке 1 значения Y. Проверяю соответствие L° и Y.
Итак разобрались:
______Y__________X_________L°________B°___
1 4175503.20 7484666.00 37.510882 55.806533
2 4175497.30 7484682.50 37.510800 55.806600
3 4176233.60 7484330.30 37.517400 55.804800
4 4177683.80 7483672.70 37.530469 55.801508
5 4180080.00 7481841.70 37.552000 55.792200
6 4178257.70 7483006.50 37.535600 55.798100
Не зная в какой СК даны геодезические координаты, преобразую их в прямоугольные (СК-42, зона 7) в двух вариантах – как L и B из WGS-84 и из СК-42.
В Транскоре для обоих вариантов прямоугольных координат СК-42 находится ключ пересчета в неизвестную МСК. Оба варианта, несмотря на плохую совместимость точек, дают решение осевого меридиана близкого к московской системе координат (той, что на элл. Бесселя) с осевым меридианом 37° 30'. Однако масштабный коэффициент пересчёта получается примерно 1.8, чего на самом деле быть не должно.
Можно предположить, что прямоугольные координаты были сняты с плана в Автокаде, где план вписали в какой-то формат, не задумываясь о масштабе и о величине сдвижки координат. Этим объясняется, почему Y и X поменялись местами, и почему у них такие "дикие" первые цифры.

Не вижу смысла далее работать с подобными материалами.

 

Давайте вернёмся к нашим вашим баранам.
Х - это У в неизвестной проекции основанной на СК-42, первая цифра (4) -номер зоны, предполагаю, что зоны 10-ти градусные
У - это Х в неизвестной проекции основанной на СК-42
Х1 - долгота WGS-84 (в долях градуса)
Y1 - широта WGS-84 (в долях градуса)
Делаю вывод, что программно можно переводить координаты точек, находящихся в одной зоне.

 

Коэффициент примерно 1.779? Если да, то это 1/cos(55.8 гр), как и должно быть в меркаторской проекции на широте Москвы. То, что мы видим в X и Y - это метры. Но в масштабе, отнесенном к экватору. Они взяли меркаторскую карту всей Земли, и ввели на ней единые метрические координаты (примерно так же, как сделал Google Maps во внутреннем представлении своих карт). Осталось понять, есть ли хотя бы одна общепринятая СК с такой структурой, а если нет - то угадать, какой моделью пользовались авторы системы.

 

эх:))) конспирологи и адепты теории заговора подтянулись.
нет, никого эксплуатировать не собирался. думал, что все просто будет.
а вон оно как вышло.
(Добавление)

В любом случае огромное спасибо. Если бы знал, что это выльется в геморрой, не стал бы никого напрягать.
(Добавление)

К сожалению инфу неоткуда раздобыть:(

 

А поменяться местами Х и Y могли по невнимательности товарищей которые при черчении в Автокаде по запаре и детской математической школьной привычке применили вместо системы плоских прямоугольных координат Гаусса -Крюгера перевернутую систему Декарта,не обратив внимание на расположение осей координат

 

А я бы еще поковырялся. Не можете кинуть известные пары из разных городов? С Москвой и Калугой все более-менее понятно, но глобальная картина пока скрыта :(

Если судить по (X,Y), то между красноярскими точками около 60 метров, а если по (X1,Y1) - то около 5 метров (если я не ошибся с подсчетом цифр). Бывает, конечно...

 

Оффтоп

Да какой тут заговор? Обычная отечественная любовь к халяве... Вместо того, чтобы обратиться с конкретной задачей к специалистам - а попробуем-ка сами... Иногда получается... Но практика показывает, что получается всё реже и реже. Чесслово, буду за вас искренне рад, если результат будет удовлетворительный. Только жаль потерянного вами времени. Держите в курсе.

ЗЫ. погуглите на тему "траблшутер"...

 

Вот ряд Уфимских координат:

X Y X1 Y1
6211856.50 7264501.60 55.802900 54.676400
6225471.60 7258650.10 55.925200 54.645900
6213290.00 7264461.50 55.815800 54.676200
6234133.00 7271887.80 56.003100 54.714800
6221538.40 7284595.00 55.889900 54.780800
6233159.20 7277866.80 55.994300 54.745900
6233074.50 7268564.40 55.993500 54.697500
6239422.90 7281817.70 56.050600 54.766400

А вот для Питера:


X Y X1 Y1
3391756.10 8345510.10 30.470876 59.922383
3391773.00 8344791.00 30.471000 59.919100
3392351.10 8343563.00 30.476220 59.913601
3393006.00 8342157.00 30.482103 59.907258
3393007.90 8342158.90 30.482120 59.907267
3393746.00 8340717.10 30.488750 59.900761
3393684.90 8340655.70 30.488200 59.900500
3388603.00 8335617.10 30.442500 59.877700
(Добавление)

:) Да дело в том, что при разработке ПО подобных задач достаточно много всплывает.
Обычно, разработчик думает, что все ответы в этой жизни можно нагуглить:)) Но не всегда это срабатывает.
Ну и плюс ко всему, большинство людей не вникают в нюансы профессий.
Я вот диву даюсь, когда мне предложили написать программу для call-центра за 50 тысяч рублей в свободное от работы время. Или просят сделать интернет-магазин типа Озона за пять тысяч рублей.

 

Еще, если можно, Иркутск, Хабаровск и Омск и/или Тюмень. Пока все получается довольно красиво (не считая Ярославля, в котором ошибка 4 км)

 

Интересно:)

Хабаровские координаты:
X Y X1 Y1
15035085.80 6166937.90 135.066475 48.496112
15035956.40 6157836.80 135.067400 48.496300
15034792.70 6166992.30 135.065214 48.496966
15035002.30 6169721.60 135.066606 48.513544
15033972.10 6163023.30 135.049500 48.527300
15035287.00 6156184.10 135.061400 48.486500

Иркутск:

X Y X1 Y1
11609031.10 6816787.80 104.284100 52.277800
11604136.20 6814300.50 104.240100 52.264100
11616210.80 6817331.30 104.348600 52.280800
11611690.40 6817973.20 104.308000 52.284300
11611822.90 6816694.80 104.309200 52.277300
11599251.30 6813624.60 104.196300 52.260300
11608520.90 6818768.70 104.279517 52.288689

Омск:

X Y X1 Y1
8168532.20 7313642.40 73.379210 54.931190
8168338.30 7326898.60 73.377500 54.999700
8158186.50 7332543.80 73.286300 55.028800
8169185.60 7316235.40 73.385100 54.944600
8157936.10 7325260.30 73.284000 54.991200
8172304.90 7324170.50 73.412644 54.985604
8169310.30 7319324.50 73.386200 54.960570
8161805.80 7323601.20 73.318800 54.982700

Ижевск:

X Y X1 Y1
5931227.60 7690971.00 53.282200 56.836000
5925299.20 7701372.40 53.229000 56.887200
5931001.10 7692140.30 53.280200 56.841800
5925157.70 7690528.70 53.227700 56.833800
5926113.30 7690651.80 53.236300 56.834400
5925638.00 7691397.60 53.232000 56.838100
5934464.70 7691979.10 53.311300 56.841000
5924668.30 7695110.60 53.223300 56.856400

 

И еще - Астрахань/Краснодар, Волгоград и Киров (для охвата по широтам).

Точности 20 метров хватит, или нужно выжимать лучше?
(Добавление)
Три хабаровских точки попали с точностью 10 м, а еще три - улетели на 6 км :)

 

Похоже всё-таки, что координаты были даны в равноугольной цилиндрической проекции Меркатора.
Создал программку в Excel для преобразования B° и L° в X и Y и пытался подобрать такие параметры эллипсоида (или сферы), чтоб пересчёт соответствовал примерам на Москву, Калугу, Ярославль, Красноярск, Владивосток. Задавал эллипсоиды и WGS, и Красовского, и произвольные. Единых параметров эллипсоида для всех регионов подобрать не удалось. Как были вычислены, и насколько точны данные не известно. Какой-то бессистемный разброд получается. Программку в Excelе прилагаю. Формулы защищены, можно менять только данные широты, долготы, размеры эллипсоида и радиус сферы.
Кто желает, может поупражняться.

 

Ellipsoidal Mercator projection


У меня получились такие формулы:

Y2=R*ln(tan(45+Y1/2)*((1-e*sin(Y1))/(1+e*sin(Y1)))^(e/2)),
X2=X1*R*pi/180

где R=6378443, e=0.082033

Y3=Y2+240*sin(X1)-4*cos(X1)-277
X3=X2-3000*(1-cos(Y1-52))-385

У всех тригонометрических функций аргументы считаются в градусах :)

Тогда (X3,Y3) будет приближением для (X,Y) с точностью до 10 метров для всех приведенных точек, за исключением 7 (для которых ошибки гораздо больше).

xls-файл (правда, без комментариев) прилагается.

 

Не знаю, как по английски.
Но не путать с равноугольной поперечно-цилиндрической - той, что на Гаусса-Крюгера похожа.