专业正规安全的配资公司炒股平台门户网
关键词不能为空

股票配资

导航

「期货参数优化」期货怎么实现资源的优化配置

未知
admin

论文推荐| 赵诣:极化SAR参数优化与光学波谱相结合的面向对象土地覆盖分类

《测绘学报》

极化SAR参数优化与光学波谱相结合的面向对象土地覆盖分类

赵诣, 蒋弥

河海大学地球科学与工程学院, 江苏 南京 211100

收稿日期:2017-12-27;修回日期:2018-07-29

基金项目:国家重点研发计划(2018YFC0407900);国家自然科学基金(41774003);江苏省自然科学基金(BK20171432);中央高校基本科研业务费专项资金资助(2018B17714)

第一作者简介:赵诣(1995-), 男, 硕士生, 研究方向为极化SAR数据处理与应用。E-mail:zyhhu@hhu.edu.cn

通信作者:蒋弥, E-mail:mijiang@hhu.edu.cn

摘要提出一种基于极化参数优化的面向对象分类方法。该方法结合光学和SAR数据,有效提高了对地物的识别能力。本文方法的关键在于:在H-A-α分解中,使用光学影像指导SAR影像选择同质点,使其更精确地估计极化参数并结合光学波谱信息作为输入特征;使用面向对象的分类方法,仅将光学影像作为分割输入,避免SAR噪声引起的分割错误。以美国Bakersfield地区的Sentinel-1/2数据为例,确定7种地物类型,对比分析不同输入与不同分类器对分类结果的影响。研究表明,优化输入参数在纹理丰富区域能够有效提高分类精度;面向对象的分类结果更加稳定并较好地维持地表几何特征;改进分类方法较传统分类方法总体精度提高了近10%,达到92.6%。

展开剩余93%

关键词合成孔径雷达极化多光谱数据融合面向对象土地覆盖分类

Integration of SAR polarimetric parameters and multi-spectral data for object-based land cover classification

ZHAO Yi, JIANG Mi

School of Earth Sciences and Engineering, Hohai University, Nanjing 211100, China

Foundation support: The National Key Research and Development Program of China (No. 2018YFC0407900); The National Natural Science Foundation of China (No. 41774003); The Natural Science Foundation of Jiangsu Province (No. BK20171432); The Fundamental Research Funds for the Central Universities (No. 2018B17714)

First author: ZHAO Yi(1995-), male, postgraduate, majors in processing and application of PolSAR dataset, E-mail:zyhhu@hhu.edu.cn.

Corresponding author: JIANG Mi, E-mail: mijiang@hhu.edu.cn.

Abstract: An object-based approach is proposed for land cover classification using optimal polarimetric parameters. The ability to identify targets is effectively enhanced by the integration of SAR and optical images. The innovation of presented method can be summarized in the following two main points:① estimating polarimetric parameters (H-A-αdecomposition) through optical image as a driver; ② a multi-resolution segmentation based on optical image only is deployed to refine classification results. The proposed method is verified by using Sentinel-1/2 datasets over Bakersfield area, California.The results are compared against those from pixel-based SVM classification using the ground truth from the National Land Cover Database (NLCD). A detailed accuracy assessment complied for seven classes of surfaces shows that the proposed method outperforms the conventional approach by around 10%, with an overall accuracy of 92.6% over regions with rich texture.

Key words: synthetic aperture radar(SAR)polarimetricmultispectraldata fusionobject-basedland-cover classification

土地覆盖分类作为遥感领域最重要的应用分支之一[1],分类结果从宏观角度提供了地物的基本信息,在环境监测、农业管理、城市规划等应用中不可或缺[2-4]。遥感技术为大范围掌握土地覆盖情况提供了契机[5],传统的光学遥感受天气等环境因素制约无法在多云雨地区作业[6]。合成孔径雷达(synthetic aperture radar, SAR)能够全天时、全天候工作,利用SAR数据进行土地覆盖分类逐渐成为研究热点[7]。自20纪90年代开始,已有学者就SAR分类能力展开研究[8],随后大量研究表明,不同的雷达参数能够提高分类精度[9-11]。此外,当获取多时相、多极化数据时,不同的地物属性得以更加完整的描述,地物识别能力得到加强[12-14]。还有许多学者使用多源数据融合光学和SAR遥感影像更加准确地分析地物属性[15-17]。

