自然资源遥感, 2021, 33(4): 98-104 doi: 10.6046/zrzyyg.2020427

技术方法

基于改进OGMRF-RC模型的SAR图像分类方法

李渊,1,2,3,4, 毋琳1,2,3,5, 戚雯雯1,2,3, 郭拯危,1,2,3, 李宁1,2,3

1.河南大学计算机与信息工程学院,开封 475004

2.河南大学河南省智能技术与应用工程技术研究中心,开封 475004

3.河南大学河南省大数据分析与处理重点实验室,开封 475004

4.商丘工学院信息与电子工程学院, 商丘 476000

5.河南大学环境与规划学院,开封 475004

A SAR image classification method based on an improved OGMRF-RC model

LI Yuan,1,2,3,4, WU Lin1,2,3,5, QI Wenwen1,2,3, GUO Zhengwei,1,2,3, LI Ning1,2,3

1. College of Computer and Information Engineering, Henan University, Kaifeng 475004, China

2. Henan Engineering Research Center of Intelligent Technology and Application, Henan University, Kaifeng 475004, China

3. Henan Key Laboratory of Big Data Analysis and Processing, Henan University, Kaifeng 475004, China

4. School of Information and Electronic Engineering, Shangqiu Institute of Technology, Shangqiu 476000, China

5. College of Environment and Planning, Henan University, Kaifeng 475004, China

通讯作者: 郭拯危(1963-),女,教授,研究生导师,研究方向为SAR遥感信息处理。Email:gzw@henu.edu.cn

责任编辑: 张仙

收稿日期: 2020-12-31   修回日期: 2021-04-23  

基金资助: 国家自然科学基金项目“基于高分辨星载SAR图像的丹江口库区生态红线保护区的水源涵养功能监测与评估”(61871175)
河南省科技攻关计划项目“众包数据辅助下基于深度学习的高分辨率土地覆盖自动分类研究”(202102210175)
“基于多极化SAR影像的小麦物候期监测研究”(212102210093)
“黄河下游游荡型河段河势雷达遥感监测研究”(212102210101)
河南省高等学校重点科研项目“河南地区冬小麦雷达极化特征及其遥感监测研究”(19A420005)
“多覆被类型地表水源涵养功能SAR反演评估研究”(21A520004)
河南省青年人才托举工程项目“基于多极化雷达影像的河南地区冬小麦长势监测研究”(2019HYTP006)

Received: 2020-12-31   Revised: 2021-04-23  

作者简介 About authors

李 渊(1996-),男,硕士研究生,主要研究方向为SAR遥感图像处理。Email: liyuan_960423@vip.henu.edu.cn

摘要

合成孔径雷达(synthetic aperture Radar, SAR)图像分类是遥感应用中的关键技术之一。针对对象高斯-马尔可夫随机场(object-based Gaussian-Markov random field,OGMRF)模型中区域类别标签对分类精度影响的问题,提出了区域类别模糊概率(regional category fuzzy probability, RCFP)标签场方法,使临界对象具有多种类别划分的可能性,避免唯一标签导致的错分类现象。该方法综合考虑区域特征与邻域特征,利用区域边缘信息和后验概率获得RCFP,并将其纳入特征场参数求解过程中,使特征场参数更加接近真实情况,从而提高SAR图像分类精度。以河南省开封市东部约1 400 km2的区域为研究区,采用Sentinel-1卫星SAR图像开展农田、建筑、水域3类地物的分类验证实验,与K-means,FCM,马尔可夫随机场和具有区域系数的OGMRF等方法相比较,所提出方法的总体分类精度达到94.16%,Kappa系数为0.895 7,在5种方法中效果最好。

关键词: SAR图像分类; 马尔可夫随机场; 特征场; 区域类别模糊概率; Sentinel-1

Abstract

