遥感原理与应用课程设计.docx
《遥感原理与应用课程设计.docx》由会员分享，可在线阅读，更多相关《遥感原理与应用课程设计.docx（16页珍藏版）》请在冰豆网上搜索。
遥感原理与应用课程设计
遥感原理与应用课程设计
—“遥感专题信息提取与专题图制作”设计报告
学院
遥感信息工程学院
班级
09031
学号
2009302590123
姓名
吴煜晖
日期
2012-2-13日至2月19日
指导教师
宋蔚林
一、课程设计的目的和意义
这门课旨在加深理解和巩固理论课上所学的有关遥感的基本原理、遥感传感器的成像机理、遥感图像的处理方法、专题信息提取以及遥感综合应用技术；锻炼熟练运用遥感软件独立分析问题、解决具体问题的实际工作能力；可以培养同学们良好的工作习惯和科学素养，为今后参加科学研究工作以及毕业设计打下良好的基础。
二、课程设计方案和流程
1、实验数据介绍
1）地点：
宜昌地区；
2）传感器类型：
TM1、3、4、5、7波段影像，空间分辨率：
30米；TM6波段影像，空间分辨率：
120米；SPOT影像，空间分辨率：
10米；
3）格式：
TIF；
4）说明文件：
Readme.txt，SPOT影像投影信息、大地坐标、分辨率。
2、硬软件介绍
本次实验使用软件ERDAS。
ERDASIMAGINE是美国ERDAS公司开发的专业遥感图像处理与地理信息系统软件。
它以模块化的方式提供给用户，使得用户可根据自己的应用程序的要求、资金情况合理地选择不同功能模块及其不同组合，对系统进行剪裁，充分利用软硬件资源，并最大限度地满足用户的专业应用要求。
在图像处理方面，直观的操作步骤使用户操作起来非常灵活方便。
本次实验主要应用到ERDAS软件主菜单中Viewer、Import、DataPrep、Interpreter、Classifier及Composer等模块。
3、总体方案
首先利用ERDAS软件将宜昌地区左右TM影像的1、3、4、5、6、7波段TIF格式单波段影像分别合成left.img和right.img多波段TM影像（不少同学先将TIF格式转换为img格式，再分别合成left.img和right.img多波段TM影像，其实这是没必要的）。
再根据宜昌地区的SPOT影像sp_yc.img以及说明文件Readme.txt，分别对left.img和right.img进行几何纠正，并将纠正后的两幅影像进行拼接，再将拼接后的TM影像进行裁剪。
最后在进行适当的图像增强处理后，对拼接后的影像进行监督分类（可以将SPOT影像和TM影像进行融合再进行分类以获得更好的分类效果），提取专题信息，制作成专题地图。
4、流程框图
三、课程设计的过程和步骤
1.图像预处理
首先需要将实习中给的tif格式的各幅图像输入转入img图像，可以通过打开tif格式的图像另存为img格式转换，或者用Import实现。
2.多波段合成
将对应的六幅影像分别作为六个波段叠加，形成宜昌地区具有六个波段的左右两幅TM影像，如图1所示。
图1（a）图1（b）
3.几何纠正
通过提供的SPOT影像对左右两幅TM影像进行几何纠正处理，纠正后的图像如图2所示。
图2（a）图2（b）
4.图像拼接
将几何纠正完后的左右两幅TM影像进行拼接处理，拼接后图像如图3所示。
图3
5.图像剪裁
将TM影像中用来制作专题图的区域剪裁出来，剪裁后图像如图4所示。
图4
6.图像分类
采取感兴趣的区域作为样本对剪裁后的TM影像进行监督分类处理。
监督分类分为5类，分别代表水体、深绿色植被、浅绿色植被、道路及裸地和居民地。
结果如图5所示。
图5
得到的结果如下：
混淆矩阵：
ERRORMATRIX
-------------
ReferenceData
--------------
Classified
DataBaldandRYoungPlanOldPlantHabitat
--------------------------------------------------
BaldandR1921012
YoungPlan4471140
OldPlant4123071
Habitat100196
Water0000
ColumnTotal201493322199
ReferenceData
--------------
Classified
DataWaterRowTotal
------------------------------
BaldandR0205
YoungPlan6495
OldPlant4328
Habitat14211
Water25732573
ColumnTotal25973812
CLASSIFICATIONACCURACYASSESSMENTREPORT
-----------------------------------------
ImageFile:
c:
/documentsandsettings/wuyuhui/桌面/2009302590123吴煜晖/2009302590123吴煜晖/wyh2009302590123/supervisedmerged.img
UserName:
WuYuhui
Date:
SunFeb2616:
37:
032012
ERRORMATRIX
-------------
ReferenceData
--------------
ClassifiedDataUnclassifiWaterHabitatOldPlant
-------------------------------------------------------
Unclassified0000
Water03100
Habitat00311
OldPlant00434
BaldandRoad0011
YoungPlant0133
ColumnTotal0323937
ReferenceData
--------------
ClassifiedDataBaldandRYoungPlanRowTotal
---------------------------------------------
Unclassified000
Water0132
Habitat0032
OldPlant0038
BaldandRoad47352
YoungPlant691102
ColumnTotal5395256
-----EndofErrorMatrix-----
ACCURACYTOTALS
----------------
ClassReferenceClassifiedNumberProducersUsers
NameTotalsTotalsCorrectAccuracyAccuracy
---------------------------------------------------
Unclassified000------
Water32323196.88%96.88%
Habitat39323179.49%96.88%
OldPlant37383491.89%89.47%
BaldandRoad53524788.68%90.38%
YoungPlant951029195.79%89.22%
Totals256256234
OverallClassificationAccuracy=91.41%
-----EndofAccuracyTotals-----
KAPPA（K^）STATISTICS
---------------------
OverallKappaStatistics=0.8860
ConditionalKappaforeachCategory.
------------------------------------
ClassNameKappa
---------------
Unclassified0.0000
Water0.9643
Habitat0.9631
OldPlant0.8770
BaldandRoad0.8787
YoungPlant0.8285
-----EndofKappaStatistics-----
7.制作专题图
用ERDAS的MapComposer工具导入裁剪好的分类图像，制作专题地图，并添加图例，比例尺和指北针等。
详见文末附录。
8.融合前和融合后图像分类精度比较*
融合后影像如图6所示，对其进行监督分类后所得结果如下图7所示：
图6图7
混淆矩阵如下：
ERRORMATRIX
-------------
ReferenceData
--------------
Classified
DataWaterOldPlantYoungPlanHabitat
--------------------------------------------------
Water621988253
OldPlant17561321294
YoungPlan42313640
Habitat10283024259
RoadandB603361337
ColumnTotal640086202136824943
ReferenceData
--------------
Classified
DataRoadandBRowTotal
------------------------------
Water862269
OldPlant36605
YoungPlan41395
Habitat325293
RoadandB7861763
ColumnTotal80497325
-----EndofErrorMatrix-----
精度评价如下：
CLASSIFICATIONACCURACYASSESSMENTREPORT
-----------------------------------------
ImageFile:
c:
/documentsandsettings/wuyuhui/桌面/2009302590123吴煜晖/2009302590123吴煜晖/wyh2009302590123/resolutionsuperised.img
UserName:
WuYuhui
Date:
SunFeb2617:
11:
562012
ERRORMATRIX
-------------
ReferenceData
--------------
ClassifiedDataUnclassifiWaterOldPlantYoungPlan
-------------------------------------------------------
Unclassified0000
Water02521
OldPlant02605
YoungPlant00010
Habitat0121
RoadandBald0021
ColumnTotal0286618
ReferenceData
--------------
ClassifiedDataHabitatRoadandBRowTotal
---------------------------------------------
Unclassified000
Water0028
OldPlant0168
YoungPlant0111
Habitat38143
RoadandBald499106
ColumnTotal42102256
-----EndofErrorMatrix-----
ACCURACYTOTALS
----------------
ClassReferenceClassifiedNumberProducersUsers
NameTotalsTotalsCorrectAccuracyAccuracy
---------------------------------------------------
Unclassified000------
Water28282589.29%89.29%
OldPlant66686090.91%88.24%
YoungPlant18111055.56%90.91%
Habitat42433890.48%88.37%
RoadandBald1021069997.06%93.40%
Totals256256232
OverallClassificationAccuracy=90.63%
-----EndofAccuracyTotals-----
KAPPA（K^）STATISTICS
---------------------
OverallKappaStatistics=0.8705
ConditionalKappaforeachCategory.
------------------------------------
ClassNameKappa
---------------
Unclassified0.0000
Water0.8797
OldPlant0.8415
YoungPlant0.9022
Habitat0.8609
RoadandBald0.8902
-----EndofKappaStatistics-----
_______________________________________________________________________________
四、课程设计的结果分析与评价
1.几何纠正
本次实习，我采用了一次多项式纠正公式，至少需要三个控制点。
为了更好地防止粗差的影响，我两次纠正都选取了7个控制点，控制点误差均在一个像素以内。
几何纠正结果精度较高，详见图3。
影响几何纠正精度的主要因素在于控制点的选取。
控制点的准确程度决定了几何纠正的精度。
所以要提高几何纠正的精度，就应该选取易于辨识的控制点，比如尖角、颜色对比度较高的点等。
2.影像拼接
我选取的影像，由于拼接线选的较好，因此整体影像并无拼接痕迹。
ERDAS只能采用基于地理坐标的影像拼接，所以做好影像拼接的前提是做好几何纠正，否则拼接出来的影像在位置或方位上将有明显的不衔接。
3、影像分类
在制作专题图的过程中，我采用了监督分类的方法，详见我在上文所附的图。
在监督分类中，要获得比较高的分类精度，必须先获得比较纯的训练样本。
检查训练样本的方法是查看训练样本的混淆矩阵，混淆矩阵的整体精度应该保证至少在90%以上。
据理论说，融合后的图像分类精度比融合前的分类精度更高，所以可以采用融合更高分辨率的影像来提高分辨率。
而我的结果却恰好相反，融合后的图像分类精度略低于融合前的。
这可能是由于我在处理融合后的图像时不如之前那么细致。
五、课程设计的总结与体会
本次课程设计从方案设计到最后专题图的输出都由我们自己完成，很大程度上地锻炼了我们实践的能力。
当老师第一天让我们独立完成遥感专题信息提取与专题图制作的总体设计方案和流程图时，我完全不知所措，但还是通过自己的学习和与同学的沟通逐渐的了解下来。
第一周的Erdas实习，我对应用型软件并不是十分感兴趣。
但在通过学习的过程中，我对Erdas的优点有了进一步的了解，这更加激励了我对亲自实现其功能的想法。
通过本次课程设计，我更好地认识了遥感，加深理解了理论课上所学的有关遥感的基本原理、遥感传感器的成像机理、遥感图像的处理方法、专题信息提取以及遥感综合应用技术；锻炼了熟练运用遥感软件独立分析问题、解决具体问题的实际工作能力；培养了我们良好的工作习惯和科学素养，为今后参加科学研究工作以及毕业设计打下了良好的基础。
最后，还是感谢这一周来老师的指导，当然我希望老师能将一些软件使用的小窍门，使我们不再一些小Bug中翻船。
附：
宜昌地区水体分布专题图
课程设计报告成绩
项目
得分
文字正确，语句通顺（10分）
内容充实、详略得当（10分）
方案正确、流程合理（10分）
过程清楚、步骤详细、实验结果正确（30分）
对实验结果分析恰当、能找出影响效果的因素和找出改进精度的方法，并叙述正确（30分）
对课程设计有一定的总结与体会（10分）
发现有拷贝、雷同、抄袭（视情况扣大于30分）
总成绩
升级会员 