荒漠化监测的理论与技术基础课件

*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,荒漠化监测,生态环境学院,水土保持与荒漠化防治学科,主讲：杨光,2009,年,9,月,荒漠化监测是对地球表面或其上的某一地区的荒漠化现象进行定期或不定期的观测、最终以某种媒介形式进行表述、评价的行为，因而荒漠化监测（特别是遥感监测）必然涉及地球形状、地球参数、地图投影系、比例尺与分幅、“3S”技术等一系列问题。因此，我们首先从与荒漠化监测相关的几个问题入手，逐渐过渡到荒漠化监测问题本身。,第二章 荒漠化监测的理论与技术基础,地球表面的复杂性使认识和表述地球的形状非常不易。现在的地球形状是由,大地水准面,决定的。所谓的大地水准面是指假设将海水引入陆地后形成的有一定水面高的面，但其形状不易为数学公式所表述，因而用,椭圆体,进行近似。这样，作为地图对象的地球形状就成为以沿椭圆短轴为转轴旋转而成的椭球面。这一形状既是地图学的基础，也是我们进行荒漠化监测的基础。,2.1 地球形状及地球参数2.1.1 地球形状,2.1.2 地球参数,在近似地确认了地球形状之后，就可以用地球参数（长半轴a，短半轴b，离心率e，扁平率f）来描述地球形状了。一般地，,世界各国都是以自己的国家为中心测定和计算地球参数的,，迄今为止的地球参数有20余种。,即使是在选定了地球椭圆体之后，也还存在一个选定的椭圆体如何与实际地球匹配、定位的问题。因此，在从事荒漠化监测，不仅要注意各国所采用的基准椭圆体（Reference Ellipseid），而且还要注意该基准椭圆体与地球的相对位置关系。,离心率定义为椭圆两焦点间的距离和长轴长度的比值。,各种地球椭圆体的形状参数表,测算者,发表年代,长半轴,（,m）,短半轴b（m）,扁平率f,离心率e,Delambre,1800,6375563.00,6356103.00,327.623998,0.07807203,Ebereslo,1810,6377304.00,6356103.00,300.802038,0.08147294,Ebereslo,1830,6377276.35,6356075.41,300.801700,0.08147298,Bessel,1841,6377397.16,6356078.96,299.152813,0.08169683,Clarke,1866,6378206.40,6356584.00,294.981427,0.08227147,Clarke,1880,6378249.15,6356515.00,293.466300,0.08248315,Harkness,1891,6377972.00,6356727.00,300.200000,0.08155303,Helmerl,1901,6378200.00,6356818.00,298.300000,0.08181366,Hayford,1906,6378388.00,6356912.00,297.000745,0.08199179,Willisch,1915,6378372.00,6356896.00,297.000000,0.08199189,国际基准椭圆,1924,6378388.00,6356912.00,297.000000,0.08199179,Heiskanen,1926,6378397.00,6356921.00,297.000000,0.08199173,Krassowskij,1940,1942,6378245.00,6356863.00,298.299738,0.08181337,Fisher,1960,6378163.00,6356781.29,298.300000,0.08181333,IAU-64,1964,6378160.00,6356775.00,298.250000,0.08181964,IAU-67,1967,6378160.00,6356774.52,298.247167,0.08182057,WGS-72,1972,6378135.00,6356750.52,298.260000,0.08181881,NWL-9D,1973,6378145.00,6356759.77,298.250000,0.08182018,SAO-SE3,1973,6378140.00,6356755.22,298.256000,0.08181936,IAU-76,1976,6378140.00,6366755.29,298.257000,0.08181922,GRS-80,1980,6378137.00,6356755.00,298.257000,0.08181406,GRS-80改订,1984,6378136.00,6356751.03,298.257000,0.08181974,WGS-84,1986,6378137.00,6356752.31,298.257223,0.08181919,在荒漠化监测、资源调查等涉及制图的工作中，一般是按,如下原则考虑和选择地球形状和参数,的：,（1）绘制世界地图、各大洲地图等以广泛地域为对象区域的小比例尺地图时，可视地球为,半径6370km的球体,；,（2）绘制100km,2,以内地域的单幅地图时，可将地表面作,平面,处理；,（3）绘制中、大比例尺地图时，一张图幅的对象地域虽然狭小，但如果需要对广大地区进行地图拼接时就得将地球作为,椭球体,考虑。对于必须保证精度的小比例尺地图，也需如此处理。,荒漠化制图方法的规范与否决定了制图质量及制图的效率，尤其是制图编码的标准化减少了数据汇总的难度。,2.2 地图投影,地图,是将呈球面的地球表面在片面上进行表示的一种图示形式。,地球是一个庞大的球体，如果只是考察地球表面上相当狭小的某一范围时，是可以近似作为平片面处理的。但当考虑面积增大到某个省、某个国家、某大陆乃至世界地图时，就象把一个乒乓球破开、压平时必然产生破裂或褶皱一样，地球的形状、大小和相互关系无法得以正确表示，必然产生许多误差。,尽可能减少误差以便将地物的形状、位置在平面上进行表示的方法，称为,地图投影,。,地图投影就是在平面上将地球上的经纬度线加以表示的一种方法。,某地经纬度,地图横纵坐标,地图投影,2.2.1 地图投影变形,在地图投影时，把球面上的经纬线网转换到投影平面上，转换后的地图上的经纬线网格必然产生变形，这种变形称为,地图投影变形,。,这种变形主要反映在三个方面，即,长度变形,、,面积变形,和,角度变形,。,长度变形,指投影后地图上不同地点和不同方位上的地球表面实际距离与相应图面距离的比值（比例）各不相同，从而无法从地图上量算和比较不同地点和不同方位的景物之间的距离的变形。,对于反映广大区域的小比例尺地图来说，在投影中心（中央经线、基准纬线、基准大圆、基准小圆）处可以获得正确的比例尺，而在其它区域就得不到正确的比例尺。离原点越远，长度变形越大；即使是同点，不同方向上的长度变形也是各异的。,面积变形,是指投影所得地图上的面积比例随地点而改变，其结果是导致不能在地图上量算和比较景物所占的面积。,角度变形,是指地图上角度同相应景物在球面上的相应角度不相等，例如在地球仪上成直角的经线和纬线，在地图上经纬线不成直角。,要完全消除投影变形是不可能的。投影时只能根据地图的应用目的，牺牲上述三个变形中某个方面的精度要求，,设计、开发或选择可保证地图应用精度要求的地图投影方法。为此，目前的地图投影法已多达千种以上。一般地，大型地理信息系统软件都支持几种常用的地图投影法，市面上也有专门进行地图投影转换的软件出售。,2.2.2 地图投影的分类,（1）根据地图投影变形性质的分类,等距离投影,等面积投影,等角投影,任意投影,对等距离、等面积、等角度投影而言，,等距离要求只能在地图上的一小部分内实现,，即使是在等角、等面积图中也可以实现。但等角、等面积条件却不可能在同一张图上实现，,等角、等面积要求互相冲突,，等面积的获取是以牺牲等角为代价；反之，等角的获得又以牺牲等面积为代价的。,（2）根据投影面与球面位置关系的分类,在地图投影中，首先将不可展的球面投影到一个可展的曲面上，然后将曲面展开成平面，从而得到不同投影。,其中，投影曲面可以是,圆柱面,、,圆锥面,及,平面,；投影面和地球的位置关系既可以是投影面同地轴,垂直,又可以是,重合,或,斜向相交,。根据投影面与地球表面的相交关系，可对地图投影进行如下分类：,方位投影,圆柱投影,圆锥投影,不同投影的投影面和球面的关系图,2.2.3 地图投影的选择,采用什么样的投影方法作图，是制图工程中必须最初决定的内容。选择投影方法时必须考虑地图的,应用目的,、,比例尺,、,对称地域的形状,、,面积,、,纬度,等因素，此外也需考虑,投影计算,、经纬线,作图难易,、展开后编辑工作,是否便利,等因素。不管是多么理想的投影方法，如果它计算繁复、经纬线形状复杂以至于影响以后的编辑工作，那么也不能采取这样的投影方法。当然，如果是应用计算机制图就另当别论。此时应主要考虑投影方法的适宜性。,（1）根据比例尺进行选择：,大、中比例尺地图：,由于作图范围是几乎可以不考虑图形变形的狭小区域，国土基本图、地形图、地势图、城市规划图、比例尺在一至二万分之一的地图、以及测量中所用的平面直角坐标等，可主要应用横轴圆柱等角投影法（横轴麦卡托投影）或多面体投影法。日本以前的地形图多用多面体投影法，现在则主要使用国际横轴麦卡托投影法。,小比例尺地图：,五十万分之一以下的地图、大区地图、省区地图、国土态势系列地图、百万分之一的国际地图等，一般应用多圆锥投影、等角圆锥投影等。,（2）根据地图使用目的进行选择,地图，一般可分为,一般地图,与,特殊地图,两类。对于一般地图，可在考虑地图比例尺、对象区域形状、纬度等条件的基础上选择投影方法；但对于特殊地图就必须选择特殊的投影方法。如：,航空图：等角圆锥投影、等角圆锥投影等；,航海图：麦卡托投影；,大圈航海图、无线方位探察图：心射投影；,统计图：各种等面积投影方法。,2.2.4 投影转换,实际从事荒漠化监测时，在下述一些情形下可能要进行投影转换：,（1）从事,跨区域的荒漠化监测工作,，而不同的区域又采用的是不同的投影系时。此时，为了综合制图，就需要把不同的投影系统一到某一种投影系上；,（2）监测中,拟使用不同时期的图形资料,，而两个或几个时期的地图投影又不同，此时也需要进行投影系转化，投影系统一后才能混合使用这些资料；,（3）当监测是,由小区域的大比例尺做起,，最后汇总到省（区）乃至全国一级时，也必须进行投影系的转换工作；,（4）GPS使用的是地心坐标系（WGS-84），使用时如果所用的,接受器没有将WGS-84坐标转换为你正在使用的坐标,，也必须进行坐标转换；,（5）虽然同为卫星定位系统，但美国的GPS（WGS-84）与前苏联的GLONASS（SGS-85）却采用的是不同的坐标系。当需要,同时使用这两种全球定位系统以提高精度时,，就需要进行投影系的转换工作。,2.2.5,我国主要应用的投影方法： 高斯,克吕格投影,高斯,克吕格投影是,等角横切圆柱投影,。,它假设一个椭圆柱面横切在椭球面的一条子午线（中央经线）上，椭圆柱的中心轴通过地球椭球中心，然后按等角条件将中央经线两侧一定正负差范围内的区域投影到椭圆柱面上，再将椭圆柱面展平而得到的投影（图a）。其投影坐标关系为，以中央经线的投影作为X轴，赤道投影作为Y轴，两轴交点O为原点，其经线和纬线的投影均为曲线（图b）。,高斯克吕格投影示意图,这种投影中的中央经线投影为直线，其长度等于球面实际长度，没有变形，其余经线为向极点收敛的弧线，距中央线越远，变形越大。为使边缘变形不致过大，采用分带投影的方法。以6投影带为例，赤道上边缘部分的最大长度变形不大于0.14%，面积最大变形不大于0.27%。我国处于中纬度地带，其变形数值更小。因为6带是从西向东每6为一带，全球共分60个带，依次编号160。我国位于13带至23带之间（如下图）。,高斯,克吕格投影分带示意图,高斯克吕格投影的,主要优点为,：,等角性质，适于系列比例尺地图的编制与使用；,经纬网与直角坐标的偏差不大，便于阅读使用；,计算工作量小，直角坐标和子午收敛角值只需计算一个带，全球可通用。,因此，我国各种大中比例尺地形图均采用该投影方法。其中万分之一地形图采用3带，二万五千分之一至五十万分之一地形图采用6带。此外，目前世界大多数国家的地形图也都使用此种投影方法。,2.3 比例尺与地图分幅,按比例尺的大小可以将地图分为大、中、小比例尺三类，具体标准如下所示：,大比例尺地图,：比例尺大于或等于1:10万的地图；,中比例尺地图,：比例尺在1:10万到1:100万的地图；,小比例尺地图,：比例尺在小于或等于1:100万的地图。,基本比例尺地形图的分幅、编号是在国际规定的1:100万地形图基础上，按经纬度进行分幅和编号的。,2.3.1 1:100万地图的分幅与编号,百万分之一国际标准图幅按纬差4、经差6划分。,即从赤道起向两极，纬差6为一列，依次用字母A、B、C、D、S表示，经度从180算起自西向东，每经差6为一行，依次用1、2、3、60表示（如下图）。编号的书写顺序是列号在先、行号在后，二者之间划一短线。比如：北京在百万分之一地图上是处于第J列、第50行，故其编号为J50。,1:100万地图的分幅与编号图,2.3.2 1:50万、1:20万、1:10万地图的分幅与编号,这三种比例尺地图的分幅、编号是在百万分之一地图的基础上按一定经差和纬差进行划分的,，划分后分别在百万分之一地图的分幅编号后加上各自的代号；,（1）1:50万地图,：,按经差3、纬差2进行划分，即一幅百万分之一的地图可划分为四幅五十万分之一的地图，后缀编号为A、B、C、D等，如J50A；,（2）1:25万地图,：,按经差1、纬差40划分，共可划分出36幅图。后缀编号是用带括号的自然数序列（1到36）从左到右、从上到下排列的。如J50（3）；,（3）1:10万地图,：,按经差30、纬差20 进行划分，共可划出144幅图。后缀编号是用不带括号的自然数序列（1到144）从左到右、从上到下排列的。如J505。,1：50万、1：20万、1：10万地形图分幅编号示例图,2.3.3 1:5万、1:2.5万、1:1万地图的分幅与编号,这三种比例尺的地图都是在1:10万地图的基础上按一定经差和纬差划分,，然后分别在该1:10万地图分幅编号的后面加上各自的分幅编号的。,（1）1:5万地图,：将一幅十万分之一的地图划分为四幅五万分之一的地图，然后在该十万分之一地图编号的后面缀以A、B、C、D等，如J505B；,（2）1:2.5万地图,：将一幅五万分之一地图划分为四幅二万五千分之一的地图，然后在该五万分之一地图编号的后面缀以1、2、3、4等，如J505B4；,（3）1:1万地图,：将一幅十万分之一地图划分为五十四幅万分之一的地图，然后在该十万万分之一地图编号的后面缀以带括号的1、2、3、54。如J505（15）。,1：5万、1：2.5万、1：1万地形图分幅编号示例图,2.4 地形图选择匹配及TM图像处理技术,2.4.1 地形图选择匹配,在图像处理的同时，必须建立一套正确的解译方法。卫星影像的目视解译，涉及植被本身和气候资源、地貌资源、地表水资源以及土地利用等各个方面。因此必须采用多学科相结合的综合性解译方法。,2.4.2 TM图像处理技术荒漠化监测遥感数据处理信息系统,（1）系统建立的目的和意义：,遥感技术以其,多平台,、,多层次,、,多波段,、,短周期,、,低价格,等特点，已广泛用于资源、环境、生态研究的众多领域，是获取实时性资料的重要手段。为了从大量遥感数据中获取我们监测荒漠化所需要的重要信息，,建立荒漠化监测遥感数据处理信息系统是必须的，也是必要的,。,（2）,荒漠化监测遥感数据处理信息系统的构成及功能,荒漠化监测遥感数据处理信息系统的构成：,硬件组成：,SGI系列工作站，CHALLENGE服务器，OCTANE1，OCTANE2等网络环境。HP750CM喷墨绘图仪、磁带机等输入输出设备。,软件组成：,ERMAPPER5.1、ERRDAS8.3等图像处理软件。,荒漠化监测遥感数据处理信息系统的功能如下图所示。目前，我们正利用该系统进行荒漠化监测遥感数据的处理工作，并将不断完善。,荒漠化监测遥感数据处理信息系统,卫星数据输入,几何校正,图象信息增强,信息特征提取,结果输出,图象自动分类,荒漠化监测遥感数据信息系统示意图,（3）,荒漠化监测的TM图像处理技术,TM图像的特点,：,TM图像是目前人们广泛应用的航天遥感图像，概括地说，它具有如下特点：较高分辨率，30米30米；波段信息大，有7个通道；覆盖面大，每景185km185km；获取速度快；几何精度高，利于计算机处理，分类能力好；价格中等。,TM图像处理技术,：,我们进行图像处理的目的主要是为了在大范围内掌握荒漠化的情况，能够全面、迅速地提供信息，输出高质量的图像，为目视判读和分类提供依据。我们的基本思路如下图。,选择信息源,几何校正,波段组合,图像增强、拉伸变换,加载地理信息,提取有关信息,分类、制图输出,TM图像处理流程图,代号,地类,TM4，5，3假彩色合成图像影像特征,11,水田,红（稻熟期）或蓝黑色（未种稻，但有水），边界清淅，色调均匀，规则块状,12,水浇地,鲜红、紫红或青灰色，条块分明，边界清淅，色调均匀，规则块状,13,旱地,多为浅红或黄红色斑块，色调不均，形状不规则,21,有林地,深红、 紫红、暗红色、内部色调较均匀，形状不规则或规则（人工林）,22,灌木林地,棕红色或红褐色，或青灰色中带“红雾”，内部色调不均，形状不规则,31,高盖度草地,红（优质草场）、浅红、棕红、等各种色调相混杂，边界不明显，形状不规则,32,中盖度草地,灰黄色，浅灰黄，呈不规则的片、块、薄雾状，边界不清淅,33,低盖度草地,浅灰、灰白、浅灰蓝色，形状不规则,41,城乡居民点及特用地,灰、淡灰、灰黑色、斑块状，边界清淅,42,独立的工矿用地,青灰，折线轮廓明显，内部色调不均，形状有规则,43,盐田,紫红色，形状规则，边界清淅,44,交通用地,灰黑色或紫褐色曲线条（直线条），平滑，色调和宽度均匀,51,河渠及湖塘、水库用地,蓝褐色，宽窄不一，椭圆，三角形等，流线状弯曲长条状轮廓线圆滑分明，水库坝址呈灰白色直线,52,滩涂,灰白色中带青、红小块，内部色调不均，与黑色水域相邻形状不规则,61,流动沙丘及沙地,浅灰黄、灰白色、呈新月形、格状、波状等，边界清淅形状不规则,62,裸戈壁,墨绿色或黑色中带绿色条带，呈明显冲积扇形，色调较均匀，边界清淅,63,盐碱地,灰白、白灰、淡蓝色，形状不规则，花斑状，较长案较均匀，边界清淅,64,裸土地,灰白色、或白色片状，内部色调较匀，边界不明显，形状不清淅,65,裸岩、石砾地,表灰色中有明显沟状、立体感强，内部色调不均，形状不规则,2.4.3 实例遥感目视解译标志,123波段合成 432波段合成,453波段合成 743波段合成,