土壤侵蚀是指人类在利用土地的过程中导致的土壤加速侵蚀，在不同土地利用条件下，土壤侵蚀强度存在较大差异；因此，在土壤侵蚀动态监测中，土地利用类型及其空间分布是重要的调查内容。随着我国高精度资源卫星的发射，高分影像已成为土地利用解译的重要数据源，大大提高了解译精度，定量评价其解译精度是准确监测土壤侵蚀动态的基础和前提。

土地利用调查是通过勘测调查手段，查清一个区域各种土地利用的类型、面积和分布状况，调查结果通常包括土地利用现状图和土地利用面积统计表。历史上，我国把土地利用调查称作“土地清丈”，所用方法主要有方田均税法、径界法、推排法和鱼鳞图册法。英国早在1930年开展了土地利用调查与制图工作，编制了全国1:25万土地利用图^[1]。目前，国内外土地利用调查方法有不同的分类方式：按工作人员是否到达现场，分为现场调查^[2]和遥感调查^[3-5]；按是否对所有地块进行全部清查，分为抽样调查和全面调查(即普查)。在实际调查中，常常将前述2种方法组合使用，如现场抽样调查和遥感普查^[4-7]。

现场抽样调查方法的调查结果与实际状况更加吻合、精度较高，但该方法是采用样本估计总体、会存在一定的误差^[5-11]。基于遥感技术的多时相、高分辨率及越来越便宜的价格，遥感普查越来越被广泛应用；但此方法存在同物异谱、异物同谱的现象以及工作人员经验参差不齐等状况，导致普查结果存在误判、错判、漏判等^[6-16]。鉴于这2种方法均有优缺点，为了分析不同方法调查结果的误差来源，评价高分影像解译的土地利用精度，笔者以怀来县为例，开展野外抽样调查和遥感普查2种方法对比分析研究，以期为全面、准确地掌握怀来县土地利用现状提供依据。

怀来县位于河北省西北部张家口市的东南，东邻北京市、西接涿鹿县、北靠赤城县。地理范围为E 115°16′~115°58′，N 40°4′~40°35′。

怀来县南北群山起伏，河谷平川居于中部，两山夹一川呈“V”型盆地，官厅水库居盆地之中。地势由盆地向南北崛起，全县平均海拔为792 m，幽州村河谷地带为最低点，海拔394 m；水口山大黑峰为海拔最高点，达1 978 m。全县总面积1 801 km²，其中河川平原、丘陵和山地面积分别占总面积的33.4%、25.0%和41.6%。怀来县是中国葡萄之乡和中国葡萄酒之乡，葡萄种植面积约占全县总面积的10%(167 km²)，年产量超过1.4万t，产业化经营率近80%，是县域经济的第一支柱产业。

怀来县地处中温带半干旱区，属温带大陆性季风气候，具有雨热同季、昼夜温差大的气候特点。年均降水396 mm，南北两山偏多，年均风速2.3 m/s。境内有永定河、桑干河、洋河、妫水河4条过境河流，汇入官厅水库；季节性二级河道14条，如东、西沙河及灵泉河、石河、东沙沟、外井沟等。

怀来县自然条件优越，山区植被类型丰富，由于受地形等多方面因素的影响，从山区到丘陵，再到河川区低洼地带，植被类型差异较大。总的分布如下：靠近村庄的区域多属草本植被，丘陵以灌丛为主要植被类型，山地则多为森林植被。森林植被主要为油松(Pinus tabulaeformis Carr)、白桦(Betula platyphylla Suk)、蒙古栎(Quercus mongoica Fisch. ex Ledeb)和辽东栎(Quercus wutaishansea Mary)等；灌丛主要以榛子(Corylus heterophylla Fisch)、胡枝子(Lespedeza bicolor Turcz)、山杏(Armeniaca sibirica (L.) Lam)和荆条(Vitex negundo L. var. heterophylla (Franch.) Rehd)为主；灌草丛植被由长芒草(Stipa bungeana Trin)、白羊草(Bothriochloa ischaemum (L.) Keng)等组成。由于人类活动日益频繁，原生植被所剩无几，平原丘陵地带尤其如此，并被灌草丛和人工植被所替代。

