Форумы GIS-Lab.info

Системная ошибка 2 — это, вроде бы, файл не найден. Убедитесь, что файл лежит в вашем текущем рабочем каталоге, или укажите полный путь к нему.

Мне подбирать параметры систем координат приходится часто, поэтому я автоматизировал процесс, написав скриптик на языке Python. Он использует метод наименьших квадратов и годится не только для проекции Гаусса—Крюгера, но и для многих других. Некоторые параметры можно зафиксировать, если они известны, и подбирать оставшиеся.

Не могу судить, насколько скрипт удобен для начинающих, так как писал изначально для себя, а до того подбирал методами, аналогичными тем, которые изложил ErnieBoyd в своей статье. Так что рекомендую пройти его руководство до конца, чтобы прочувствовать процесс и приобрести некоторый опыт. Для использования моего скрипта необходимо понимание параметров проекции в PROJ.4 и умение устанавливать питоновские модули (под Windows это не всегда тривиально).

НЕ АКТУАЛЬНО
Всем спасибо, вроде справился с задачей, вот моё решение: Теперь возникла новая проблема, при помощи сервиса http://cs2cs.mygeodata.eu/ всё работает, координаты преобразуются, однако при попытке что-то сделать в своей командной строке возвращается вот такой результат: Из командной строки proj открываю вот так: пары координат просто через пробел ввожу Что делаю не так?

НЕ АКТУАЛЬНО
Понял в чем ошибка, надо было так вызывать: Теперь другой вопрос) как их вывести в таком формате 71.1234 51.1234 , а не 71d25'5.215"E 51d10'20.475"N 0.000 ?

man cs2cs
-f format

Всё, разобрался, Извиняюсь за флуд.
Вот окончательное решение для меня: Если есть замечания по нему - пожалуйста скажите. Спасибо большое всем кто содействовал решению задачи.

А почему вы решили, что у вас проекция Ламберта (lcc)? Это маловероятно на просторах СНГ, тем более что ваши данные довольно хорошо ложатся на поперечную проекцию Меркатора (tmerc) с круглыми значениями параметров.

Ariki писал(а):А почему вы решили, что у вас проекция Ламберта (lcc)? Это маловероятно на просторах СНГ, тем более что ваши данные довольно хорошо ложатся на поперечную проекцию Меркатора (tmerc) с круглыми значениями параметров.

На самом деле я это не решил. Просто методом проб и ошибок нащупал наиболее приемлемое решение) Попробовал tmerc, с высоты моего понимания выглядит не менее хорошо , но если вы считаете что так будет правильнее, то наверное оставлю tmerc.

Просто местные системы координат, как правило, построены на основе системы координат Гаусса—Крюгера (она же поперечная проекция Меркатора), и смещения по осям в них часто представляют собой круглые числа в метрах (либо смещение на целое число минут плюс круглое число в метрах). То же касается долготы осевого меридиана: там обычно целое число минут или даже градусов. Бывают, конечно, исключения, когда система координат строится от произвольной точки с произвольной ориентацией осей, но это реже.

В данном случае подбор долготы осевого меридиана сразу дал число, близкое к 70 градусам: Я его смело округлил до этой величины, а затем заново определил смещения, которые оказались кратными полукилометру плюс-минус несколько десятков метров: Эти величины тоже просятся, чтобы их округлили.

Параметры +towgs84, характеризующие взаимное положение эллипсоидов (а также часто используемые для компенсации деформаций геодезической сети), здесь взяты из российского ГОСТа и могут отличаться от реальных на данную территорию, но это может дать ошибку в худшем случае в пару десятков метров.

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