基于高分二号的城市黑臭水体动态监测

doi:10.6046/gtzyyg.2020167

基于高分二号的城市黑臭水体动态监测

胡国庆^,¹, 陈冬花¹^,², 刘聪芳³, 谢以梅¹, 刘赛赛², 李虎^,¹

1.安徽师范大学地理与旅游学院,芜湖 241000

2.滁州学院计算机与信息工程学院,滁州 239000

3.新疆师范大学地理科学与旅游学院,乌鲁木齐 830001

Dynamic monitoring of urban black-odor water bodies based on GF-2 image

HU Guoqing^,¹, CHEN Donghua¹^,², LIU Congfang³, XIE Yimei¹, LIU Saisai², LI Hu^,¹

1. College of Geography and Tourism, Anhui Normal University, Wuhu 241000, China

2. College of Computer and Information Engineering, Chuzhou University, Chuzhou 239000, China

3. College of Geography Science and Tourism,Xinjiang Normal University, Urumqi 830001, China

通讯作者: 李虎(1962-),男,教授,主要研究方向为国产卫星在资源环境领域的研究与应用。Email:lihu2881@aliyun.com。

责任编辑: 陈理

收稿日期: 2020-06-11 修回日期: 2020-07-13 网络出版日期: 2021-03-15

基金资助:

安徽省高校协同创新项目“国产高分辨率对地观测系统安徽区域综合应用示范”.  GXXT-2019-047
芜湖市重点研发项目“长江芜湖段黑臭水体遥感监测”.  2020ms1-3
安徽省科技重大专项项目“现代农业遥感监测系统构建与产业化应用”.  202003a06020002
滁州市院士助滁行动项目“基于国产卫星数据的河(湖)长制遥感监测”共同资助.  2019 1/S002

Received: 2020-06-11 Revised: 2020-07-13 Online: 2021-03-15

作者简介 About authors

胡国庆(1995-),男,硕士研究生,主要研究方向为水环境遥感。Email: 2221001119@qq.com。

摘要

目前,城市黑臭水体的遥感识别处于算法初探阶段,由于受到水深、阴影等因素的影响,在实际应用中精度较低,并且对于长时间的黑臭水体动态监测研究较少。以芜湖鸠江区为研究区,分析黑臭水体的成因及表观特征,针对单波段阈值法、波段差值法、归一化指数法和斜率指数法,基于GF-2影像进行阈值修正并评价精度,同时联合黑臭水体目视解译标志进行动态监测。结果表明: ①水体发生黑臭通常伴有颜色异常、河道淤塞、次生环境问题等特征; ②单一识别算法中差值法识别效果最好,总精度为87.5%; ③GF-2的高空间分辨率特性提升了目视解译的效率和精度,可以有效减少水深和建筑物阴影对其遥感识别的干扰,相比使用单一算法,进一步提高了动态监测的识别精度和可靠性; ④利用2014—2020年4景GF-2影像提取鸠江区主城区黑臭水体面积分别为0.313 km²,0.152 km²,0.069 km²和0.008 km²。结果显示,鸠江区水体黑臭现象正逐步得到改善,但神山公园水系黑臭现象依然严重。

关键词： 高分二号 ; 黑臭水体 ; 动态监测 ; 芜湖

Abstract

At present, the remote sensing identification of urban black-odor water bodies is in the preliminary stage of algorithm; due to the influence of water depth, shadow and other factors, the accuracy is low in practical applications, and there is little research on the long-term dynamic monitoring of black-odor water bodies. In this study, the Jiujiang District of Wuhu was chosen as a research area to analyze the causes and apparent characteristics of black-odor water bodies. For single-band threshold method, band difference method, normalized index method and slope index method, threshold correction was performed based on GF-2 images, the accuracy was evaluated, combined with the visual interpretation of the black-odor water bodies for dynamic monitoring at the same time. The results are as follows: ① The occurrence of black and odor in the water body is usually accompanied by features such as color abnormality, river siltation, and secondary environmental problems. ②The band difference method has the best recognition effect in the single recognition algorithm, and the total accuracy is 87.5%. ③ The high spatial resolution feature of GF-2 improves the efficiency and accuracy of visual interpretation, which can effectively reduce the interference of water depth and building shadows on its remote sensing recognition; compared with the use of a single algorithm, it further improves the recognition accuracy and reliability of dynamic monitoring. ④ The four GF-2 images from 2014 to 2020 were used to extract the areas of black-odor water bodies in the main urban area of Jiujiang District, which are 0.313 km², 0.152 km², 0.069 km², and 0.008 km² respectively. The results show that the black and odor phenomenon in the water body of Jiujiang District has been gradually improved, but the black and odor phenomenon in the water system of Shenshan Park is still serious.