土地利用遥感普查选用的影像有高分1号和资源3号卫星数据。由于数据质量直观上是资源3号卫星较好，所以优先选用资源3号卫星数据，但是受时相和云量的影响，有些数据缺乏，在没有资源3号卫星数据的情况下，采用高分1号的数据进行补充。卫星数据包括全色影像和多光谱影像，其中多光谱影像为蓝、绿、红和近红外4个波段的组合，成像时间为2015年7—9月。高分1号的全色影像分辨率为2 m，多光谱影像分辨率为8 m/16 m^[7-18]；资源3号的全色影像分辨率为2.1 m，多光谱影像分辨率为5.8 m^[8]。

土地利用通过遥感解译配合野外实地调查和验证获取，遥感解译采用人机交互切割多边形的方法，以保证解译成果中图斑不重叠、无缝隙和无其他类型拓扑错误。在解译过程中，控制解译精度的比例尺为1:5 000~1:1万，最小成图图斑面积为4 mm²，条状图斑短边长度不小于1 mm^[9]。

土地利用采用第1次全国水利普查土壤侵蚀普查土地利用分类。该分类是在GB/T 21010—2007《土地利用现状分类》 ^[10]基础上，根据土壤侵蚀特点，保留了耕地、园地、林地和草地的分类特征，对其他一级类别进行适当归并，具体分类见表 1。

表 1(Table 1)

表 1 土地利用分类体系 Table 1 Classification system of land use 一级 First class 二级 Second class 含义 Meaning 耕地 Cultivated land 旱地 Dry land 指无灌溉设施，主要靠天然降水种植旱生农作物的耕地 It refers to the xerophytic crops-cultivated land with no irrigation facilities, mainly by natural precipitation. 水浇地 Irrigated cropland 指有水源保证和灌溉设施，在一般年景能正常灌溉，种植旱生农作物的耕地 It refers to the cultivated land with water supply and irrigation facilities guaranteed, which can be irrigated normally in the normal years, and the cultivated land with xerophytic crops. 园地 Garden plot 指种植果树的园地 It refers to the garden plot where fruit trees are planted. 林地 Forest land 有林地 Closed forest land 指树木郁闭度≥0.2的乔木林地 It refers to the arbor forest land with a canopy density ≥ 0.2. 其他林地 Other forest land 包括疏林地、未成林地、迹地、苗圃等林地 Including the sparse woodland, young afforested land, cut-over land, tree nursery, etc. 灌木林地 Shrub land 指灌木覆盖度≥40%的林地 It refers to the forest land with shrub coverage ≥40%. 草地 Grass land 天然牧草地 Natural grazing land 指生长草本植物为主的土地 It refers to the land of mainly growing herbs. 指以天然草本植物为主，用于放牧或割草的草地 It refers to the grass land that is mainly natural herbs used for grazing or harvesting grass. 其他草地 Other grass land 表层为土质，生长草本植物为主，不用于畜牧业的草地 It refers to the land where the surface is soil and the growth of herbaceous plants is not used for animal husbandry grassland. 其它 Other 城镇居民点、商服及公共用地 Urban settlements，commercial business，and public land 指城镇用于生活居住的房屋和用于商业、服务业以及公共设施等用地的集合 It refers to the collection land for living in the town and the commercial, service and public facilities. 独立工矿用地 Independent industrial and mining land 指主要用于工业生产、物资存放场所的土地 It refers to the land used mainly for industrial production and storage of goods and materials. 交通运输用地 Transportation land 指用于运输通行的地面路线、场站等的土地 It refers to the land used for transportation of the ground route, field station, etc. 农村居民点 Rural settlement land 指农村用于生活居住的宅基地 It refers to the homestead used for living in the countryside. 水域及水利设施用地 Water and water conservancy facilities land 指河流水面、湖泊水面、水库水面、坑塘水面等用地 It refers to water surface, lake surface, reservoir surface, ponds surface and other land. 特殊用地 Special uses land 指用于军事设施、涉外、宗教、殡葬等的土地 It refers to the land used for military facilities, foreign affairs, religion, graveyards, etc. 其他土地 Other land 指其他分类中未提到的各种土地利用，如裸土、裸地 It refers to the various land uses not mentioned in above categories, such as barren soil and barren land.

表 1 土地利用分类体系 Table 1 Classification system of land use

