遥感图像处理课程设计长沙市植被覆盖率变化.docx
- 文档编号:27991742
- 上传时间:2023-07-07
- 格式:DOCX
- 页数:24
- 大小:7.82MB
遥感图像处理课程设计长沙市植被覆盖率变化.docx
《遥感图像处理课程设计长沙市植被覆盖率变化.docx》由会员分享,可在线阅读,更多相关《遥感图像处理课程设计长沙市植被覆盖率变化.docx(24页珍藏版)》请在冰豆网上搜索。
遥感图像处理课程设计长沙市植被覆盖率变化
实习序号及题目
长沙市3年的植被等的变化
学生姓名
专业班级及编号
地信1班
任课教师姓名
实习指导教师姓名
实习地点
实习日期时间
实习目的
1、本课程设计是综合运用ENVI进行遥感图像处理
2、掌握使用ENVI进行多幅图像镶嵌的方法
3、掌握使用ENVI进行图像裁剪的方法
4、掌握使用ENVI进行图像分类来分析问题的方法
5、掌握使用ENVI进行遥感动态监测对比的方法
实习内容
1、数据下载:
课程设计所需遥感影像数据的下载
2、图像镶嵌:
将下载的多幅遥感影像进行拼接
3、图像裁剪:
裁剪出长沙市的遥感影像
4、图像分类:
对长沙市的主要地物进行分类
5、遥感动态监测:
对长沙是2006-2009年分类后的遥感影像进行动态监测,分析本地区植被的地物的变化
基本原理
1、NDVI(归一化植被指数)
NDVI=(NIR-R)/(NIR+R),或两个波段反射率的计算。
它的取值范围是-1—1。
对于陆地表面覆盖来说,云、雨、雪在可见光比近红外波段有较高的反射作用所以NDVI为负值,岩石、裸土的NDVI一般为0,有植被覆盖的地方一般大于0。
2、图像镶嵌(mosaic image)
是具有地理参考的若干相邻的图像合并成一幅图像或一组图像,需要拼接的图像必须含有地图投影也就是说图像必须经过几何校正处理,虽然所有的输入图像可以具有不同的投影类型,不同的象元大小,但必须有相同的波段数。
在进行图像拼接时需要确定一幅参考影像,参考图像作为图像拼接的基准,决定输出图像的地图投影和象元大小和数据类型。
3、图像裁剪
目的是将研究之外的区域去除。
常用的方法是按照行政区划边界或自然区划边界进行图像裁剪。
在基础数据生产中,还经常要进行标准分幅裁剪,按照ENVI的图像裁剪过程,可分为规则裁剪和不规则裁剪。
4、遥感图像分类
就是通过对遥感图像中地物的光谱信息和空间信息进行分析,选择特征,将图像中每个象元按照某种规则或算法划分为不同的类别,然后获得遥感图像与实际地物的对应信息,从而实现遥感图像的分类。
一般的分类方法可分为两类:
监督分类和非监督分类。
将多源数据应用于图像分类中,发展成基于专家知识的决策树分类。
5、遥感动态监测
指利用多时相获取的覆盖同一地表区域的遥感影像及其它辅助数据来确定和分析地表变化。
它利用计算机图像处理系统,对不同时段目标或现象状态的变化进行识别、分析;它能确定一定时间间隔内地物或现象的变化,并提供地物的空间分布及其变化的定性与定量信息。
6、NDVI计算方法
NDVI植被指数1植被指数(VegetationIn-dex),又称光谱植被指数,是指由遥感传感器获取的多光谱数据,经线性和非线性组合而构成的对植被有一定指示意义的各种数值。
研究表明,植被指数与植被覆盖度具有密切的相关关系,其中归一化植被指数(NDVI)是植被生长状态及植被覆盖度的最佳指示因子。
TM影像中NDVI的计算公式为:
NDVI=TM4-TM3TM4+TM3(3-1)
式中,TM4和TM3分别为LandsatTM影像的第四近红外)和第三(红)波段亮度值131212植被覆盖度1像元二分模型[7-9]假设像元只由两部分组成:
植被覆盖地表和无植被覆盖地表,遥感传感器观测到的光谱信息也只由这两个组分因子线性加权合成,它们各自的面积在像元中所占的比率即为因子的权重,其中植被覆盖地表占像元的百分比即为该像元的植被覆盖度,因而可以使用此模型来估算植被覆盖度1其计算公式为:
fc=(S-Ssoil)/(Sveg-Ssoil)(3-2)
其中Ssoil与Sveg是像元二分模型的两个参数,因此只要知道这两个参数就可以根据以上公式利用遥感信息来估算植被覆盖度1由于归一化植被指数(NDVI)也是一种由遥感传感器所接收的地物光谱信息推算而得的反映地表植被状况的定量值1根据像元二分模型,一个像元的NDVI值可以表达为由绿色植物部分所贡献的信息NDVIveg与由裸土部分所贡献的信息NDVIsoil这两部分组成,同样满足公式(2-2)的条件,因此可以将NDVI代入公式(3-2),得:
fc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)(3-3)
其中,NDVIsoil为无植被覆盖区域(裸土)的NDVI值,NDVIveg为完全植被覆盖区域的NDVI值1NDVIsoil与NDVIveg的计算公式为:
NDVIsoil=(fcmax#NDVImin-NDVImax)/(fcmax-fcmin)
NDVIveg=[(1-fcminNDVImax-(1-fcmax)NDVImin]/(fcmax-fcmin)(3-4)
fcmax、fccmin为区域中像元植被覆盖度可能的最大值与最小值,NDVImax、NDVImin为图像中NDVI在给定置信度的置信区间的最大值与最小值,由于长沙地处亚热带地区,通常情况下植被覆盖度的最大值都可以取100%,植被覆盖度的最小值一般都取0,将fcmax=100%,fcmin=0代入公式(3-4)得:
NDVIsoil=NDVImin
NDVIveg=NDVImax(3-5)
数据准备
1、2006年长沙TM影像_L4_30米4幅
2、2009年长沙TM影像_L4_30米4幅
3、长沙Shp矢量文件
4、长沙DEM影像_L4_30米7幅
操作方法及过程
1、实习概述
本次研究,结合研究目标和实际情况,综合权衡各种因素后,主要采用的卫星遥感影像为长沙市2个时相的LandsatTM影像,成像时间分别为2006年和2009年,影像分辨率为30m,影像已经过大气校正、几何精校正及裁减处理。
数据来源:
TM影像数据从中国科学院遥感与数字地球研究所(
下图为影像查询结果
DEM影像数据从地理空间数据云(
下图为影像查询结果
Shp矢量数据从GIS帝国论坛分享中下载。
2、长沙矢量文件提取
利用ArcGIS10提取要素来导出长沙市的矢量Shp数据文件。
(1)用ArcGIS10打开中国地级市矢量图
(2)选择要素-长沙,使用筛选工具导出长沙市的矢量Shp数据文件。
导出结果后第二图
3、图像镶嵌
(1)打开影像文件2006年
(2)设置参数
(3)输出参数设置
(4)图像镶嵌结果
2009年的图像镶嵌按此操作
4、图像裁剪
(1)打开长沙矢量Shp文件
(2)设置文件打开参数
(4)将矢量文件转换为ROI
(5)输出图像裁剪后的结果
2009年的图像裁剪按此操作
5、图像分类
(1)按分类要求绘制感兴趣区域
(2)检查感兴趣区域的可分离度
(3)设置输出参数
(4)分类后结果输出
(5)图像分类后处理-majority分析
2009年的图像分类操作按此进行
6、遥感动态监测
(1)选择2006年和2009年分类后的遥感时相图
(2)定义变化分析类别
(3)设置输出参数
(4)得出分析结果,分析数据报告导出
PixelCounts
Unclassified植被[Green]3550points建筑[Red]15187points水系[Blue]4773pointsRowTotalClassTotal
Unclassified1018084410151598099952421128224411282244
植被[Green]15849points7036975271680847087105668235206823520
建筑[Red]53838points35591326792761688509706247307604730760
水系[Blue]4817points758964288360479990163778163778
ClassTotal11316350897040326201999335000
ClassChanges113550636987239316901336000
ImageDifference-34106-214688321105617042800
Percentages
Unclassified植被[Green]3550points建筑[Red]15187points水系[Blue]4773pointsRowTotalClassTotal
Unclassified89.96611.3173.0915.615100.000100.000
植被[Green]15849points6.21858.76732.3291.131100.000100.000
建筑[Red]53838points3.14529.86864.4427.565100.000100.000
水系[Blue]4817points0.6710.0480.13885.688100.000100.000
ClassTotal100.000100.000100.000100.0000.0000.000
ClassChanges10.03441.23335.55814.3120.0000.000
ImageDifference-0.301-23.93380.55075.4450.0000.000
Area(SquareMeters)
Unclassified植被[Green]3550points建筑[Red]15187points水系[Blue]4773pointsRowTotalClassTotal
Unclassified9162759600.00913643100.0072899100.004717800.0010154019600.0010154019600.00
植被[Green]15849points633327300.004744512000.00762378300.00950400.006141168000.006141168000.00
建筑[Red]53838points320321700.002411348400.001519658100.006355800.004257684000.004257684000.00
水系[Blue]4817points68306400.003859200.003243600.0071991000.00147400200.00147400200.00
ClassTotal10184715000.008073362700.002358179100.0084015000.000.000.00
ClassChanges1021955400.003328850700.00838521000.0012024000.000.000.00
ImageDifference-30695400.00-1932194700.001899504900.0063385200.000.000.00
Reference
InitialStateImage:
E:
\学习资料\ENVI\课程设计\2006_jdm.img
Dimensions:
(1:
7726,1:
2977),Band1,PixelSize:
(30.00,30.00)Meters
FinalStateImage:
E:
\学习资料\ENVI\课程设计\2009_jdm.img
Dimensions:
(1:
7726,1:
2977),Band1,PixelSize:
(30.00,30.00)Meters
EquivalentClassPairings
Unclassified<==>Unclassified
植被[Green]3550points<==>植被[Green]15849points
建筑[Red]15187points<==>建筑[Red]53838points
水系[Blue]4773points<==>水系[Blue]4817points
ReportProducedon:
SunMay1016:
16:
002015
结果与分析
一、几何精纠正控制点位置及其误差
地面收集的控制点见下表
其中控制点的误差RMS在1左右不超过1.5
二、影像处理结果分析
1、图像裁剪前后影像比较
(1)图像裁剪前的多光谱影像
(2)图像裁剪后的影像
【结果与分析】按照行政区划边界或自然区划边界进行图像裁剪,通过对图像的裁剪,来进行特定区域的地物问题分析
2、长沙两年的遥感影像比较
(1)2006年的多光谱影像
(2)2009年的影像
【结果与分析】从上两图明显可以看出来的是眼和湘江流域,整个居民区在不断扩张,人类的活动影响十分明显。
3、长沙两年分类后的遥感影像比较
(1)2006年的多光谱影像
(2)2009年的多光谱影像
【结果与分析】从上两图清楚的看出长沙是的植被减少了很多,取而代之的是连绵的居民区特别是眼水系的区域,居民区以水系为中心在不断向外扩张。
4、长沙两年分类处理后的遥感影像比较
(1)2006年的多光谱影像
(1)2006年的多光谱影像
【结果与分析】右边的平原地区居民建筑明显增加,右方高山去的居民建筑变化不大,而沿水系的居民建筑则增加迅速,这是由于城市化进程带来的城市的发展。
5、动态监测报告
【结果与分析】从上表中的数据可以清楚的看出长沙是的植被减少了很多,取而代之的是连绵的居民区,水系也在缩减。
很明显人类的活动下,城市化的扩张,都不可避免的占用植被来满足人类生活的需要。
存在问题与解决办法
1、部分遥感影像的云量较多,怎么减少其对影像结果的影响?
分析:
遥感影像的云量可以从数据源上进行解决,当遥感影像的云量较多难以避免时,通过计算ndvi值,分析云值范围,对该段数据进行忽略,尽量减少云量对实验结果的影响。
2、导入ENVI进行查看,发现未拼接重投影之前影像数据范围为-3000—0.99880,拼接重投影之后影像数据范围在-3000—9988是数据值发生了改变吗?
分析:
ENVI—BasicTool—Preprocessing—Data-SpecificUtilities—
ViewHDFAttribute在envi里面查看NDVI波段参数可以发现它的有效值范围为-2000—10000,因此数据值-3000是无效值,其他的值是乘以了10000这个系数scale_factor,因此在拼接重投影过程中ndvi值并没有发生变化。
还有些值如32767等表示有云。
3、按上面所述,它的根本值没有发生改变,但是毕竟一部分数值已经扩大了,怎么处理比较合理?
分析ENVI—BasicTools—BandMath波段运算,把负值去掉,在弹出的对话框BandMath中Enteranexpression(b1lt0)*0+(b1ge0)*(b1*0.0001)。
这个公式意思就是要是值小于0就乘以0使其变为0同时值大于等于0的话就乘以0.0001这个系数。
这样就OK了。
关于B1代表所选的那个影像。
总结
通过此次实验,能够对遥感图像的处理方法有更加深刻的了解,为以后的学习打好了基础,同时使我对ENVI软件操作有了更熟练的应用!
本次实习分别通过长沙这个实例实现了变化监测的相关内容的学习和操作,变化监测一般分为三个大的步骤:
数据与处理,变化信息检测、变化信息提取,通过实验操作之后,对变化检测的基本原理有了深入的理解,对于影像直接比较法和分类后比较法变化检测的基本方法和步骤有了熟练的掌握。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 遥感 图像 处理 课程设计 长沙市 植被 覆盖率 变化