高分辨率遥感影像在土地整治项目可研阶段的应用

2020-04-16耿黎娜朱兰艳徐燕黄志坚

江苏农业科学 2020年4期

耿黎娜朱兰艳徐燕黄志坚

摘要：依托于云南省昆明市寻甸县某地土地整治项目，利用高分辨率遥感数据与第二次全国土地调查数据相结合的方法为土地整治项目可行性研究阶段进行了应用研究。通过遥感影像对二调数据进行更新，获取符合项目区实际情况的基础数据，对比分析二调数据与经影像更新后数据，分析土地利用现状及限制性因素，深挖土地整治的潜力;且在可研阶段有效控制建设规模，避免可研阶段和规划设计阶段出现较大调整。通过高分辨率遥感影像对项目区耕地图斑进行分割并引入景观格局指标，分析发现，项目区耕地细碎化程度严重，且耕地分割成地块产生的非有效耕地面积较大，造成了耕地资源的浪费，降低了耕地的利用率、产出率，严重制约着农业现代化发展。结果表明，将高分辨率遥感影像应用于土地整治项目可行性研究阶段，能够为分析结果的准确性提供现实性较强的基础数据保障，从而为土地整治项目的可行性提供切实可靠的依据。

关键词：土地整治;可行性研究阶段;遥感;土地利用现状;景观格局;耕地细碎化

中图分类号： TP79;F321.1;S127文献标志码： A

文章编号：1002-1302（2020）04-0240-06

收稿日期：2018-12-16

基金项目：国家自然科学基金（编号：41261043）。

作者简介：耿黎娜（1989—），女，云南宣威人，硕士研究生，主要从事遥感监测、土地资源管理等研究。E-mail：371637891@qq.com。

通信作者：朱兰艳，硕士，副教授，主要从事测量数据处理及误差分析、GPS和GIS的质量评估、土地资源管理等研究。E-mail：kmoy2222@sina.com。

近年来，随着我国建设用地占用耕地问题的日益突出，耕地后备资源愈加不足，耕地资源安全问题受到各级政府和广大学者的高度关注。国务院在关于《全国土地整治规划（2016—2020年）》的批复中提到，要大力推进土地整治和高标准农田建设，夯实农业现代化基础[1-3]。土地整治是指對未利用、不合理利用或低效利用的土地进行整理、复垦和开发，以提高土地利用效率，是补充耕地数量、提高耕地质量、实现耕地占补平衡的重要手段[4-5]。土地整治过程包括可行性研究、规划设计、工程施工和验收等阶段。其中可行性研究阶段（以下简称可研阶段）是项目前期工作的重要内容，它从项目土地利用现状、限制性因素、整治潜力等方面分析项目的可行性，论证项目建设必要性以及建设方案的合理性，为项目的投资决策提供科学依据[6]。一般情况下，在土地整治项目中规划设计阶段不能改变可研阶段已经批复的内容，但在实际工作中，可研阶段所用的土地利用现状图等基础图件数据往往更新不及时，现势性较差，而项目区内实际土地利用方式已发生很大的变化，例如村庄的增加，导致不动工面积随之增加，使得在规划设计阶段实测所得的建设规模与可研阶段相比发生变化。随着遥感技术的普及、影像分辨率的逐渐提高[5]，高分辨率遥感影像数据越来越多地被应用在土地整治项目中，尤其是在可研阶段，为土地整治可行性研究提供现势性的土地利用数据，有利于研究人员深入挖掘土地整治潜力，从而分析项目的可行性，在项目前期对土地整治进行合理有效地管控。目前，有不少学者提出了用遥感影像辅助土地整治项目，主要将遥感影像应用于监测、监管土地整治项目实施情况方面。杨清华等利用遥感影像提取路、渠、井等相关信息，结合规划图，监测土地整治规划的执行情况[7]。李少帅等以多时相高分辨率遥感数据为数据源，提取土地整治工程地物，构建进度评价模型，评价项目进展情况，监测管理项目实施[5]。

