成功案例

第三次全国土地调查内业信息提取 39 分包项目

2020-06-29 中勘 5

1 项目概述

1.1 项目目标

中国土地勘测规划院对第三次全国土地调查内业信息提取进行公开招标,北京中勘迈普科技有限公司中标该项目39分包,承担了第三次全国土地调查内业信息提取39分包项目。项目39分包目标是在国家要求的时间节点内,完成56个县的内业不一致信息提取工作,具体目标是河北省55个县区,面积3.7万平方公里;新疆维吾尔自治区若羌县,面积19.8万平方公里。

任务目标是使用甲方提供的高分辨率数字正射影像图(DOM)和土地利用调查数据库,以县区为单位,按照土地利用现状分类标准和统一的图斑提取原则,采用以人机交互为主的信息发现和提取方法逐个图斑对比分析数据库地类与DOM地物特征的一致性,根据对比核实结果,依据影像特征勾绘两者类别不一致的地类不一致、偏移和在建三类图斑,制作系列成果文件。

1.2 工作任务

1)任务范围

39分包计划任务是完成56个县的内业不一致信息提取工作,包括河北省55个县区,面积3.7万平方公里,河北省具体县区名称如下:藁城市、 故城县、 磁县、 深州市、永年县、晋州市、鸡泽县、肥乡县、桃城区、成安县、井陉矿区、巨鹿县、峰峰矿区、 武安市、阜城县、深泽县、临城县、高邑县、武邑县、任县、、裕华区、隆尧县、广宗县、辛集市、沙河市、安平县、南宫市、邯郸县、新华区、枣强县、饶阳县、内丘县、冀州市、邢台县、平山县、武强县、平乡县、井陉县、元氏县、邱县、丛台区、桥东区、景县、赞皇县、涉县、复兴区、桥西区、栾城县、长安区、威 县、邯山区、曲周县、南和县、赵县、桥西区。

新疆维吾尔自治区若羌县,面积19.8万平方公里。

2)任务内容

接收到基础资料后,首先对基础资料进行检查,主要包括以下内容:

1)土地调查数据库检查

检查土地调查数据库内容是否包含DLTB、XZQ、XZDW图层,检查数据是否可读等。

2)影像检查与处理

A.DOM检查

检查DOM影像是否存在影响内业信息提取的问题,例如,曝光、增益、掉线、偏色、纹理丢失、亮度偏暗等,若存在问题,需要按照检查反馈要求,及时将问题反馈给甲方。

B.后批次影像处理

若区县影像第一批次下发未覆盖完整,对后批次影像进行检查与处理。

3)镶嵌检查与处理

A.检查镶嵌块与土地调查数据库中行政区范围的一致性。对于不一致的情况,需要对镶嵌块范围进行处理,使其范围与数据库中行政区范围一致。

B.检查影像实际覆盖范围与镶嵌块信息的一致性。对于不一致的情况,需要对镶嵌块属性信息进行处理,使镶嵌块信息反映的情况与影像实际覆盖情况一致。

4)行政界线处理

行政界线通过土地调查数据库的XZQ矢量层合并形成。

5)基础资料问题反馈

A.对于基础资料内容缺失、损坏、基础资料内容与任务不一致等问题。

B.对于影像存在严重影响内业信息提取的问题,如严重纹理丢失、偏色等无法判读地类的情况,洪水淹没等部分区域无法判读地类的情况;对于后批次影像与前批次影像基本一致,重复下发影像等问题。

遇到以上及其他影响内业信息提取工作开展的问题,需要在接收到数据24小时内向甲方反馈。

6)DLTB层纠正

加载数据库地类图斑文件和最新DOM,目视判断数据库矢量与最新DOM是否精确配准,若整体普遍存在偏差,影响内业信息提取时,则对DLTB层进行整体几何纠正;若局部区域存在偏差,局部影响内业信息提取时,则将DLTB层中局部偏差部分图斑筛选出来,单独导出后进行局部区域纠正。

7)综合信息查询

为了准确的进行遥感影像解译,提高内业信息提取图斑的正确性,需要对该地区的综合信息进行查询,了解当地的土地利用状况及特色产业等内容。查询内容主要包括:坐落分布、地形地貌、气候条件、水文特征、规划政策、主要产业结构、农业产业结构、生态保护情况、当季自然状况等。