The classification of synthetic aperture Radar (SAR) images is one of the key technologies in the field of remote sensing applications. To address the problem that regional class labels affect the classification accuracy in the object-based Markov random field (OMRF) model, this paper proposes the concept of regional category fuzzy probability (RCFP) label field, which can effectively avoid the misclassification caused by wrong class labels by fully considering the possible categories of a single object. The RCFP of every region can be obtained using the regional edge information and posterior probability according to the features of the region and its adjacent regions. Then it is included in the calculation of feature field parameters to make the feature field parameters highly close to the real conditions of objects. The study area is located in the eastern part of Kaifeng City, Henan Province, covering an area of about 1 400 km2. Sentinel-1 SAR images were used for the classification experiment of farmlands, buildings, and water in the study area, and the performance of the improved method in this study was compared with that of the method of K-means, fuzzy C-means (FCM), MRF, and OGMRF-RC. The experimental results show that the overall accuracy (OA) and the Kappa coefficient of the proposed method are 94.16% and 0.8957 respectively, which are higher than those of other methods.

Keywords: classification of synthetic aperture Radar (SAR) images; Markov random field (MRF); feature filed; regional category fuzzy probability (RCFP); Sentinel-1

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

本文引用格式

李渊, 毋琳, 戚雯雯, 郭拯危, 李宁. 基于改进OGMRF-RC模型的SAR图像分类方法[J]. 自然资源遥感, 2021, 33(4): 98-104 doi:10.6046/zrzyyg.2020427

LI Yuan, WU Lin, QI Wenwen, GUO Zhengwei, LI Ning. A SAR image classification method based on an improved OGMRF-RC model[J]. Remote Sensing For Natural Resources, 2021, 33(4): 98-104 doi:10.6046/zrzyyg.2020427

0 引言

合成孔径雷达(synthetic aperture Radar, SAR)图像分类是遥感应用中的关键技术之一,其分类精度直接决定着SAR图像应用的质量。常用的SAR图像分类方法有K-均值(K-means)方法[1,2],模糊C均值(fuzzy C-means, FCM)方法[3] ,支持向量机(support vector machine, SVM)方法[4] ,以及马尔可夫随机场(Markov random field, MRF)模型等。其中,MRF模型充分考虑了空间上下文信息,具有较强的抗噪能力,在近些年的研究中得到了广泛的应用[5,6,7,8,9]

MRF是一种概率图模型,分别使用特征场和标签场表示图像特征和类别。面向对象的MRF(object-based Markov random field, OMRF)模型采用超像素分割处理,按像素特征将图像划分为同质区域,作为分类的基本单元,以区域邻接图(region adjacency graph, RAG)[10] 表示各区域概率关系,并以此进行OMRF建模,将图像分类转化为区域类别最优标记问题。与传统MRF模型相比,OMRF方法用于SAR图像分类可有效利用区域对象大小、形状和边缘等信息,降低相干斑噪声对分类结果的干扰。

近年来,基于OMRF模型的图像分类方法不断涌现各种改进形式。Chen等[11] 建立了统一的MRF(unified Markov random field, UMRF)模型,将像素和区域的似然函数相乘,同时获取像素和区域信息,使得图像分类精度得到提高; Zheng等[12] 提出了一种具有区域系数的对象高斯-马尔可夫随机场模型(object-based Gaussian-Markov random field model with region coefficients, OGMRF-RC),以边缘信息和区域大小为系数,使用基于对象的线性回归方程(object-based linear regression equation, OLRE)获取相邻区域的线性关系,用于特征场建模; Xu等[13] 提出了一种极化辅助场OMRF(OMRF with polarimetric auxiliary fields, OMRF-PA)模型,采用极化指数评估OMRF模型迭代过程中的极化信息损失,结合改进的条件概率分布,有效提高了PolSAR图像分类精度。

然而,上述改进模型中,迭代计算的区域对象标签只有一个类别,这极有可能影响图像中临界区域的类型界定。本文针对该问题,提出了区域类别模糊概率(regional category fuzzy probability, RCFP)的方法,利用边缘信息和后验概率共同获取区域对象对所有可能类别的RCFP,作为OGMRF-RC模型的标签场信息,同时在特征场参数计算中考虑RCFP的影响,并将该改进模型应用于SAR图像分类中。经实验验证,与K-means,FCM和其他MRF模型等方法相比,本文方法在分类过程中以概率形式赋予各区域多类别标签,充分利用了区域对象的图像特征与空间关系信息,使得分类精度有了明显的提升。

1 OGMRF-RC模型

1.1 OLRE邻域关系

设图像I超像素分割后区域集合为 R={r1,r2,,rn},n为区域个数,各区域像素数目为 |R|={|r1|,|r2|,,|rn|}。特征场 Y={y1,y2,,yn}表示各区域特征信息,标签场 X={x1,x2,,xn}表示各区域类别标签。

