(2)在Toolbox中,选择/Regionsoflnterest/SubsetDatafromROIs,选择裁剪图像。
(3)在SpatialSubsetviaROIparameters中,设置以下参数:
在
在"MaskpixelsoutsideofROI?
”项中,选择yes。
裁
(4)剪选择输出路径及文件名14_rad_tianshui.dat,单击OK裁剪图像,结果图如
图5所示:
o
2.3大气校正
在Toolbox中,选择/RadiometricCorrection/FLAASHAtmosphericCorrection,打开FLAASH工具将相关参数填入shp。
(1)文件输入与输出信息项目
单击VnputRadianeeImage按钮,选择上一步准备好的辐射亮度数据的裁剪结果a2014_rad_tianshui.dat”。
在RadianceScaleFactors对话框中选择
Usesinglescalefactorforallbands(Singlescalefactor:
1,在定标中已经进行了单位转换,这里直接输入:
。
0。
单击OutputReflectaneeFile按钮选择输出文件名和路径。
(2)传感器与图像目标信息
Lat:
55.81904406°Lon:
139.02248813°(ENVI5中自动获取)
SensorType:
Landsat-8OLI
FlightTime:
03:
32:
04(从元数据文件中获取或者从下载
GFIightDate:
2014-04-07r网站中获取)。
(o
3(u
4(n
5(d
大7)多光谱设置(MultispectralSettings)
气ED
模e
He如果内存大于4G,则UseTiledProcessing可设置为No;否则设置TileSize(MB)为
模/a00M。
其余按照默认设置。
(9)u单击Apply按钮,执行FLAASH
sMidL-LatitudeSummer
iR
3.1ND0。
计算
直接使用ENVI中的NDVI计算工具计算。
(1)在Toolbox中,选择Spectral/Vegetation/NDVI。
NDVICalculationInputFile面板中,选
择14rad_tianshui_ref.dat图像,点击"OK”。
⑵在NDVI^CalculationParameters对话框中,单击"InputFileType”下拉菜单,选择LandsatOLI,手动填入“Red”为4,“NearIR”为5。
提示:
因为目前NDVI计算工具只能识别原始Landsat8OLI数据,我们这里输入的是经过
处理的结果,1暂时只能手动输入,希望以后改进。
⑶在"OutputDataType”下拉菜单选择输出浮点型(FloatingPoint)。
(4)选择输出路径及文件名,点击“OK”。
t
a
图6
NDVI
由于大气校正后的结果有部分像元为负值,主要集中在阴影地区,这部分区域计算得到的NDVI在[-1,1]之外,为了便于后面的分析,我们这里统一将这部分像元进行处理,即值大于1的变为1,小于-1的变成-1。
去除异常值:
bandmath表达式:
-1>b1<1B1:
NDVI
得到去除异常值文件:
NDVI_去除异常值.dat
3.2掩膜文件制作
在ENVI中打开classimage.dat,这个数据与NDVI数据严格空间配准。
总共有5类土地覆盖类型:
林地:
1
农业用地:
2
城市用地:
3
水体:
4
其他:
5
图8土地覆盖分类图
我们利用这五种土地覆盖类型分别为NDVI生成5个掩膜文件。
下面以林地为例子介绍
这个过程:
(1)打开NDVI结果“NDVI_去除异常值.dat”。
(2)在Toolbox中选择/RasterManagement/Masking/ApplyMask,在SelectInputFile面板中
选择“NDVI_去除异常值.dat”NDVI图像文件,点击MaskOptions->BuildMask。
⑶在MaskDefinition面板中,单击Options->lmportDataRange,选择土地覆盖分类图
classimage.dat。
在最大、最小值中同时填入1。
(4)选择输出路径输出,完成掩膜文件的生成。
不关闭MaskDefinition面板,分别输入其他土地覆盖类型的像元值生成其他四类掩膜文件。
3.3获取阈值
这一步就是求解NDVImax和NDVImin,使用上一步获取的掩膜文件分别对NDVI图像文
件进行统计,在一定置信度范围内获取每个掩膜文件(也就是土地覆盖类型)对应的最大和
最小NDVI值。
下面以林地为例介绍这个过程:
(1)在Toolbox中,选择/Statistics/ComputeStatistics,如下图所示,分别选择ndvi文件
“NDVI_去除异常值.dat”作为统计文件,以及掩膜文件。
(2)在统计类型面板中,勾选直方图(Histogram)选项。
(3)得到统计结果。
地物覆盖类型的NDVI阈值,如下表:
表
2014年
2018年
土地覆盖类型
NDVImin(NDVIsoil)
NDVImax(NDVIve
g)
NDVImin(NDVIsoil)
NDVImax(NDVIveg)
林地
0.121569
0.678431
0.32765
0.71365
农业用地
0.147411
0.214603
0.163286
0.275404
城市用地
-0.035294
0.270588
0.184314
0.505882
水体
0.08948
0.165195
0.053726
0.433725
其他
0.406863
0.427602
0.129412
0.615686
3.4生成参数文件
这一步是根据上面得到的NDVI阈值分别生成NDVIsoil和NDVIveg参数文件,主要使用bandmath工具,表达式如下:
NDVIsoil:
b1*0.121569+b2*0.147411+b3*(-0.035294)+b4*0.08948+b5*0.406863
其中b1:
林地掩膜文件b2:
农业用地掩膜文件b3:
城市用地掩膜文件
b4:
水体掩膜文件b5:
其他用地掩膜文件
NDVIveg:
b1*0.678431+b2*0.214603+b3*0.270588+b4*0.165195+b5*0.427602
其中b1:
林地掩膜文件b2:
农业用地掩膜文件b3:
城市用地掩膜文件
b4:
水体掩膜文件b5:
其他用地掩膜文件此生成的参数文件为2014年,结果图如图10下:
图10
这一步是根据上面得到的NDVI阈值分别生成NDVIsoil和NDVIveg参数文件,主要使用bandmath工具,表达式如下:
NDVIsoil:
b1*0.32765+b2*0.163286+b3*(0.184314)+b4*0.053726+b5*0.129412
其中b1:
林地掩膜文