以上研究已经较好地利用雷达数据的优点进行地物分类,但仍存在不足之处。第一,上述成果主要通过增加特征数量达到提高分类精度的目的,鲜有文献考虑提取特征对地物的识别能力及其估计的准确性。文献[18]表明较分类器而言,输入特征对分类精度影响更加显著;文献[19]使用不同滤波算法抑制噪声,研究表明滤波后影像分类精度提高约5%;文献[20]通过检验假设更加准确地估计相干性,分类结果保持良好的几何特征,总体精度提高10%。以上研究均表明优化特征参数能获取精度更高的分类结果。第二,常用的基于像素分类方法难以保持地物几何特征,椒盐现象严重,部分研究考虑分割雷达影像[21-22],然而SAR影像中大量噪声导致分割结果不可靠。最后,随着新一代Sentinel系列卫星发射,已有学者验证其分类能力[23-25],但Sentinel-1双极化特性相关研究仍不成熟。

针对上述问题,本文对输入特征与分类方法分别进行优化以期提高分类精度。首先,改进极化参数提取方法,利用光学影像良好的波谱特性指导SAR影像选取同质点,取代原有正方形窗口选点,估计极化参数;其次,使用面向对象的分类器,将光学影像分割结果应用于极化SAR分类;最后,评估Sentinel-1/2融合地表分类精度,并与现行基于像元的分类方法进行对比,论证本文方法的可行性。

1 优化的输入参数及分类方法1.1 H-A-α目标分解

极化SAR通过收发不同极化状态电磁波增加观测量,能够完整描述地物散射特性,近年来广泛应用于地物分类[26-28],这些特性往往需要通过提取极化参数得以展现[29]。文献[30-31]提出H-A-α分解方法(式(1)-式(3)),从协方差矩阵中提取散射熵H和平均散射角α,结合反熵A识别不同地物的散射机制。

(1)

(2)

(3)

为协方差矩阵特征值;αi为协方差矩阵特征向量对应的散射角。

1.2 基于光学影像的同质像素选择

为了抑制噪声,目标分解通常使用正方形滑动窗口平均邻域像素。假设中心像素和邻域像素属于同一地物,然而实际应用中滑动窗口内包含了大量异质点,直接平均影响地物统计特性,导致极化参数估计存在误差且易出现分辨率损失,地物边缘模糊等现象,在纹理丰富区域如城市、农田区域这一问题更加严重。就如何选取同质点,国内外学者从像素相似性和统计意义等不同角度提出大量算法[32-33]。这些方法计算了像素分布特性,通过假设检验从SAR影像中直接选取同质点。相比SAR影像,地物在光学影像中几何特征和纹理信息更加显著,因此本文从光学影像选取同质点指导极化参数估计。

不同地物在灰度直方图中具有不同的峰值,文献[34]首先提出一种基于灰度直方图的图像分割方法,通过迭代计算阈值区分地物,随后不同学者提出了大量的改进方法[35-36]。然而阈值分割仅通过灰度直方图将像素分为不同类型,不能直接用于中心像素的同质点选择。由于同质像素各类属性相似,因此本文使用欧氏距离作为判别标准。首先计算中心像素与邻域像素的距离,将同质点选择问题转换为阈值分割问题。为了避免“同谱异物”和“同物异谱”现象的干扰[37],试验选择小尺寸滑动窗口,在小窗口内地物种类较少,通过阈值能够准确区分。

(1) 首先在光学和SAR影像相同位置选择同样大小滑动窗口,计算中心像素与窗口内各像素间的欧氏距离

(4)

式中,dp, c表示像素p与中心像素c的距离;xpi、xci分别表示像素p和中心像素c的第i个波段观测值;N为影像总波段数量。

(2) 假设初始阈值T=T0为窗口内所有像素距离均值。

(3) 距离小于等于T的像素集合为P1,其余像素集合为P2,P1、P2t1、t2