区域 ri与其邻域集合 Nri的相关性参数可以表示为:

θi,j=ei,jei×rjrjNrirj,

式中: θi,j为区域 ri与其邻域集合 Nri的相关性参数; i, j为区域的编号; ei,j为区域 rirj相交边界的像素数; ei=rjNriei,j为区域 ri边界处的像素数。

若区域 ri的类别为 h, 1hk(k为总分类数目),区域 ri与其邻域集合 Nri的特征关系可用OLRE获取,表示为:

yi-μh=rjNriθi,j(yj-μh)+gh,

式中: 特征参数 μh表示类别 h的区域均值; gh为类别 h均值为零的高斯白噪声。

1.2 OGMRF-RC模型

OGMRF-RC模型是将图像分类问题转化为给定特征场条件下,求解区域类别标签最大后验概率的问题。根据Bayesian公式,区域类别标签 x^的计算公式为:

x ^=ArgmaxP(Y=y|X=x)·P(X=x),

式中: P(Y=y|X=x)为特征场中给定类别标签 x的条件概率分布; P(X=x)为标签场中给定类别标签 x的先验概率分布。

将图像分类转化为求解最小化能量的问题,则式(3)可转化为:

x*=Argmin{-log[P(Y=y|X=x)]-log[P(X=x)]}  =Argmin{Ef+El},

式中: x*为新的类别标签; Ef为特征场能量; El为标签场能量。

设待分类图像满足高斯分布,则OGMRF-RC模型的目标函数为:

x*=Argmin12log(2πσh2)+gh22σh+rjNriV(xi,xj),

式中: 特征参数 σh2为类别 h的区域方差; V(xi,xj)为势能函数,其定义为:

V(xi,xj)=-β×ei,j (xi=xj) β×ei,j (xixj),

式中 β为势能参数。

2 改进的OGMRF-RC模型

根据地理学第一定律[14] ,空间距离越近的地物,其相互影响越大。区域对象的边缘像素特征与该区域类别一致,两相邻区域类别相同的概率可由二者边缘的像素特征的差异性反映,差异越小,类别相同概率越高。例如,对如图1(a)所示的SAR图像进行超像素划分,可得如图1(b)所示的区域对象分割结果。其中,区域2,3,4属于区域1的邻域集合,区域1,3,4为相同类别,其边缘像素在特征上差异不大; 而区域1与区域2为不同类别,其边缘像素则有明显的特征差异。

图1

图1   超像素分割结果

Fig.1   Result of super-pixel segmentation


因而,图像分类时,除了区域对象本身的纹理特征外,将其与邻域区域对象的边缘特征作为分类参数之一,对区域类别的正确划分具有较高的参考价值。

OGMRF-RC模型中,标签场采用对象唯一类别标签的形式。对于不同类别交界的区域对象,一旦出现错误的类别标签将会直接影响特征场参数的计算和最终分类结果的准确性。针对此问题,本文提出依据区域对象自身的后验概率和邻域相关边缘信息,获得该区域对所有可能类别的RCFP,作为标签场,并参与特征场参数的迭代计算,使特征参数的更新考虑了图像的全局信息,利用图像空间特征不断优化并纠正每个区域对象的类别归属,获得更好的分类效果。

2.1 标签场模糊概率

基于边缘信息获取的区域 ri的类别概率可以表示为:

P1(h-ri)=rjNri,xj=hei,jrjNriei,j,

式中 P1(h-ri)为区域 ri属于类别 h的概率。

后验概率 P2(h|yi)的计算公式为:

P2(h|yi)=G(yi;θh)P(h|xNri)P(yi),

式中: P2(h|yi)为区域 ri为类别 h的后验概率; G(yi;θh)为特征场参数 θh=(μh,σh2)的高斯分布函数; P(h|xNri)为先验概率; P(yi)为常数。 G(yi;θh)P(h|xNri)可以定义为:

G(yi;θh)=12πσhexp-(yi-μh)22σh,
P(h|xNri)=1Z×exp-rjNriV(xi,xj),

式中Z为归一化常数。

区域 ri属于类别 h的模糊概率 RCFPh-ri可以表示为:

