[GRASS/Советы] Пакетный импорт растров

[KolesovDmitry]

В каталоге ~/tmp/rastr/ хранятся растры GeoTiff, нужно экспортировать все эти растры в текущий набор GRASS.
Названия файлов содержат суфикс "_longlat.tif". Пример названия файла: "A2007091_longlat.tif", этот файл нужно импортировать под названием "A2007091"

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

GRASS> for file in ~/tmp/rastr/*tif; do r.in.gdal input=$file output=$(basename $file '_longlat.tif'); done

В результате каждый канал растра будет импортирован отдельно: например, двух-канальный растр "A2007091" даст в итоге два растра в GRASS - "A2007091.1" и "A2007091.2"

[Максим Дубинин]

не лишний ли пробел между $file и 'longlat' в этой части скрипта?

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

output=$(basename $file '_longlat.tif')

[KolesovDmitry]

не лишний ли пробел между $file и 'longlat' в этой части скрипта?

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

output=$(basename $file '_longlat.tif')

С одной стороны, пробел не лишний, но, с другой стороны, эту часть скрипта можно и модифицировать под конкретные нужды.

Команда basename (стандартная unix-совская команда, входит в пакет coreutils) выводит последний компонент пути к файлу. Может применяться в двух вариантах -- с суфиксом или без.

Параметры:

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

basename NAME [SUFFIX]

Применение поясню на примерах (привожу команду и результат выполнения):

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

basename /tmp/rastr/abcd_longlat.tif
abcd_longlat.tif

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

basename /tmp/rastr/abcd_longlat.tif .tif
abcd_longlat

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

basename /tmp/rastr/abcd_longlat.tif _longlat.tif
abcd

[Максим Дубинин]

спасибо, понял, убирает подстроку, а без пробела наоборот, будет добавлять

хороший рецепт, подтверждаю - работает