8)数据库地类合并

考虑到基础数据库图斑多因权属、年度变更等因素而将相同利用类型的图斑分割的比较细碎,导致数据库图斑数量增多,为了提高内业信息对比分析的效率,将数据库相邻同地类图斑合并到二级类,在ArcToolbox工具箱中,选择Data Management Tools-》Generalization-》Dissolve工具对相邻且地类编码相同的地类图斑进行合并。

9)建成区绘制

将影像特征为成规模的城市、建制镇的较为密集连片的建(构)筑物范围勾绘成建成区,勾绘时参考数据库中图斑类型为201(城市)、202(建制镇)的图斑边缘线,一般可参考与城镇相邻的环城(镇)路、护城(镇)河、护城(镇)林内侧,以及相对密集的建(构)筑物边缘等,原则上集中连片建成区范围内不能过大包含非建设用地。

10)不一致图斑勾绘

通过对比数据库与影像特征,通过复制、切割、重新勾绘等方法,提取不一致图斑。对于数据库整块图斑地类与影像特征不一致,但图斑边界与影像特征反映边界情况相符合,可沿用数据库图斑边界,通过复制方法提取不一致图斑。如果数据库图斑地类与影像特征不一致,并且无法沿用数据库图斑边界时,则采用重新勾绘的方式提取不一致图斑,重新勾绘图斑边界需严格按不一致地物边界勾绘,不得扩大、缩小勾绘范围。对于复杂的变化类型,采用综合处理方法勾绘。若数据库的多个图斑类型与影像特征存在不一致,发生多对一(即数据库为多种地类,影像特征为单一地类)或多对多(即数据库为多种地类,影像特征为多种地类)的变化,或数据库图斑错乱情况较多,则可以根据影像特征重新矢量化。

11)属性信息表达

所有内业提取图斑编号以县级行政辖区为单元,按照自上而下、从左到右的采用阿拉伯数字方式从“1”开始统一编号。基于GIS软件对提取图斑建立拓扑关系,按实施方案要求生成图斑属性表

1.3 任务完成情况

截止2019315日,我公司已提交完成56个县的批次成果、全覆盖成果和最终成果。中国土地勘测规划院已整理完39包提交的的56个县的最终成果,成果内容有:

1)正射影像图

以县级行政辖区为单位制作的遥感正射影像图共56幅。文件命名为:“县级行政代码”+数据源DOM,数据格式为:imgtif

2)影像镶嵌块

以县级行政辖区为单位制作的影像镶嵌块文件共56个,文件命名为:“县级行政代码”xq,数据格式为:shp

3)县级行政界线

以县级行政辖区为单位制作的县级行政界线共56个文件,文件命名为:“县级行政代码”xzjx,数据格式为:shp

4)内业信息提取图斑

以县级行政辖区为单位制作的内业信息提取图斑数据共56个文件,文件命名为:“县级行政代码”nytb,数据格式为:shp

5)集中连片建成区范围线

以县级行政辖区为单位制作的集中连片建成区范围线数据共56个文件,文件命名为:“县级行政代码”jcq,数据格式为:shp

6)纠正后的数据库DLTB

39分包任务范围内无需要纠正的数据库DLTB层,故无需提交纠正后的数据库DLTB层。

7)区县综合信息查询记录

以县级行政辖区为单位制作的区县综合信息查询记录文件共56个,文件命名为:××区县综合信息查询记录,数据格式为:doc

8)成果修改说明

以县级行政辖区为单位制作的成果修改说明文件共56个,文件命名为:××区县成果修改说明,数据格式为:doc

9)补充说明文件

部分县区有以县级行政辖区为单位制作的补充说明文件共43个,文件命名为:××区县补充说明文件,数据格式为:doc


1.4  主要成果

1.4.1 批次成果

截止2018年11月13日,我公司已提交完成56个县的批次成果,其中河北省邢台县、新疆若羌县涉及第二批次的成果提交,其他县区均为第一批次提交完全县全覆盖的成果,本次共涉及提交批次初检成果58次,复检成果56次。

1.4.2 批次全覆盖成果

截止2018年12月31日,我公司已提交完成56个县的全覆盖成果。

1.4.3 最终成果

截止2019年3月15日,我公司已提交完成56个县的最终成果。

最终成果资料以县(区、市)为单位按提交,成果上交前做了如下整理:

1) 最终内业全覆盖图斑整理

合并各批次全覆盖成果,形成覆盖全县的最终全覆盖成果。采用组织实施单位下发各区县单位提交的最终批次成果和检查记录。首先核实下发的最终批次成果是否正确,对于准确无误的批次成果,进行最终成果整理。整理过程:内业全覆盖图斑仅需要将原各批次全覆盖图斑层合并,保持原批次全覆盖图斑数量、位置、形状和图斑编号等内容不变,批次图斑之间无重叠、无缝隙、不重新编图斑号。

2) 最终镶嵌块文件整理

镶嵌块属性信息是后期影像整理的唯一依据,确保每块内业图斑提取使用的影像信息正确。最终镶嵌块为覆盖全县的最终镶嵌块成果。最终镶嵌块文件是在批次成果基础上的累积和更新,理论上应为最后一批次的镶嵌块。对存在以下问题的做了修正处理并提交。

3)上交成果命名与格式

最终成果以区县为单位,成果文件夹命名为“县级行政代码+县级辖区名称”,最终成果文件内容及命名如下表所示。

1-6 三调内业信息提取最终成果文件夹内容、命名与格式

内容

文件命名

格式

影像镶嵌块

县级行政代码xq

Shapefile

县级行政界线

县级行政代码xzjx

Shapefile

内业全覆盖图斑层

县级行政代码nytb

Shapefile

集中连片建成区范围

县级行政代码jcq

Shapefile

4)最终成果上交

分三级,一级目录为“省级行政名称”,二级目录为“县级行政代码+县级辖区名称”,三级为成果文件夹中各成果文件,如下图所示。

图片关键词图片关键词

1-4 提交成果文件截图


2.技术流程与方法

2.1 技术依据

2.1.1技术依据:

1)《土地调查条例》(国务院令第518号);

2)《土地调查条例实施办法》(国土资源部令第45号);

3)《国务院关于开展第三次全国土地调查的通知》(国发〔201748号);

4)《第三次全国土地调查技术规程》;

5)《第三次全国土地调查总体方案》;

6)《第三次全国土地调查实施方案》;

7)《土地基本术语》(GB/T 192312003);

8)《土地利用现状分类》(GB/T 210102017)。

2.1.2技术指标:

1)最小上图面积

内业预判为明确的固定建(构)筑物(20)、建筑JZ)、机场用地(1007)港口码头用地(1008)水工建筑用地1109)的图斑,最小上图面积400平方米,其余图斑最小上图面积1500平方米。

2)宽度指标

道路、河流、沟渠称为线性地物,线性地物宽度为上图指标,且满足最小上图面积。影像与数据库套合,新增线性地物(除铁路外)上图指标8米;对于拓宽的线地物(铁路外)拓宽后总宽度大于等于8米的,拓宽侧的上图指标大于等于3米。铁路不设宽度指标,铁路及其附属用地按其实际情况提取。

3)勾绘精度

图斑勾绘边界相对于DOM上明显同名地物的位移不得大于2倍采样间隔;受地物阴影情况影响时,勾绘边界精度根据实际情况可以适当放宽,难以分辨清楚的地物可不提取

2.2  技术流程

2.2.1 项目总体技术流程图

图片关键词

2-1 第三次全国土地调查内业信息提取工作流程图


2.2.2 总体技术流程

根据本项目工作内容及相关技术依据,本项目技术流程可以划分为以下几个阶段:准备阶段、基础资料处理阶段、内业信息提取阶段、全面质量检查阶段、成果整理与提交阶段。

1.准备工作

项目准备工作包括组建项目组、项目实施方案制定、软硬件准备、人员培训。

1)项目组组建

根据本项目工作内容及工作量,抽调满足工作需要的专业技术人员,组成第三次全国地图调查内业信息提取项目组。

图片关键词

2-2 项目组织结构图

2)实施方案制定

组织技术人员,结合项目工作内容和要求,明确工作任务,细化工作流程,制定工作计划,制定详细可行的项目实施方案,并组织专家进行评审论证,为项目的工作开展提供依据和指导。

3)软硬件准备

根据工作需要,提前调配完成项目工作所需要的台式机、笔记本、服务器、打印机、绘图仪、移动硬盘、建立局域网的相关设备等硬件,以及Erdas、Envi、ArcGIS、PhotoShop、自主研发的坐标转换软件等处理遥感影像、矢量数据的软件,为项目的顺利开展保障。