RCFPh-ri=(P1(h-ri)+P2(h|yi))/2

2.2 特征场参数求解

用RCFP作为区域的类别标签,会使区域 ri属于每个类别 h都有一个 RCFPh-ri,将其参与特征场参数计算时,则是考虑了图像的全局信息。假设已知所有区域属于类别 hRCFPh-ri,类别 h的特征场参数 (μh,σh2)的计算定义为:

μh=riRCFPh-ri·yiriRCFPh-ri,
σh=riRCFPh-ri(yi-μh)T(yi-μh)riRCFPh-ri

3 改进模型流程

改进的SAR图像分类算法流程如图2所示。首先采用像素级分类方法与图像超像素分割方法对SAR图像进行处理,用于获取区域对象分割结果及初始类别标签。然后应用改进的OGMRF-RC模型对区域对象特征进行迭代计算,不断更新标签场信息,直至满足结束条件完成对SAR图像分类,获得分类结果。算法执行结束条件为: 达到最大迭代次数T,算法结束。

图2

图2   本文方法流程图

Fig.2   Flow chart of method in this paper


4 实验与结果

4.1 研究区域及数据

经实地考察,选取河南省开封市东部约1 400 km2的区域为研究区,其位置如图3(a)所示。该区域为典型的乡村地物分布,具有大片的农田,其间分布着房屋建筑,北部有黄河主河道流经,即研究区内具有3种地物类型,分别为农田,建筑和水。实验数据见图3(b)和表1

图3

图3   研究区域与样本分布

Fig.3   Study area and sample distribution


表1   实验使用的Sentinel-1卫星SAR图像参数

Tab.1  Parameters of Sentinel-1 Satellite SAR image used in experiment

参数SAR图像
产品类型SLC
成像模式干涉宽幅
波段C
成像日期2020年2月27日
分辨率4 m×20 m
图像尺寸1 900像素×1 900像素
极化方式VV极化

新窗口打开| 下载CSV


实验数据为一景Sentinel-1卫星SAR图像,在图像覆盖区域内,选取若干农田、建筑、水域等不同类别地物的样本,进行真实地物标记,样本分布如图3(c)所示。

4.2 实验结果

为验证本文方法的有效性,选择基于像素的K-means,FCM和MRF模型,以及基于区域的OGMRF-RC模型等方法进行对比实验。实验中本文改进的OGMRF-RC方法选用MRF模型作为像素级分类方法,采用简单线性迭代聚类(simple linear iterative clustering, SLIC)[15] 方法对图像进行超像素分割,依据经验设置势能参数 β值为0.06。

图4为5种对比方法的分类结果,其中图4(a)为谷歌地图光学图像。对照真实地物分布情况,中部区域不存在水域,但图4(c)和图4(d)在该区域中存在大量农田错分为水域的情况。左上矩形区域内地物多为农田,仅有少量的建筑和水域,但图4(b)—(e)中该区域存在一些地物错分的情况。直观对比不同方法的分类结果,图4(f)所示的本文方法分类结果与真实地物图最为接近。为了更清晰地展示算法改进后的分类效果,选取了图4左上矩形区域进行放大展示,结果如图5所示。其中图5(a)为所选区域的谷歌地图光学图像,图5(b)—(f)为所选区域SAR图像的分类结果。结果显示图5(b)—(e)的分类结果中,农田和建筑存在一些错误分类,本文方法(图5(f))对选定区域的分类结果最接近地物真实情况。

图4-1

图4-1   5种对比方法的SAR图像分类结果

Fig.4-1   SAR image classification results of five methods


图4-2

图4-2   5种对比方法的SAR图像分类结果

Fig.4-2   SAR image classification results of five methods


图5

图5   选定区域的SAR图像分类结果

Fig.5   SAR image classification results of selected area


4.3 定量分析

总体精度(overall accuracy,OA)、Kappa系数和用户精度(user accuracy,UA)常用作分类精度评价指标,其计算方法分别为:

OA=MN,
UA=EF,
Kappa=OA-P1-P,
P=1N2i=1kD-(w,:)·D-(:,w),

式中: M为图像中正确分类的像素数目; N为图像总的像素数目; E为某一类别中被正确分类的像素数目; F为分类结果中该类别的像素数目; D为混淆矩阵; k为类别数; D-(w,:)为第 w行元素的和; D-(:,w)为第 w列元素的和。

