ArcGISArcMap中修改栅格图的像元值.pdf

-在处理图象数据时，我们经常会碰到要求修改栅格图象象元值的问题，比如说 DEM 图的部分数据错误，我们要进行修改；再比如说栅格图象中有些与周围均匀色彩不一致的错误斑块要更正等等， 那我们如何来处理这一类问题呢？现我以一 DEM 栅格图（名字为 eldodem）为例， 现在我要修改它的部分象元值， 总结出以下三种方法，大家可以参考一下。第一种方法不大实用，但可借鉴，第二三种方法针对的条件不一致，大家可以在具体情况下进行选择。一、直接运用转换，思路简一、直接运用转换，思路简单，易操作，但实用性单，易操作，但实用性1、栅格ascii 文件栅格这种方法是先将栅格图用工具直接转成 ascii 文件，然后在 ascii 文本文件中直接修改需要修改的象元的值， 修改好后又用工具-.可修编 .-转换成栅格图。这种方法可行，但是不实用，因为我们要搜索到指定的象元好像不是那么简单。那么有没有别的比较好的方法呢？二、在栅格二、在栅格计算器中操作，方法灵活，可操作性强，实用性强计算器中操作，方法灵活，可操作性强，实用性强准备工作准备工作先要在 option 中设置保留的栅格围，通常情况下默认的为相交后的部分，这里我们要保留整个DEM，所以要改为以下设置：2、 通过点的位置修改点象元值 21* 问题一 、要是我想修改图象中的指定的行列的点的象元值 （比如说把第 100 行，200 列的点的值修改为 0），那该如何操作呢？这时候我们可以在栅格计算器中输入以下公式：1.con($rowmap = 99 & $colmap = 199),0,eldodem.img)(栅格图的编号是从 0 行 0 列开始的)-.可修编 .-22* 问题二、那若是要修改指定行列围的栅格的象元值，比如说把第101 行，251 列到第401 行，301 列的部分的象元值改为0，又改如何进行呢？1.con($rowmap 100 & $colmap 250),0,eldodem.img)结果如下-.可修编 .-23 * 问题三、这些行和列是给定的，但是实际情况往往不是这样的，那若看到一些要修改的点， 我如何通过以上的方法来修改它的属性值呢， 这时要解决的首要问题是如何获知这个点的行列值， 我们可以通过以下方法来解决：（1） 准备工作， 设置 spatial analyst 中的 option，设置所需围，象素以及捕捉-.可修编 .-这个确定后再次打开设定捕捉，使生成的栅格和原栅格象元重合-.可修编 .-（2）在栅格计算器中分别做出对应的行列栅格图行栅格图的生成：在栅格计算器中输入以下公式1.r=$rowmap1.c=$colmap（在这里 r-行，c-列）（3）然后可以用 identify 工具， 对要修改的点的行列号进行查询-.可修编 .-如图，所要修改的点位于第 166 行，129 列，若要修改此点象元值为6，可以在栅格计算器中输入以下公式：1.con($rowmap = 165 & $colmap = 128),6,eldodem.img)3、通过点面文件来修改象元值*问题、现在我们有一幅栅格图，由于某种需要要修改其中的部分点单元格的象元值或者是修改部分斑块的象元值，这时候该如何实现呢？本人觉得以下方法可行。思路主要是（1）先建一点 shp 或面 shp 图层，（2）然后在点或面图层上数字化需要修改的点或面， （3）接着将点或面转换成栅格图，（4）在栅格计算器中进行替换运算。现我仍以eldodem 为例，假设其中有部分象元的值是错误的，要对其进行修改， （同时进行点面转换操作）具体步骤如下：（1） 新建一点文件point_11,然后在上面数字化出需要修改的点。 新建一面文件eldodem，数字化出需要修改的面。-.可修编 .-（2）将 point_11 和 eldodem 在 feature to raster 工具中分别转换成栅格文件poin1 和和 eldo1.（3）在栅格计算器中进行运算假设现在要把个别点象元值修改为0，我们在栅格计算器中填入以下的公式：1.con(IsNull(poin1),eldodem.img,0)-.可修编 .-然后点右下角的 Evaluate 运行因为点象元值很小， 放大可以看到我们修改后的微小的黑点对于面状象元值的修改，相似的，输入以下公式：1.con(IsNull(eldo1),eldodem.img,0)-.可修编 .-然后运行，得到的图如下（大家很可能得到的图与此图有异， 要在属性的 symbol 里设置一下， 使其为拉伸样式）可以看到，修改过的部分变为黑色。-.可修编 .