Keywords： GF-2 ; black-odor water bodies ; dynamic monitoring ; Wuhu

PDF (5462KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

胡国庆, 陈冬花, 刘聪芳, 谢以梅, 刘赛赛, 李虎. 基于高分二号的城市黑臭水体动态监测. 国土资源遥感[J], 2021, 33(1): 30-37 doi:10.6046/gtzyyg.2020167

HU Guoqing, CHEN Donghua, LIU Congfang, XIE Yimei, LIU Saisai, LI Hu. Dynamic monitoring of urban black-odor water bodies based on GF-2 image. REMOTE SENSING FOR LAND & RESOURCES[J], 2021, 33(1): 30-37 doi:10.6046/gtzyyg.2020167

0 引言

随着城市化的进程不断加快,居民人数的迅速增加,大量的工业废水、生活污水等排入城市河道,导致水体黑臭现象频发,不仅对城市生态环境造成恶劣影响,还影响着市容市貌以及居民的身体健康。2015年4月2日国务院颁发的《水污染防治行动计划》(“水十条”)^[1]明确提出,到2020年,地级及以上城市建成区黑臭水体均控制在10%以内,到2030年,城市建成区黑臭水体总体得到消除。卫星遥感技术具有的效率高、成本低、可长时间周期监测等优势,为环境保护^[2]和黑臭水体的动态监测提供了新思路和解决方案。需要注意的是,利用遥感手段监测水体时,只能依靠光学信号监测水体的“黑”而无法监测水体的“臭”,对于只发出臭味而颜色正常的水体,光学遥感手段是无法监测到的,而大多情况下,水体受到污染时,常常同时出现刺鼻气味以及颜色异常等现象,因此,以下统称为黑臭水体的遥感监测。

早期的黑臭水体遥感识别研究集中在对类似的“黑水现象”或“黑水团”的研究,Nichol^[3]运用Landsat TM影像成功识别了新加坡一个河口的黑水区域; Zhao等^[4]利用MODIS和SeaWiFS数据,基于443 nm的离水辐射率和CDOM吸收系数提取了2002年和2012年2年的佛罗里达礁群黑水。随着国内外学者对水体光学特性和水质分析的研究不断深入,涌现了各种遥感识别算法,靳海霞等^[5]利用高分二号(GF-2)卫星融合影像对北京市城镇水体进行水体岸线提取和水质参数反演,利用水质参数浓度构建综合营养状态指数来判定河流黑臭情况; 温爽等^[6,7]根据黑臭水体光学特性提出了波段比值法,并验证了其在南京市黑臭水体识别中的可靠性; 姚月等^[8,9]提出用瑞利散射校正反射率代替遥感反射率,在波段比值法的基础上提出了一种改进后的归一化比值模型(black and odorous water index,BOI); 七珂珂等^[10,11]提出了一种城市黑臭水体遥感分级指数(black and odorous water classification index,BOCI),用于重度黑臭水体与轻度黑臭水体的分级识别,并在此基础上得到了2015—2018年沈阳市黑臭水体动态变化过程。

以上学者的研究表明目前城市黑臭水体的遥感识别尚处于算法初探阶段,并且集中在对城市主要干流黑臭水体的提取,而对黑臭现象经常出现的细小河道沟渠的研究较少,不可避免出现漏判现象。此外,大多学者利用单一算法监测黑臭水体动态变化,由于容易受到影像质量、拍摄时间等因素干扰,无法验证其在历史影像中的精度,在实际应用中可靠性较低。

针对当前黑臭水体遥感识别及其动态监测中的不足之处,本文结合黑臭水体的光学特征和表观特征,利用GF-2卫星的高空间分辨率特性,使得提取城市细小河道、黑臭水体的目视解译更为准确,并以长江流域芜湖段鸠江区为研究区,根据2016年2月5日政府公布的《全国地级及以上城市黑臭水体名单》^[12](以下简称《名单》),获取同时期GF-2影像,针对单波段阈值法、波段差值法、归一化指数法和斜率指数法进行阈值修正和精度评价得到最优算法,并在此基础上利用2014—2020年GF-2影像,综合利用目视解译标志得到黑臭水体的分布及变化,探索基于光谱特征和表观特征相结合的识别方法用于长时间黑臭水体动态监测的可靠性。

1 研究区概况与数据源

1.1 研究区概况

研究区位于安徽省东南部的芜湖市,地理坐标介于E117°40'~118°44',N30°19'~31°34'之间,流经芜湖市的河流及芜湖市境内河流,包括青弋江、漳河、水阳江、裕溪河和扁担河等,以长江芜湖段为主干构成了一个较为完整的水系。芜湖市下辖鸠江区跨长江两岸,水系发达、湖塘密布,经济发展迅速。根据《名单》和《安徽省水污染防治工作方案》^[13]公布的黑臭水体清单,芜湖市黑臭河段79条,其中鸠江区有弋江站主沟,上、下新塘水系,火石埂东和西沟等25条黑臭河段,占全市的31.6%,是黑臭水体分布的典型区域。

1.2 高分影像数据及预处理

GF-2卫星^[14]是我国自主研制的首颗亚米级民用光学卫星,回归周期69 d,搭载2台PMS(全色/多光谱),左右相机幅宽为45 km,侧摆能力±35°,可获取的波谱范围如表1所示。

表1 GF-2 载荷参数

Tab.1 GF-2 load parameters

载荷	谱段号	谱段范围/ $μm$	空间分辨率/m
全色多光谱相机	1	0.45~0.90	1
	2	0.45~0.52	4
	3	0.52~0.59
	4	0.63~0.69
	5	0.77~0.89

新窗口打开| 下载CSV

获取2景2016年3月2日GF-2影像,用于算法阈值修正和精度评价; 获取2014年12月14日、2016年11月3日、2018年1月1日、2020年3月19日各1景,用于鸠江区主城区黑臭水体动态监测,影像质量较好。基于ENVI5.3软件平台完成影像数据的预处理工作,主要包括正射校正、辐射定标、大气校正、影像融合、影像镶嵌等。首先,利用天地图无偏移影像,对全色和多光谱影像分别添加控制点进行正射校正; 之后进行辐射定标和大气校正,大气校正采用ENVI软件自带的FLAASH大气校正模块; 然后,利用校正后的GF-2影像4 m多光谱和1 m全色数据,采用NNDiffuse Pan Sharpening方法进行融合,融合后的影像空间分辨率为1 m,既具有多光谱特征,又有较高的空间分辨率,如图1所示。并且校正后的影像可以真实反映地物信息及水体遥感反射率信息。采用归一化水体指数(normalized difference water index,NDWI)提取水体,掩模得到水体分布,并剔除混合像元,对主要的细小河流进行补全。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 融合前后遥感影像对比

Fig.1 Comparision of remote sensing images before and after fusion

1.3 黑臭水体目视解译标志

在GF-2影像上黑臭水体的表观特征与一般水体有明显差别(图2),根据前人研究^[15],利用以下黑臭水体的目视解译标志作为辅助判别手段可以有效提高其遥感识别的精度。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 黑臭水体与一般水体特征对比

Fig.2 Comparison of characteristics between black and odorous water body and general water body

1)水体颜色。由于生活污水、工业废水的大量排放,城市河道有机物增多,导致水中耗氧速率大于复氧速率,造成缺氧环境,厌氧微生物分解有机质产生大量的恶臭水体,使水体发黑发臭。因此一般水体与黑臭的颜色差别是判别水体类型的首要依据。如图2(a)和(b)所示,黑臭水体呈现黑色、墨绿色,而一般水体呈现绿色、浅绿色。