5种方法的OA和Kappa定量分析结果如图6所示,对不同类别地物的UA如表2所示。图6中,3种MRF模型分类结果的OA和Kappa系数都高于FCM方法和K-means方法,这也进一步验证了MRF模型对于SAR图像分类的适用性。而本文方法由于考虑了对象多类别属性概率,获得了最高的分类精度,OA为94.16%,Kappa系数为0.895 7。由表2的UA结果可知,5种方法对农田的分类结果差别不大; 对水域的分类结果,由于FCM和MRF方法中存在将农田错分为水域的情况,因此这两种方法的UA精度较低,而其他方法都取得了较高的精度; 对于建筑的分类结果,其他方法UA处于80%左右,而本文方法UA为89.05%。综合考虑本文方法的分类效果最好。综上所述,与其他方法相比,本文方法的分类结果,对于水域、农田和建筑物等都具有较好的效果,量化指标OA和Kappa结果同样优于对比算法,证明了本文方法的有效性。

图6

图6   5种对比方法的OA和Kappa结果

Fig.6   OA and Kappa Values of five methods


表2   5种对比方法的UA结果

Tab.2  UA value of five methods

地物类型K-meansFCMMRFOGMRF-RC本文方法
农田0.934 80.953 70.952 00.954 40.959 4
建筑0.754 20.798 30.824 40.829 00.890 5
0.969 00.509 40.838 60.954 30.955 9

新窗口打开| 下载CSV


5 结论

本文提出了一种基于改进OGMRF-RC模型的SAR图像分类方法,该方法中充分考虑区域对象自身与其邻域的特征信息,提出RCFP的方法。利用区域后验概率与边缘信息共同计算RCFP,作为区域对象标签场,并纳入特征场参数计算过程,让参数的计算考虑了图像的全局信息,使SAR图像分类精度得到了有效的提高。

为验证本文方法的有效性,采用K-means,FCM,MRF,OGMRF-RC等方法开展对比实验。结果表明,本文方法对水域、农田、建筑等地物的总体分类精度最高可达0.941 6,Kappa系数最高为0.895 7,优于其他对比方法,验证了本文方法的有效性。

本文在进行特征场参数计算时,只利用了图像的灰度特征和纹理特征。在后期的研究中,可以在分类过程中增加极化等特征,从而进一步提高SAR图像的分类精度。

参考文献

范文兵, 孙志远.

基于小波域广义高斯分布的SAR图像分割算法

[J]. 计算机工程与应用, 2020,56(5):222-226.

[本文引用: 1]

Fan W B, Sun Z Y.

SAR image segmentation algorithm based on generalized gauss distribution in wavelet domain

[J]. Computer Engineering and Applications, 2020,56(5):222-226.

[本文引用: 1]

来全, 张继权, 席宇阳, .

基于三聚类中心K-Means算法的SAR船只检测方法

[J]. 东北师大学报(自然科学版), 2017,49(3):61-67.

[本文引用: 1]

Lai Q, Zhang J Q, Xi Y Y, et al.

A ship detection method for SAR imagery based on three-centroid K-Means clustering algorithm

[J]. Journal of Northeast Normal University(Nature Science Edition), 2017,49(3):61-67.

[本文引用: 1]

韩子硕, 王春平.

基于改进FCM与MRF的SAR图像分割

[J]. 系统工程与电子技术, 2019,41(8):1726-1734.

[本文引用: 1]

Hang Z S, Wang C P.

SAR image segmentation based on improved FCM and MRF

[J]. Systems Engineering and Electronics, 2019,41(8):1726-1734.

[本文引用: 1]

国贤玉, 李坤, 王志勇, .

基于SVM+SFS策略的多时相紧致极化SAR水稻精细分类

[J]. 国土资源遥感, 2018,30(4):20-27.doi: 10.6046/gtzyyg.2018.04.04.

[本文引用: 1]

Guo X Y, Li K, Wang Z Y, et al.

Fine classification of rice with multi-temporal compact polarimetric SAR based on SVM+SFS strategy

[J]. Remote Sensing for Land and Resources, 2018,30(4):20-27.doi: 10.6046/gtzyyg.2018.04.04.

