из WKT в proj4

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

Спасибо, буду пробывать...
просто есть много определений МСК в формате PRJ файлов ArcGIS, хотел перегнать в формат proj.4...,

Тогда однозначно надо напустить на них gdalsrsinfo в цикле на каком-либо скриптовом языке.

[Denis Rykov]

Вот пример такого скрипта:

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

import os, fnmatch



def main():

  root_path = os.path.abspath("C:/personal/gislab")

  for root, dirs, files in os.walk(root_path):

    prjfiles = fnmatch.filter(files, '*.prj')

    if (len(prjfiles) != 0):

      for filename in prjfiles:

        os.system('gdalsrsinfo -o "%s" "%s" > "%s.proj4"' % ('proj4', os.path.join(root,filename), os.path.join(root, os.path.splitext(filename)[0])))



  return 0



if __name__ == '__main__':

  main()

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

На всякий случай вариант в линуксовом шелле:

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

find . -type f -name "*.prj" | while read PRJ; do gdalsrsinfo -o proj4 $PRJ > $(dirname $PRJ)/$(basename $PRJ .prj).proj4; done

[Донецков]

Пытался запустить размещенный выше скрипт, увы, пока ничего не выходит...
пытался 2 способами (из под установленного NextGIS, так и установленного в чистую с помощью osgeo4w)
в скрипте путь отредактировал на каталог где разместил несколько файлов prj...
но похоже "Питон" не загружается..., может что-то не так делаю?

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

Пытался запустить размещенный выше скрипт, увы, пока ничего не выходит...
пытался 2 способами (из под установленного NextGIS, так и установленного в чистую с помощью osgeo4w)

К сожалению счастью, винды под рукой нет, чтоб проверить... Но, судя по скриншотам из архива, вы запускаете скрипты (wkt2proj и proj2wkt) без аргументов (типа "+proj=utm +zone=38 +ellps=WGS84 +units=m +no_defs"), а надо бы с ними. И можно без предваряющего "python" -- система с установленным питоном сама поймёт, что это питоновский скрипт.

[Denis Rykov]

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

[Донецков]

На другой (рабочей) машине получилось (W7x64)....
Запускал OSGeo4W.bat из установленного NextGIS 1.9.0, который сделал в виде Portable...
столкнулся с тем, что файлы prj должны быть в названиях без пробелов (иначе появлялась ошибка), переименовал файлы в 1, 2, 3, 4 и т.д и получил такие же с расширением proj.4 ...
Но почему на первой машине не запускается, т.к. NextGIS 1.9.0 Portable один и тот же и размещен на одном и том же диске (D)...

[Denis Rykov]

Обновил текст скрипта, теперь должны корректно обрабатываться имена файлов с пробелами.

[Denis Rykov]

2Донецков
Получилось разобраться почему не работает на одной машине?

[Донецков]

Нет, с помощью OSGeo4W.bat никак не получается...
а из консоли "питона" с помощью измененного скрипта все работает ...

[Донецков]

Взможно, у меня не запускается из-за того, что ранее был установлен Python25, как наследие от ArcGIS 9.3.1...
при копировании скрипта в этот каталог и его запуска proj4 создается, но пустой (размер 0).
Думаю, не снести ли этот Python25..., только не повлияетля это на ArcGIS?

[Донецков]

Интересная программа:
p4Vis
из proj.4 в WKT точно сохраняет. кроме этого еще геокалькулятор на основе proj.4