Форумы GIS-Lab.info

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

#!/bin/bash
# Файл создан TraCE, с сайта gis-lab.info, особая благодарность потльзователям
# IRC-канала gis-lab за своевременную помощь.
# Предназначен для трансформации и обрезки файлов из одной проекции в другую.
# Для правильной работы скрипта необходим gdal-dev из пакета OSGeo4W и шейп
# сохраненый в нужной проекции. Шейп-файл взят с сайта gis-lab.info. Растры взяты
# с сайта poehali.org (для других растров нужна трансформация имен).
# Запрашиваемый в начале путь, указывается в конце без "/"
# Например /c/temp, а не /c/temp/. иначе выдаст ошибку.
# Запускается с помощью программы Msys и команды sh /путь/до/файла/map2tif.sh
echo Укажите входные параметры для начала работы скрипта
echo Например /c/temp
echo ВНИМАНИЕ!!! Если есть данные в папке с конечным результатом
echo и данные имеют одинаковые название, то все данные будут заменены
echo без подтверждения!
read -p "Введите путь до папки с файлами gif и map -> " DIR_GIF
read -p "Введите путь до папки в которую надо сохранить конечный результат ->" DIR_TIF
read -p "Введите путь до папки с шейп-файлом ->" DIR_SHP
read -p "Введите название шейп-файлас расширением shp (Имя_файла.shp) ->" NAME_SHP
echo Введиде код проекции. Например 28418 - 18 Зона, 28417 - 17 зона
read -p "EPSG: -->" EPSG_NUM
echo Трансформация файлов из gif в geotif
cd $DIR_GIF
    for a in $( ls *.gif); do 
		/c/OSGeo4W/apps/gdal-dev/bin/gdal_translate -of GTiff $a ${a%%.gif}".tif" > null
		/c/OSGeo4W/apps/gdal-dev/bin/gdalwarp -t_srs "EPSG:"$EPSG_NUM -of GTiff ${a%%.gif}".tif" $DIR_TIF"/"${a%%.gif}".tif" > null
		echo $a, готов
	done
echo Удаление файлов с расширением .tif из папки gif
cd $DIR_GIF
rm *.tif

echo Обрезка файлов по вектору и дополнительному полю
cd $DIR_TIF
    for b in $( ls *.tif ); do
	     c="$(echo "${b%%.tif}" | sed 's/200k--//' | tr '[[:lower:]]' '[[:upper:]]' | sed 's/_/,/' | sed 's/-/-0/' | sed 's/,/,0/')"
		 gdal_rasterize -i -b 1 -burn 0 -l ${NAME_SHP%%.shp} $DIR_SHP"/"$NAME_SHP -where 'indexfull="'$c'"' $b > null
		 echo "$c", готов
    done
echo Трансформация закончена.
exit 0