[本文引用: 1]

Voisin A, Krylov V A, Moser G, et al.

Classification of very-high resolution SAR images of urban areas using copulas and texture in a hierarchical markov random field model

[J]. IEEE Geoscience and Remote Sensing Letters, 2013,10(1):96-100.

DOI:10.1109/LGRS.2012.2193869      URL     [本文引用: 1]

Delia C, Ruscino S, Abbate M, et al.

SAR image classification through information-theoretic textural features,MRF segmentation,and object-oriented learning vector quantization

[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014,7(4):1116-1126.

DOI:10.1109/JSTARS.4609443      URL     [本文引用: 1]

Shen Y, Chen J Y, Xiao L, et al.

A bilevel contextual MRF model for supervised classification of high spatial resolution remote sensing images

[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2019,12(12):360-5372.

[本文引用: 1]

Bi H X, Xu L, Cao X, et al.

Polarimetric SAR image semantic segmentation with 3D discrete wavelet transform and Markov random field

[J]. IEEE Transactions on Image Processing, 2020,29:6601-6614.

DOI:10.1109/TIP.2020.2992177      URL     [本文引用: 1]

Jiang X F, Zhang Y, Li Y, et al.

Hyperspectral image classification with transfer learning and Markov random fields

[J]. IEEE Geoscience and Remote Sensing Letters, 2020,17(3):544-548.

DOI:10.1109/LGRS.8859      URL     [本文引用: 1]

Lee S H.

Hierarchical classification of land-cover types using RAG-based merging

[C]//2007 IEEE International Geoscience and Remote Sensing Symposium, 2007:2060-2062.

[本文引用: 1]

Chen X H, Zheng C, Yao H T, et al.

Image segmentation using a unified Markov random field model

[J]. The Institution of Engineering and Technology, 2017,11(10):860-869.

[本文引用: 1]

Zheng C, Yao H T.

Segmentation for remote-sensing imagery using the object-based Gaussian-Markov random field model with region coefficients

[J]. International Journal of Remote Sensing, 2019,40(11):4441-4472.

DOI:10.1080/01431161.2018.1563841      URL     [本文引用: 1]

Xu X F, Zou B, Zhang L M.

PolSAR image classification based on object-based Markov random field with polarimetric auxiliary label field

[J]. IEEE Geoscience and Remote Sensing Letters, 2020,17(9):1558-1562.

DOI:10.1109/LGRS.8859      URL     [本文引用: 1]

朱阿兴, 闾国年, 周成虎, .

地理相似性:地理学的第三定律?

[J]. 地球信息科学学报, 2020,22(4):673-679

DOI:10.12082/dqxxkx.2020.200069      [本文引用: 1]

定律对知识的理解和传播具有极大的推动作用。本文通过对比地理学中的定律和物理学中的定律的异同之处,认为地理学的定律虽然在形式上和提出的方式上与经典的力学定律有所区别,但仍符合人们对“定律”的定义。在共识地理学定律存在的前提下,本文阐述了地理学的第一定律和第二定律中所体现出的地理学定律的特征——描述性和归纳性。基于这些对地理学定律所具有的特征的认识,本文首先介绍了Zhu等从空间推测的角度将所提出的地理学第三定律,即“地理环境越相似,地理特征越相近”,或称地理相似性定律。本文从普遍性、独立性和应用性3个方面探讨了地理相似性是否可以作为地理学的一个定律,认为地理相似性与地理学第一和第二定律相比,有着同样的普遍性,且与第一和第二定律所涉及的规律有着本质上的区别,在应用上可能解决第一和第二定律所面临的挑战,因此应可作为地理学第三定律。

Zhu A X, Lv G N, Zhou C H, et al.

Geographic similarity:Third law of geography?

[J]. Journal of Geo-information Science, 2020,22(4):673-679.

[本文引用: 1]

李智, 曲长文, 周强, .

基于SLIC超像素分割的SAR图像海陆分割算法

[J]. 雷达科学与技术, 2017,15(4):354-358

[本文引用: 1]

Li Z, Qu C W, Zhou Q, et al.

A sea-land segmentation algorithm of SAR image based on the SLIC superpixel division

[J]. Radar Science and Technology, 2017,15(4):354-358.

[本文引用: 1]

/

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编:100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发