本研究基于高分辨率遥感影像，以云南省昆明市寻甸县某土地整治项目为例，利用遥感影像对二调数据进行更新，获取符合项目区实际现状的土地利用数据，分析项目的整治潜力;提取耕地图斑，并利用遥感影像对耕地图斑进行分割，以获得更为精细的耕地地块数据，从景观格局角度定量评价耕地细碎化问题，探讨耕地细碎化对耕地利用的影响，分析通过土地整治项目优化耕地利用结构和提高耕地质量，旨在为项目区土地整治可研分析提供切实可靠的依据。

1 研究区概况和数据

1.1 研究区概况

本研究中的土地整治项目位于云南省昆明市寻甸县某乡镇，丘陵地貌，海拔为1 960～2 040 m，地形坡度为0°～15°，平均坡度为5°，地势相对较平缓，便于组织机械和人工施工。项目区自然条件良好，利于粮食作物、经济作物的种植。该项目为中、低产田地改造项目，其重心是改造集中连片、水利条件较好、耕地增加潜力较大的中、低产田地，从单纯的耕地数量增加向质量提高转变，从而达到既增加耕地有效利用面积，又提高耕地质量的目的。

1.2 数据

本研究以项目区第二次全国土地调查数据（2015年变更）、2017年4月份的空间分辨率为0.5 m 的World View-1遥感影像数据为主要数据源。从图1可以看出，该遥感影像清晰，无大面积噪声和条纹，云量少，便于对土地利用信息进行提取，使用遥感影像数据前先对其进行正射纠正、配准融合、镶嵌与裁剪等预处理，以消除误差，改善影像质量。此外，本研究还收集了与项目区土地整治工作有关的土地利用总体规划、水文地质、地理国情监测等数据资料。

2 研究方法

在ArcGIS软件平台上，利用高分辨率遥感影像采用人机交互方式对项目区二调数据进行更新，获取符合项目区实际情况的土地利用现状数据;对更新前后的数据进行叠加，分析土地利用的实际情况，研究项目区建设规模和各地类数量的变化，结合土地利用程度指标，探讨土地整治的潜力及限制性因素;利用ArcGIS软件的空间分析功能从更新后的土地利用数据中提取耕地图斑数据，并进行影像分割，从影像的光谱差异性、田埂纹理特征、是否统一耕作、耕作方向一致性等多个方面将耕地图斑分割成耕地地块，即在耕地图斑的内部再分割出更为精细的符合项目区农业生产实际情况的耕地地块，然后利用景观格局的方法，构建耕地细碎化程度评价指标，分析耕地利用限制性因素和整治潜力。

2.1 土地利用程度指标

通过对项目区土地资源的种类、数量、结构和分布进行分析，明确项目区内土地资源的整体优势和劣势，挖掘土地整治潜力，可为土地整治的方向和重点以及制定合理规划方案提供科学依据。本研究选取土地利用率和土地垦殖率2个指标进行分析。

土地利用率：U=LA×100%。（1）

式中：U为土地利用率，%;L为项目区已利用土地面积，hm2;A为项目总面积，hm2。

土地垦殖率：R=CA×100%。（2）

式中：R为土地垦殖率，%;C为项目区耕地面积，hm2;A为项目总面积，hm2。

2.2 耕地地块细碎化评价指标

传统意义上的耕地图斑是根据沟渠、道路等线状地物以及地类界线划分而成的基本单元，但在实际农业生产生活中，由于农村土地施行家庭联产承包责任制等原因，使得耕地更加分散、破碎[8-9]。耕地细碎化的土地利用格局制约着耕地规模化经营，降低了耕地的利用率，提高了农业经营成本，限制了农业生产生活规模化和现代化的推进[10-11]。很多学者采用多种方法，从多种视角对耕地细碎化问题展开了探讨和研究。薛剑等基于地理信息系统（GIS）平台和景观指数法分析现代农业和传统农业的土地利用景观格局差异，旨在为土地整治项目布局提供依据[12]。张超等基于高分辨率遥感影像，采用面向对象分类法识别田块后，利用景观格局的方法以土地整治实例项目为试验区，通过客观评价得出，经过土地整治，研究区田块形态得到了较大的提升[8]。因此，可利用景观格局分析法对土地整治项目可研阶段的耕地细碎化问题进行分析，以此考量土地整治对耕地细碎化所带来的影响，为今后土地整治项目的规划设计提供依据。

