Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
r.rescale.eq ?bolotoved писал(а):Можно: http://grass.osgeo.org/grass70/manuals/r.rescale.htmlrhot писал(а):В GRASS никак нельзя преобразовать 16-бит в 8-бит?
Но нужно разбираться как правильно настроить гистограмму.
i.landsat.rgb тут не поможет?rhot писал(а): Хочется добиться такой же насыщенности и контраста как и после атмосферной коррекции
756, Метод Brovey -- странный результат...rhot писал(а):Привожу результаты паншерпенинга. Хочется добиться такой же насыщенности и контраста как и после атмосферной коррекции, пока не выходит.
Распишите подробнее, пож-та, как именно вы делали.Strix писал(а): У меня в той же комбинации каналов без атмосферной коррекции преобразование Бровея в GDAL (в связке gdal_calc.py + gdal_contrast_stretch) дает такое изображение:
, тем не менее, с помощью форума gis-lab и гугла я написал нечто, что дает приемлемый для моих целей результат. Скрипт запускаю в Ubuntu в каталоге с обрезанными по нужному региону каналами сцены Landsat 8 (на полной сцене не пробовал -- боюсь, что не хватит оперативки в моем компьютере). В системе должен быть установлен dans-gdal-scripts и, как писал выше, gdal_calc.py должен быть допилен по рецепту Voltrona. Да и еще, это в общем-то вариант скрипта для одной комбинации каналов; в расширенном его варианте я получаю все комбинации каналов сразу.
Код: Выделить всё
#!/bin/bash
mkdir ./mult
for i in *B5.TIF *B6.TIF *B7.TIF; do gdalwarp -tr 15 15 -r bilinear $i ./mult/$i; done
find . -name "*B8.TIF" -exec cp {} ./mult \;
cd ./mult
gdal_calc.py -A *7.TIF -B *5.TIF -C *6.TIF -D *8.TIF --outfile=7.tif --calc="(D*A)/(A+B+C)"
gdal_calc.py -A *7.TIF -B *5.TIF -C *6.TIF -D *8.TIF --outfile=5.tif --calc="(D*B)/(A+B+C)"
gdal_calc.py -A *7.TIF -B *5.TIF -C *6.TIF -D *8.TIF --outfile=6.tif --calc="(D*C)/(A+B+C)"
mkdir ./8bit; for i in *.tif; do gdal_contrast_stretch -histeq 100 $i ./8bit/$i; done
cd ./8bit
gdal_merge.py -separate 7.tif 5.tif 6.tif -o 756Lt.tifКод: Выделить всё
gdal_calc.py -A *7.TIF -B *5.TIF -C *6.TIF -D *8.TIF --outfile=7.tif --calc="(D*A)/(A+B+C)"По вашей ссылке выше по теме я рецептов не нашёл (ну и блог Voltrona закрыт для простых смертных). В чём заключается эта доделка? Выложите, пож-та, допиленную версию <gdal_calc.py>.Strix писал(а):gdal_calc.py должен быть допилен по рецепту Voltrona.
Потому как в QGIS теперь мощные инструменты настройки гистограммы, rhot - хочет увидеть такую же картинку в GRASS, какую Strix автоматом получает в QGIS.Strix писал(а):756, Метод Brovey -- странный результат...rhot писал(а):Привожу результаты паншерпенинга. Хочется добиться такой же насыщенности и контраста как и после атмосферной коррекции, пока не выходит.
У меня в той же комбинации каналов без атмосферной коррекции преобразование Бровея в GDAL (в связке gdal_calc.py + gdal_contrast_stretch) дает такое изображение:
Я гистограмму в QGIS не настраивал, т.к. у меня изображение, полученное с помощью данного приема, выглядит одинаково где угодно: в любом графическом вьюере, редакторе etc.bolotoved писал(а):Потому как в QGIS теперь мощные инструменты настройки гистограммы, rhot - хочет увидеть такую же картинку в GRASS, какую Strix автоматом получает в QGIS.
Скрипт я писал изначально для обработки каналов одной сцены Landsat 8, находящихся в одном каталоге. Ситуацию с несколькими сценами в одном каталоге я не предусматривал.KolesovDmitry писал(а):Немного смущает строка такого вида:у меня есть подозрение, что она будет работать, только если в каталоге лежит только один файл, оканчивающийся на ...7.TIF. Вы не проверяли, что случится, если таких файлов будет несколько?Код: Выделить всё
gdal_calc.py -A *7.TIF -B *5.TIF -C *6.TIF -D *8.TIF --outfile=7.tif --calc="(D*A)/(A+B+C)"
Там по ссылке нужно было пролистать вниз до сообщения "Подводные камни gdal_calc" (прошу прощения что не уточнил). Привожу сам рецепт:Александр Мурый писал(а):По вашей ссылке выше по теме я рецептов не нашёл (ну и блог Voltrona закрыт для простых смертных). В чём заключается эта доделка? Выложите, пож-та, допиленную версию <gdal_calc.py>.Strix писал(а):gdal_calc.py должен быть допилен по рецепту Voltrona.
Код: Выделить всё
myval=BandReadAsArray(myFiles[i].GetRasterBand(myBands[i]),
xoff=myX, yoff=myY,
win_xsize=nXValid, win_ysize=nYValid)Код: Выделить всё
myval=BandReadAsArray(myFiles[i].GetRasterBand(myBands[i]),
xoff=myX, yoff=myY,
win_xsize=nXValid, win_ysize=nYValid).astype(numpy.float32)Код: Выделить всё
r.info -r map=band1
r.recode input=band1 output=band1_recoded rules=- << EOF
0.:1.65006412008:0:255
EOF
r.colors map=band1_recoded color=grey