(4) 重复迭代步骤(3)直至T不再变化,最终所有距离小于T的像素被选为同质点。

图 1为同质点选择(a)和直方图阈值分割结果(b),其中符号*和?分别代表中心像素与选中的同质点,多次模拟试验表明迭代阈值同质点误选率低于0.1。

图 1 同质点选择与直方图分割结果Fig. 1Results of homogeneous pixels selection and threshold segmentation

虽然上述理论分析能够更加准确地提取极化参数,但噪声影响、分辨率损失依然存在,若直接作为输入变量参与分类,特别是基于像素的分类方法,分类精度会因此衰减[38]。使用面向对象的分类方法能够抑制噪声影响,分类结果能够保持良好的几何特征[39]。在众多分割方法中,多尺度分割[40]根据影像灰度、纹理特征及像素间背景关系计算异质性参数f,将性质相同的像素合并为对象参与分类,分割效果优越,如式(5)、式(6)所示

(5)

(6)

式中,Δhcol、Δhsha、Δhcom、Δhsmo分别表示灰度、形状、精致度、光滑度的特征参数值;ωcol、ωsha、ωcom、ωsmo分别表示具体描述对应特征的权重值,且ωcol+ωsha=1,ωcom+ωsmo=1。相比SAR影像,光学影像灰度和纹理特征更加清晰,因此,本文仅将光学影像作为分割输入特征以避免污染纹理,SAR影像采用相同分割,避免相干斑引起的分割错误。

2 试验分析2.1 研究区及数据源

本文选择Bakersfield作为研究区域,该区域(35°15′N-35°42′N, 118°29′W-119°20′W)位于美国加州南部沙漠,是加州主要城市之一。该地区地物类型多样,以城市建筑用地与农业用地为主,纹理复杂,适合用于检验本文方法的分类能力。试验区域地理位置如图 2所示,其中红色阴影为研究区域的范围,红色线框为影像覆盖范围。

图 2 研究区域地理位置Fig. 2Geographic location of study area and data set

雷达数据采用Sentinel-1B干涉宽轨模式单视复数双极化影像,极化方式为VH+VV组合,获取于2017年7月13日,影像方位向和距离向分辨率分别为5 m和20 m。为了避免影像获取时间不同对分类结果的影响,光学数据选取成像于2017年7月7日的Sentinel-2A1C级别影像。本文试验仅使用4个10 m分辨率波段(B2(蓝)、B3(绿)、B4(红)、B8(近红外))。

在数据验证方面,本文采用2011年美国地质勘测局(USGS)发布的30 m分辨率全美土地覆盖类型数据集(National Land Cover Database,NLCD)作为真实地物类型的参考[41]。考虑到NLCD数据与试验数据时间间隔长达6年,与真实地物分布情况仍有差异,所以在训练样本和验证样本选择过程中还参考了高分辨率GoogleEarth光学影像,以确保样本获取的质量。

地物类型的定义根据美国2011年NLCD数据[41],共划分为7类:水体、湿地、灌木、草地、农田、高密度人工建筑及中低密度人工建筑。参照光学影像和NLCD数据随机层析抽样,每类地物各选取15个样方作为训练样本(图 3(a)),检验样本由每类10个样本构成(图 3(b))。

图 3 训练样本和检验样本(底图为Sentinel-2A光学影像R、G、B:B4、B3、B2)Fig. 3Training sample and validation data(base map: Sentinel-2A optical image R, G, B:B4, B3, B2)

本文使用IDAN滤波算法抑制噪声[42-43],图 4(b)中仍存在大量噪声,图 4(c)、(d)比较了使用正方形窗口和利用光学影像选取同质点两种方法的平滑效果,窗口大小均为11×11。本文提出的方法能够有效减少边界模糊现象,部分地物边缘得到加强。

图 4 两种方法平滑影像效果对比(VV)Fig. 4Comparison of images averaged using two methods (VV)