4)人员培训

项目开始前,组织项目组全体人员参加项目实施方案、相关技术依据的培训。实施过程中,提前建立影像样本库并不断更新,根据作业过程中发现的问题及时安排集中专题培训,重点解决漏提、误提图斑的问题,不断提高作业人员提取信息的质量和效率。

2.基础资料处理

项目的基础资料包括高分辨率数字正射影像图(DOM)、土地利用调查数据库和区县行政界线。基础资料检查与处理主要包括数据检查、几何纠正、数据库部分地类合并。

1)资料接收与登记

由项目组安全管理员接收甲方提供的的基础资料,检查数据硬盘是否可以正常打开,并检查资料的可用性,并填写资料接收登记表。

2)基础资料检查与处理

1)对高分辨率数字正射影像图(DOM)、土地利用调查数据库和区县行政界线进行全面的分析检查,及时发现问题并汇总反馈给甲方,对有问题的数据及时反馈给甲方做进一步处理,以便于后续工作的顺利开展。

2)针对土地利用调查数据库与DOM坐标系不一致无法套合的问题,开展非2000国家大地坐标系数据的坐标转换处理工作。

3)数据库图斑数量多、类型复杂,为了提高内业信息对比分析的效率,对相邻同类型图斑以及部分相同一级类的二级类图斑进行合并,合并规则如下表所示。

2-1 数据库合并规则

地类编码

备注

数据库原地类代码

合并后地类代码


011、012、013

01


021、022、023

02


031、032、033

03


041、042、043

04


101、102、104、105、106、107


全部保留二级类

111、112、113、114、

115、116、117、118、119


全部保留二级类

121、124、125、126、127

12

122保留二级类

201、202、203

20

204、205保留二级类

3.内业信息提取

内业信息提取工作主要包括理解技术规范,图斑勾绘、图斑属性表生成。

1)明确不一致信息提取的标准

根据甲方提供的影像图和作业区地形地貌、经济发展、土地利用情况等,建立内业信息提取类型、必须提取、禁止提取以及应该提取的图斑的样本库,以图文并茂的方式体现,明确信息提取的标准。

2)图斑勾绘

根据第三次全国土地调查工作要求,内业信息提取以下三类图斑:

1)偏移图斑

当数据库图斑与影像上地物套合偏移大于(含)8米时,提取为偏移图斑,图斑类别标注为“PY”。

2)在建图斑

在建图斑指在建道路。当数据库图斑地类为非建设用地,影像特征为动土在建道路,道路作业面宽度大于(含)8m时,提取为在建图斑,图斑类别标注为“ZJ”。

3)地类不一致图斑

除以上偏移图斑和在建图斑外,数据库地类与影像特征反映的土地利用类型不一致情况提取为地类不一致图斑,图斑类别标注为“BYZ”。

3)图斑属性表生成

基于GIS软件对提取图斑建立拓扑关系,所有内业提取图斑编号以县级行政辖区为单元,按照自上而下、从左到右的采用阿拉伯数字方式从“1”开始统一编号,统计图斑面积,填写属性信息生成图斑属性表。

4.成果整理与提交

内业信息提取完成后,以县级辖区为单位,按照内业信息提取成果管理文件夹的形式,与区县综合信息查询记录、成果修改说明、补充说明文件一起提交监理,并按照监理意见整改提交。成果提交分为批次成果提交和最终成果提交。

1)内业信息提取成果管理文件夹

以县级辖区为单位,以文件夹形式统一管理内业信息提取成果。批次任务成果管理文件夹命名采用“县级行政代码+任务批次数编号(两位数字,不足前面补0)+县级辖区名称”,如北京市朝阳区内业信息提取第一批次任务成果管理文件夹命名为“13010201长安区”。最终成果管理文件夹命名,在批次成果命名基础上取消“任务批次数编号”,所辖文件内容及命名如所示。

2-2 内业信息提取成果管理文件夹内容、命名与格式

内容

文件命名

格式

备注

正射影像图

“县级行政代码”+数据源DOM

img

注1

影像镶嵌块

“县级行政代码”xq

shp


县级行政界线

“县级行政代码”xzjx

shp


