image-processing - Gdal:如何有条件地为光栅图像的像素分配一个新值？

Question

给定地形 GIS 栅格 一国crop.tif :

# download:
curl -o ETOPO1.zip 'http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO1/data/ice_surface/grid_registered/georeferenced_tiff/ETOPO1_Ice_g_geotiff.zip'
# unzip:
unzip ETOPO1.zip
# crop:
gdal_translate -projwin -005.50 051.30 10.00 041.00 ETOPO1_Ice_g_geotiff.tif crop.tif

给定海拔阈值 n = 50(米)

如何将 z >= n 的所有像素的值设置为 50 ？

相关: https://stackoverflow.com/questions/18300210/ , Gdal: How to get a pixel's value from raster image?

帮助？: gdallocationinfo 可以得到一个点的值(无需迭代整个图像):

$ gdallocationinfo crop.tif 1 1 -valonly
> 73

Answer 1

使用 gdal_calc.py .举例:

gdal_calc.py -A crop.tif --outfile=level0100.tif --calc="100*(A>100)"     --NoDataValue=0

读作:`给定输入-A(crop.tif)，输出level0100.tif，当像素具有第一个波段的值A>100时，将值设置为100。否则，设置为0。

内 --calc="" ，您可以设置其他条件。如果要将 100 < v < 200 等值设置为 100，则:

gdal_calc.py -A crop.tif --outfile=level0100.tif --calc="100*(100<A<200)"     --NoDataValue=0

应该管用。 (最后一个没有测试)