样本点缓冲区的的光谱特征值
其实提取的最好办法就是把所有的影像给合并成一个大的图像,这样的话不论是加载图像或者是加载点,底图的话都会很快的,你要是单独加载进去的话,十来G的影像是要加载很长时间的(一般的电脑)。

步骤:

一、首先把所有样本点的影像给合并在一个图像中

数据管理工具–>栅格–>栅格数据集–>镶嵌至新栅格

光谱特征提取python代码 光谱数据特征提取_其他

输入栅格:选择你要合成图片的dat文件

输出位置,只能选择文件夹名,所以说提前建立一个文件夹

具有拓展名的栅格数据集名称:我们这里输入名字,以及后缀,我这里生成的 名字.dat文件

像素类型大小:要根据波段的数据值来选择,一般的话选择一下这个就够了

光谱特征提取python代码 光谱数据特征提取_光谱特征提取python代码_02


波段数:根据你要合并的影像的波段数据来看,我这里是8个波段

其他的就可以默认不设置

光谱特征提取python代码 光谱数据特征提取_数据_03


光谱特征提取python代码 光谱数据特征提取_数据_04

ps:图像很大,波段很多的话一般是要跑一个半小时左右,而且内存很大,最好准备充足的内存,我这次跑了8个波段,大约60G

二、给年份shp数据样本点数据建立缓冲区

光谱特征提取python代码 光谱数据特征提取_加载_05


三、加载合并的影像,一个个波段的加载进去,分别提取八个波段的样本点缓冲区内的光谱特征值

光谱特征提取python代码 光谱数据特征提取_其他_06


第一个填写的是缓冲区文件,区域字段的话最好选择建立缓冲区后生成的一个序号文件ORIG_FID ,这样回头好进行连接到缓冲区文件,这样可以进行点的匹配,知道对应的点位置了

光谱特征提取python代码 光谱数据特征提取_加载_07


生成的数据可以右击点击查看,会发现只有数据值

光谱特征提取python代码 光谱数据特征提取_光谱特征提取python代码_08


光谱特征提取python代码 光谱数据特征提取_经验分享_09


连接的时候,我们选择刚刚的ORIG_FID文件连接,要连接的表选择缓冲区,他已经包含了样本点的所有内容数据,还有新建立的序号字段

光谱特征提取python代码 光谱数据特征提取_数据_10


到这一步其实已经结束了,但是我们还可以把生成的结果表给输出出来

通过转换工具就可以实现

光谱特征提取python代码 光谱数据特征提取_数据_11


注意:使用的是哨兵数据,因为没有ENVI3.4以上的版本,所以处理数据的步骤很繁琐,先是下载哨兵数据,接着使用官方发布的Sen2cor工具进行大气校正,接着使用SNAP软件进行重采样,输出各个波段,接着使用ENVI软件进行波段合成,我们使用前8个波段,所以,8个波段给合成在一块了(注意,合成的时候要注意波段的顺序,是可以调整的。链接:ENVI波段顺序的调整 ),最后在Arcgis中打开,进行样本点光谱特征值的提取(多值提取至点和样本点缓冲区光谱特征值提取两类(做的只有这两类,具体的话不知道有多少方法))