进一步分析本文方法对解译地物属性的影响。图 5为两种不同方法所得极化参数。传统方法所得结果无法保持地物边缘,白色窗口内线状地物甚至无法识别,本文提出的方法能够有效避免这一现象。在保持地物边界同时,由于减弱了异质点的影响,散射机制逐渐统一,随机性得到抑制,因此极化参数H显著降低。优化输入参数对分类结果的影响,将结合分类精度定量分析。

图 5 两种方法提取极化参数对比Fig. 5Comparison of parameters extracted using two methods

由于试验区范围较大且地物形状较为规则,经过多次试验比较最终确定分割尺度参数为300,形状因子和紧致度参数分别设为0.6和0.5。原始影像中含有1120万个像素,经过影像分割后仅剩3853个对象,运算量大幅减少。局部分割结果如图 6所示,农田、人工建筑等地物均能被准确分割。

图 6 影像分割局部结果Fig. 6The results of segmentation at local scale

为确保各类要素在试验数据中可分,对各类地物观测值进行特征统计,其光谱特征与极化参数统计如图 7(a)-(d)所示。可以看出,在光学影像中农作物、草地、灌木及湿地几类自然植被波谱特征相近,与人工建筑和水体差异显著;相反,几类自然地物由于物理结构差异显著因此通过后向散射系数与极化参数特征能够轻易区分,而部分人工建筑与草地、湿地散射机制相似,以二次散射为主,仅通过SAR参数不易区分。这种光谱特征和散射机制差异说明单个传感器不能将所有特征进行区分,而将两者相结合,能够增加特征对比度,进而提高分类精度。

图 7 样本观测值特征统计Fig. 7Statistics of all classes samples

为了获取稳定的分类结果,避免“维数灾难”,本文选择支持向量机分类器(support vector machine, SVM)[44-45],使用交叉检验优化RBF核函数,惩罚因子C与核参数g分别为1和0.125,优化精度为91%(图 8)。在分类之前采用Z转换将所有特征进行归一化处理,减小不同特征层数值范围差异对分类精度的影响。

图 8 SVM参数优化结果Fig. 8SVM parameter optimization

本文试验结果分辨率为10 m,相比30 m分辨率NLCD数据,地物细节特征更加丰富,能够准确识别尺寸较小及不连续分布的地物。NLCD数据与试验影像获取时间相隔6年,城市范围扩张、植被面积减少等变化对比光学影像能够明显辨别。为强调本文分类方法的优势,选取纹理丰富的城市区域(图 9(a))和湿地区域(图 9(b))为测试区。

图 9 不同输入特征和不同分类方法组合的分类结果比较Fig. 9Land cover classification maps obtained with different input variables

定性分析,在两个子区域,面向对象的分类结果更加接近地物在光学影像中的分布,相比基于像素的分类结果,地物几何特征更加规则。考虑输入参数对分类结果的影响,从图 9中可以看出,两种数据源所得结果并无明显差异,仅有部分中低密度人工建筑和湿地分类存在分歧。对比光学影像,改进参数所得结果更加合理。这一现象可以归因于影像分割时,已将同质像素合并作为整体分析,噪声得到抑制,因此未能体现参数质量的重要性。对极化参数改进的优势在基于像素分类时得以展现,高密度人工建筑范围估计更加准确;湿地区域内,草地和人工建筑误分显著降低,不规则形状的灌木和湿地纹理信息得以保持。

定量分析,本文选择制图精度、用户精度及Kappa系数作为评价指标,基于对象和像素的分类结果评定分别如表 1、表 2所示。

表 1基于对象分类方法精度评估Tab. 1Assessment of classification accuracy obtained with two kinds of parameters using object-based classifier

类别 传统方法提取参数 本文方法优化参数
制图精度/(%) 用户精度/(%) 制图精度/(%) 用户精度/(%)
100 59.1 100 69.1
84.7 100 85.7 100
水体 81.0 66.2 82.7 84.3
农田 100 91.7 100 93.6
草地 90.1 90.1 100 90.8
灌木 88.0 100 89.1 100
湿地 79.2 100 72.9 100
总体精度/(%) 90.1 92.6
Kappa系数 0.89 0.91

表 2基于像素分类方法精度评估Tab. 2Assessment of classification accuracy obtained with two kinds of parameters using pixel-based classifier

