Страница 1 из 1

Преобразовать локальную проекцию (wkt) в формат +proj (решено)

Добавлено: 06 апр 2018, 20:56
romanshuvalov
Добрый день. Переношу из этой ветки конкретный вопрос про координаты.

Дано: Spatial Reference вот в таком виде:

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

PROJCS["Local Conditional",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",27478.878],PARAMETER["False_Northing",19463.17],PARAMETER["Scale_Factor",1.0],PARAMETER["Azimuth",-1.272472],PARAMETER["Longitude_Of_Center",49.3933174],PARAMETER["Latitude_Of_Center",53.4992761],UNIT["Meter",1.0]] 
При попытке использовать его в python-скрипте или в утилите gdalwarp я получаю

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

ERROR 6: No translation for Local to PROJ.4 format is known.
Как преобразовать эту строку в формат "+proj..."?

Цель: преобразовать картинку так, как это делается вот в этой статье.

Глобальная цель: нарезать "OSM-совместимые" тайлы из изображений, которые можно загрузить только в определенной проекции.

Re: Преобразовать локальную проекцию (wkt) в формат +proj

Добавлено: 06 апр 2018, 21:23
romanshuvalov
Помог Павел с telegram-канала:

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

+proj=omerc +lat_0=53.4992761 +lonc=49.3933174 +alpha=-1.272472 +k=1 +x_0=27478.878 +y_0=19463.17 +gamma=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Вроде всё правильно. Вопрос закрыт.