基于RS与GIS的土地利用变化研究.docVIP

PAGE  PAGE 7 基于RS与GIS的土地利用变化研究 　　摘要：基于遥感数据，提取NDVI和MNDWI，利用决策树自动阈值分类等方法研究了长汀县2001～2013年间的土地利用类型变化，结果表明：自2001～2013年，长汀县林地整体大约增加了345km2，约占福建省总面积的11.1%；裸地面积大约减少了403km2，约占福建省总面积的13%，是2001年裸地面积的58.2%；道路面积增加了24.3km2，是2001年的1.1倍多；居民用地增加了约7km2，是2001年的117%。 　　关键词：土地利用；自动阈值决策树；动态监测；长汀县 　　中图分类号：N945.11 　　文献标识码：A文章编号：1674-9944（2016）22-0129-05 　　1引言 　　土地是人类赖以生存和发展的物质基础，是宝贵的不可再生的自然资源，也是可持续利用的生产资料。土地的利用与变化，不仅直接影响社会经济的可持续发展，而且还影响到全球环境的变化。土地利用动态变化逐渐成为研究热门和前沿问题。有研究表明：在将来的数十年间，由土地利用与土地覆被变化所带来的全球影响，要等于甚至大于潜在气候变化所带来的影响[1]。利用遥感影像对长汀县的土地利用覆盖变化进行动态监测，探讨土地利用格局空间演变机制及其生态环境效应，为当地政府部门制定生态安全条件下的土地利用规划提供科学依据，对于区域全面、协调与可持续发展具有重要的现实意义[2]。 　　2数据处理 　　本次研究所使用的主要资料为从美国航空航天局（USGS）网站下载的2001年10月、2005年3月、2009年10月和2013年10月4个时间的Landsat/TM影像数据，空间分辨率都为30m×30m，前3个时间的影像光谱为7个波段的多光谱数据。2013年的影像数据为Landsat8采集，含有11个波段，遥感信息更加丰富，但其波长范围更窄，各个波段的位置也有所变化。此次用于研究的四幅影像比较清晰，云量基本为0%，影像成像质量?^好，确保研究数据的准确性。研究平台以ENVI5.0及Arcgis10.0为主，其他辅助资料及平台有并结合长汀县数字高程（DEM）图、行政区域图、谷歌地球等。 　　2.1数据预处理 　　由于下载的影像为单波段的tif影像，需要对影像进行影像合成，运用ENVI5.0的BasicTools模块中的LayerStack工具将各个时单波段影像合成多光谱影像。 　　辐射矫正：目的是消除大气（水蒸气、氧气、二氧化碳、甲烷和臭氧等）和光照等因素对地物反射的影响，使得影像真实反映对太阳光的辐射情况，以获取地物真实反射率数据，准确地获取研究区的土地利用信息，减小地物识别误差，从而使得土地利用类型分类更加准确可靠。先做辐射定标，然后再做大气校正，且采用直方图法大气校正对影像进行辐射校正。首先对影像进行一次直方图均衡，然后以2001年的影像为参照影像，分别对另外三期影像与2001年影像进行直方图匹配，最终使得四幅影像的辐射特征基本一致[3]。 　　几何纠正：为了消除由于非系统性因素（传感器姿态、高度、地球曲率、地形等因素）造成影像发生的畸变，以提高研究结果的准确性，在ENVI5.0中，以长汀县已通过几何精校正的影像作为参考图像，均匀的选取多个地面相同控制点，采用最近邻法（nearestneighbor）对四幅遥感影像进行几何精校正，且控制纠正误差小于0.5个像元[4～8]。 　　2.2数据处理与分析 　　2.2.1NDVI空间分布 　　在多光谱遥感影像中，利用近红外波段的反射值减去红光波段的反射值，再比上它们两者之和，即可得到归一化植被指数（NDVI）。即： 　　NDVI=（NIR-R）/（NIR+R） 　　其中，NIR为近红外波段的反射值；R为红光波段的反射值。 　　研究使用的数据中，2001年、2005年及2009年的Landsat/TM影像直接使用第四波段与第三波段之差比上两者之和即可得到长汀当年的NDVI值，即： 　　NDVI（2001、2005、2009）=（Band4-Band3）/（Band4+Band3） 　　2013年10月的影像为Landsat8拍摄，与之早期的Landsat卫星相比，相同光谱成的影像位于不同的波段位置，计算NDVI时，其第五波段为近红外波段（NIR），第四波段为红外波段（RED），计算公式为： 　　NDVI2013=（Band5-Band4）/（Band5+Band4） 　　通过上述方法分别计算得到长汀县2001年、2005年、2009年及2013年的NDVI值空间分布图，统计NDVI的最大值、最小值、加权平均值和标准差，统计结果见表1。 　　结合NDVI值空间分布图表得出结论：自2001～2013年期间