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：如何从光栅图像中获取像素值？

帮助？：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），输出级别为0100.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

应该起作用。（最后一个没有测试）