现我以一DEM栅格图(名字为eldodem)为例,现在我要修改它的部分象元值,总结出以下三种方法,大家可以参考一下。第一种方法不大实用,但可借鉴,第二三种方法针对的条件不一致,大家可以在具体情况下进行选择。
一、直接运用转换,思路简单,易操作,但实用性
1、栅格——ascii文件——栅格这种方法是先将栅格图用工具
直接转成ascii文件,然后在ascii文本文件中直接修改需要修改的象元的值,修改好后又用工具
转换成栅格图。这种方法可行,但是不实用,因为我们要搜索到指定的象元好像不是那么简单。那么有没有别的比较好的方法呢?
二、在栅格计算器中操作,方法灵活,可操作性强,实用性强准备工作
先要在option中设置保留的栅格范围,通常情况下默认的为相交后的部分,这里我们要保留整个DEM,所以要改为以下设置:
2、通过点的位置修改点象元值
2.1*问题一、要是我想修改图象中的指定的行列的点的象元值(比如说把第100行,200列的点的值修改为0),那该如何操作呢?
这时候我们可以在栅格计算器中输入以下公式:
con(($$rowmap == 99 &$$colmap == 199),0,[eldodem.img])
(栅格图的编号是从0行0列开始的)
2.2*问题二、那若是要修改指定行列范围内的栅格的象元值,比如说把第101行,251列到第401行,301列的部分的象元值改为0,又改如何进行呢?
con(($$rowmap <400 &$$rowmap >100 &$$colmap <300 &$$colmap >250),0,[eldodem.img])
结果如下
2.3 *问题三、这些行和列是给定的,但是实际情况往往不是这样的,那若看到一些要修改的点,我如何通过以上的方法来修改它的属性值呢,这时要解决的首要问题是如何获知这个点的行列值,我们可以通过以下方法来解决:
(1)准备工作,设置spatial analyst中的option,设置所需范围,象素以及捕捉
这个确定后再次打开设定捕捉,使生成的栅格和原栅格象元重合
(2)在栅格计算器中分别做出对应的行列栅格图
行栅格图的生成:在栅格计算器中输入以下公式——
r=$$rowmap
c=$$colmap
(在这里r---行,c------列)
(3)然后可以用identify工具,对要修改的点的行列号进行查询
如图,所要修改的点位于第166行,129列,若要修改此点象元值为6,可以在栅格计算器中输入以下公式:
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,我们在栅格计算器中填入以下的公式:
con(IsNull([poin1]),[eldodem.img],0)
然后点右下角的Evaluate运行
因为点象元值很小,放大可以看到我们修改后的微小的黑点
对于面状象元值的修改,相似的,输入以下公式:
con(IsNull([eldo1]),[eldodem.img],0)
然后运行,得到的图如下
(大家很可能得到的图与此图有异,要在属性的symbol里设置一下,使其为拉伸样式)
可以看到,修改过的部分变为黑色。
假设有两组栅格数据,一组代表2019年中国每月降雨量,一组代表2019年中国每月植被叶面积指数(LAI)。想要得到中国月降水量与LAI的相关性分布,那么需要对两组栅格数据对应的栅格点进行逐栅格的相关性分析。
将降水数据导入栅格栈中,这个过程可以理解为将降水数据按时间顺序从上到下堆叠。同理,按相同的时间顺序将LAI数据堆叠。值得一提的是,stack()函数在堆叠栅格数据时是按文件名拼音和数字大小顺序自动堆叠的,具体规则可以亲自尝试。最后,将这两个栅格栈合并成一个。
对相关性分析函数稍作改变。
以上方法是可以推广的,线性回归函数lm()和相关性分析函数cor()的输入都可以是向量,因此只要函数支持向量输入,理论上讲都可以类比上述过程实现。但是如果函数只支持数据框输入,如gbm包中的函数gbm(),那就只能另辟蹊径了。