内业信息提取图斑

“县级行政代码”nytb

shp


集中连片建成区范围线

“县级行政代码”jcq

shp


纠正后的数据库DLTB层

“县级行政代码”dltb

shp

注2

注1:行政代码与数据库的行政代码一致

注2 :纠正后的数据库DLTB层只有自行纠正后的区县需要提交,没有纠正的区县无需提交。

2)区县综合信息查询记录

每个区县作业前,作业员需查询和分析作业区县的坐落分布、地形地貌特征、主要产业结构、生态保护情况、自然、社会、经济等综合信息等,便于初步宏观掌握该区县因各种因素影响导致的自然及土地利用结构与分布特征,格式为doc。

3)成果修改说明

每个县级行政辖区的批次成果,监理检查要求修改的,每次提交修改后成果时,均需提交成果修改说明,简要说明存在问题、修改情况、截图及注意事项。成果为电子格式,随批次成果提交。

4)补充说明文件

1)对于影像存在整体或局部质量较差,例如偏暗、曝光等,影响内业信息图斑提取的,对问题区域简要说明存在问题、截图;

2)对于区县综合信息查询记录中,对内业信息提取具有重要参考价值的信息,进行简要的描述;

3)对于成果提取中,发现的大面积的不一致图斑(将导致土地利用类型流向产生非常大的变化),例如耕地、园地之间等,或对于其他特殊图斑等情况,通过认真查阅和分析后,对图斑如何合理处理的,进行简要的说明、截图。

4)在内业信息提取中需要注意的其他问题。

对以上内容提交补充说明文件(格式为doc)。

5)成果提交

批次成果提交:成果资料以县(区、市)为单位按批次提交,每批次任务原则上自影像图提供之日起一周内完成。具体按任务书时间要求完成并通过内部自检后,按提交成果内容和要求在监理现场提交成果和成果交接单。每批次成果仅提交该批次新增部分,前批次已提交成果无需在后批次重复提交。

最终成果提交:各批次成果经监理检查通过后整理合并,制作平面投影坐标系、经纬度坐标系最终成果各一套提交甲方。

5.全面质量检查

质量检查工作贯穿项目的全过程。开展作业人员的自检、互检,作业组的过程检查及单位质检机构的最终检查,重点检查基础资料的处理、内业信息的提取及成果整理与提交。对质检不合格的一律返工,通过全面质量检查后的成果提交监理,确保成果质量。

2.3 基础资料检查与处理 

2.3.1 技术流程

项目的基础资料包括高分辨率数字正射影像图(DOM)、土地利用调查数据库和区县行政界线。基础资料检查与处理主要包括数据检查、数据库部分地类合并。

1)资料接收与登记

由项目组安全管理员接收甲方提供的的基础资料,检查数据硬盘是否可以正常打开,并检查资料的可用性,并填写资料接收登记表。

2)基础资料检查与处理

1)对高分辨率数字正射影像图(DOM)、土地利用调查数据库和区县行政界线进行全面的分析检查,及时发现问题并汇总反馈给甲方,对有问题的数据及时反馈给甲方做进一步处理,以便于后续工作的顺利开展。

2)针对土地利用调查数据库与DOM坐标系不一致无法套合的问题,开展非2000国家大地坐标系数据的坐标转换处理工作。

3)数据库图斑数量多、类型复杂,为了提高内业信息对比分析的效率,对相邻同类型图斑以及部分相同一级类的二级类图斑进行合并。

2.3.2 技术方法

1.数字正射影像图检查与处理

根据《第三次全国土地调查技术规程》(试行)对调查精度的要求:农村土地利用现状调查采用优于1 米分辨率覆盖全国的遥感影像资料;城镇内部土地利用现状调查,采用优于0.2 米分辨率的航空遥感影像资料。

本项目下发的数字正射影像图的主要遥感数据源有SuperView-1、高分二号、北京二号、GeoEye-1、QuickBird、WorldView、Pleiades、IKONOS、KOMPSAT2/3等优于1米分辨率的卫星影像。

质量检查采用软件自动检查、人机交互检查、人工校对检查等方法。视具体检查内容,确定采用一种或多种方法。    


2.4 数据库几何纠正

