Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
Для координат WGS84 четыре знака после запятой думаю будет маловато.Dob писал(а):Во вложении исходные данные, буду признателен, если направите на путь истинный.
Я так понимаю, у вас наверно проблема с определением систем координат который в Эксел файле вы написали как "неизвестная СК", не так ли? Если так, то могу вам советовать не использовать "Долгота/Широта(WGS84)_EPSG_4326", это вам не понадобится. Вам надо преобразовывать координаты "GK Зона 15 (Пулково 1942)_EPSG_28415" на геодезическую (Долгата\Широта Пулково 1942) а потом с помощью программу который я предлагал Ахе (Вычисление параметров проекции), определит параметры проекции. Вот и всё.Dob писал(а): Доброго времени суток!
Во вложении исходные данные, буду признателен, если направите на путь истинный.
bingeomap писал(а):Вам надо преобразовывать координаты "GK Зона 15 (Пулково 1942)_EPSG_28415" на геодезическую (Долгата\Широта Пулково 1942)
bingeomap писал(а):потом с помощью программу который я предлагал Ахе (Вычисление параметров проекции), определит параметры проекции.
Не совсем. На самом деле параметры обычно более круглые, чем получается из решения. В данном случае подойдет такой набор:bingeomap писал(а):Вот и всё.
Спасибо, за оперативность!bingeomap писал(а): Я так понимаю, у вас наверно проблема с определением систем координат который в Эксел файле вы написали как "неизвестная СК", не так ли? Если так, то могу вам советовать не использовать "Долгота/Широта(WGS84)_EPSG_4326", это вам не понадобится. Вам надо преобразовывать координаты "GK Зона 15 (Пулково 1942)_EPSG_28415" на геодезическую (Долгата\Широта Пулково 1942) а потом с помощью программу который я предлагал Ахе (Вычисление параметров проекции), определит параметры проекции. Вот и всё.
), а в WGS 84 выложил, потому что gamm просил именно в них.ну, раз просил, тоже отпишусь - хотя использовал ГК15. Файл прилагается, результаты те же (запись проекции в формате стрки PROJ4). Да, X - смотрит на Восток, Y - смотрит на СеверDob писал(а): в WGS 84 выложил, потому что gamm просил именно в них.
Код: Выделить всё
Исходное состояние после грубой подгонки
Параметры (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
Долгота осевого меридиана и смещение Efalse имеют коэффициент корреляции очень близкий к единице:Параметры (Lat0 Lon0 k dx dy) и СКО:
0.0000000000 83.9998294663 1.0000000000 19981.3297116726 -5480002.8357249955 0.025619
в программке, которая это считает, при определении осевого меридиана не используется ни Efalse, ни Nfalse. Там дисперсия минимизируется сначала (она не зависит от сдвига), если точек много и координаты с некоторой ошибкой, то работает лучше, чем аналитические методы, основанные на точных вычислениях (сближение, и т.д.). А потом, при фиксированных остальных параметрах, определяется сдвиг. С разными весовыми схемами, в том числе устойчивыми к выбросам (грубым ошибкам).stout писал(а):На практике долготу осевого меридиана находят из формул, в которых нет ни Efalse, ни Nfalse.