из WKT в proj4

[Филиппов Владислав]

Доброго всем!
как имея описание СК в формате proj4 (использую QGIS)
получить его WKT? и наоборот в WKT из proj4

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

PROJCS["WGS 84 / UTM zone 35S",
  GEOGCS["WGS 84",
    DATUM["World Geodetic System 1984",
      SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],
      AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
    UNIT["degree", 0.017453292519943295],
    AXIS["Geodetic latitude", NORTH],
    AXIS["Geodetic longitude", EAST],
    AUTHORITY["EPSG","4326"]],
  PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
  PARAMETER["central_meridian", 27.0],
  PARAMETER["latitude_of_origin", 0.0],
  PARAMETER["scale_factor", 0.9996],
  PARAMETER["false_easting", 500000.0],
  PARAMETER["false_northing", 10000000.0],
  UNIT["m", 1.0],
  AXIS["Easting", EAST],
  AXIS["Northing", NORTH],

[rhot]

http://spatialnotes.blogspot.com/2010/1 ... proj4.html

[Филиппов Владислав]

интересно, интересно....
спасибо за ссылку!

[vlalexey]

Подскажите, а кроме того кода по ссылке есть какая-то готовая утилита (под Win), или онлайн сервис для конвертации Proj4<->WKT?

[Slinger]

Если проекция общеизвестная, то:

http://www.spatialreference.org/

Если своя, то лично я не видел

[Филиппов Владислав]

вот именно что "своя"

[Донецков]

Имеется ли перевод

Cartographic Projection Procedures for the UNIX
Environment—A User’s Manual
ftp://ftp.remotesensing.org/proj/OF90-284.pdf

на русский язык?

[Александр Мурый]

Имеется ли перевод

Cartographic Projection Procedures for the UNIX
Environment—A User’s Manual
ftp://ftp.remotesensing.org/proj/OF90-284.pdf

на русский язык?

Нет, конечно. Всё ждёт своих героев

[Донецков]

Возможность перевода WKT <-> proj.4 встречал в онлайновском Coordinate Converter Сканэкса...
но формат WKT программа использует "старую" версию..., т.к. импорт prj из набора ArcGIS 9.3.1 происходит с ошибкой..., формат файла не тот...
Может кто знает более удобную программа по созданию пользовательской СК в формате proj.4?

[Дмитрий Барышников]

OGRSpatialReference позволяет, но это надо программировать.
http://www.gdal.org/ogr/classOGRSpatialReference.html

[Донецков]

Как понимаю на spatialnotes.blogspot.com/2010/11/converting-wkt-projection-info-to-proj4.html приведены 2 макроса, вопрос, как им воспользоваться...

[Дмитрий Барышников]

Так там же написано:
wkt2proj.py [WKT Projection Text]
proj2wkt.py [Proj4 Projection Text]

Чего не понятно?

[Донецков]

да запустить как? никогда с *.py не связывался до этого...
понятно что Питон...
можно ли их запустить из консоли Питона в составе QGIS?
и каким образом? составлять командную строку?

[Александр Мурый]

да запустить как? никогда с *.py не связывался до этого...
можно ли их запустить из консоли Питона в составе QGIS?
и каким образом? составлять командную строку?

wkt2proj_win.py

(295 байт) 1167 скачиваний

proj2wkt_win.py

(297 байт) 1020 скачиваний

Сохраняем скрипты куда-нибудь на диск --> открываем командную строку --> переходим в папку со скриптами --> запускаем скрипты из командной строки с аргументами в виде строк WKT или PROJ4 соответственно.

Пример с UTM зоной 38:

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

wkt2proj.py "PROJCS["UTM_Zone_38_Northern_Hemisphere",GEOGCS["GCS_unnamed ellipse",DATUM["unknown",SPHEROID["Unknown",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",45],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]"

+proj=utm +zone=38 +ellps=WGS84 +units=m +no_defs

И наоборот:

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

proj2wkt.py "+proj=utm +zone=38 +ellps=WGS84 +units=m +no_defs"

PROJCS["UTM Zone 38, Northern Hemisphere",GEOGCS["WGS 84",DATUM["unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",45],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]

-----------------------------------------------------------------------------------
Кстати, с некоторых пор в составе GDAL появилась спец. утилита gdalsrsinfo, которая отлично работает без этих "костыльных" скриптов.

Утилита выводит информацию о СК/проекции из геоданных, файлов *.prj и по номеру EPSG в форматах:

default proj4 and wkt (default option)
all all options available
wkt_all all wkt options available
proj4 PROJ.4 string
wkt OGC WKT format (full)
wkt_simple OGC WKT (simplified)
wkt_noct OGC WKT (without OGC CT params)
wkt_esri ESRI WKT format
mapinfo Mapinfo style CoordSys format
xml XML format (GML based)

[Донецков]

Спасибо, буду пробывать...
просто есть много определений МСК в формате PRJ файлов ArcGIS, хотел перегнать в формат proj.4...,
попробывал СаанЭкс калькулятором, он не обрабатывает новые PRJ...
конечно, как вариант, можно думаю загружать например пустые данные в формате shp поочередности с разными PRJ в QGIS. он должен сохранять их в пользовательские СК и уже в нем можно видеть строку в формате proj.4