在以上述数据的基础上，本研究结合项目区耕地情况，选取耕地地块的平均面积、密度指数、形状指数以及斑块边界密度指数等4个景观格局指数指标来分析项目区内耕地细碎化程度[8，12-13]，具体过程见图2。

2.2.1 耕地地块平均面积耕地地块平均面积是最基本的空间特征，地块平均面积越大，细碎化程度越低，计算公式如下：

Ma=An。（3）

式中：Ma為地块平均面积，hm2;n为地块数，个。

2.2.2 耕地地块密度指数耕地地块密度指数是指耕地被分割成地块的破碎程度，在一定程度上反映了农业耕作习惯以及农民拥有零散地块对耕地的干扰程度，其值>0，且值越大，表示细碎化程度越高，计算公式如下：

C=nA。（4）

式中：C为地块密度指数，个/hm2。

2.2.3 耕地地块形状指数在土地整治的相关研究中，多采用耕地地块形状指数来度量耕地形状规整度，本研究中以正方形（形状指数为1）为几何形状参照，计算公式如下：

Si=0.25PiAi。（5）

式中：Si为研究区内第i个地块形状指数;Pi为第i个地块周长，m;Ai为第i个地块面积，hm2。

2.2.4 耕地地块边界密度指数耕地地块边界密度指数是指项目区内单位面积耕地地块的边界度，边界长度越大，分割程度、破碎化程度越高，同时耕地地块边界密度也反映了因地块边界的存在而造成的耕地浪费程度，计算公式如下：

B=LA。（6）

式中：B为地块边界密度指数，km/hm2;L为研究区内地块边界总长度，km。

3 结果与分析

3.1 土地利用现状分析

依据GB/T 21010—2007《土地利用现状分类》，研究区主要由耕地、园地、林地、草地、交通运输用地、水域及水利设施用地、其他用地、城镇村及工矿用地8种土地利用类型构成（表1）。

3.1.1 主要土地利用转化及现状特征分析

由表1可知，项目区总面积为576.513 1 hm2。无论是从更新前还是更新后的土地利用数据均可看出，耕地都是项目区内所占比例最高的土地利用类型，且集中连片分布。对更新前后的土地利用数据进行叠加分析发现，耕地和林地面积变化较为明显，耕地面积由496.064 2 hm2减少为443.256 8hm2，土地垦殖率随着耕地面积的减少而降低，由86.05%降低到76.89%，田坎面积也随着耕地面积减少而减少。减少的耕地主要变成了林地、草地、村庄等，这些变化多集中在项目区地势较高的西部、水库沟渠周围以及邻近村庄的区域。林地面积由0.204 6 hm2 增加到50.659 7 hm2，增加面积主要来自耕地的转化，此外还来自裸地、特殊用地等的转化。项目区总体地势虽然较为平缓，但存在零散的局部微地貌，地势凸起不利于耕作以及区内道路、水域周围的耕地疏于管理等，使得部分耕地变成草地，因此草地面积由10.634 0 hm2增加到19.932 3 hm2。草地、裸地面积的变化使得项目区土地利用率由97.88%减小为96.54%。此外道路、水域的伸缩变化、灭失、新生等使其面积也发生了变化。建设规模是指项目区内动工面积，根据二调数据统计出的建设规模为558.376 8 hm2，而经过影像更新后，项目区实际建设规模为508.209 8 hm2，建设规模减少了50.167 0 hm2，项目建设规模的减小将影响项目的规划方案以及项目总投资等。

3.1.2 耕地内部转化及现状特征分析