2)河道淤塞。河道狭窄、排水不畅,或河道封堵,水动力不足形成断头浜,如图2(c)所示,也是导致水体黑臭的重要原因,因此,河道淤塞也是黑臭水体的重要特征和判据。

3)次生环境。当含有大量氮、磷元素的生活污水、农业、工业废水排入河道,会使水体富营养化,导致水华、浮萍泛滥等现象出现,造成次生环境问题,从而引发黑臭,并伴随大量挺水植被的出现,如图2(d)所示。

4)岸边垃圾堆放。当河道两岸堆积大量生活垃圾和建筑垃圾时,若无法得到及时清理,经发酵后散发异味并产生其他污染物,在雨水冲刷等外力作用下进入河道,是造成河流黑臭的主要来源之一。如图2(e)所示,方框内为垃圾堆放处,灰黑色水体为黑臭水体。

1.4 样点分布

由于缺少2016年鸠江区水系实测数据,本文根据官方公布的黑臭水体清单定位黑臭河段,结合目视解译标志判断水体类别,在弋江站主沟,上、下新塘水系,火石埂东和西沟等地共判别并标定了24个黑臭水体样点。同时在青弋江、裕溪河、下坝支渠、青山河、清水河选取了24个一般水体样点。样点分布如图3所示。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 样点分布

