Програмка/скрипт для пересчета координат

[doomhammer]

Здраствуйте,
помогите пож-ста разобраться новичку в проблеме, а точнее понять шаги и что делать.
Дан векторный слой в shp формате,
Нужно взять какой-то базовый слой и наложить данный векторный слой на него.
Мои действия,
нашел програмку ogr2ogr для перевода shp -> geojson, ниже скрипт который использовал.
ogr2ogr --config SHAPE_ENCODING "UTF-8" -f geoJSON output.geojson input.shp
далее при помощи библиотеки leafletjs слепил базовый слой openstreetmap и сформированный geojson слой.
Слой соответственно не спроецировался куда надо, а лег куда-то в антактику. Ну а как без этого.
И тут выяснилось что координаты в данном векторном слое не в стандартной мировой системе координат которую понимает leaflet, а в какой-то непонятной местной системе координат.

Далее нашел в интернетах что можно в скрипт огра добавить примерно вот такое.
-s_srs EPSG:26917 -t_srs EPSG:4326
но как узнать какие именно цифры подходят?
как я понял .prj файл может в этом как то помоч.
И вообще в правильном ли я направлении. Возможно какие действия вам покажутся странными, если что-то не так направьте кто в этом понимает, просто карты для меня пока в новинку.

[trir]

а в какой-то непонятной местной системе координат

параметры этого непонятного есть?

http://gis-lab.info/qa/local-cs.html

В данном случае, EPSG вам не поможет

[doomhammer]

спасибо за разъяснение,
на счет параметров есть только данные в .prj файле

GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Local"],PARAMETER["False_Easting",7096.2717],PARAMETER["False_Northing",8623.8408],PARAMETER["Scale_Factor",1.00011880188265],PARAMETER["Azimuth",0.0791562083333333],PARAMETER["Longitude_Of_Center",77.0407772555556],PARAMETER["Latitude_Of_Center",43.3363979194444],UNIT["Meter",1.0]

Как я понял там взято за центр точка указанная в долготе широте, и от нее указывают координаты в метрах.
а как их использовать, или какой проге скормить, чтоб координаты пересчитать или как нибуть leafletу показать.
Есть какие нить идеи?

[doomhammer]

невошедший кусок

Спойлер

.

GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],
PROJECTION["Local"],
PARAMETER["False_Easting",7096.2717],
PARAMETER["False_Northing",8623.8408],
PARAMETER["Scale_Factor",1.00011880188265],
PARAMETER["Azimuth",0.0791562083333333],
PARAMETER["Longitude_Of_Center",77.0407772555556],
PARAMETER["Latitude_Of_Center",43.3363979194444],UNIT["Meter",1.0]

[trir]

PROJECTION["Local"],

Серъёзно? Откуда такие параметры? Это где то работает?

Читать:
http://gis-lab.info/qa/proj-sk-faq.html#6
https://ru.wikipedia.org/wiki/%D0%94%D0 ... 1%83%D0%BC
https://ru.wikipedia.org/wiki/%D0%9A%D0 ... 0%B8%D1%8F
http://gis-lab.info/qa/mapinfo_to_wkt_proj4.html

Нашёл, это только ArcGIS поддерживает

Вариантов 2:
1. В ArcGIS сохранить в EPSG:4326
2. Расчитать параметры для Gauss-Kruger

[Донецков]

Пересчитывать в ArcGIS в нужную СК...

Не секрет на какую территорию такая МСК?

[trir]

Алматы