类别 传统方法提取参数 本文方法优化参数
制图精度/(%) 用户精度/(%) 制图精度/(%) 用户精度/(%)
60.2 65.2 85.3 86.2
85.3 53.4 80.5 91.3
水体 64.3 67.5 85.2 70.0
农田 50.4 71.1 96.3 71.6
草地 90.7 90.3 86.1 96.4
灌木 50.2 80.9 50.6 90.7
湿地 70.4 100 88.6 90.8
总体精度/(%) 79.8 86.6
Kappa系数 0.74 0.84

使用改进输入参数的面向对象分类方法获得最高总体精度92.6%。就分类方法而言:使用面向对象分类方法时,两种输入参数所得结果总体精度均提高8%左右。在各类地物中人工建筑和灌木纹理复杂,基于像素分类受异质性影响严重,因此两类地物精度最低。影像分割后,减少了异质点的影响,对地物散射特性分析更加准确,因此分类精度在人工建筑和灌木区域精度提高更加显著,这一趋势与目视解译结果相吻合,也体现了面向对象分类方法的优势。就数据源而言,如前分析,优化输入参数对面向对象分类精度没有明显改变,使用基于像素分类方法时,各地物分类精度均有所提高,在纹理丰富的城市和湿地区域,精度提高达到8%。

综合各数据的分类结果及精度评价,可以看出,使用本文方法提取的输入参数在纹理丰富区域优势显著;面向对象SVM分类方法能够有效区分离散分布的地物,提高分类的精度,保持地物几何特性。利用此法,绘制的研究区整体土地覆盖类型分布如图 10所示。

图 10 使用改进参数的面向对象分类结果Fig. 10Object-based classification result using refined parameters

相干斑噪声和有限的空间分辨率制约了合成孔径雷达获得可靠的地物分类结果。本文提出了一种利用光学影像优化极化参数的面向对象分类方法,同时结合两种数据源优点改善分类效果。本文方法的优势在于:

(1) 利用光学影像较高的光谱分辨率指导SAR影像选取同质像素,避免了噪声的干扰,较传统的正方形窗口,剔除异质点后极化参数估计更加准确,提高了对地物特性的分析能力。

(2) 结合极化特征和光谱信息,相似地物得以区分,弥补了光学和SAR影像各自的不足。同时利用光学影像分割结果,地物几何特征得以保留,在纹理复杂的城市和湿地区域分类表现卓越。

使用Sentinel 1/2数据进行验证,分类精度提高约10%,试验结果展现出了Sentinel-1 SAR和Sentinel-2光学影像在大尺度范围内土地覆盖分类方面的巨大潜力。虽然本文方法能够有效提高分类精度,保持地物边缘信息,但仅使用了一个时相的数据,未结合时变的物候信息,因而造成部分农作物与其他地物误分现象。结合多时相数据和多源数据,揭示地物在不同季节下的特性将成为今后的研究主题。

致谢:欧洲空间局通过哥白尼计划提供Sentinel系列卫星数据;美国地质调查局USGS提供全美土地类型数据

【引文格式】赵诣, 蒋弥. 极化SAR参数优化与光学波谱相结合的面向对象土地覆盖分类. 测绘学报,2019,48(5):609-617. DOI: 10.11947/j.AGCS.2019.20170746

期货怎么实现资源的优化配置

我觉得期货主要的作用是发现价格,资源的优化配置这个光期货实现不了,但在资源配置时参考期货价格对于优化应该是有帮助的。

如何解决在程序化交易中参数优化

程序化交易是反应你的交易思想的工具,首先你要建立自己的交易系统,如开平仓的条件,资金的投入量多少,风险控制标准等,有了基本的交易系统以后,设计程序量化交易条件,建立交易模型用历史数据测试优化参数,最后才能确定自己交易思想的程序化交易。

期货MACD优化参数

参数你可以自己改的,修改到自己感觉合适的好了

更适合中国期货市场的macd的参数为

谈不上适合“中国期货市场”的MACD参数。 国内期货市场50多个交易品种,每个品种的情况跟每个品种的都不一样,有对其中某个品种特别专注的,才去调。而且MACD调整也主要调快慢线,本质跟均线差不多。

