Поиск сдвига координат

[NBAH1990]

Добрый день, уважаемые знатоки! Отчаявшись найти решение самостоятельно зарегался здесь. Есть у меня целая куча координат в неведомой системе координат, но так же я знаю к какому адресу относится каждая пара координат, основываясь на этом написал скрипт, который по адресу находит координаты на Яндекс карте, таким образом для некоторых пар моих координат есть соответствующая пара координат в системе Яндекса(например, 306,0 ; -200,9 и 71,425363;51,167798). Теперь к вопросу, как мне найти закономерность этих смещений? Как вывести формулу для всех остальных точек, чтобы я мог все свои корявые координаты преобразовать в реальные (яндексовские, например)? поставил себе триалку ArcGis, пытаюсь найти решение в ArcMap, но пока безуспешно. Надеюсь на вашу поддержку, спасибо.

[trir]

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

[Донецков]

306,0 ; -200,9

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

[NBAH1990]

можете выложить пример данных и расказать откуда они взялись

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

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

[trir]

очень местная

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

Прямоугольные координаты и географические, что выдаёт Яндекс по адресам, не относятся к одним и тем же точкам. Расхождения довольно велики с точки зрения картографа. Возможно, Вас вполне устроит определение положения с точностью до дома. Но и в этом случае для надёжности крайне желательно использовать все имеющиеся пары. Вы выложили всё?

[NBAH1990]

Вы выложили всё?

Всего более 200 000 точек, и они довольно часто меняются, я скинул несколько пар точек, по которым я нашел соответствие с реальными координатами. Точности до дома более чем достаточно

[NBAH1990]

очень местная

Читаю по ссылке, там для вычисления реальных координат задаётся 4 параметра, а у меня их только 2 последних. Или в первые 2 можно подставить координаты Яндекса?

[trir]

Их надо перевести в СК42, а заодно проверить растояние между точками, оно долно быть, хотя бы, пропорционально

[NBAH1990]

Их надо перевести в СК42, а заодно проверить растояние между точками, оно долно быть, хотя бы, пропорционально

http://clip2net.com/s/3j7q1Xx вот смог вывести в arcmap адреса по моим исходным данным(те что в непонятной СК);
http://c2n.me/3j7qkwO вот вывел эти же адреса, вручную найденные в координатах Яндекса.
Не уверен что тут можно говорить о пропорциональности, но сходство в отображении определенно есть, на текущем этапе большая точность не принципиальна.

"Их надо перевести в СК42", не понимаю, как это нужно сделать? Или если есть возможность как-то перевести в СК42, почему не перевести сразу в координаты Яндекса?

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

Переводить в СК 42 ничего не надо. Первое действие по статье — переход к географическим координатам, а тут они уже в готовом виде.

[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.

[NBAH1990]

В грубом приближении параметры вашей местной системы координат такие:

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

+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 точек) из одного населенного пункта, и похоже что в каждом населенном пункте эти параметры будут разные.

[NBAH1990]

Насколько я понял меняться будет только эта часть:

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

+lat_0=0 +lon_0=70 +k=1 +x_0=-99500 +y_0=-5672000

[NBAH1990]

Переводить в СК 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