加载数据库地类图斑文件和最新DOM,目视判断数据库矢量与最新DOM是否精确配准,若整体普遍存在偏差,影响内业信息提取时,则对DLTB层进行整体几何纠正;若局部区域存在偏差,局部影响内业信息提取时,则将DLTB层中局部偏差部分图斑筛选出来,单独导出后进行局部区域纠正。


2.5 集中连片建成区绘制

将影像特征为成规模的城市、建制镇的较为密集连片的建(构)筑物范围勾绘成建成区,勾绘时参考数据库中图斑类型为201(城市)、202(建制镇)的图斑边缘线,一般可参考与城镇相邻的环城(镇)路、护城(镇)河、护城(镇)林内侧,以及相对密集的建(构)筑物边缘等,原则上集中连片建成区范围内不能过大包含非建设用地。


2.6 内业信息提取方法

2.6.1 技术流程

在进行内业信息提取工作前,对区县综合信息进行分析,根据情况有针对性的建立区县影像样本库,明确地类解译标志。将最新遥感数据与土地利用调查数据库进行叠加对比,按照统一的图斑提取原则,采用自动化提取与人机交互结合的方式提取信息,并分析偏移图斑、在建图斑、地类不一致图斑的影像特征,做好图斑类型界定并完善图斑属性表的属性信息,最终做好自查以利于成果整理与提交。

技术流程图如下图所示。


图片关键词

2-24 内业信息提取技术流程图


2.6.2 技术方法

通过对比数据库与影像特征,通过复制、切割、重新勾绘等方法,提取不一致图斑。对于数据库整块图斑地类与影像特征不一致但图斑边界与影像特征反映边界情况相符合,可沿用数据库图斑边界,通过复制方法提取不一致图斑。如果数据库图斑地类与影像特征不一致,并且无法沿用数据库图斑边界时,则采用重新勾绘的方式提取不一致图斑,重新勾绘图斑边界需严格按不一致地物边界勾绘,不得扩大缩小勾绘范围。对于复杂的变化类型,采用综合处理方法勾绘。若数据库个图斑类型影像特征存在不一致,发生多对一(即数据库为多种地类,影像特征为单一地类)或多对多(即数据库为多种地类影像特征为多种地类)的变化,或数据库图斑错乱情况较多,则可以根据影像特征重新矢量化。

2.6.3 本环节技术特色、创新及应用效果

根据内业信息提取成果质量的要求,为了保证内业信息提取的准确性,以人机交互解译提取信息为主,结合遥感影像处理软件(Erdas、ENVI等)的监督分类等方式自动提取地类,通过遥感图像解译、地类判读来降低图斑提取的漏提率,提高内业信息提取的准确性。

监督分类,又称训练分类法,用被确认类别的样本像元去识别其他未知类别像元的过程。它就是在分类之前通过目视判读和野外调查,对遥感图像上某些样区中影像地物的类别属性有了先验知识,对每一种类别选取一定数量的训练样本,计算机计算每种训练样区的统计或其他信息,同时用这些种子类别对判决函数进行训练,使其符合于对各种子类别分类的要求,随后用训练好的判决函数去对其他待分数据进行分类。使每个像元和训练样本作比较,按不同的规则将其划分到和其最相似的样本类,以此完成对整个图像的分类。

监督分类中的样本选择和分类器的选择比较关键。在样本选择时,为了更加清楚的查看地物类型,可以适当的对图像做一些增强处理,如主成分分析、最小噪声变换、波段组合等操作,便于样本的选择;分类器的选择需要根据数据源和影像的质量来选择。


3 数据分析

3.1 监测图斑统计与分析

北京中勘迈普科技有限公司承担了的39分包任务共56个县,分别为河北省55个县区,面积3.7万平方公里;新疆维吾尔自治区若羌县,面积19.8万平方公里,共计23.5万平方公里。共提取不一致图斑420039个,面积869.32万亩。

提取的不一致图斑中,占比较大的几类为01/03、SM、04。其中耕地类85232个,占比28.75%,林地类65768个,占比21.69%。水面SM占比17.06%,其他草地04类占比12.47%。

3.2 其他相关分析

以河北省平山县为例,平山县面积2643.7平方公里(3965540.6亩)。提取A类图斑16638个,B类图斑38723个,提取A类图斑个数占总图斑个数的30%。A类图斑面积184758.1亩,B类图斑3780782.5亩,提取的A类图斑占总面积的4.7%。