参考文献背景方法概述遥感.doc

1 基于 DMSP-OLS-夜间灯光数据的长江三角洲地区城市化格局与过程研究 数据源：本文使用 AcrGis Desktop 作为地理信息处理软件(数据源包括美国国家地球物理数据中心提供的DMSP-OLS夜间灯光数据以及通过国家基础地理信息系统网站(htttp://nfgis.nsdi.gov.cn) 下载的全国1:400万数据库中的县级以上境界面状数据和主要交通线状数据(二者的数据都是地理坐标)，椭球体为 WGS-84，可以相互叠加(为了计算面积，需对二类数据的坐标加以转换，本文将其转换为横轴墨卡托( UTM) 投影，椭球体仍为WGS-84。在县级以上境界面状数据中提取长三角洲地区各市和市辖区边界，用以裁剪全球灯光数据的计算。原始数据是消除了云( 通过热红外通道阈值法) 及火光等偶然噪声影响后对全年VNIR通道灰度值直接平均化处理后得到的，数据灰度值范围1~63，饱和灯光灰度值为63。 研究方法： 构建城市化指数。 通过分析了城市灯光数据（区域灯光面积属性指标顿、平均灯光强度属性指标、城市灯光亮度总值、复合夜间灯光数据）与长三角地区各城市经济（全市GDP、全市人均GDP、第二和第三产业产值比重）、人口（全市人口密度、市辖区人口密度、非农业人口占总人口比重）和土地利用（市辖区建设用地面积、市辖区建成区面积）的相关性。采用平均灯光强度属性指标作为基于DMSP-OLS的城市化水平指数。 （2）灯光阈值的确定。 阈值的确定方法包括经验法、突变检测方法和基于辅助资料的空间比较法等。本文采用空间比较法中的二分法来确定灯光阈值（阈值为DN时的面积与建设用地面积之差绝对值最小) 具体步骤： 图中，DNT为夜间灯光数据的阈值：S(DNT)为统计阈值内的灯光斑块面积；△S(DNT)为S(DNT)与各市市辖区建设用地面积（Area）的差值绝对值；DNmax、DNmin分别为阈值的最大、最小值；△S(DNT-1）、 △S(DNT+1）分别为阈值DNT-1、DNT+1时的面积与建设用地面积之差的绝对值；IF[△S(DNT-1）≥△S(DNT）≤△S(DNT+1）表示阈值为DNT时的面积与建设用地面积之差绝对值最小。 结果分析：1城市群城市化空间模式。城市化过程可以理解为城市用地在空间上通过占有非城市用地而实现扩展的过程。 2城市群城市化过程分析。通过区域内平均灯光强度均值的变化分析城市群整体城市化的过程，平均灯光强度的变差系数来反应城市群内部城市化水平的差异。 3城市群用地扩展类型。基于DMSP-OLS数据特点，像元灯光强度随时间的变化在一定程度上可以反映其所在区域开发利用强度的变化过程，因此可以根据城市像元灯光强度的时间变化特征来识别不同的城市用地主导扩展类型。 2基于DMSP-OLS灯光数据新疆天山北坡经济带城市化研究 研究理论背景：用夜间灯光数据提取城市主要是基于城市夜间灯光强度与乡村黑暗背景的强烈反差特征，采用阈值法进行提取，即设定一个阈值 DN0，如果象元的 DN 值大于或等于 DN0将其划为城市类型，如果小于 DN0，则认为该象元代表非城市用地类型。由于 DMSP/OLS 具有不同的数据产品格式，而且不同的研究者分析思考问题的方式、角度不一，因此阈值的确定方法也不一样。目前确定阈值常见的方法主要有三种： A ) 经验法。即根据以往研究的经验确定阈值。Croft[4]、Kramer[12]等曾先后用它来绘制城市分布图。 B ) 突变理论。Imhoff 认为[13]：在一段时间的连续观测中，像元发出的灯光被探测 到的频率越高，它属于城市类型的概率越大；正确提取城市的阈值应该保持市中心的完 整(不破脆)。基于这种假设，可以通过逐渐增加阈值来减小噪声带来的影响。当将阈值 向 100 逐渐增加时，城市多边形逐渐收缩，当阈值增加到某一个点时，明亮的城市多边 形不再沿着边缘变小，而是开始在内部出现破脆，即明亮的多边形内部开始析出黑暗的 小岛，此时所研究的城市多边形的周长会突然增加，那么使周长突然增加开始析出黑暗 的小岛，此时所研究的城市多边形的周长会突然增加，那么使周长突然增加的阈值点即 为所寻点。 C ) 空间比较法。就是将 DMSP 夜间灯光图像与其他遥感图像做空间比较分析以确 定城市化边界。宫鹏等人基于 DMSP /OLS 拼合的两种数据产品：稳定灯光图像，辐射 标定灯光图像，确定了具有不同城市化水平和经济发展水平的三个城市—美国旧金山、 中国北京、中国拉萨的边界。并将这些边界一一与 TM 提取出的城市边界比较，以此调 整并确定不同发展水平下城市的阈值。 数据：非辐射定标的夜间平均灯光强度数据。 方法：主要依据