Страница 1 из 2
Поиск сдвига координат
Добавлено: 10 июн 2015, 13:15
NBAH1990
Добрый день, уважаемые знатоки! Отчаявшись найти решение самостоятельно зарегался здесь. Есть у меня целая куча координат в неведомой системе координат, но так же я знаю к какому адресу относится каждая пара координат, основываясь на этом написал скрипт, который по адресу находит координаты на Яндекс карте, таким образом для некоторых пар моих координат есть соответствующая пара координат в системе Яндекса(например, 306,0 ; -200,9 и 71,425363;51,167798). Теперь к вопросу, как мне найти закономерность этих смещений? Как вывести формулу для всех остальных точек, чтобы я мог все свои корявые координаты преобразовать в реальные (яндексовские, например)? поставил себе триалку ArcGis, пытаюсь найти решение в ArcMap, но пока безуспешно. Надеюсь на вашу поддержку, спасибо.
Re: Поиск сдвига координат
Добавлено: 10 июн 2015, 13:29
trir
берём учебник по картографии...
или можете выложить пример данных и расказать откуда они взялись
Re: Поиск сдвига координат
Добавлено: 10 июн 2015, 13:43
Донецков
306,0 ; -200,9
- что-то слишком условное..., точно не градусы (не бывает таких), вряд ли метры (футы и т.п.), возможно км (мили), либо это просто безразмерные единицы...
Re: Поиск сдвига координат
Добавлено: 10 июн 2015, 14:30
NBAH1990
trir писал(а):можете выложить пример данных и расказать откуда они взялись
Это точки покрытия города одной организацией, от неё и взялись. Заведомо известно что координаты все верные, но они каким-то образом зашифрованы, и если привести их к человеческому виду, то можно спроецировать их на нормальную карту.
306,0 -200,9 ; 71,425363 51,167798
2176,0 1941,5 ; 71,454486 51,187973
-1313,7 -430,6 ; 71,401674 51,165697
-1853,5 4079,6 ; 71,396024 51,207631
359,0 931,8 ; 71,426234 51,178392
356,9 952,8 ; 71,426234 51,178392
-286,5 986,8 ; 71,417565 51,178352
Re: Поиск сдвига координат
Добавлено: 10 июн 2015, 14:33
trir
Re: Поиск сдвига координат
Добавлено: 10 июн 2015, 18:46
Игорь Белов
Прямоугольные координаты и географические, что выдаёт Яндекс по адресам, не относятся к одним и тем же точкам. Расхождения довольно велики с точки зрения картографа. Возможно, Вас вполне устроит определение положения с точностью до дома. Но и в этом случае для надёжности крайне желательно использовать все имеющиеся пары. Вы выложили всё?
Re: Поиск сдвига координат
Добавлено: 11 июн 2015, 06:06
NBAH1990
ErnieBoyd писал(а):Вы выложили всё?
Всего более 200 000 точек, и они довольно часто меняются, я скинул несколько пар точек, по которым я нашел соответствие с реальными координатами. Точности до дома более чем достаточно
Re: Поиск сдвига координат
Добавлено: 11 июн 2015, 06:09
NBAH1990
Читаю по ссылке, там для вычисления реальных координат задаётся 4 параметра, а у меня их только 2 последних. Или в первые 2 можно подставить координаты Яндекса?
Re: Поиск сдвига координат
Добавлено: 11 июн 2015, 07:28
trir
Их надо перевести в СК42, а заодно проверить растояние между точками, оно долно быть, хотя бы, пропорционально
Re: Поиск сдвига координат
Добавлено: 11 июн 2015, 09:40
NBAH1990
trir писал(а):Их надо перевести в СК42, а заодно проверить растояние между точками, оно долно быть, хотя бы, пропорционально
http://clip2net.com/s/3j7q1Xx вот смог вывести в arcmap адреса по моим исходным данным(те что в непонятной СК);
http://c2n.me/3j7qkwO вот вывел эти же адреса, вручную найденные в координатах Яндекса.
Не уверен что тут можно говорить о пропорциональности, но сходство в отображении определенно есть, на текущем этапе большая точность не принципиальна.
"Их надо перевести в СК42", не понимаю, как это нужно сделать? Или если есть возможность как-то перевести в СК42, почему не перевести сразу в координаты Яндекса?
Re: Поиск сдвига координат
Добавлено: 11 июн 2015, 10:37
Игорь Белов
Переводить в СК 42 ничего не надо. Первое действие по статье — переход к географическим координатам, а тут они уже в готовом виде.
Re: Поиск сдвига координат
Добавлено: 11 июн 2015, 12:50
Ariki
В грубом приближении параметры вашей местной системы координат такие:
Код: Выделить всё
+proj=tmerc +lat_0=0 +lon_0=70 +k=1 +x_0=-99500 +y_0=-5672000 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,-0.35,-0.79,-0.22 +units=m +no_defs
Подбирал по семи точкам, приведённым вами. Остаточная погрешность на исходных точках до 104 метров. Вы можете использовать эти параметры для пересчёта координат в WGS84, например, с помощью QGIS или с помощью
онлайн-версии cs2cs. В последнем случае в левом окошке вводите параметры вашей СК, в правом выбираете WGS84.
Re: Поиск сдвига координат
Добавлено: 12 июн 2015, 07:59
NBAH1990
Ariki писал(а):В грубом приближении параметры вашей местной системы координат такие:
Код: Выделить всё
+proj=tmerc +lat_0=0 +lon_0=70 +k=1 +x_0=-99500 +y_0=-5672000 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,-0.35,-0.79,-0.22 +units=m +no_defs
Подбирал по семи точкам, приведённым вами. Остаточная погрешность на исходных точках до 104 метров. Вы можете использовать эти параметры для пересчёта координат в WGS84, например, с помощью QGIS или с помощью
онлайн-версии cs2cs. В последнем случае в левом окошке вводите параметры вашей СК, в правом выбираете WGS84.
Спасибо большое, это именно то что нужно!
Если можно, расскажите пожалуйста, как Вы вычислили эти параметры для преобразования? Чтобы я смог их просчитать по всем точкам, какие смогу найти, для большей точности. К тому же эти точки (7 точек) из одного населенного пункта, и похоже что в каждом населенном пункте эти параметры будут разные.
Re: Поиск сдвига координат
Добавлено: 12 июн 2015, 08:26
NBAH1990
Насколько я понял меняться будет только эта часть:
Код: Выделить всё
+lat_0=0 +lon_0=70 +k=1 +x_0=-99500 +y_0=-5672000
Re: Поиск сдвига координат
Добавлено: 12 июн 2015, 15:00
NBAH1990
ErnieBoyd писал(а):Переводить в СК 42 ничего не надо. Первое действие по статье — переход к географическим координатам, а тут они уже в готовом виде.
Пытаюсь идти по шагам по ссылке, подставляя свои данные. Застрял на шаге:
Код: Выделить всё
proj -f "%.4f" +proj=tmerc +lat_0=71.425363 +lon_0=51.167798 +k=1 +x_0=306.0 +y_0=-200.9 +ellps=krass pt_longlat.dat > pt_tmerc.dat
Возвращается следующая ошибка:
Код: Выделить всё
Rel. 4.9.1, 04 March 2015
<proj>:
Sys errno: 2: <system mess. texts unavail.>
pt_longlat.dat
содержимое файла pt_longlat.dat :
Код: Выделить всё
71.425363 51.167798
71.419560 51.167939
71.454486 51.187973
71.401674 51.165697
71.396024 51.207631
71.426234 51.178392
71.426234 51.178392
71.417565 51.178352
71.422003 51.185353
71.422003 51.185353
71.422003 51.185353
71.436062 51.166894