Fig.3 Sample distribution

2 黑臭水体遥感识别模型修正与精度分析

2.1 模型修正

目前对于黑臭水体的遥感识别,主要是通过波段的线性组合凸显黑臭水体与一般水体的光谱特征和颜色差异,进而构建黑臭水体遥感识别模型,实现黑臭水体的识别与监测。本文尝试利用单波段阈值法、波段差值法、归一化指数法和斜率指数法分别构建识别模型。通过1.4节选取的样点,利用校正后GF-2影像水体的遥感反射率进行建模和阈值修正,对比分析各算法精度,从而选择最优算法。不同算法建模结果如图4所示,算法公式及修正前后阈值如表2。建模结果表明,4种算法均能够区分黑臭水体与一般水体,但结果各有差异。从图4(a)可以看出,对于校正后的GF-2绿光波段反射率,黑臭水体与一般水体均存在较高或较低的样点,两者值域重叠导致基于峰值反射率的单波段阈值法无法有效辨别水体类型; 从图4(b)—(d)可以看出,基于光谱曲线变化趋势的3种算法提升了对这2种水体的区分效果,使得识别精度进一步提高。对比修正前后的阈值可以看出,由于各地区黑臭水体光谱特征和轻重程度不同,同一算法在不同地区的建模结果有明显的差别,不同算法在不同地区的适用性不同,需要根据实际情况,综合评价算法精度,以得到最优算法。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 算法建模结果

Fig.4 Algorithm modeling results

表2 公式及阈值选取

Tab.2 Formula and threshold selection

算法	公式^①	阈值
算法	公式^①	修正前	修正后
单波段阈值法	$I_{1} = Rrs (G) \leq T_{1}$	0.019 0 $s r^{- 1}$	0.144 3 $s r^{- 1}$
波段差值法	$I_{2} = Rrs (G) - Rrs (B) \leq T_{2}$	0.003 6 $s r^{- 1}$	0.002 2 $s r^{- 1}$
归一化指数法	$I_{3} = \frac{Rrs (G) - Rrs (R)}{Rrs (B) + Rrs (G) + Rrs (R)} \leq T_{3}$	0.065	0.064 2
斜率指数法	$I_{4} = \frac{\| Rrs (G) - Rrs (B) \|}{Δ λ_{1}} \cdot \frac{\| Rrs (G) - Rrs (R) \|}{Δ λ_{2}} \leq T_{4}$	0.005 $s r^{- 2}$	0.011 1 $s r^{- 2}$

①: $Rrs (B), Rrs (G) 和 Rrs (R)$ 分别为经过大气校正后GF-2数据的蓝光波段、绿光波段和红光波段反射率; $Δ λ 1 和 Δ λ 2$ 分别为绿光波段与蓝光波段和红光波段波长之差; $T 1, T 2, T 3, T 4$ 为阈值。

新窗口打开| 下载CSV

2.2 精度评价

