龙口市海岸带土地利用动态变化分析_吴泉源.pdf
《龙口市海岸带土地利用动态变化分析_吴泉源.pdf》由会员分享,可在线阅读,更多相关《龙口市海岸带土地利用动态变化分析_吴泉源.pdf(10页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/第25卷 第5期2006年9月地 理 研 究GEOGRAPHICALRESEARCHVol125,No15Sept1,2006 收稿日期:2006205206;修订日期:2006207231 基金项目:国家自然科学基金资助项目(40271001,40471122)作者简介:吴泉源(19592)男,山东龙口市人,教授,硕士生导师。主要从事遥感教学和研究工作。E2mail:wqy6420582 1631com龙口市海岸带土地
2、利用动态变化分析吴泉源1,2,侯志华3,于竹洲4,姜春玲2,邹 敏2,杨圣军2,李轶平2,韩丛丛2(11 山东科技大学地球科学与工程学院,山东 泰安271019;21 山东师范大学人口资源与环境学院,济南250014;31 太原师范学院城市与旅游学院,山西 太原030012;41 龙口市金泰规划建筑设计院,山东 龙口265701)摘要:利用多期遥感数据,采用目视解译和人机交互计算机分类技术提取1984至2004年间龙口市海岸带土地利用信息,从土地利用总量变化、土地利用变化速度、土地利用类型之间的相互转化、土地利用类型变化的海岸区位效应等方面分析龙口市海岸带动态变化特点。研究结果表明,20年间龙
3、口市海岸带一半以上的土地类型发生了变化,变化速度逐年增加,后期明显大于前期;转化方向上主要表现在耕地向园地的转化,其次是耕地、园地、林地向建设用地的转化;变化的区位特征是距海越近,海岸带特征地类的变化越明显,反之,距海越远,内陆土地利用变化的特点表现得越突出。关 键 词:遥感;海岸带;土地利用;龙口市文章编号:100020585(2006)05209212101引言 海岸带是海洋系统与陆地系统相连接、复合和交叉的地理单元,既是地球表面最为活跃的自然区域,也是资源环境条件最为优越的人文活动区域,同时又是全球变化影响最为敏感的地带。对它的研究体现了陆 海 人三大系统交互作用的动力机制,落实于调整“
4、人地关系”这一地理学的中心问题,是地理学研究深化的一个突破口和切入点,成为全世界广为关注和研究的问题之一13。近年来,随着社会经济的发展,人们在海岸带进行了大规模的土地资源开发活动,使得海岸带面临的压力越来越大,土地资源和环境发生了前所未有的变化,出现许多有碍于可持续发展的问题。因此,人们迫切希望能及时、准确、高效地掌握海岸带土地变化信息,以实现海岸带土地资源的科学管理和持续利用。在科技迅速发展的今天,遥感技术在这方面的应用脱颖而出,它对地观测具有宏观、快速、综合、高频、动态和节省投资等突出优势,为海岸带土地资源的动态监测研究提供了一种方便和重要的现代化方法49。本文利用多源、多时相、多分辨率
5、卫星遥感影像,在遥感技术和GIS技术支持下,对1984年以来龙口市海岸带土地利用信息实现快速、准确的提取,并对其演变状况进行分析研究,为龙口市政府有关机构、行业及应用部门在海岸带土地利用、海岸工程、环境和资源评价等方面提供动态的、科学的、及时有效的信息,对促进龙口市海岸带资源的可持续发展,具有十分重要的意义。1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/922地 理 研 究25卷2 研究区概况 研究区位于山东半岛北部,地理坐标:1201242E至1203542E,37
6、3154N至图1 研究区位置图Fig11The geographical locationof the study area374553N(图1)。东与蓬莱市接壤,南与招远市毗邻,西、北濒临渤海,面积约187km2。区内人口稠密,交通便利,资源丰富,社会经济优越,其中龙口港是国家一级对外开放口岸,也是国内最大的地方港口。研究区地势平坦,地面标高在20m以下。海岸线长7614km,其中人工岸线长1617km。气候属温带半湿润大陆性季风气候,年平均气温1117,平均年降水量62419mm10。西部沿海为不正规半日潮,北部沿海属正规半日潮。自东向西有黄水河、泳汶河、北马河、八里沙河四条河流注入渤海。
7、31 数据源与土地信息提取311 数据源 选用了物候相近的1984年至2004年20年间5期TM/ETM+、SPOT5数字遥感图像。地面高程数据选用了1984、1996、2004年110000的等高线和分辨率为1米的DEM。并对这些数据进行了几何校正、配准等项预处理工作1113。312 分类系统确定 考虑龙口市海岸带的典型覆盖类型,参照全国农业区划委员会1984年土地利用现状调查规程,在目视判读和实地考察基础上,确定研究区土地利用类型为7大类:即耕地-1、园地-2、林地-3、居民点及工矿用地-4、水域-5、沙地-6和滩涂-7。313 土地信息提取 除滩涂外的地类,即高潮线以上的各种地类,运用B
8、P人工神经网络(BPNN)分类器,将遥感影像的光谱信息、纹理信息、区位指数3种特征结合起来,将其分类识别1416,可从影像图上直接获取。首先采用最佳波段组合指数法对各期遥感图像进行计算,选定三个最佳的合成波段17;其次抽取纹理信息,运用离散分形布朗运动(DFBM)统计模型对全色波段数据进行分维值计算18。19,1998年以前的TM没有全色波段,可对6个多光谱波段作主成分分析,用第一主分量计算分维值。由于分维值是介于23之间的实数,为了与影像光谱值保持一致,将分维数映射到0255之间并进行离散化(取整),得到与原影像同样大小的二维分维值矩阵,即“分形图像”或纹理图像,作为纹理特征向量,以逻辑波段
9、的方式加入光谱数据中,相当于“第四波段”进行人工神经网络的分类训练;再次以距海远近为区位指数,采用类似DEM的做法,用等距线代替等高线,通过一个灰度图来表征距海的远近。即以低潮线为基线,200m为间距,向陆地依次缓冲,并赋予距离属性,然后将这些缓冲线离散化,生成一个0255范围的灰度图,像元灰度值大小表示该像元距海的远近,此图像即为“区位图像”,作为“第五波段”配合光谱信息和纹理信息共同参与神经网络训练。最后经过样本选取、网络结构设计、网络训练、分类结果的优 1994-2010 China Academic Journal Electronic Publishing House.All rig
10、hts reserved.http:/5期吴泉源 等:龙口市海岸带土地利用动态变化分析923化处理、分类精度分析及目视修正与矢量化等工序,获取土地类型分类结果。滩涂是低潮线和高潮线所围部分。其中,高潮线,即一般高潮线,通过建立潮滩解译标志,采用目视解译可以较方便地从影像图上提取。而低潮线,绝大多数情况均淹没于海水以下,没有办法从图像上直接获取,本文以水深遥感理论为指导,对海水深度与影像灰度进行相关回归分析,通过建立遥感水深模型,反演水深,再结合一定潮汐资料实现低潮线的提取。建模所用的实测水深是经过潮汐改正到卫星过境时的瞬时水深,某地瞬时水深可用下式计算:某地某时瞬时水深=当地实测水深+该地该时
11、潮高+(实测深度基准面-潮高基准面)据资料可知,研究区实测水深的实测深度基准面与潮汐表的潮高基准面相同,只需将实测水深值加上当地某一时间的潮高,即可得到该地该时的瞬时水深。对于已有验潮站的地区,直接利用我国海洋测绘部门提供的卫星过境时的潮汐预报资料进行潮位校正;没有验潮站的地区,依据该区潮位分布规律,用已有的潮位资料内插获取各个断面的瞬时潮值。本文所用的实际水深值都是经过潮汐校正后的与遥感图像同步的瞬时水深,也就是说,本文所提取的各期滩涂面积已经过潮时订正处理,具有一定可比性。314 信息叠合 考虑到矢量数据较栅格数据更加便于对点、线、面各种要素的集成操作和图形 属性的查询管理,所以确定矢量数
12、据作为本次研究的主要数据格式。首先将BP神经网络下提取的各期海岸带土地利用类型栅格图转为矢量图;然后,运用Arc/Info空间分析功能,将已完成的5期土地利用矢量数据进行“Intersect”叠加,并用“Eliminate”模块去除因空间匹配上的误差而在叠合过程中产生的一些破碎单元,完成本次研究所需的海岸带土地利用数据库。4 土地利用动态变化分析411 土地利用类型总量变化分析 分别对数据库中的5期土地利用图形数据统计分析,结果表明,龙口市海岸带土地利用类型以耕地、园地、建设用地和林地为主,其次为水域、沙地和滩涂。该地区20年来土地利用方式发生了显著变化,3种地类面积明显减少,3种地类面积明显
13、增加。在面积净减的地类中,耕地的减少量居首,达2540195hm2,林地减少335168hm2,沙地减少938180hm2;面积显著增加的地类有园地、建设用地和水域,分别增加1797119 hm2、1785199 hm2和390156 hm2。滩涂面积总量变化不大。土地利用结构的变化表现在,1984年以耕地为主,占总面积的30146%,其次为林地、建设用地和园地,分别占总面积的19159%、18168%和18165%,发展到2004年转变为以园地和建设用地为主,比例分别增加到32137%和32131%,林地和耕地次之,下降到16173%和10133%(见图版3,图2)。412 土地利用变化速度
14、分析 土地利用动态度可定量描述土地利用变化速度,分为单一土地利用动态度和综合土地利用动态度。单一土地利用动态度可表达区域一定时间范围内,某种土地利用类型变化速度的情况,公式表达为20:K=Ub-UaUa1T100%式中,Ub、Ua分别为研究期初和研究期末某一种土地利用类型的数量;T为研究步 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/924地 理 研 究25卷长。当T设定为年时,K为研究时段内某种土地利用类型的年变化率。综合土地利用动态度可描述研究时段的区域土地利用
15、变化速度,用公式表示为20:LC=ni=1L Ui2ni=1LUi1T100%式中,L Ui为监测起始时间i类土地利用类型面积,L Ui为监测时段i类土地利用类型转为非i类土地利用类型面积的绝对值,T为监测步长。当T设定为年时,L C的值就是该研究区的土地利用年变化率。表119842004年间龙口市海岸带的土地利用动态度表(%)Tab11Annual gradient of Longkou citys coastal zone land use between 1984 and 2004研究时段耕地园地林地建设用地水域沙地滩涂LC19841992-01700133-01120196-0123-
16、01170168012619921996-51636192-21137120-0181-81423130218119962000-101255189-1144414031197-8130-0195310420002004-51102115013411868195-16162-2179114719842004-31283179-016731769170-412701101153 计算研究区7种土地利用类型,4个研究时段的土地利用动态度(表1),可以看出,20年来龙口市海岸带土地利用变化速度较快,土地利用动态度为1153%,其中水域的变化速度最大,达到9170%,滩涂只有0110%。从研究时段来看
17、,具有初期和末期慢、中期快的特点,19962000年间的土地利用变化速度最快,动态度达3104%。不同地类在不同研究时段也表现出不同的变化速度,耕地和水体在19962000年间土地利用动态度最大,园地、林地、建设用地和滩涂在1992年-1996年间动态度最大,沙地则在20002004年间具有最大的土地利用动态度。413 土地利用类型相互转化分析 在Arc/info系统内,对建立的海岸带土地利用数据库统计分析,结果表明,20年来累计有55165%的土地利用方式发生了转变,达一半以上。为了更清楚地认识土地利用类型之间的相互转化情况,把握土地变化主要特征,求得19841996年和19962004年两
18、个时期的土地利用转化矩阵(表2和表3),并对主要土地利用类型转化进行空间定位,分析主要用地类型转化的来龙去脉(见图版3,图3)。结果表明:(1)两个时期耕地和沙地都是明显转出大于转入,是以转出为主的变化类型。其中耕地主要流向是园地,而且空间分布几乎遍及整个研究区,1996年之前主要发生在北部沿海地区,之后扩展到西部沿海地区。其次是流向建设用地,主要是龙口港扩建占用大量耕地;沙地在19841996年间绝大部分流向的是林地,是由于海岸的后退导致沿海防护林沙地化。在19962004年间除了32%流向林地外,还有向建设用地、园地的明显流入,主要是南山集团扩建造成的。(2)园地、建设用地和水域的转入明显
19、大于转出,是以转入为主的变化类型,其中园地的主要补给来源为耕地,耕地转向园地的面积占整个园地面积的30%以上,后期林地和建设用地对园地的补给较前期略有增大。建设用地的补给来源较广,除大量的耕地转入外,还有林地、园地以及新增用地的补给,主要是 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/5期吴泉源 等:龙口市海岸带土地利用动态变化分析925表2 龙口市海岸带19841996年土地利用转换矩阵Tab12Transition matrix of Longkou citys
20、 coastal zone land use between 1984 and 19961996年1984年耕地园地林地建设用地水域沙地滩涂消失合计转出耕地225619099917767122452100913485197010001001614130园地16816818911346214317618551536512601000100478176林地34713710417816781142251591715211314521880107811164建设用地44129431042317822531243143613301220100121109水域21637133121107180153147
21、161001193011047188沙地1114159112411169781741194432191971107109667109滩涂0100010001143174010001012591924012244110新增010101180100941740101010701660100合计转入5741381214121577137103914637177287108102178表3 龙口市海岸带19962004年土地利用转换矩阵Tab13Transition matrix of Longkou citys coastal zone land use between 1996 and 20042
22、004年1996年耕地园地林地建设用地水域沙地滩涂消失合计转出耕地10301291114149137138515154331590100010001001800199园地1701892390161152197265114125192010201000100714195林地1316232412914081382761352021752216371480100847113建设用地941232191251911652738133341391117531100100554137水域019512148301721210513015941370108010060165沙地20127106114230191
23、121127641641171634219916122602144滩涂010001002110391610103417925513360183107137新增0100010401001921130100010011060100合计转入299196177616974517314221094611334315754172龙口港和南山集团向周围扩展占用大量其他用地而形成的,龙口港向海开发也使建设用地得到一定补充。水域在19841996年间,只有林地的少量补给,而在19962004年间,其他地类补给占到总面积的78%,主要是煤矿塌陷形成的水面和沿海沙地养殖水面的增加。(3)林地的转入、转出在整个研究期
24、间都很活跃,19841996年间主要流向耕地和建设用地,19962004年间的流向主要是园地、建设用地和水域,而补给源主要是沙地,后期耕地、园地、建设用地的补给较前期有一定增加。在空间分布上,林地与建设用地的转化主要发生在屺姆岛、南山集团以及龙口港周围,而向水域的流转集中分布于煤矿塌陷区。(4)滩涂在19841996年间转入大于转出,主要是高潮线后退,使一些沙地转为滩涂,而在19962004年间则转出大于转入,主要是龙口港扩建占用大量的滩涂。414 土地利用类型变化的海岸区位分析 从“海岸区位”的角度分析海岸带各种土地利用类型转换的规律性。以研究区1984年低潮线为基线,500m为间距,依次向
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 龙口市 海岸 土地利用 动态 变化 分析 泉源
限制150内