Проекция в километрах

[RMSea]

Вопрос, может, глупый, но до меня решение не доходит.

Проблема такая: мое все исходное в UTM координатах, зона 48 (окрестности Байкала). Появились старые данные, которые сделаны, как я подозреваю, в Surfer'е.. Там есть и гриды, и шейп-файлы (мне нужны только последние пока). Все в координатах Гаусса-Крюгера, НО.. в километрах. Казалось бы, переконвертнуть из ГК в UTM проще простого, однако, ничего не выходит.

Пыталась выложить все длинное описание того, что я предпринимала, но по-моему, это все никак не поможет.

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

[trir]

Может так:

+proj=tmerc +lat_0=0 +lon_0=105 +k=1 +x_0=18500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=km +no_defs

[dab]

В... Все в координатах Гаусса-Крюгера, НО.. в километрах. Казалось бы, переконвертнуть из ГК в UTM проще простого, однако, ничего не выходит. ..

Согласен с trir
Если у Вас не получится так решить проблему, предлагаю пример данных в километрах опубликовать.
Для проверки идей на практике

[RMSea]

Если не ошибаюсь, то я так уже пробовала.. Сейчас еще раз попробую на всякий случай.

[ericsson]

Наверняка, дело в какой-нибудь ошибке или чем-то еще таком.
Например, в самих данных используется проекция в километрах, а в prj-файлах (или что там еще есть) - метры. Или наоборот. Либо false easting задан с номером зоны. Либо еще что-то в этом духе. Вот потому пример данных и был бы полезен.

[RMSea]

Сейчас восстановила все данные, как были изначально, как они мне достались. Там в принципе не было файлов .prj. Я вручную указывала проекцию.

У меня в NextGIS-QGIS 1.8.0 не позволяет изменить в командной строке метры на километры, я меняла ручками через блокнот. В итоге файл так и остается с дурацкими координатами, никаким образом не соотносится с моими данными.

[Максим Дубинин]

в какой еще командной строке?

то что вы в блокноте меняете описание системы координат никоим образом не меняет сами координаты
http://gis-lab.info/qa/proj-sk-faq.html#22

[RMSea]

Вот, например, данные. Как были изначально, без .prj

[RMSea]

Максим Дубинин, неправильно выразилась и про командную строку и про координаты. Я понимаю, о чем вы пишете. Но мне надо, чтобы отображалось у меня в метрах, даже если и гаусс-крюгер, то пусть это будут координаты 185ххххх, а не 185хх.ххх, как это на данный момент выглядит.

Командной строкой я назвала ту самую строку:

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

+proj=tmerc +lat_0=0 +lon_0=105 +k=1 +x_0=18500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs

я не могу в ней сразу изменить метры на километры, не дает ее редактировать. Говорю, может, не в том месте это делаю?

[Максим Дубинин]

Как перепроецировать данные:

http://gis-lab.info/qa/ogr2ogr-examples ... 0.B8.D0.B5

Сделайте так:

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

ogr2ogr -s_srs "+proj=tmerc +lat_0=0 +lon_0=105 +k=1 +x_0=18500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=km +no_defs" -t_srs "+proj=tmerc +lat_0=0 +lon_0=105 +k=1 +x_0=18500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs" -f "ESRI Shapefile" ves2.shp ves1.shp

[dab]

Сделайте пользовательскую систему координат

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

+proj=tmerc +lat_0=0 +lon_0=105 +k=1 +x_0=18500000 +y_0=0 +ellps=krass towgs84=24.47,-130.89,-81.56,-0,-0,0.13,-0.22 +units=km +no_defs

Укажите её для SHP файла
В свойствах проекта включите преобразования координат "на лету" и укажите стандартную Pulkovo 1942 / Gauss-Kruger zone 18, то есть с параметрами

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

+proj=tmerc +lat_0=0 +lon_0=105 +k=1 +x_0=18500000 +y_0=0 +ellps=krass towgs84=24.47,-130.89,-81.56,-0,-0,0.13,-0.22 +units=m +no_defs