Arcgis修改栅格数据像素值

Arcgis修改栅格数据像素值

遇到对DEM部分像素的高程值进行修改，记录一下过程。

方法一：栅格计算器

若要对一个范围内的像元值进行统一的修改，可以：

1. 新建shp
2. Feature to raster ，注意设置environment：设置extent与输入DEM一致，开启捕捉，以及像素大小
3. Raster Calcultor 键入：Con(IsNull(a),b,c)。a:第二步生成的 b：dem c:统一修改的高程值。
   网上有具体的。

方法二：利用Raster To ASCII逐一修改像元值

1. 利用Raster To ASCII将DEM转化为.txt格式
2. 在txt文件里根据行列号修改像元值
3. 利用ASCII To Raster将.txt格式转回DEM（.tif）格式
4. 生成的DEM没有空间参考，设置空间参考：Data Management Tools/Projections and Transformation/Raster/Define Projection 给DEM设置空间参考。

问题

方法一，适合将部分像素统一修改为同一值，法二可以逐像素修改。但法二需要知道待修改像元的行列号，本篇遇到修改地方DEM范围小，进行裁剪后才进行了该步，行数列数只有几十行*几十行，数了下修改的，确定栅格数据的行列号本篇没有研究，网上提供的利用Raster Calcultor输入r=$$rowmap只适用于arcgis9.x版本。10.x版本Raster Calculator支持用Python语句找到几篇解决该问题的解答，还未尝试成功，附在下面：
https://community.esri.com/t5/arcgis-spatial-analyst-questions/ncols-rowmap/m-p/272290#236639
https://gis.stackexchange.com/questions/138754/getting-latitude-and-longitude-raster-from-dem-using-arcgis-desktop?answertab=votes#tab-top