在开展怀来县土地利用调查时，采用野外抽样调查和遥感普查2种方法时，都遵守上述分类体系。本文仅仅针对一级土地使用类型进行比较分析研究。

野外调查采取抽样调查的方法。抽样调查单元采用第1次全国水利普查土壤侵蚀普查野外调查单元布局方法，即分层不等概系统抽样方法^[11]，在怀来县按4%密度布设野外调查单元。首先确定调查单元所在1:1万地形图中心点坐标，然后在对应1:1万地形图中心位置确定调查单元具体位置和边界。调查单元是指进行野外调查的空间范围，平原区为1 km×1 km网格，山丘区为0.2~3 km²的小流域。全县共布设70个野外调查单元(如图 1)，其中7个单元因无1:1万地形图不进行调查，因此实际调查的单元共计63个。

图 1(Figure 1)

图 1 怀来县野外抽样调查单元分布 Figure 1 Distribution map of field sampling survey in Huailai county

野外调查于2016年7月9—12日开展，历时4 d完成所有调查单元的调查工作。调查具体内容包括勾绘调查单元内的地块、将每个地块信息填入调查表、拍摄地块和景观照片。地块信息包括地块编号、土地利用类型和水土保持措施类型^[12]。

野外调查的基本步骤为：先拍摄标识照片和小流域远景照片，然后规划调查路线。按照规划线路行走，勾绘地块边界，按顺序给地块编号，标注在调查底图上，给该地块拍照，照片编号也同时标注在调查底图对应的地块上，且编号书写方向或箭头标识方向即为拍摄方向。如果地块有水土保持措施，照片应体现水土保持措施近景，最后将地块信息填写至调查表。最终完成调查表填写。

野外调查单元的总面积为25 km²，调查得到各类土地利用面积从大到小依次为：灌木林地5.8 km²、其他5.5 km²、耕地4.6 km²、园地3.9 km²、草地3.0 km²、林地2.2 km²。

遥感解译得到全县各类土地利用面积从大到小依次为：园地392.4 km²、草地332.1 km²、林地298.2 km²、耕地287.7 km²、灌木林地262.1 km²、其他228.5 km²。2种调查方法所得土地利用面积如表 2。

表 2(Table 2)

表 2 2种调查方法所得土地利用面积 Table 2 Land use area by two survey methods km2 一级 First class 二级 Second class 抽样调查面积 Area by field sampling survey 遥感解译面积 Area by remote sensing survey 耕地 Cultivated land 旱地 Dry land 3.3 183.7 水浇地 Irrigated cropland 1.3 104.0 小计 Subtotal 4.6 287.7 园地 Garden plot 小计 Subtotal 3.9 392.4 林地 Forest land 有林地 Closed forest land 1.3 292.7 其他林地 Other forest land 0.9 5.5 小计 Subtotal 2.2 298.2 灌木林地 Shrub land 小计 Subtotal 5.8 262.1 草地 Grass land 小计 Subtotal 3.0 332.1 其他 Other 城镇居民点、商服及公共用地 Urban settlements, commercial business, and public land 1.3 82.0 独立工矿用地 Independent industrial and mining land 0.7 24.4 交通运输用地 Transportation land 0.2 31.6 农村居民点 Rural settlement land 0.2 16.8 水域及水利设施用地 Water and water conservancy facilities land 2.2 72.1 其他特殊用地 Other special uses land 0.9 1.6 小计 Subtotal 5.5 228.5 总计 Total 25.0 1 801.0

表 2 2种调查方法所得土地利用面积 Table 2 Land use area by two survey methods

对于抽样方法来说，按照如下方法计算怀来县土地利用结构：首先，计算每个调查单元不同土地利用占调查单元的面积比例；其次，按照耕地、园地、林地、灌木林地、草地和其他等6个一级类，分别进行面积比例求和；再除以调查单元总数63，得到调查单元的土地利用结构。某种土地利用比例的计算公式如下：

式中：L是某种土地利用结构的比例；L_i为U_j中第i种土地利用类型的面积，km²；i为耕地、园地、林地、灌木林地、草地和其他等6种之一；U_j为第j个调查单元的面积，km²；63为调查单元总数。