做期货要看哪些技术指标

技术指标最好自己做,最好能了解技术指标的设计原理,技术指标实际就是用数学语言描述价格走势形态,大概常见技术指标也不到200个,先学会编辑技术指标,回测技术指标,优化技术指标,建立风险控制,资金管理,这时候你从事的就是量化交易了。

对做期货来说,什么指标是经常用的

传你三个指标,是我期货炒作的法宝,最重要的三个指标! 第一第二个指标RSI+SKD这两个指标是最重要的,判断价格走势的法宝! RSI>50为涨势区,在这个区域,往上弹一点给价格很大的上涨力量,可能形成顶背离,就是价格越走越高,RSI形成一顶比一顶低的走势,可能形成一-六重顶,这是我的经验! RSI<50为跌势区,在这个区域,往下弹一点给价格很大的下跌力量,可能形成底背离,就是价格越走越低,RSI形成一底比一底低的走势,可能形成一-六重底,这是我的经验! 第二个指标SKD是KD指标优化的指标,博易大师有,D线可以预测RSI的走向,RSI形成金叉D线走出很陡的走势上攻,是很好的买入信号,想判断最低点,可以看RSI在跌势区的哪个区域运行,一般刚进入最低区域不会创出最低价,一般是次几重低的区域创出最低价,这个时候买入最好!反之K线也很重要,K线的走势一般比较断,最短只有两个1 2个周期。RSI形成死叉D线走出很陡的走势下跌,是很好的卖出信号,想判断最高点,也是看RSI在涨势区的哪个区域运行,一般刚进入最高区域是不会创出最高价,最高价一般在次数重最高区域创出,这是RSI的特性决定的! 第三个指标MACD,这个指标白黄线如果不是很陡走势一般不大!这个指标会滞后而且经常纠缠在一起很多次,看你的经验判断了! 记住,价格走势由RSI左右,SKD又能左右RSI的大概走势!

程序化交易,换合约参数怎设置,用旧合约参数或优化新参数,新参数用哪个数据优化,指数,主连还是新合约

文华的可以用TRADE_OTHER('AUTO'),加载到主连合约,实现自动换月

做期货的一般看什么指标?

期货商品做@单刷@单BOLL和KDJ共振稳赚操作方法,详细看视频讲解#

期货的作用和重要性

期货的作用一、规避风险功能——套期保值1.利用期货市场和现货市场价格走势相同的原理,同事在两个市场上交易,用一个市场上生产的收益弥补另一个市场的亏损,从而锁定生产利润和生产成本,这就是企业利用期货市场进行套期保值。 2.通过期货交易,套期保值者可以把现货市场上的价格波动风险转移给寻求风险利润的投机者。 3.期货市场采用“公平、公正、公开”的集中竞价交易,产生能反映未来市场供求状况的期货价格。在成熟的期货市场作用下,期货企业往往采用期货价格作为现货交易的报价。 二、期货市场的作用 期货在宏观经济中的作用: 1.提供分散、转移价格风险的工具,有助于稳定国民经济。 2.为政府制定宏观经济政策提供参考依据。 3.有助于本国夺曲国际定价权。 4.促进本国经济的国际化。 5.有助于市场经济体系的建立和完善。 期货在微观经济中的作用: 1.锁定生产成本,实现预期利润。 2.利用期货价格信号,组织现货生产。 3.期货市场拓展现货销售和采购的渠道。 4.期货市场促进企业关注产品质量问题。期货的重要性期货的重要性主要体现在两方面,第一是发现价格,因为在期货市场中买方和买方对未来价格的预期通过提前交易达成一致,这个提前的交易综合了市场上对未来商品的供求关系,提高了资源配置效率;第二个作用是为企业规避风险提供一个手段,比如一个企业在未来要卖出一批产品,但是根据公司预判由于各种原因导致未来的商品价格会下跌,导致收入下降,那么这时企业则可以选择以高于预期价格的价格卖出期货,以达到为未来商品保值的目的。
相关阅读

期货参数优化

猜您喜欢

更多与期货参数优化相关文章