统计样点判别结果如表3所示,本文利用用户精度、产品精度和总精度评价模型精度,并利用Kappa系数检验模型识别一致性,公式如表4所示,计算结果如表5所示。从表5可以看出,波段差值法的用户精度、产品精度以及总精度在4种算法中处于中上等,总精度最高表明其对一般水体和黑臭水体均有较好的识别效果,斜率指数法次之,归一化指数法最差,单波段阈值法虽然对黑臭水体识别效果最好,但对一般水体识别效果较差。并且波段差值法的Kappa系数为0.75说明其识别结果具有高度一致性,明显高于其他算法。算法精度评价结果表明,波段差值法用于识别芜湖鸠江区黑臭水体的精度最高,将其用于黑臭水体的识别与监测可以更好地反映黑臭水体的动态变化。

表3 样点判别统计

Tab.3 Sample discrimination statistics

水质类型	单波段阈值法		波段差值法		归一化指数法		斜率指数法
水质类型	一般水体	黑臭水体	一般水体	黑臭水体	一般水体	黑臭水体	一般水体	黑臭水体
一般水体	12	0	21	3	12	3	17	3
黑臭水体	12	24	3	21	12	21	7	21

新窗口打开| 下载CSV

表4 精度评价公式

Tab.4 Accuracy evaluation formula

指标	公式^①
用户精度	$U = \frac{n_{正确识别数}}{n_{总识别数}} \times 100 %$
产品精度	$P = \frac{n_{正确识别数}}{n_{实际识别数}} \times 100 %$
总精度	$P_{0} = \frac{n_{总正确识别数}}{n_{总数}} \times 100 %$
Kappa系数	$Kappa = \frac{P_{0} - P_{e}}{1 - P_{e}}$

①: $P e$ 为黑臭水体与一般水体分别对应的“实际与预测数量的乘积”之和,再除以“样本总数的平方”。

新窗口打开| 下载CSV

表5 精度统计

Tab.5 Accuracy statistics

指标	单波段阈值法	波段差值法	归一化指数法	斜率指数法
用户精度/%	100	87.5	87.5	87.5
产品精度/%	67	87.5	64	75
总精度/%	75	87.5	68.7	79
Kappa系数	0.5	0.75	0.375	0.58

新窗口打开| 下载CSV

3 黑臭水体动态监测

鸠江区黑臭水体主要分布在长江以东,扁担河以西的主城区,该区域人口密集,河道沟渠密布,全面掌握其黑臭水体的分布情况和消长变化对于治理和监测水体环境具有重要意义。

利用获取的2014年12月14日、2016年11月3日、2018年1月1日以及2020年3月19日的GF-2影像,使用波段差值法初步提取黑臭水体,如图5所示,并联合目视解译标志判断水体是否黑臭,验证该方法在动态监测中的可靠性,并分析黑臭水体的变化情况和原因,为相关部门评价黑臭水体的治理效果提供参考。图5表明,2014年鸠江区主城区水体黑臭情况已相当严重,监测区域内不同河段、大小沟渠都出现了不同程度的黑臭现象,主要集中在保兴垾、大阳垾和弋江站主沟; 2016年末,大公沟黑臭现象消失,大阳垾黑臭水体显著减少,新增黑臭河段保兴垾九华北路支沟; 到2018年初,水体环境整体有了明显的改善,不过大阳垾、大公沟出现了水体反黑现象; 截至2020年3月,鸠江区主城区水体除神山公园、保兴垾鸠兹家苑支沟和上、下新塘水系仍有黑臭现象外,其他河段黑臭水体已全部消失。2014—2020年黑臭水体逐渐减少,这与政府的治理和监管有直接的关系,并且从GF-2影像上可以得知,火石埂及东河是通过整改河道治理黑臭(图6),而大阳垾等河道无明显变化,是通过物化或生物等方法治理和修复的。

图5-1

新窗口打开| 下载原图ZIP| 生成PPT

图5-1 单一算法识别结果

Fig.5-1 Single algorithm recognition results

图5-2

新窗口打开| 下载原图ZIP| 生成PPT

图5-2 单一算法识别结果

Fig.5-2 Single algorithm recognition results

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 东河整改前后对比

Fig.6 Comparison before and after Dong River rectification