抽样调查方法得到的怀来县土地利用结构分别为：耕地17.1%，园地10.7%，林地8.9%，灌木林地33.9%，草地15.5%，其他13.9%。遥感解译得到的怀来县土地利用结构分别为：耕地16%，园地21.8%，林地16.5%，灌木林地14.6%，草地18.4%，其他12.7%。

对于整个怀来县，遥感普查与野外抽样调查的土地利用结构如图 2。与野外抽样调查相比，遥感普查耕地比例偏小1.1%，园地比例偏大11.1%，林地比例偏大7.6%，灌木林地比例偏小19.3%，草地比例偏大2.9%，其他用地比例偏小1.2%。总体来看，遥感普查方法的园地比例和林地比例明显偏大，灌木林地比例明显偏小，耕地、草地和其他用地的比例相当。

图 2(Figure 2)

图 2 怀来县遥感解译与野外抽样调查土地利用比例 Figure 2 Land use percentage by remote sensing survey and field sampling survey in Huailai county

在比较分析上述2种方法的调查结果时，将野外抽样调查土地利用数据作为比较基准，比较遥感普查土地利用结果的偏差，分析这种偏差在土地利用类型及其空间分布、数量等方面的表现。

遥感普查与野外抽样调查土地利用数据的偏差，可分为3种情况：一是土地利用类型正确、面积接近，如图 3(a)；二是土地利用类型正确、面积偏差，如图 3(b)；三是土地利用类型偏差(其中也包含面积偏差)，如图 3(c)。在全部63个野外调查单元中，第1种共有30个，第2种共有15个，第3种共有18个。由此可知，对于63个野外抽样调查单元，遥感普查的土地利用解译结果基本可靠，但仍有个别存在一定偏差。在第3种情形中，遥感普查将灌木林地误判为草地的较多，如13、19、27、28、46、55、56、57、58、59、68、69号共12个调查单元。

图 3(Figure 3)

图 3 遥感普查与野外抽样调查两种方法土地利用类型和面积偏差的典型示例图 Figure 3 Typical examples of deviation in land use type and area by remote sensing survey and field sampling survey

遥感普查与野外抽样调查土地利用结果的偏差，可通过分析遥感普查土地利用图上对应于野外抽样调查单元范围内的土地利用状况得到，包括土地利用的类型、面积及分布等。

对应于野外抽样调查单元的遥感普查图上土地利用的数据采集方法为：首先，利用具有63个野外抽样调查单元边界的图层1裁切遥感普查土地利用图，得到对应于63个野外抽样调查单元的遥感普查土地利用图层2；然后，将图层1和图层2在ArcGIS中进行交集处理，生成相应的表格(表 3)；最后，以野外抽样调查土地利用数据为基准，筛选出土地利用有差异的图斑，并统计其面积及比例。

表 3(Table 3)

表 3 遥感普查解译对应野外抽样调查单元的土地利用转移结果 Table 3 Land use transfer by remote sensing survey corresponding to unit in field sampling survey 野外调查土地利用 Land use by field samping survey 遥感普查土地利用 Land use by remote sensing survey 面积比例 Area ratio/% 耕地 Cultivated land 耕地 Cultivated land 85.3 园地 Garden plot 2.0 林地 Forest land 0.1 灌木林地 Shrub land 0.5 草地 Grass land 3.7 其他土地 Other land 8.4 园地 Garden plot 耕地 Cultivated land 0.4 园地 Garden plot 91.0 林地 Forest land 1.4 灌木林地 Shrub land 1.2 草地 Grass land 2.2 其他土地 Other land 3.8 林地 Forest land 耕地 Cultivated land 2.0 园地 Garden plot 0.7 林地 Forest land 78.6 灌木林地 Shrub land 0 草地 Grass land 15.4 其他土地 Other land 3.3 灌木林地 Shrub land 耕地 Cultivated land 0.2 园地 Garden plot 0.7 林地 Forest land 0 灌木林地 Shrub land 68.2 草地 Grass land 30.5 其他土地 Other land 0.4 草地 Grass land 耕地 Cultivated land 0.06 园地 Garden plot 2.9 林地 Forest land 1.54 灌木林地 Shrub land 3.9 草地 Grass land 89.0 其他土地 Other land 2.6 其他 Other 耕地 Cultivated land 3.0 园地 Garden plot 13.4 林地 Forest land 0.2 灌木林地 Shrub land 2.5 草地 Grass land 16.7 其他土地 Other land 64.2

