edvi操作课堂作业3

《edvi操作课堂作业3》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、马钰琪201211191003实习三实习三主要内容用直方图，bandmath等方法对图像进行分割。虽然都是根据光谱影像，但分割和分类有一定的区别，可以说，最后的目的还是为了分类。分割的功能大致相当于arcgis中的矢量化，把需要的部分用函数关系等方式提取出来，而分类是为了把地物都区别开。1.直方图的使用1）打开图像菜单中Enhance-interactivestretching，直方图只能选择一个波段的，所以查看时要选择查看的波段。2）要在屏幕上使草地消失可以打开G波段，将最小值设为150，所有草地部分就都变黑了。继续试验可以发现，当RGB三个波段的值拉伸最小值150以

2、上时，留下没有在拉伸范围外的就只有天空部分了，所以考虑用这个方法将天空和其他区域分开。图12.用bandmath去除天空因为天空部分三个值都大于150，可以把天空分割出来，所以考虑用算式1-((b1-150)*(b2-150)*(b3-150))得到下图10马钰琪201211191003实习三图2然后将这个新得到的图与原图相乘，选出青山的部分。得到最后图图33.彩色图像的分割3.1提取兰花1.如图观察直方图可得知，兰花蓝色的值比其他波段大，所以考虑用这个条件选出兰花。输入算式(b1-b2)*(b1-b3)，其中b1对应蓝波波段，b2和b3对应红波和绿波波段。10马钰琪2

3、01211191003实习三图42.输入算式(b1-b2)*(b1-b3)，其中b1对应蓝波波段，b2和b3对应红波和绿波波段。选出兰花部分后，与原图层相乘，得到只有兰花部分的图像。如下图，从左到右依次为：原图，选出的兰花部分和最后相乘得到的兰花图。图53.2去除娃娃的噪声1.分别打开三个波段的图，可以看出在背景以外的区域，R图像上的点比对应G上的点更亮，而在背景区域，R上的点亮度小于G，用这个特点我们可以把背景提取出来。10马钰琪201211191003实习三图6该图的两个波段2.用bandmath输入公式float(b1)/float(b2)，其中吧b1为R波段，b

4、2为G波段。得到新图如下，可以查看它的直方图。如图可知，0.98左右是一个分界值，0.98以上是背景部分，以下是有人的部分。图73.用bandmath输入公式b1gt0.98,产生图像m2，得到背景。10马钰琪201211191003实习三图84.使用原始图像的RGB分别作为b1，使用m2作为b2，进行代数运算b1*b2,产生新的图像m3，进行彩色显示。图95.提取背景，进行运算（1-b1）*255，得到背景图像m4，然后将背景与原图像中的R图层叠加，b1+b2，得到m5图1010马钰琪201211191003实习三6.将m5和前一步得到的新G和B图层叠加，得到最后图像

5、。图11可以看到，娃娃的头和手有一点瑕疵，大体上提取出了娃娃。3.3提取水体信息1.查看直方图和光谱信息。用波段4,3,2进行标准假彩色显示，在scoll窗口上右键，打开Zprofile，可以得到波段的信息。Enhance-InteractiveStretching可以打开直方图查看。用Dataparameter还可以设置范围输出。图122.查看空间剖面的光谱变化。当在图层中画的线变化，即剖面变化时，图上的值也相应变化。10马钰琪201211191003实习三图133.拉伸突出水体。打开波段5灰度显示，enhance-interactivestretching,打开直方

6、图，设置拉伸最大值为25，得到图前后对比，除水体以外的部分被拉伸。图144.合成图像。10马钰琪201211191003实习三使用波段5,4,2合成，再将波段5按照上一步拉伸，查看结果。图15（不是很明白这个的意义是什么）5.查看不同类型地物光谱差异。在Image窗口中，移动光标位置到长江中，然后在SpectralPixelEditor窗口中，点击Edit->Extract8PixelAverage。如图所示，依次得到江中（蓝），玄武湖（紫），山地（红，建筑（白）和农田（黄）的值。图166.比较地物光谱差异。观察上图可以看到，蓝色的线和紫色的线和黄色的线差别较小，说明农

7、田与他们的差异较小，而红色差距较大，说明山地和水体差别大，而两种水体之间由于各种性质差异，也有很大的区别。10马钰琪201211191003实习三7.提取水体。输入公式b2-b5，b2对应波段2，b5对应波段5。得到图如下。可以看出，水体被清晰得提取出来了。图188.直方图拉伸确定地物界限。图19打开这个图的直方图，分别用40,25,10作为最大值进行拉伸，得到结果如下。可以看到，40可以分割出江水和玄武湖，25只能分割出玄武湖，而10图上的一切除了边缘都没有差别。9.密度分割。打开tools-colormapping-densitysl