根据单一算法识别结果,新塘及其周边河段自2014年一直被识别为黑臭水体,但通过GF-2目视解译发现,新塘主体水域面积较大,水面平整无次生环境问题,无明显黑臭水体特征(图7(a)),因此判断为一般水体。查阅相关资料^[9]获知,新塘水质较好,水体较深,在可见光波段表现为强吸收,因此在遥感影像呈现暗像元的特征,导致单一算法错判。此外,新塘周边河段2014年和2016年水面浮萍较多、水体呈现黑色,确为黑臭水体,2018年之后,浮萍消失,水面整洁,已变为一般水体,但由于此处河道狭窄,两岸高楼林立,和保兴垾鸠兹家苑支沟相同,均是由于阴影的覆盖,导致其被错判为黑臭水体(图7(b))。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 水深及阴影的干扰

Fig.7 Water depth and shadow interference

通过以上分析最终得到2014—2020年黑臭水体动态监测结果如表6,2014年、2016年、2018年、2020年黑臭水体面积分别为0.313 km²,0.152 km²,0.069 km²和0.008 km²。

表6 黑臭水体动态监测结果

Tab.6 Dynamic monitoring results of black and odorous water

编号	名称	水质情况				编号	名称	水质情况
编号	名称	2014年	2016年	2018年	2020年	编号	名称	2014年	2016年	2018年	2020年
1	水岸星城排水渠	黑臭	黑臭	一般	一般	8	保兴垾赤铸山路支沟	黑臭	黑臭	黑臭	一般
2	东河	黑臭	黑臭	黑臭	一般	9	保兴垾鸠兹家苑支沟	黑臭	黑臭	一般	一般
3	上、下新塘水系	黑臭	黑臭	一般	一般	10	保兴垾九华北路支沟	一般	黑臭	黑臭	一般
4	火石埂东沟	黑臭	黑臭	一般	一般	11	弋江站主沟	黑臭	黑臭	轻微	一般
5	火石埂西沟	黑臭	黑臭	一般	一般	12	大阳垾湿地公园	黑臭	少量黑臭	部分水体反黑	一般
6	旭日天都水系	黑臭	黑臭	黑臭	一般	13	河清路明渠	黑臭	黑臭	一般	一般
7	大公沟	黑臭	一般	部分水体反黑	一般	14	神山公园水系	黑臭	黑臭	黑臭	黑臭

新窗口打开| 下载CSV

4 结论

根据2016年2月5日公布的芜湖市鸠江区黑臭水体清单,利用单波段阈值法、波段差值法等构建遥感识别模型,联合目视解译标志用于城市黑臭水体长时间动态监测,得到以下结论:

1)单一算法中,波段差值法用于识别鸠江区黑臭水体精度最高,总精度为87.5%,Kappa系数为0.75,说明其识别结果具有高度一致性。

2)结合黑臭水体光谱特征和表观特征,将目视解译用于遥感识别算法的补充,可有效提高识别精度,消除水深、阴影的影响,进一步提高解译精度,并可获知河段黑臭成因,有利于相关部门及时制定治理方案。

3)将本文发展的黑臭水体识别方法用于2014—2019年4期GF-2遥感影像,提取黑臭水体面积从0.313 km²减少为0.008 km²,表明鸠江区主城区黑臭水体逐年减少,相关部门整治效果显著,但应注意水体的反黑反臭等问题,以达到城市水体长治久清的目的。

4)由于本文选择的样点是典型的重度黑臭水体,其判定阈值较低,在识别重度黑臭水体上效果较好,不过对于轻度黑臭水体识别效果较差,未来的研究应该考虑将黑臭水体分级研究,分别确定识别方法及阈值,进一步提高识别精度。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

国务院.

国务院关于印发水污染防治行动计划的通知

[EB/OL].(2015-04-16)[2020-06-06]. http://www.gov.cn/zhengce/content/2015-04/16/content_9613.htm.

URL [本文引用: 1]

The State Council.

Action plan for water pollution control

[EB/OL].(2015-04-16)[2020-06-06]. http://www.gov.cn/zhengce/content/2015-04/16/content_9613.htm.

URL [本文引用: 1]

[2]

赵少华, 王桥, 游代安, 等.

高分辨率卫星在环境保护领域中的应用

[J]. 国土资源遥感, 2015,27(4):1-7.doi: 10.6046/gtzyyg.2015.04.01.