表 3 遥感普查解译对应野外抽样调查单元的土地利用转移结果 Table 3 Land use transfer by remote sensing survey corresponding to unit in field sampling survey

为定量说明遥感普查与野外抽样调查土地利用面积的偏差，对63个调查单元的耕地、园地、林地、灌木林地、草地和其他等6个一级土地利用类型在遥感普查中的转移进行统计。分析表 3，便可得到遥感普查解译不同于野外抽样调查土地利用的面积及比例，得出遥感普查判读正确率以及误判的土地利用类型及其面积、比例。

结果显示：耕地的解译正确率达85.3%，耕地最易错判为其他土地、错误率为8.4%；园地的正确率达91.0%，同样较易判读为其他土地、错误率为3.8%；林地的正确率为78.6%，最易混淆为草地、错误率达15.4%；灌木林地的正确率为68.2%，主要是灌木林地与草地难以区分，将灌木林地误判为草地的比率达30.5%；草地的正确率为89.0%，草地中只有3.9%的被误判为灌木林地，但错误率远低于前者；其他土地的解译正确率最低、为64.2%，最易误判为草地，错误率达16.7%。各类土地利用解译精度如图 4。

图 4(Figure 4)

图 4 各类土地利用解译精度 Figure 4 Interpretation accuracy of all land use types

总体来说，遥感普查土地利用解译的耕地、园地和草地精度较高，林地的精度次之，灌木林地的精度更低，其他土地的精度最低。在遥感解译中，其他土地主要被误判为园地和草地，而耕地、园地、林地、草地也均有被误判为其他土地，因此其他土地转移比例不大；灌木林地被误判为草地比例较大，导致灌木林地比例偏小较多。

本研究将野外抽样调查的土地利用数据作为基准(即：以野外抽样调查的土地利用结果作为真值)，对比分析了怀来县遥感普查土地利用调查结果的偏差，分析了这种偏差在土地利用类型及其空间分布、数量等方面的表现。

1) 对于怀来县而言，无论从全域的层面(全县)、还是在调查单元范围的层面，遥感普查与野外抽样调查2种方法所得的土地利用结构具有良好的一致性，尤其是对于耕地、园地和草地等3种土地利用类型，这表明两种方法相互印证了对方可以用来开展县域土地利用调查。

2) 对于抽样调查单元而言，野外抽样调查与遥感普查两者的结果对比分析，不仅可以揭示遥感普查所得土地利用面积与结构的偏差，而且可以反映土地利用类型、分布的差异。2种方法对于63个调查单元的土地利用调查结果中，具有良好一致性(土地利用类型正确、面积接近)和较好一致性(土地利用类型正确、面积偏差)的占71.4%。这不仅表明野外抽样调查与遥感普查2种方法在县域土地利用调查方面具有较好的适用性，而且说明高分遥感影像在大比例尺土地利用调查与制图中具有较高的使用价值。

3) 遥感普查对应于野外抽样调查单元的土地利用转移结果分析表明，遥感普查灌木林地的解译正确率偏低，主要是因为灌木林地被误判为草地的比例较大。导致这种偏差的主要原因可能包括：一是灌木林地和草地存在“异物同谱”现象，影像的色彩比较接近，易导致混淆和误判^[13-14]；二是在山丘区，当灌木林地和草地混交且地块较小时，边界不是很明显，会增加判读难度^[14]。而怀来县山丘区面积比例较大，导致2种类型之间误判较多。

为提高这2类土地利用遥感解译的正确率，可借鉴如下几种方式和方法进行弥补和修正：1)在解译时，可仔细分析2种土地利用类型的纹理与明暗的差异，如当影像的多光谱组合为RGB时，灌木林地呈现一簇簇颗粒状纹理、亮度较低，而草地的纹理则均匀平滑，亮度较高^[15-17]；2)增加野外实地调查，建立较多的各种土地利用类型解译标志，尤其是对影像特征易混淆的土地利用类型，保证提高土地利用解译精度^[16-20]；3)可通过多光谱影像处理技术尽可能消除异物同谱现象，提高土地利用解译精度^[17-26]。