Определить грид по имеющимся двум парам координат

[lam]

Во вложении исходные данные, буду признателен, если направите на путь истинный.

Для координат WGS84 четыре знака после запятой думаю будет маловато.

[bingeomap]

Доброго времени суток!
Во вложении исходные данные, буду признателен, если направите на путь истинный.

Я так понимаю, у вас наверно проблема с определением систем координат который в Эксел файле вы написали как "неизвестная СК", не так ли? Если так, то могу вам советовать не использовать "Долгота/Широта(WGS84)_EPSG_4326", это вам не понадобится. Вам надо преобразовывать координаты "GK Зона 15 (Пулково 1942)_EPSG_28415" на геодезическую (Долгата\Широта Пулково 1942) а потом с помощью программу который я предлагал Ахе (Вычисление параметров проекции), определит параметры проекции. Вот и всё.

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

Вам надо преобразовывать координаты "GK Зона 15 (Пулково 1942)_EPSG_28415" на геодезическую (Долгата\Широта Пулково 1942)

GeoCalculator.png (20.14 КБ) 9337 просмотров

потом с помощью программу который я предлагал Ахе (Вычисление параметров проекции), определит параметры проекции.

ComputeDatumParams.png (29.09 КБ) 9337 просмотров

Вот и всё.

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

Долгота осевого меридиана: 84 00 00
Сдвиг вдоль оси X, м: -5480000
Сдвиг вдоль оси Y, м: 20000

[Slinger]

Подтверждаю. Для MapInfo проекция.

"Условная система", 8, 1001, 7, 84, 0, 1, 20000, -5480000

image.JPG (39.38 КБ) 9291 просмотр

[Dob]

Я так понимаю, у вас наверно проблема с определением систем координат который в Эксел файле вы написали как "неизвестная СК", не так ли? Если так, то могу вам советовать не использовать "Долгота/Широта(WGS84)_EPSG_4326", это вам не понадобится. Вам надо преобразовывать координаты "GK Зона 15 (Пулково 1942)_EPSG_28415" на геодезическую (Долгата\Широта Пулково 1942) а потом с помощью программу который я предлагал Ахе (Вычисление параметров проекции), определит параметры проекции. Вот и всё.

Спасибо, за оперативность!
поразбираюсь!
P.S. кстати, в Долгата\Широта Пулково 1942 тоже переводил - программа постоянно выдавала ошибку (видимо, сам косячил ), а в WGS 84 выложил, потому что gamm просил именно в них.
еще раз спасибо!

[stout]

Подтверждаю

[Axe]

Спасибо всем большое! Большое дело сделано.
Мы сейчас попробуем применить найденные параметры смещения для координат всех имеющихся устьев в пределах 15 зоны.
И, главное, будем дальше разбираться с параметрами условных прямоугольных координат, которые использовались по другим рудникам.

[gamm]

в WGS 84 выложил, потому что gamm просил именно в них.

ну, раз просил, тоже отпишусь - хотя использовал ГК15. Файл прилагается, результаты те же (запись проекции в формате стрки PROJ4). Да, X - смотрит на Восток, Y - смотрит на Север

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

Исходное состояние после грубой подгонки
Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 84.0000000000 1.0000000000 19993.6162142626 -5480002.7681545150 0.020517

Начальное приближение, эвристический метод
Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 83.9998293978 1.0000000000 19981.3266250160 -5480002.8205153728 0.020517

Квази-ньютоновский метод, среднеквадратическая ошибка
 ID	    dX	    dY	  Dist
 -----------------------------
  1	 -0.01	 -0.02	  0.02
  2	 -0.00	 -0.02	  0.02
  3	  0.01	  0.01	  0.02
  4	 -0.01	  0.03	  0.03

Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 83.9998293978 1.0000000000 19981.3266250160 -5480002.8205153728 0.020517
+proj=tmerc +lat_0=0.000000 +lon_0=83.99982939778 +k=1.000000 +x_0=19981.326625 +y_0=-5480002.820515 +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,1.1 +units=m +no_defs

Квази-ньютоновский метод, абсолютная ошибка, без весов
 ID	    dX	    dY	  Dist
 -----------------------------
  1	  0.00	 -0.02	  0.02
  2	  0.00	 -0.02	  0.02
  3	  0.02	  0.01	  0.02
  4	 -0.00	  0.02	  0.02

Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 83.9998294503 1.0000000000 19981.3243402268 -5480002.8170133363 0.021676
+proj=tmerc +lat_0=0.000000 +lon_0=83.99982945027 +k=1.000000 +x_0=19981.324340 +y_0=-5480002.817013 +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,1.1 +units=m +no_defs

Квази-ньютоновский метод, средняя ошибка, без весов
 ID	    dX	    dY	  Dist
 -----------------------------
  1	 -0.00	 -0.01	  0.01
  2	  0.00	 -0.00	  0.00
  3	  0.01	  0.02	  0.03
  4	 -0.01	  0.04	  0.04

Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 83.9998291942 1.0000000000 19981.3104021836 -5480002.8349128393 0.025077
+proj=tmerc +lat_0=0.000000 +lon_0=83.99982919423 +k=1.000000 +x_0=19981.310402 +y_0=-5480002.834913 +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,1.1 +units=m +no_defs

Квази-ньютоновский метод, среднеквадратическая ошибка, с весами
 ID	    dX	    dY	  Dist
 -----------------------------
  1	 -0.00	 -0.00	  0.01
  2	  0.00	 -0.00	  0.00
  3	  0.01	  0.02	  0.03
  4	 -0.00	  0.04	  0.04

Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 83.9998294663 1.0000000000 19981.3297116726 -5480002.8357249955 0.025619
+proj=tmerc +lat_0=0.000000 +lon_0=83.99982946633 +k=1.000000 +x_0=19981.329712 +y_0=-5480002.835725 +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,1.1 +units=m +no_defs

[stout]

Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 83.9998294663 1.0000000000 19981.3297116726 -5480002.8357249955 0.025619

Долгота осевого меридиана и смещение Efalse имеют коэффициент корреляции очень близкий к единице:
(84°-83.9998294663°)×(π/180°)*6378245 м = 18.984 м
20000 м - 19981.3297 м = 18.670 м
(Величина большой полуоси эллипсоида большой роли не играет)
Находить их совместно из одной системы уравнений очень сложно.
Любая ошибка в долготе осевого меридиана легко компенсируется изменением в Efalse. И наоборот.
На практике долготу осевого меридиана находят из формул, в которых нет ни Efalse, ни Nfalse.
Это может быть формула для расстояния между точками или, например, сближения меридианов.
В программке, которая на скриншоте выше, я использовал приращения координат. Так меньше сказываются ошибки в координатах исходных пунктов на результат.

[gamm]

На практике долготу осевого меридиана находят из формул, в которых нет ни Efalse, ни Nfalse.

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