项目区内耕地的种类和空间分布格局也发生了明显变化，如图3所示，项目区内中部的部分旱地转变成了水田和水浇地，中部偏东和东南区域旱地大片转化成水浇地。由表1可知，项目区水田和水浇地面积分别增加了35.046 4、100.489 6 hm2，旱地面积减少了 188.343 4 hm2，除部分旱地转化为其他地类外，其余135.536 0 hm2 旱地转化成水田和水浇地，且旱地和水浇地之间的变化相当明显，更新后主要耕地类型也由旱地变成了水浇地。旱地转化成水田和水浇地的主要原因是项目区本身具有良好的地势、水热条件，当地农民通过自发开挖箐沟和修建临时喷灌设施来引水灌溉，但缺乏系统的规划和建设，虽然使得部分旱地转化成水浇地，但转化不稳定，且灌溉面积小。一般而言，水田和水浇地的单位面积产量和经济效益高于旱地，从耕地内部的变化可以发现，项目区耕地整治潜力大，通过开展土地整治项目进行合理的规划设计建设，项目区耕地的利用率能够大大提高，同时耕地质量可明显提高。

3.2 耕地细碎化分析

利用高分辨率遥感影像对耕地图斑进行识别切割后得到，耕地地块数为10 195个（图4）。

如图5所示，项目区耕地地块面积在0.004 9～0.492 0 hm2之间，平均面积为0.046 9 hm2，标准差为0.037 2 hm2，差异系数为79.4%，面积大小分布不均;面积在0.300 0 hm2之上的耕地地块数量仅占0.2%，这严重制约着机械化耕作的发展。耕地地块密度指数为21.3，即项目区内平均每1 hm2耕地面积上分割有21个耕地地块。地块面积越大，地块数量越少，本研究中面积较大的地块分布在项目区中部道路条件较好的区域，而在项目区沟渠等周围分布的耕地地块面积都普遍偏小。

耕地地块形状指数在0.9～2.8之间，平均形状指数为1.2，差异差系数为16.7%，其中形状指数近似为1.0的地块数占24.9%;随着形状指数的增大，地块数量比例呈下降趋势。由于项目区总体地势较为平缓，因此大部分地块形状较为规整，在对项目区进行实地踏勘时发现，项目區存在很多地势凸起的微地貌以及沟渠分布较为密集和杂乱的区域，它们周围的耕地地块形状均不规整。

耕地图斑分割程度越高，地块越破碎，边界密度越大。项目区地块边界长度为558.07km 耕地地块边界密度指数为1.26 km/hm2，且据实地踏勘发现，项目区内边界宽度（田埂）均在0.4 m左右，因此耕地在被分割成地块时产生了22.322 9 hm2的耕地资源浪费，即1 hm2耕地因被分割成地块而会产生非有效耕地面积0.05 hm2，因此可以看出，耕地越细碎，对耕地资源的浪费越严重，这会在一定程度上降低耕地的利用率。

在项目区东边4 km处有一片经过土地整治项目且取得良好效果的耕地，对比分析项目区和该区域耕地地块，结果如表2所示。

从表2可以看出，经过土地整治后，地块平均面积为1.292 2 hm2，是未经整治的27.55倍，按照云南省《土地开发整理工程建设标准》，根据耕作机械工作效率、田块平整度、渠道布设以及地势等因素所确定的丘陵地区地块面积标准在1 hm2以上，因此经过整治后的耕地地块适合机械化耕作，有利于提高耕地的产出率和经济效益。经过土地整治的耕地地块密度指数为0.8，地块密度指数越大，表示耕地地块细碎化程度越高，本研究中项目区的地块密度指数为21.3，说明其线状耕地地块分布十分零碎，而整治后耕地地块密度指数可以提高到0.8，说明项目区内耕地整治潜力大。经过土地平整之后，地块形状基本上都很规整。整治后的耕地地块边界密度指数为0.2 km/hm2，按边界宽度为0.4 m 来计算，则1 hm2耕地因耕地分割成地块只产生非有效耕地面积0.01 hm2，相比未整治的区域，1 hm2耕地减少了0.04 hm2的非有效耕地浪费。对比2个研究区可以看出，经过整治后的耕地地块空间形态得到十分明显的改善，地块破碎程度有所好转。

