《第五章 城市遥感ppt课件.ppt》由会员分享，可在线阅读，更多相关《第五章 城市遥感ppt课件.ppt（37页珍藏版）》请在三一办公上搜索。
1、第五章 城市遥感,城市与环境学院王细元xiyuan80_http:/ 遥感地学分析，黄家柱教授,遥感地学分析,内容提要,5.1 城市遥感概述5.2 城市扩张遥感监测5.3 城市热岛效应遥感监测5.4 城市人口密度遥感监测,5.1 城市遥感概述,城市是人口集中、集约经济活动及不同生活方式并存的复杂社会。城市又是人类活动的缩影，并且不断地经历着迅速变化的过程，需要及时地进行监测与分析。城市规划和城市建设者面临的重大任务之，就是获取与分析那些能有效地用于城市规划、建设和管理的资料。城市遥感的任务就是为城市规划、建设和管理提供多方面的基础地理信息和其他与城市发展有关的资料，诸如城市土地利用现状、城市演
2、变、城市及区域的自然状况、城市人门及其分布情况、城市道路与交通状况、,城市热岛、通讯受地理限制的因素等。城市遥感与传统的城市相关资料调查相比，既省时，又省钱，而且效率很高，因而具有广阔的应用和发展前景。,图中的目的、内容和要求出调查的项目需要所决定；遥感图像包括两类，一是航天遥感图像(主要用于区域性和市域性的宏观调查)和航空遥感图像(主要用于建成区和城市局部地区的较微观调查)；地形图主要用于划分调查的空间层次、地理单元和影像解译时的参考以及作为遥感调查的基础底图。,典型地物在Landsat TM影像上的光谱曲线,地物光谱特征,5.2 城市扩张遥感监测,一、归一化建筑用地指数 NDBI 在TM4
3、和TM5波段之间除了城镇用地DN值走高之外，其他地物DN值都变小。 因此图像上NDBI值大于0的地物则认为是城镇用地。,原始图,NDBI结果图,NDBI二值化结果图,二、改进的归一化裸露指数 MNDBI,NDBI主要反映的是城镇和裸露地信息，所以将NDBI和（1-NDVI）相加就可以更加突出居民地信息。因此将之称为改进的归一化裸露指数（Modified Normalized Difference Barren Index, MNDBI），即,MNDBI是对NDBI的改进突出了城市信息，使其与周围地物的反差增大，有利于提取城市信息。,原始图,NDBI结果图,NDVI结果图,MNDBI二值化,MN
4、DBI结果图,NDBI结果图,三、城镇用地指数 ULI,在NDBI基础上引入NDVI，通过对图像的二值化求交（相乘）运算，利用这两个指数的各自的优势来提取城镇用地，把这两个指数的结合用于提取城镇用地称为城镇用地指数（Urban Land-use Index，ULI）。,典型地物二值化后的像元值,NDBI结果图,反NDVI结果图,ULI结果图,$n1_ndbi * $n2_no,三种基于谱间特征分析的城镇用地提取方法比较,NDBI,MNDBI,ULI,5.2.2 遥感监测城市扩张的实例,选用的是Landsat TM影像，研究区覆盖浙江省金华地区，成像时间是1996年9月6日，TM6波段是热红外波
5、段，仅选用了15波段和7波段，其空间分辨率为30m。其他数据包括金华市1996年的土地利用现状图以及1:5万地形图，用于几何校正和精度检验。,NDBI法提取的城镇用地,（a）大气校正后的TM影像(543合成),(b) NDBI图,(c) NDBI二值化图,ULI提取的城镇用地,(a) NDVI图 (b) NDVI二值化图 (c)求交运算后的二值化图,NDBI提取城镇用地精度分析,ULI提取城镇用地精度分析,5.3 城市热岛效应遥感监测,5.3.1遥感监测城市热岛效应的原理与方法,城镇扩展对大气环境的另外一个显著影响就是城市的“热岛效应”。快速城市化进程改变了地表下垫面的理化性质。原本是土壤、草
6、地和水体等比热大的自然表面被水泥、沥青等比热小的表面代替，这不仅改变了反射和吸收面的性质，还改变了近地面层的热交换和地面的粗糙度，使大气的物理状况受到影响。大量的观测对比和分析研究确认，城市热岛是城市气候中最普遍存在的气候分布特征。这种城市气温高于四周郊区气温的现象称为城市“热岛效应”，有时也统称为城市热岛。城市气温与郊区同期（瞬时、日平均、月平均、年平均等）气温差值大小，则称为城市热岛强度。,热岛效应示意图,第一种方法需要已知地表比辐射率；主要包括大气校正法和单窗算法。适用于只包含一个热红外波段的遥感影像，如Landsat TM影像。第二种方法根据分裂窗通道对水汽吸收的差异将地表比辐射率作为
7、输入变量进行大气和地表比辐射率订正。适用于包含2个或多个热红外波段的遥感影像，如NOAA /AVHRR和MODIS遥感影像。第三种方法：白天/夜间MODIS LST方法，该方法利用MODIS 的7个热红外通道的白天/夜间资料同时反演地表温度和通道平均比辐射率，而不需要高精度的大气温度和水汽廓线。,二、地表温度遥感反演的方法,1、适用于TM遥感影像数据的算法 由于Landsat TM影像只有第6波段是热红外波段，因此，在反演地表温度的时候只能用单一热红外通道方法进行地面温度的演算。算法具体演算过程如下： 对于Landsat 5，求算卫星高度的像元亮度温度： 式中 为TM6的像元DN值，0 255
8、，T6 为亮度温度值，单位为K。由于大气辐射和地表热特性的影响，卫星高度的亮度温度与实际地表温度有较大差距。对于要求精度较高的地面热量空间分析，有必要进行较为精确的地表温度演算。,式中Ts单位为K；a6=-67.355351，b6 =0.458606，C6和D6为中间变量，分别用以下两式表示 ：,因此，只要知道参数 即可用上述单窗算法推算任何像元的实际地表温度。对 、 进行估计时可通过一些简单易行的方法：大气平均作用温度 主要是根据当地的地面气象观测数据（地面附近的气温和水分含量）进行估计，如果不能获得当地地面气象观测数据，则可以利用适当的方法进行模拟替代；大气透射率 主要是根据大气水分含量来
9、估计。当参数估计没有误差时,该方法的地表温度演算精度达到0.4,在参数估计有适度误差时,演算精度仍达1.1。从而能满足大多数应用的精度要求。,2、适用于MODIS遥感影像数据的算法,1）分裂窗方法 劈窗方法是到目前为止应用最广泛的地表面温度反演方法，尤其是用来分析NOAA/AVHRR数据。Becker等通过研究认为，地面温度可以表示为两个热红外波段在大气顶层亮温的线性组合即:,其中，Ts为地表温度，P、M为常数，且,T4和T5分别为NOAA卫星AVHRR第4和第5波段的亮温值 ， 和 分别为NOAA卫星第4和第5波段发射率， 为二者之差。,2）白天/夜间MODIS LST算法,该算法适用于包括
10、具有变化的或未知比辐射率的地表在内的所有地表。在MODIS通道j测量到的辐射率可表示为：,其中各项均为通道平均值, 是通道比辐射率， 是地表温度为Ts的黑体辐射率， 散射太阳辐射， 大气层顶的太阳辐射率， 和 分别为到达地表的通道平均太阳漫射辐射和大气向下的热辐射。 为通道有效传递函数。该方程可由统计回归法或最小二乘法解出。Ts即为所要求算的地表温度。,地表温度反演单窗算法（TM6）,亮度温度：是遥感器质卫星高度所观测到的热辐射强度相对应的温度。从TM6数据中求算亮度温度：DN值转化为相应辐射亮度；根据辐射亮度推算相应亮度温度。,DN值转化为相应辐射亮度,L为TM所接收到的辐射强度Lmin、L
11、max为TM遥感器接收到的最大最小辐射强度，可在元文件中查找DN为像素值，DNmax为最大像素值=255。,辐射亮度推算相应亮度温度,亮度温度可以用Planck辐射函数计算，或者用以下近似式求算：T6为像元亮度温度（K），K1和K2为发射前预设的常量，对于Landsat5的TM数据K1=60.776 mWcm-2sr-1um-1K2=1260.56 K,卫星高度所接收到的热辐射强度 为地表温度， 为TM6的亮度温度， 为大气透射率。 为遥感器接收到的热辐射强度， 为地表在TM6波段区间内的实际热辐射强度，直接取决于地表温度， 分别为大气中TM6波段区间内的向上和向下热辐射强度。,式中Ts单位为
12、K；a6=-67.355351，b6 =0.458606，C6和D6为中间变量，分别用以下两式表示 ：,因此，只要知道参数 即可用上述单窗算法推算任何像元的实际地表温度。对 、 进行估计时可通过一些简单易行的方法：大气平均作用温度 主要是根据当地的地面气象观测数据（地面附近的气温和水分含量）进行估计，如果不能获得当地地面气象观测数据，则可以利用适当的方法进行模拟替代；大气透射率 主要是根据大气水分含量来估计。当参数估计没有误差时,该方法的地表温度演算精度达到0.4,在参数估计有适度误差时,演算精度仍达1.1。从而能满足大多数应用的精度要求。,5.3.2 遥感监测城市热岛效应实例,以西安地区为例
13、，介绍一个遥感监测城市“热岛效应”的例子。选用了1988年8月23日和2003年5月29日上午10点的TM6数据进行地表亮温反演，通过对比两期的地表亮温的变化情况，分析15年因城镇扩展所带来的城市“热岛效应”变化的影响。,(a) 1998年 (b) 2003年,研究区TM6热红外伪彩色图像,(a) 1998年 (b) 2003年,研究区TM6热红外三维伪彩色图像,(a) 1988年 (b) 2003年,研究区亮温等级分布图,5.4 城市人口密度遥感监测,5.4.1遥感监测城市人口密度的原理与方法 利用土地密度法估算人口密度的方法应用较多。假设城市某区域内的同一类型住宅的密度一致，通过抽样调查住
14、宅面积和人口，求得每一住宅类型的代表人口密度D1，D2，Dn，利用航空像片或卫星影像判读出各种住宅类型并求出其面积S1，S2，Sn，则可用下式估算该区人口总数P估：,P估S1* D1+ S2* D2+ Sn* D 此时在已知该区的统计人口数量P统，则可求出估算误差P：P|P统P估| / P统当误差在一定小的范围内时，可以进行全面估算。,5.4.2 遥感监测城市人口密度的实例,设某城市区域有n种居住用地类型，每一住宅类型的人口密度为Di（i1，2，n）；又设该区域被分为m个行政街道，且已知每一街道的统计人口数量为Pj（j1，2，m）；航空遥感调查获得的每个街道内的各种住宅用地面积为Sij（i1，2，n；j1，2，m），则可建立线性方程组：,在此方程组中mn时，即街道个数多于居住类型数时，可采用最小二乘法原理，将方程组）转化成为nn正规方程组。解正规方程组，即可求得该区域内统计人口数总误差为最小的各类住宅密度估计值。,本章结束,
