对地遥感图像地图投影程序的设计与开发-毕业论.docx

济南大学毕业设计毕业设计题 目 对地遥感图像地图投影程序 的设计与开发 学 院 信息科学与工程学院 专 业 通信工程 班 级 通信1203 学 生 毕* 学 号 20121221* 指导教师 刘* 二一六 年 六 月 六 日- 2 -毕业设计摘 要地图投影技术是将地球椭球面投影到地图平面的技术，高斯-克吕格地图投影方法是其中使用最广泛的一种投影方法，也是我国地图制图所使用的投影方法。因此，研究该投影技术具有实际意义。本设计的主要内容包括以下几个方面：首先理解地图投影的原理、分类方法及各类方法的适用条件，其次，掌握高斯克吕格地图投影的理论基础，研究该投影方法的电算公式，推导基于各种地球椭球模型的高斯克吕格地图投影的电算公式，利用程序开发语言编写电算公式的实现程序并选择MATLAB软件开发界面设计开发一个电算公式换算界面，本文细致描述了基于MATLAB GUI 图形界面开发技术的电算公式计算软件的制作过程，并进行了仿真测试，无论在理论还是结果上都取得了满意的结果。关键词：地图投影；遥感图像；高斯-克吕格；MATLAB图形界面ABSTRACTMap projection technology is the earth ellipsoid which is used to map projection plane technology,Gauss - Kruger projection method is using a projection method and it is the most widely, projection method which also as Chinas cartography are used.So,the study of this projection technology has practical significance The main contents of this design include the following: first understand the principle of map projection, classification method and the applicable conditions of the method, secondly, grasp the theoretical foundation of Gauss Kruger projection on the projection method, the formula and the inversion formula, calculation formula is derived based on the Gauss earth ellipsoid model Kriging map projection, using programming language program and select the MATLAB software interface design and development of a computing formula calculation formula of interface, make a detailed description of the process of making MATLAB GUI graphical interface, and use the MATLAB software to make the simulation of success, satisfactory results are obtained. Whether in theory or in the. Key words：Map projection；Remote sensing image；Gauss - Krueger；MATLAB GUI目 录摘要. .IABSTRACT.II1 前言.11.1选题背景与意义.11.2 研究现状.11.3 设计方案.22 地图投影.32.1 地图投影概述.32.2 地图投影的分类.32.3 高斯克吕格正反算公式及其推导.52.4地球椭球模型.83 电算公式软件开发. 10 3.1 MATLAB简介.103.2软件的界面设计.10 3.2.1 GUI的启动.10 3.2.2界面的设计.12 3.2.3 各控件属性设置.174 软件测试.224.1电算公式界面开发结果.224.2 地图投影结果. 24结论.25参考文献.26致谢.27附录. . . . .28- 31 -1 前言1.1 选题背景与意义地图投影技术是将地球椭球面投影到地图平面的技术，高斯-克吕格地图投影方法是其中使用最广泛的一种投影方法，也是我国地图制图所使用的投影方法。随着遥感技术的发展，遥感图像的地理编码也必须使用地图投影技术，因此，研究高斯克吕格投影方法对遥感图像的应用具有实际意义。1.2 研究现状目前，国内外众多专家学者对这项技术进行了深入的研究，例如徐州师范大学的周兴东在文献一中利用遥感图像技术对我国土地的分类和利用状况进行了研究，对遥感图像中的一些值比如光谱特征进行了分析，这样就可以识别用地信息，这类方法可以在农村及城市的土地问题研究中发挥应有的作用，具有很强的现实指导意义。重庆大学通信工程学院马大玮在文献二中介绍了他所设计的遥感图像传输系统，这个系统主要是将压缩电路进行了规范化的设计处理，利用此系统，遥感飞机和地面以及卫星之间就可以实现遥感图像的传输，这类技术在军用和民用领域得到了广泛的认可和重视。湖南科技大学林剑在文献三的光谱学与光谱分析中利用对地遥感图像技术对我国土地覆盖多光谱特征进行了研究分析，过程中使用到了空间距离的分析方法，从而可以准确的反映研究区域土地的多光谱特征，对土地分析有着重要的作用。北京大学遥感与地理信息系统研究所的陈超在文献四中利用高分辨率的遥感图像分析光谱特征进而进行水上桥梁图像的提取，实验表明，这种方法能够适应不同的复杂环境，克服周围复杂状况的影响，并且 能够达到一个比较高的精度。西安电子科技大学的刘德连，在文献中对大气的遥感图像和大气效应的校正原理进行了分析，使用亮度聚类探测算法作为基础，在此之上经过深入研究得到了纹理分割的目标探测算法，可以实现复杂地面背景中目标的探测。中南大学的高光明教授利用遥感图像进而对地壳运动进行分析，在对我国许多地区的图像进行分析后，使用对地遥感技术总结出了地壳运动对地震的影响，并且总结出一系列的方法，可以在以后我国及世界地震的预测方面提供所需的帮助，为人类应对自然灾害从而减小损失做出应有的贡献。内蒙古大学的刘鹏涛使用遥感技术对中国北方地区的植被进行了研究，其中特别研究了草地这种植被，利用高分辨率的遥感图像，可以实现区域草地覆盖率的分析和研究工作，利用这种分析方法，可以全面的分析我国的植被覆盖情况，分析我国的植被所存在的问题，进而对我国植被的保护和绿化工作做出贡献。北京航空航天大学教授李树楷利用GPS系统研究了对地的动态定位工作，在航空航天学中，对地面目标的准确定位一直是一项重要而且非常关键性的一个步骤，它是利用一种投影变换的公式，将测量到的目标点的数据投影到已建好的坐标系中，完成坐标的转换，进而实现对地面的动态定位工作。北京化工大学的朱炜，深入研究了对地遥感图像中提取各种图像特征的方法，这种方法可以从已有的遥感图像中按需并且准确的提取出人们所需要的各种要素，这样不仅使最终的结果变得精确，而且使图像特征提取的整个过程和方法变得简便易行。1.3 设计方案本设计的主要内容包括以下几个方面：首先应该了解地图投影对于遥感图像地理编码的重要作用，理解地图投影的原理、分类方法及各类方法的适用条件，其次，掌握高斯克吕格地图投影的理论基础，研究该投影方法的电算公式，推导基于各种地球椭球模型的高斯克吕格地图投影的电算公式，然后，利用程序开发语言编写电算公式的实现程序并选择一种界面开发软件设计开发一个电算公式换算界面，最后对所开发的程序和界面进行测试和验证。2 地图投影2.1 地图投影概述所谓的地图投影，简单地说就是将椭球模型上的距离，方位和坐标等元素，按照一定的数学规则将它投影到平面上，这种专门研究地图投影问题的学科称之为地图投影学。这里所用的数学法则，总的来说可以用以下的两个方程式加以概括： x=F1(L,B)y=F2(L,B) (2.1)该式中，L和B是椭球模型上的大地坐标，x和y则是投影平面上的直角坐标。该式表示了椭球面上一点通过特定的投影方式将其投影到相应的直角坐标系上的坐标换算关系式，F1和F2叫做投影函数。由此我们可以看出，地图投影的主要研究内容是探讨研究者所用的研究方法与所用的椭球模型上的元素和投影平面坐标相应元素之间关系的解析式。投影的方法以及种类有很多，不同的方法其投影的基本原理是相同的，只是他们的投影函数即F1和F2不相同，下面我们将介绍不同的地图投影方法。2.2 地图投影的分类地图投影的分类方法有很多种，这里只对其中常用的几种作简要介绍。1.按投影变形的性质分类1）等角投影等角投影，顾名思义，这种投影方法就是能够保证投影前后图像的角度不变形，即，在进行等角投影的过程中，微分圆的投影结果仍是微分圆，投影前和投影后具有微小圆形上的相似性，投影的长度之比和它的方向没有关系，即投影图像的某点的长度比是一个不变的定值。正因为如此，也把等角投影称作正形投影。2）等积投影由名字可知，这种投影就是要求在投影前后的图像面积保持不变。3）任意投影这种类型的投影是最后的结果是既不等角，也不等积，这种投影方法比较多，在实际中应用也比较广泛。其中，在某一个方向上保持它的长度比是1，这种就称其为等距离投影。2.按经纬网的投影形状分类在这个分类方法中，是按照正轴投影图像的经纬网形状来划分，采用投影面的名称来命名。1）方位投影取一个与椭球极点相切的平面，将椭球极点附近范围内的区域投影到这个与极点相切的平面上，纬线在进行投影之后是圆形，它们都以极点为圆心，经线在投影之后是它的向径，并且经线的交角不变。方位投影的直角坐标公式为 (2.2)式中：是球面极坐标。 方位投影的变形公式为垂直圈长度比 (2.3)等高圈长度比 (2.4)面积比 (2.5)最大角度变形 (2.6)由上面的公式可以看出不同方位投影之间的区别。2）圆锥投影取一个圆锥面，使其和椭球的某条纬线相切，将纬圈附近的这片区域投影到圆锥面上，然后再将投影后的圆锥面沿着它的某一条经线剪开，使之成为平面。在这个投影中，纬线在投影后会成为同心圆，经线投影后则成为这些同心圆的半径，并且经差与经线交角的比例是一定的。圆锥投影的直角坐标公式为 (2.7)式中：圆锥投影的变形公式为 (2.8)在圆锥投影中最小纬线的长度比为 (2.9)而最小纬线的纬线投影半径为 (2.10)3）圆柱投影取一个与椭球赤道相切的圆柱，将赤道附近的这片区域投影到与赤道相切的圆柱面上，之后将这个圆柱面的沿某条线剪开使之成为平面，在这个投影中，纬线在投影后成为一组相对于赤道对称的平行线；经线在投影后也会成为一组平行线，且与纬线投影后所成的平行线是垂直的。圆柱投影的直角坐标公式为 (2.11)圆柱投影的变形公式为 (2.12) (2.13) 2.3高斯-克吕格投影正反算公式及其推导图2.1高斯投影假想图高斯克吕格(Gauss-Kruger)投影是等角横切椭圆柱投影。如上图所示，假想在地球椭球体外面横套一个椭圆柱，并且使这个椭圆柱与某一条子午线相切，这条子午线就称之为中央子午线或轴子午线，椭球体的中心在椭圆柱的中心线上，然后使用高斯克吕格投影的方法，以中央子午线为中心，将两侧一定经差范围内的部分区域投影到外套的椭圆柱面上，再将这个圆柱面展开就成为了我们想要的投影面。 我国规定按照经差6和3两种规格进行投影分带,在为大比例尺进行测图和进行工程测量时通常采用3带投影，在一些特定的情况下，工程测量时也会采用1.5带或者任意带。一般情况下，为了测量以及实验结果的通用性，需要和国家的6带和3带联系起来。高斯投影6带，从0度子午线开始算起，每隔6经差自西向东分带，1，2，3依次编号，我国使用投影6带的高斯投影的中央子午线经度，由69开始每隔6至135，一共有12带，带号使用n表示，L0用来表示中央子午线的经度，两者的关系为 L0=6n-3。如图所示高斯投影3带，3带投影是以6带为基础形成的，它的中央子午线的其中一部分带（其中的单数带）和6带的中央子午线相重合，而另外一部分带（其中的偶数带）和6带的分界子午线相重合。3带的带号用n表示，3带的中央子午线经度用L表示，两者之间的关系是L=3n，如图所示图2.2 6度带与3度带在这个投影面上，赤道和中央子午线在投影之后成为直线，坐标原点O是中央子午线和赤道的交点，纵坐标轴是中央子午线的投影，横坐标轴是赤道的投影，这样就形成了高斯平面直角坐标系。高斯投影是正形投影，这样就可以保证投影角度前后保持不变，前后图形保持较高的相似度，任意一点各方向上的长度比具有同一性。高斯投影使用了相同法则的分带投影，这样不仅限制了长度的变形，还可以使得在计算变形引起的各项改正数时公式和表达式尽可能简洁。种种优点使得高斯投影得到了广泛的推广。高斯-克吕格投影需要满足下面三个条件：1、经纬线投影后是相对于中央经线对称的曲线；2、投影前后没有角度变形；3、投影前后中央经线的长度保持不变。高斯投影的公式： 已知椭球面上一点的大地坐标是(L,B),其中，L表示经度，B表示纬度，在进行高斯投影的时候，为了减小长度上的变形，投影是在分带内进行的，这里的分带就是我们前面说的6度带或3度带，所有投影都是在每一个分带内相对于此分带的中央子午线和赤道线进行的，其中每个分带的中央子午线叫做投影带的中央子午线，如果它的经度是，那么在正算公式计算过程中使用的经度值是相对于的经度差，即 （2.14)如果以6度带为例，那么它的计算方法是 (2.15)上式中的意思是取整，即取出其中的整数部分，为了方便计算，以上公式中使用的角度全部是以弧度为单位，因此6度角和半个6度角需要先转成弧度。 由于是在不同的分带内进行投影，如果把最后的计算结果放到统一的高斯直角坐标系中，则需要将计算出的结果加上相应投影带中央子午线的高斯直角坐标的值，即 (2.16)式中为地球椭球模型的长半轴。考虑到上面所讲的转换关系，那么由和计算统一高斯直角坐标的公式是： (2.17) (2.18)这里要说明的是，对应的是南北方向，对应的是东西方向，而且和的单位都使用弧度数。还有，计算式中的是赤道到纬度值处的子午线的弧长，是卯酉圈曲率半径，两者的计算公式会在之后加以说明，另外两个参数和与地球参数和纬度的关系式为： (2.19) (2.20)式中为地球第二偏心率，可以通过地球椭球模型半长轴和扁率计算出来（具体的关系在后面的地球椭球模型上加以介绍），即, (2.21)其中为地球椭球模型的短半轴。 X是赤道至纬度处的子午线弧长，其计算公式是： (2.22)式中， (2.23) 并且 (2.24) (2.25) (2.26) (2.27) (2.28) (2.29) 卯酉圈曲率半径的计算公式为： (2.30) 2.4地球椭球模型表2.1各地区及其对应的椭球模型适用区域 椭球名称 年代 长半轴 扁率北美 克拉克 1866 6378206 1:294.98南美及欧洲 海福特 1910 388 297.00前苏联 克拉索夫斯基 1940 245 298.30印度 埃弗瑞斯 1830 6377276 300.80日本 白赛尔 1841 397 299.15澳大利亚 1967年参考系 1967 6378160 298.247 克拉克 1880 249 293.47非洲 第16届国际推荐值 1975 140 298.257 第17届国际推荐值 1980 137 298.257从上表可以看出，椭球模型一般是给定它所对应的长半轴值和扁率的值，此外，在计算过程中还可能会使用到其他的一些中间量，均可由这两个计算出来：第一偏心率 (2.31) 第二偏心率 (2.32) 短半轴 (2.33) 长半轴 (2.34) 扁率 (2.35) 3投影电算公式软件开发3.1MATLAB简介MATLAB是美国mathworks公司研发出来的用于商业用途的数学软件，可以用于很多比较高级的计算语言和不同的交互式环境，例如：算法开发，数据可视化和各种数据分析，这个软件中主要包含MATLAB和simulink这两个部分。MATLAB的全称是matrix&laboratory，是矩阵实验室的意思。它将许多强大的功能（例如：数值的分析计算，矩阵的计算，动态数据系统的建模以及仿真等）聚集在一个方便用户使用的视窗环境中，可极大的方便科学研究工作者在科研，系统设计和一些科学实验数据的计算，为众多科技领域的发展做出了应有的贡献。MATLAB 将计算、可视化和编程集于一身。在MATLAB中，问题的提出，计算的过程和最终结果的呈现都使用用人们比较习惯的数学方法来描述，免去了运用传统编程语言的处理环节。这一系列特点让MATLAB成为研究学者们进行数学分析、算法开发以及所需应用程序开发的相对良好的环境。它是Math Works产品家族中一系列衍生产品的基础。在通信技术领域，MATLAB也具有明显的优势，因为系统性能问题是通信领域经常研究的一个问题，构建实验系统进行测量的这种传统方法需要大量的资金花费以及很长的构建周期，实现起来比较麻烦，而使用MATLAB进行系统性能的测试，在计算机上就能够实现一系列操作，避免了系统构建过程中的麻烦，只要在MATLAB中键入不同的参数数据就可以得到不同情况下系统的性能，种种优势使得MATLAB在通信仿真方面得到了广泛的应用。3.2软件的界面设计3.2.1 GUI的启动MATLAB GUI，即MATLAB图形用户界面，其中GUI是Graphic User Interface的简写,它经常用于通信系统的模拟，使编程更加人性化，易于用户操作，图形用户界面，它不需要写命令，对初学者来说相对容易。图形用户界面在人与计算机的交互操作中占有重要的地位，它的界面既生动形象，又能够让用户方便快捷的操作。图形用户界面被用于当前绝大多数开发环境和应用程序，也可以用于比较流行的开发工具。图形用户界面有以下的几个部分构成：布局区域（Layout Area），对象对齐工具（Align Object），顺序切换编辑器（Tab Order Editor），菜单编辑器（Menu Editor），M文件编辑器（M-file Editor），属性检查器（Property Inspector），激活按钮（Run），图形对象浏览器（Object Brower），重复或者取消上次操作的按钮（Undo和Redo），组件面板（Component Pallette）。图3.1 新建GUI界面在MATLAB的主窗口中，单击File菜单，之后点击其中的New选项，再选择其中的GUI选项，系统就会弹出图形用户界面的设计窗口。图3.2 新建GUI的编辑界面MATLAB中的GUI设计一共有四种模板，它们分别是Blank GUI（这是系统默认的选项），GUI with Uicontrols（这是带有控件对象的模板），GUI with Axes and Menu（这是附带有坐标系和菜单项的GUI模板）和Modal Question Dialog（这是带有模式问答对话框的模板），当用户单击模板对应的按钮时，在GUI的模板设计界面的右侧会出现用户所选择的相应的GUI模板。在GUI的设计模板类型中选择其中一个模板类型，点击OK按钮，系统就会显示出GUI的设计窗口。图3.3 GUI编辑界面在GUI的布局区中，为拖拽任何控件时，其下方的状态栏会显示当前当前鼠标所在点的窗口所对应的坐标（单位为像素）以及当前GUI对象的位置和大小（单位为像素）。3.2.2界面的设计界面布局设计的步骤主要包括以下几个：1）通过鼠标的拖拽将控件面板中所需的控件放到界面设计区域，2）使用对象对齐工具（Align Objects）对所拖拽控件的布局进行调整，3）在使用Tab顺序编辑器（Tab Order Editor）对界面设计区的各个控件的Tab顺序进行相应的调整设置，4）如果用户所设计的界面需要菜单，还需要使用菜单编辑器（Menu Editor）对界面进行菜单的设计，5）添加完所需的控件之后，用户可以在对象浏览器（Object Browser）中看到所有刚添加的图形对象，到此，就完成了整个界面的布局设计。所设计的界面中，所拖拽的各控件及其位置如图3.4所示，界面上还包括高斯克吕格投影的电算公式以及在电算公式计算过程中所用到的部分参数的值图3.4 所设计的界面 图3.5所示为MATLAB的图像对齐工具，可以利用此工具对界面上的控件进行比较整齐的排列。图3.5 图像对齐工具下面图3.6和图3.7所示是MATLAB中切换顺序编辑器，此图显示了界面中所拖拽控件的前后顺序，利用它，可以设置用户计算机键盘上的Tab键时界面上控件被选中的先后顺序，在Tools菜单下单击Tab Order Editor项，即可以打开下图所示的顺序编辑器。图3.6 顺序编辑器下图为部分编辑框所对应的顺序编辑器界面图3.7 顺序编辑器使用菜单编辑器对界面上所需的菜单进行编辑，如下图所示，编辑完成之后，它就会显示在所设计界面的上方。图3.8 菜单编辑器对象浏览器如图3.9和图3.10所示，打开它之后，就可以看到所设计界面上所有控件均显示在其中。图3.9 对象浏览器图3.10 对象浏览器如下图3.11所示为GUI选项界面图3.11 GUI选项3.2.3各控件属性设置每个控件在最开始不可能是都符合所设计界面的要求，在设置和查看界面中各控件的属性时，都需要用到属性查看器（Property Inspector），属性查看器里面显示的列表中包括用户所能够设置的控件属性以及控件属性的当前值，双击布局编辑器的其中一个控件或右键选择Property Inspector即可打开属性检查器，如图所示图3.12 属性查看器打开属性检查器后，就要对每个控件的属性进行设置，其中最重要的是String和Tag属性，String属性就是界面中的标签文字，Tag属性是图形对象的标识，此属性是全部图形对象的最重要属性，对象的回调函数名字是自动以Tag属性命名的。Tag属性的主要用途是来构造相应控件的响应函数的名字还有 handles结构数组的名字，由此可得，在用户设置Tag属性时，必须为其选择一个有效的变量名。此外，还有Style属性，这是控件的类型；Visible属性，这是对控件的可见性进行设置，可见on/不可见off；Enable属性，有效on/无效off/不活动inactive；Fontsize属性。在为一些比较特别的控件设置属性时，需要注意一些特定属性的设置，例如：在对可编辑文本框设置属性时，必须要对Max和Min的值进行设置，在未设置时，系统默认的值分别是1和0；当Max-Min的值大于1时，此控件的编辑文本框可允许多行输入，本设计在控件属性设置时就要对Max的值进行修改，将其改为2.0，这样就可以进行多行输入；当Max-Min的值小于等于1时，编辑文本框只能允许单行输入。图3.13 下拉菜单属性设置单击按钮即出现多行输入对话框，如图3.14所示，所输入的内容就会显示在下拉框中。图3.14 多行输入框图3.15 按钮属性设置在对界面上每个控件的属性设置完成后，这时就可以生成fig文件和m文件，在tools选项中单击GUI Options，系统就会出现GUI选项这个对话框，在这个对话框中选择”Generate FIG-file and M-file”,如图所示，这就表明GUIDE要建立fig文件和m文件。图3.16 GUI选项界面点击OK之后，GUI界面正式建立完成，并且生成了对应的m文件，m文件的界面如图所示图3.17 m文件界面在打开m文件的界面上点击按钮，就会出现GUI界面上每个控件所对应的函数如图所示图3.18回调函数图3.19回调函数单击某一个回调函数，界面就会自动跳转到相应的回调函数的位置，然后将所要实现功能所对应的函数写入callback函数后面，到此就完成了界面的设计任务。4软件测试4.1电算公式界面开发结果按照上述过程所完成的界面设计，在点击按钮之后会出现如下结果，由于种类较多，这里只选择其中的几个进行展示。界面上方为地图区域，其下方为这个区域所对应的模型，生成的年代，这个模型的长半轴和扁率；再往下则是显示了高斯克吕格投影的电算公式，电算公式下方为在投影计算过程中所使用到的参数；界面的最下方是三个按钮，其一为“测试数据”，单击之后会显示遥感原图，其二为“结果图形”，单击之后会显示该地区遥感原图在进行高斯克吕格投影之后的投影图像，其三为“关闭”，单击之后就会关闭窗口。图4.1 北美地区参数南美及欧洲地区对应椭球模型及计算过程中各参数如图4.2所示图4.2 南美及欧洲参数澳大利亚地区对应的模型参数及计算过程中的参数如图4.3所示图4.3 澳大利亚地区参数4.2 地图投影结果单击界面上的“测试数据”按钮后，系统运行相应程序，之后会弹出显示框，即显示北美地区的卫星遥感图像原图，如图4.4所示图4.4 遥感原图在单击“结果图形”按钮后，系统运行相应的投影程序，之后会弹出显示框，即显示北美地区在经过高斯克吕格投影之后的投影图形，如图4.5所示图4.5 投影图像结 论通过本次毕业设计，让我有机会能够将大学本科期间所学到的专业知识进行了比较系统的总结，并且进行了一系列实践。从最初论文研究方向的确定一直到最后运用MATLAB程序实现该课题的最初要求，都是在老师的指导以后，独立思考认真完成的。本篇论文所研究的课题是我之前从未接触过的领域，这个领域所涉及的多方面的专业知识也是我一无所知的，所以，我在本校图书馆和网络上查阅了这方面很多的专业资料和文献，以此来弥补自己专业知识上的不足。经过对此课题所涉及专业知识长时间不间断的学习和实践，对高斯地图投影技术有了一定的理解。本论文在一开始首先介绍所研究课题的目的、背景和国内外前沿人士的研究现状，之后讲述了地图投影技术在遥感图像地理编码方面的重要作用，分析研究了高斯投影的电算公式，进一步推导出了不同地球模型下的电算公式，然后使用相应的程序在MATLAB图形用户界面上成功实现了该算法，文末详细描述了MATLAB GUI图形界面的设计与建立过程。本论文重点研究了高斯-克吕格投影方法，知道了它的公式及原理，并使用MATLAB程序成功实现了该算法，从研究理论和实验结果两方面看，都取得了比较理想的效果。通过对本课题的研究，我知道了自己多方面的不足与欠缺，尤其是在投影方法和原理的专业知识方面，地图投影技术在当今社会的许多领域都有着比较广泛的应用，本论文只是其应用的冰山一角，我深知自己需要更多专业知识的学习与增强，在以后的学习中须将所学理论知识与实践相结合，不断地完善自我。参 考 文 献1 周兴东. 利用遥感图像进行土地利用分类方法的研究J. 煤炭学报, 2007, 2(28):4814842 马大玮. 机-星-地遥感图像传输系统的设计J. 重庆大学学报, 2003, 9(26):1151183 林剑. 遥感图像土地利用/覆盖多光谱特征分布空间距离分析方法J. 光谱学与光谱分析, 2009, 2(29):4374394 陈超, 秦其明. 基于地物光谱特征分析的高分辨率遥感图像水上桥梁提取J. 专题论述, 2013, 2(27):7957975 刘德连, 张建奇. 对地遥感图像的目标探测技术研究J. 大学教育, 2006, 4(25) :78856 张剑清, 张勇, 程莹. 基于新模型的高分辨率遥感影像光束法区域网平差J. 武汉大学学报, 2001, 5(23):55587 张剑清, 张祖勋. 高分辨率遥感影像基于仿射变换的严格几何模型J. 张家口职业技术学院学报, 2008,3(21):59618 袁修孝, 曹金山, 姚娜. 顾及扫描侧视角变化的高分辨率卫星遥感影像严格几何模型J. 测绘学报, 2012, 5(23):1731799 李立钢. 星载遥感影像几何精校正方法研究及系统设计J周口师范学院学报, 2010, 5(27):697810 刘景正. 无人飞行器影像场景配准与目标监视技术研究J. 测绘学报, 2001,2(23):11111811 林霖, 杨丰, 张志德. 基于Matlab 的信号与系统课程演示软件的应用J. 电气电子教学学报, 2009, 5(31):9810012 陶冀. 遥感图像在森林资源二类调查中的应用J. 中南林业科技大学学报, 2003,2(24):181913汤毅,万建伟. 地物空间分布特性的高光谱遥感图像解混算法J. 红外与毫米波学报, 2011,3(6):25125414 Li, J., Jarvis, C.H., and Brunsdon, C. (2010). The use of immersive real-time 3D computer graphics for visualisation of dilution of precision in virtual environments. International Journal of Geographical Information Science, 24(4), 591605.15 Savas Ko, Zafer Aydomus a mat lab/gui based fault simulation tool For power system education Mathematical and Computational Applications, 2009, 5(23):207217致 谢在此首先需要感谢在整个毕业设计过程中对我提供帮助的老师和同学。这次的毕业设计是在老师的悉心指导之下完成的。由于论文所研究的课题是我之前从未接触过的一个领域，涉及许多方面的知识，为此我浏览阅读了大量的资料和文献，以此来弥补我在该领域专业知识的不足。从最初设计的选取，到开题报告的撰写，图形界面的布局到最后结果的实现，刘老师都对我进行了详细的指导。通过本次毕业论文，使我能够把大学本科期间所学的专业知识进行了比较系统的实践和运用。从论文开题对本课题的分析研究直到最终的在MATLAB上运用程序实现算法，都是经过老师的指导后自己独立完成的。在毕业论文进行的整个过程中，刘老师定期询问论文的进程，时刻关注论文撰写的进展情况，并且给我有效的建议，为我指点迷津，帮助我开拓研究思路，精心点拨、热忱鼓励。刘老师一丝不苟的治学态度，严谨求实的工作态度，不仅授我以文，而且教我做人，在我一点一滴进步的道路上都凝结了刘老师孜孜不倦的指导与教诲，从刘老师的身上，我不仅学到了应该具备的理论知识，还学到了踏实严谨的态度和遇到困难迎面而上的无畏精神，这些无形的价值在我今后学习或工作中都将会是受益终身的宝贵财富。在此，对刘老师悉心的指导致以最真挚的感谢。在临近毕业之际衷心感谢各位老师，是你们让我在大学本科的四年时间内不仅学到了许多宝贵的理论知识和专业技能，更让我明白了许多为人处世的道理。对以后的学习工作和生活均会有积极的影响。最后感谢各位老师在百忙之中抽出时间对本文进行评审和指导。附 录部分代码如下：fid=fopen(D:workdata-bl.dat,r); m=1425;n=597; x=zeros(m*3,n);x=fread(fid,m*3,n,f