遥感ENVI实验报告.docx
- 文档编号:11518934
- 上传时间:2023-06-01
- 格式:DOCX
- 页数:37
- 大小:2.11MB
遥感ENVI实验报告.docx
《遥感ENVI实验报告.docx》由会员分享,可在线阅读,更多相关《遥感ENVI实验报告.docx(37页珍藏版)》请在冰点文库上搜索。
遥感ENVI实验报告
序言………………………………………………………………………………3
一、试验目………………………………………………………………3
二、试验内容………………………………………………………………3
三、试验时间………………………………………………………………3
四、组织人员………………………………………………………………3
1.专题概述……………………………………………………………………4
2.处理步骤介绍……………………………………………………………4
2.1图像获取……………………………………………………………4
2.2数据读取和定标…………………………………………………4
2.3图像配准……………………………………………………………5
2.4大气校正……………………………………………………………5
2.5反演模型构建及模型应用………………………………………5
2.6植被改变……………………………………………………………6
3.具体处理过程……………………………………………………………7
3.1数据预处理…………………………………………………………7
3.1.1安装环境小卫星数据处理补丁…………………………………7
3.1.2数据处理和定标…………………………………………………7
3.1.3工程区裁剪………………………………………………………9
3.1.4图像配准…………………………………………………………14
3.1.5大气校正…………………………………………………………17
3.1.6裁剪浑善达克区…………………………………………………23
3.2植被覆盖度反演……………………………………………………27
3.2.1计算归一化植被指数……………………………………………27
3.2.2计算植被覆盖度…………………………………………………28
3.3植被改变监测………………………………………………………29
3.3.1植被覆盖区提取…………………………………………………29
3.3.2植被改变检测……………………………………………………31
3.4结果后期处理与应用……………………………………………32
3.4.1植被改变区域图背景值处理…………………………………32
3.4.2植被改变区域制图………………………………………………33
试验心得………………………………………………………………………36
序言
一、试验目
1、掌握ENVI软件基础操作。
2、掌握卫星影像预处理基础步骤。
3、经过实习,学会自己去处理部分问题。
4、深入提升学生分析问题、处理问题能力,增强实践技能,并培养学生勇于动手、勤于动手、热爱本专业思想。
5、深刻地了解和巩固基础理论知识, 掌握基础技能和动手操作能力, 提升综合观察分析问题能力
二、实习内容
1、了解ENVI基础操作。
2、实现影像图像几何校正、融合、镶嵌及剪裁。
3、掌握ENVI对影像信息提取
4、了解ENVI部分应用分析
专题:
基于环境小卫星草原荒漠化监测
一、专题概述
浑善达克地域位于内蒙古草原阴山北麓锡林郭勒高原中部,是亚洲草原荒漠化土地东部边缘区关键组成部分,经纬度在东经114°55’~116°38’,北纬41°46’~43°07’之间,平均海拔高度在1100米左右。
退化区属温带半干旱区、中温带干旱大陆性季风气候。
浑善达克退化土地多为沙地,植被稀少,尤其是春季地面回暖解冻,地表裸露,多细沙土,狂风起时沙尘弥漫,形成沙尘天气。
多年来频频发生在京津地域春季沙尘天气与该地域生态环境恶化,人地关系严重失调相关。
据统计,京津地域沙尘暴70%沙源来自于这个区域。
依据该区域不一样时期植被覆盖数据能够实现该区域植被改变监测。
本专题具体介绍了利用环境小卫星CCD-1A图像反演植被覆盖图完整步骤,专题包含环境小卫星数据读取、辐射定标、大气校正、植被覆盖反演模型建立、遥感反演过程、植被覆盖改变监测等内容。
除了使用ENVI主模块功效外,还需要用到FLAASH大气校正扩展模块、IDL开发环境小卫星数据读取补丁、波段运算等功效。
二、处理步骤介绍
依据环境小卫星CCD数据特点及草原植被改变监测要求,采取技术路线为:
先对环境小卫星CCD数据进行数据预处理:
数据读取、辐射定标、大气校正、研究区裁剪,建立反演模型,利用波段运算工具,反演出整个浑善达克地域归一化植被指数、植被覆盖度,依据不一样时期植被覆盖数据,实现草原土地退化遥感监测。
步骤说明:
2.1图像获取
本文基础数据包含年土地覆盖图数据、环境小卫星CCD-1A数据。
环境小卫星CCD-1A数据能够在环境保护部卫星环境应用中心无偿下载获取,土地覆盖图数据制作时间是年8月份,本专题选择了影像质量良好年8月11日环境星数据作为数据源。
环境小卫星CCD数据下载地址:
.cn
2.2数据读取和定标
网上无偿获取HJ-1A卫星CCD分发格式为Geotiff,每一个波段为一个Geotiff文件,并提供一个元数听说明(.XML)。
能够依次打开每个波段,并用元数听说明中提供定标参数用BandMath工具进行手动定标,再用Layerstcaking功效将4个波段图像合成为一个多波段数据文件。
本专题中用是环境小卫星读取补丁,直接读取CCD数据,输出结果为一个多波段ENVI标准栅格文件,并带有中心波长等信息,可直接输出定标结果(辐射亮度)。
补丁下载地址为:
2.3图像配准
获取年土地利用分类图已经过正确地理定位,以该景作为基准影像,对环境小卫星数据进行配准,便于后续改变监测。
2.4大气校正
对于环境小卫星数据,提供光谱响应函数,可用FLAASH大气校正模块中进行大气校正,波谱响应函数下载地址:
。
2.5、反演模型构建及模型应用
2.6植被改变监测
依据年8月和年8月植被覆盖数据,采取波段运算,实现草原植被改变遥感监测。
三、具体处理过程
本专题数据存放在“19-草原植被改变遥感监测”文件夹内。
3.1数据预处理
3.1.1安装环境小卫星数据处理补丁
将ENVI_HJ1A1B_Tools.sav补丁放在:
home\ITT\IDL\IDL80\products\envi48\save_add目录下。
3.1.2数据读取和定标
主菜单->File->OpenExternalFile->HJ-1A/1BTools,打开环境小卫星数据读取补丁。
在HJ-1A/1BToolsV3.0面板中,选择CCD,点击InputFile输入“\1-环境小卫星数据\HJ1A-CCD2-2-64-0811-L0154793\154793”文件夹中.xml文件,点击OutputPath设置数据输出路径,勾选“Calibration”“LayerStacking”两个选项,
单击Apply按钮。
3.1.3工程区裁剪
因为整景影像数据范围非常大,本专题工程区只是其中一小部分,在进行大气校正之前,先将浑善达克以及周围区域裁剪出来。
(1)打开上一步生成文件:
HJ1A-CCD2-2-64-0811-L0154793_Calbrated_LayerStacking.img,
(2)主菜单->File->SaveFileAs->ENVIStandard,弹出NewFileBuilder面板;
(3)在NewFileBuilder面板中,单击ImportFile,弹出CreateNewFileInput
—File面板;
(4)在CreateNewFileInputFile面板中,选中SelectInputFile列表中裁剪
数据,单击SpatialSubset按钮;
(5)在SelectSpatialSubset面板中,单击Image,弹出SubsetbyImage对话框,
(6)在SubsetbyImage对话框中,按住鼠标左键拖动图像中红色矩形框确定裁剪区
域,裁剪出包含浑善达克区域一部分,单击OK;
(7)在SelectSpatialSubset面板中,能够看到裁剪区域信息,单击OK;
(8)在CreateNewFileInputFile对话框中,单击OK;
(9)在NewFileBuilder,单击Choose设置输出文件名0811-Cal-sub.img及路
径,单击OK。
3.1.4图像配准
下面以土地利用图作为基准影像对环境小卫星图像进行图像配准。
(1)分别打开和显示基准影像“\19-基于环境小卫星草原植被改变监测\浑善达克
年8月土地利用分类图.img”。
(2)主菜单->Map->Registration->SelectGCPs:
ImagetoImage,打开几何校正模块。
(3)选择显示年土地利用分类图文件Display为基准影像(BaseImage),
显示环境星文件Display为待校正影像(WarpImage),点击OK进入采集地面
控制点。
(4)打开Tools->Link->Geograficlink,将两个窗口都选择为on,单击确定,找到
定位大致区域后,再Tools->Link->Geograficlink,改为off,关闭链接。
(5)在两个Display中找到相同区域,在Zoom窗口中,点击左小下角第三个按钮,
打开定位十字光标,将十字光标到相同点上,点击GroundControlPoints
selection上AddPoint按钮,将目前找到点加入控制点列表。
(6)用一样方法继续寻求其它点,当选择控制点数量达成3时,RMS被自动计
算。
GroundControlPointsSelection上Predict按钮可用,选择
Options->AutoPredict,打开自动估计功效。
这时在BaseImage上面定位点,
WarpImage上会自动估计区域。
(7)完成控制点选择,RMS值小于1个像素,点击GroundControlPoints
Selection上File->SaveCoefficientstoASCII,将控制点保留。
(8)在GroundControlPointsSelection上,选择Options->WarpFile(asImage
Map),选择校正文件(HJ数据文件)。
(9)在校正参数面板中(图专5-7),默认投影参数和像元大小与基准影像一致,30米。
(10)重采样选择NearestNeighor,背景值(Background)为0.
(11)OutputImageExtent:
默认是依据基准图像大小计算,能够做合适调整。
(12)选择输出路径和文件名,单击Ok按钮。
3.1.5大气校正
环境小卫星提供了波谱响应函数,以文本形式提供,第一列表示波长(nm),后面四列分别表示4个波段对应波长波谱响应值。
需要制作波谱曲线来描述波谱响应函数,用于大气校正。
制作波谱曲线:
(1)主菜单Window->StartNewPlotWindow,打开ENVIPlotWindow面板,在
波谱绘制窗口中,选择File->InputData->ASCII,导入“681_HJ1ACCD2.txt”
文本文件,如图,自动将第一列作为x轴,后面4列作为y轴,波长单位
选择Nanometers,单击OK。
(2)如图,在绘制窗口生成了4条曲线,选择Edit->DataParameters,编辑每条线名称为b1,b2,b3,b4,便于区分。
(3)选择File->SavePlotAs->SpectralLibrary,在OutputPlotstoSpectral
Library面板中,单击SelectAllItems,单击OK。
(4)在OutputSpectralLibrary面板中,有输出曲线相关参数设置,按默认,
选择输出路径和文件名,单击OK,将波谱曲线保留为波谱库文件:
环境1A星
CCD2光谱响应.sli。
FLAASH大气校正:
第一步:
数据准备
FLAASH对图像文件有以下多个要求:
(1)数据是经过定标后辐射亮度(辐射率)数据,单位是:
(μW)/(cm2*nm*sr)。
(2)数据带有中心波长(wavelenth)值,假如是高光谱还必需有波段宽度(FWHM),这两
个参数都能够经过编辑头文件信息输入(EditHeader)。
(3)数据类型:
支持四种数据类型:
浮点型(floating)、长整型(longinteger)、整型
(integer)和无符号整型(unsignedint)。
数据存放类型:
ENVI标准栅格格式文件,
且是BIP或者BIL。
4)波谱范围:
400-2500nm。
此次用环境小卫星经过以上处理,已经定标为W*m^(-2)*sr^(-1)*um^(-1)单位、浮点型辐射率数据,有中心波长信息,下面将BSQ格式转成BIL格式。
选择主菜单BasicTools->ConvertData(BSQ、BIL、BIP),选择已经经过定标和配准数据0811-cal-jz.img,在ConvertFileParameters中,OutputInterleave选择BIL,选择ConvertInPlace:
yes,单击OK。
第二步:
设置参数进行FLAASH大气校正
(1)主菜单Spectral->FLAASH打开FLAASH大气校正模块;
(2)点击InputRadianceImage,选择BIL格式环境小卫星数据
0811-cal-jz.img,在RadianceScaleFactors面板中选择Usesinglescale
factorforallbands,因为定标辐射量数据与FLAASH辐射亮度单位相差
10倍,所以在此Singlescalefactor选择默认:
10,单击OK;
(3)设置输出文件及路径设置
(4)传感器基础信息设置:
成像中心点经纬度、传感器高度、成像区域平均高度、成像
时间设置,这些都能够从数据头文件中读取:
HJ1A-CCD2-2-64-0811-L0154793.XML。
(5)大气模型,选择MLS,气溶胶模型,选择Rural,气溶胶反演方法选择None,能见
度设置为40km。
(6)单击MultispectralSetting按钮,在FilterFunctionFile导入之前做好光
谱响应曲线“环境1A星CCD2光谱响应.sli”,单击OK;
(7)单击AdvancedSettings,在高级设置中,TileSize默认是Cashsize大小,
手动改为100Mb,单击OK;
(8)设置好后,在大气校正模块面板中,单击Apply。
(9)大气校正完成后,检验大气校正结果,分别加载校正前后图像,将两幅影像进行地理链接,移动到植被区域(植被波谱曲线比较特殊),在影像上右键,选择ZProfile(Spectrum)打开光谱曲线窗口,显示两幅图像同一位置光谱曲线图。
3.1.6裁剪浑善达克区
(1)显示大气校正后图像0811-cal-jz-FLAASH.img,在image窗口选择
Overlay->Vectors,打开VectorParameters面板,选择File->OpenVectorFile,
打开hunshandake.evf;
(2)在AvailableVectorsList面板中选择该矢量文件,点击LoadSelected,选择
显示图像Display,单击OK,矢量叠加在影像上;
(3)在AvailableVectorsList面板中,选择File->ExportLayerstoROI,在Select
DataFiletoAssociatewithnewROI面板中,选择0811-cal-jz-FLAASH.img,单击OK,在ExportEVFLayerstoROI中,选择ConvertallrecordsofanEVFlayertooneROI,单击OK,将矢量转为一个ROI;
(4)在图像窗口,选择Overlay->RegionofInterest,打开ROI面板,浑善达克ROI
显示在图像上,在ROITool面板中,选择File->SubsetDataviaROIs,在Select
InputFiletoSubsetviaROI面板中,选择0811-cal-jz-FLAASH.img,单
击OK;
(5)在SpatialSubsetviaROIParameters面板中选择浑善达克ROI,Maskpixels
outsideofROI选择Yes;
(6)设置输出路径及文件名0811-yanjiuqu.img,单击OK。
到此,数据预处理工作已经全部完成,下面介绍怎样应用ENVI进行植被覆盖度反演和植被改变监测。
3.2植被覆盖度反演
3.2.1计算归一化植被指数
(1)选择主菜单->File->Transform->NDVI,打开NDVI计算模块。
(2)选择裁剪后数据:
0811-yanjiuqu.img。
(3)在NDVICalculationParameters面板中进行对应设置,环境星波段参数与TM数
据相同,应用被植被强吸收红光波段(环境星第3波段)和被植被强反射近
红外波段(环境星第四波段)计算归一化植被指数。
数据类型默认为浮点型。
(4)选择文件保留名(0811-yanjiuqu-NDVI)和路径,单击OK实施。
3.2.2计算植被覆盖度
计算植被覆盖度FC采取是混合像元分解法,将整景影像地类大致分为水体、植被和建筑,具体计算公式以下:
FC=[(NDVI-NDVISoil)/(NDVIVeg-NDVISoil)]
(1)
其中,NDVI为归一化植被指数,NDVISoil为完全是裸土或无植被覆盖区域NDVI值,NDVIVeg则代表完全被植被所覆盖像元NDVI值,即纯植被像元NDVI值。
取经验值NDVIVeg=0.70和NDVISoil=0.00,且有,当某个像元NDVI大于0.70时,FC取值为1;当NDVI小于0.00,FC取值为0。
利用ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入:
(b1gt0.7)*1+(b1lt0.)*0+(b1ge0andb1le0.7)*((b1-0.0)/(0.7-0.0))
b1:
选择NDVI图像
3.3植被改变监测
3.3.1植被覆盖区提取
(1)8月植被覆盖区提取
利用ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入:
(b1le0.3)*0+(b1gt0.3)*1
b1:
选择HJ1A-CCD2-0811-cal-sub-flaash-ref-hunshandake-VCI图像,0.3为经验值。
选择文件保留名(年8月植被覆盖区)和路径,单击OK实施,得到年8月植被盖度区。
(2)8月植被覆盖区提取
获取8月标准土地利用分类图,分类代码与地物类型以下表所表示,耕地、林地、草地,属于植被覆盖区域。
ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入:
(b1ge1andb1le3)*1+(b1lt1)*0+(b1gt3)*0
b1:
选择“浑善达克年8月土地利用分类图”。
选择文件保留名(08植被覆盖图.img)和路径,单击OK实施,得到年8月植被盖度图。
3.3.2植被改变检测
利用ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入:
b1-b2
b1:
选择年8月植被覆盖区图像。
b2:
选择年8月植被覆盖区图像;
选择文件保留名(浑善达克植被覆盖改变(-))和路径,单击OK实施,得到浑善达克年至年植被盖度改变区域图像。
3.4结果后期处理与应用
3.4.1植被改变区域图背景值处理
(1)显示上面得到-植被改变区域图,主菜单Basic
Tools->Masking->Applymask,选择--植被改变区.img,单击SelectMask
Band,选择掩膜图像,单击OK;
(2)在ApplyMaskParameters面板,设置MaskValue为5,设置输出路径与文件名,
单击OK;
3.4.2植被改变区域制图
将掩膜背景值后改变图像--植被覆盖改变区-maskback.img打开,在主图像窗口Overlay->DensitySlice,对-1、0、1值以及背景值分别进行密度分割。
点击Apply查看密度分割结果。
在DenstySlice面板,选择File->OutputRangetoClassImage,能够将分割结果输出为ENVI分类格式,
打开Overlay->Classification面板,选择输出分类图像,叠加各类显示,制图
实习心得
经过此
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 遥感 ENVI 实验 报告