1011植被覆盖度提取的后续处理.docx

植被覆盖度提取的后续制图处理1、我们打开植被覆盖提取完成的影像文件名，打开影像如下：2、再打开之前裁剪完成后的543波段影像对照着看：我的文件目录，打开。对比图：对比可以发现，影像上树越密，FVC影像上越白。2、想知道FVC上每一个象元的值，可以在影像上右键找到找到中间的，或者双击鼠标，得到。3、我们选择后，可以让两幅影像同时移动，方便观察。此时，出现两个Date值。我们接下来的后续制图操作就是将Date值中大于1的情况处理掉。（12月22日FVC影像的背景值也大于0。）（植被覆盖完后的后续数据处理，因为0=FVC=1，而我们目前的影像中有大于1或者小于0的情况，必须处理这些象元。）4、在影像中右键找到此时可以关闭543波段影像了。打开，得到了关于本影像的统计结果。上面一半是图像，有几个加号，总体上反映了影像大部分值集中分布的地方，当然这部分只是一个定性的描述，下面的才是定量的描述，DN值中文意思是灰度值，应该是位于0~1，小于0的情况应该没有，大于1的全是不符要求的，当然这是光谱原因，所以，我们要通过某种方法将大于1的值变成1.（12月22日做的影像没有大于1的情况，只是边界都大于0，所以接下来应该是将边界外的批量处理为0，所以此次跳过7、8步骤，直接用掩膜操作.）图中“+”指大部分数据所在的范围，但这只是一个大概的定性描述，定量描述在下面的表格里。此处普及几个概念：①：灰度值（0——1）②：某个DN值象元的数量③：DN值百分比的和④：DN值总共的数量⑤：百分比我们接下来的操作就是把DN值中大于1的象元批量处理为1。5、用ENVI 的Band Math进行处理，并且处理完后要用“Build Mask”将黄石市边界以外的象元处理为0。打开波段计算器：，6、公式：处理FVC值大于1或者小于0的象元：(b1 lt 0)*0+(b1 gt 1)*1+(b1 le 1 and b1 ge 0)*b1关于逻辑符号的说明：ge大于或等于，gt大于,le小于或等于，lt小于and 并且or 或者b1代表FVC 影像中植被覆盖情况这个公式表示b1小于0时作为0，b1大于1时作为1，而b1小于等于1并且大于等于0时保留b1.注意：此处的字母大小写都可以，但是要保持统一。不能用搜狗输入法，只能用微软输入法写出来。7、将写好的公式，复制到波段计算器中，保存为。完成后弹出的对话框如下：并选中我们目前要操作的影像完成后，命名为。在中会生成下面的一个文件：打开此文件，注意观察当光标在有植被覆盖的地方和没有植被覆盖的地方Data是不同的，有0和1的变化。8、改进计算后要将之前的计算结果覆盖掉，否则还是有问题。将公式（如果公式考虑的情况不够全面，会出现的问题是大量在0到1之间的象元都被赋值为0了。所以，我们的改进公式就是将象元值小于0的弄成0，大于1的弄成1，确确实实是介于0到1之间的我们将它保留，也就是说还是弄成b1。遇到问题时一定要学会检查问题。）复制到波段计算器中，，记住，在此时弹出的对话框中应该选中的是之前的植被覆盖文件，而不是刚才用部分公式形成的文件。。完成之后打开处理图像：9、掩膜处理（1）打开影像，对其进行建立掩膜然后，再打开黄石市边界，弹出命名保存。点击“应用”，就会在中生成文件名。打开这个掩膜，再应用掩膜，命名为。12月22日接下来的背景就是合理的可以观察到黄石市边界范围之外的地方DN值都为0。10、植被覆盖分级图（1）进行植被覆盖度分级计算公式：b1:表示处理后的植被覆盖度影像，b2：表示掩模影像（b1 ge 0 and b1 le 0.1 and b2 eq 1）*1+(b1 gt 0.1 and b1 le 0.3)*2+(b1 gt 0.3 and b1 le 0.5)*3+(b1 gt 0.5 and b1 le 0.7)*4+(b1 gt 0.7 and b1 le 1)*5+(b2 eq 0)*0公式表示把植被覆盖度影像与掩模影像的象元分为5个等级，但是边界以外的为0的也是一种情况，要考虑，所以一分级要将掩膜影像考虑其中。