综上，项目区虽然地势较为平缓，但区内耕地处于比较严重的细碎化状态，不利于机械化耕作，这严重制约着耕地的利用率和产出率，耕地分割产生的非有效耕地造成了耕地资源的浪费，且项目区内耕地有较大的整治潜力。

4 讨论与结论

将高分辨率遥感影像应用到土地整治项目可研阶段，更新现势性较差的土地利用数据，分析发现，各地类面积和位置均发生了变化，其中耕地和林地面积变化较为明显;区内土地疏于管理，在一定程度上造成了耕地资源的浪费;区内未利用地已发生变化，不宜开发，对这些有良好环境水热条件的中、低产田地进行改造，可在保护耕地资源的同时提高耕地质量，挖掘有限耕地资源的产出潜力。项目区内村庄和道路周围的居民点数量有所增加，同时林地面积大幅度增加，使得项目区建设规模由二调数据统计的558.376 8 hm2变成了符合实际情况的508.209 8 hm2，建设规模大小影响着项目的建设方案和总投资等。

耕地内部发生了明显的变化，部分旱地变成了水田和水浇地，说明耕地整治潜力大。通过影像切割获取耕地地块单元，构建耕地细碎化评价指标，结果发现，在项目区耕地虽然集中连片分布，但细碎化程度十分严重，地块面积较小，耕地分割成地块产生的非有效耕地面积较大，造成了耕地资源的浪费，严重制约耕地的机械化耕作，降低耕地利用率，较为乐观的是由于项目区地势整体平缓，因此耕地地块形状较为规整，可通过土地整治来改善耕地地块空间格局，从而提高耕地的有效利用率和耕地质量。

将高分辨率影像应用到土地整治项目可研阶段获得现势性土地利用现状数据，可为土地整治项目可研阶段提供真实可靠的数据和分析依据，促进项目立项;还可以在可研阶段对项目范围、建设规模进行把控，从而有效避免规划设计阶段与原可行性研究的相关方案发生较大调整，有利于提高项目工程建设的经济性、合理性和可行性。

遥感技术在土地整治项目领域的应用还可以进一步拓宽到项目立项审批、范围划定、审核项目立项时基础调查工作的真实性等方面;同时可结合数字高程模型（DEM）数据辅助规划设计，监测项目施工进度和工程验收，对比土地整治项目实施前后的遥感影像，评价项目建设的成效[14]。

参考文献：

[1]杜国明，李昀，于凤荣，等. 基于遥感的2000—2009年三江平原北部耕地变化特征分析[J]. 农业工程学报，2012，28（1）：225-229.

[2]王万茂，董祚继，王群，等. 土地利用规划学[M]. 北京：科学出版社，1996.

[3]国土资源部信息中心. 中国国土资源安全状况分析报告：2003—2004[M]. 北京：中国大地出版社，2005：93-208.

[4]张兵，崔希民，赵彦博，等. 基于High-1卫星影像的土地整治遥感监测方法研究与实践[J]. 农业工程学报，2015，31（20）：225-233.

[5]李少帅，陈原，张超，等. 基于多时相高分辨率遥感影像的土

地整治项目进展评价模型研究[J]. 中国土地科学，2014，28（3）：83-88.

[6]邹连敏. 土地开发整理项目规划设计实用技术[M]. 北京：中国水利水电出版社，2011：7-8.

[7]杨清华，李景华，韩旭，等. QuickBird遥感数据在土地整理工作中的应用[J]. 国土资源遥感，2007（4）：72-75.

[8]张超，刘佳佳，陈英义，等. 土地整治区田块空间形态变化遥感监测与评价[J]. 农业机械学报，2017，48（1）：104-111.