CN110069868B - Gnss测站非线性运动建模方法与装置 - Google Patents
Gnss测站非线性运动建模方法与装置 Download PDFInfo
- Publication number
- CN110069868B CN110069868B CN201910345177.8A CN201910345177A CN110069868B CN 110069868 B CN110069868 B CN 110069868B CN 201910345177 A CN201910345177 A CN 201910345177A CN 110069868 B CN110069868 B CN 110069868B
- Authority
- CN
- China
- Prior art keywords
- frequency part
- reconstruction
- low
- layer
- singular spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 59
- 238000001228 spectrum Methods 0.000 claims abstract description 38
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 7
- 238000010183 spectrum analysis Methods 0.000 claims description 10
- 230000001932 seasonal effect Effects 0.000 abstract description 7
- 239000011159 matrix material Substances 0.000 description 25
- 230000000737 periodic effect Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 238000012935 Averaging Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 125000003275 alpha amino acid group Chemical group 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 238000000700 time series analysis Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/12—Timing analysis or timing optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种GNSS测站非线性运动建模方法与装置,属于全球定位导航系统技术领域。该方法包括:获取GNSS测站的原始坐标时间序列;对原始坐标时间序列进行小波分解和重构,得到低频部分和多层的高频部分;对所述低频部分和每层高频部分分别进行奇异谱分解和重建;将奇异谱分解和重建后的低频部分和每层高频部分进行合成得到原始坐标时间序列的拟合序列;对所述拟合序列进行建模。该方法先利用小波分解和重构将原始坐标时间序列分解成频率成分上相对单一、平滑的低频部分和各层高频部分,然后针对各部分分别进行奇异谱分解和重建,在一定程度上降低了部分周期项如季节周期项、月周期项被当作噪声剔除的概率,从而提高了建模精度。
Description
技术领域
本发明涉及一种GNSS测站非线性运动建模方法与装置,属于全球定位导航系统技术领域。
背景技术
随着GNSS(Global Navigation Satellite System,全球卫星导航系统)观测数据的不断丰富和数据处理技术的不断发展,分析测站非线性运动规律可以获得各种地球物理现象和季节性影响,如果对这部分残差规律进行建模拟合从而消除其影响,那么可以进一步提高GPS台站坐标的精度。
尽管目前国际地球参考框架ITRF(International Terrestrial ReferenceFrame)中基准站的历元坐标和速度场已达到了毫米级,但几乎所有IGS(InternationalGPS Service)站坐标时间序列(尤其是垂直方向)都呈现显著的非线性运动趋势,且振幅可达1~2cm。因此,地面点相对于参考框架的运动无法用线性速度进行完整描述,仅基于线性速率的模型来维持参考框架会有一定的局限性,而且其精度只能达到厘米级。一般情况下,要精确表达站点坐标,通常采用两种方式描述非线性运动:①不考虑引起测站坐标非线性变化的各种物理机制,只根据坐标时间序列本身的运动趋势建模;②从产生形变的物理机理入手,分析各因素影响从而对测站坐标进行改正。目前,国内外对IGS站坐标时间序列的分析已有不少相关研究,也取得较大进展。张鹏,蒋志浩等在2007年武汉大学学报(信息科学版)第32卷第3期第251-254页,公开了一篇名称为《我国GPS跟踪站数据处理与时间序列特征分析》的文献,并在该文献中对我国IGS站坐标时间序列进行频谱分析,发现高程方向周年特性表现明显;傅彦博,孙付平等在2018年测绘学报第47卷第10期第1337-1345页,公开了一篇名称为《全球GPS测站垂向周年变化统计改正模型的建立》的文献,并在该文献中对全球461个GPS测站垂向时间序列所包含的各类周期项进行了统计,发现周年项在全球具有普遍性,其次是半周年项。另有学者证实了大气负荷、水文负载、非潮汐海洋负载等负载变化是引起GPS测站垂向位移的主要因素。
奇异谱分析(Singular Spectrum Analysis,SSA)是时间序列分析中一种常用的方法,在气候学、测量学和海洋学等领域应用广泛,它的优越性在于无需任何先验信息和假设条件,能稳定识别和强化周期信号,把最可预报的分量聚集到若干个时间序列中,以此选择若干有意义的分量重建序列,降低噪声影响,因此特别适合于分析有周期震荡的时间序列数据。然而,只根据贡献率较大的前几项来对时间序列进行重构,仅能有效地将半年及半年以上的周期项提取出来,而短周期项如季节周期项、月周期项则容易当作噪声被忽略掉。因此,需要一定方法削弱这种影响。
发明内容
本发明的目的是提供一种GNSS测站非线性运动建模方法,用以解决现有技术会忽略短周期项导致建模精度低的问题;本发明还提供一种GNSS测站非线性运动建模装置,用以解决现有技术会忽略短周期项导致建模精度低的问题。
为实现上述目的,本发明提供了一种GNSS测站非线性运动建模方法,包括以下步骤:
获取GNSS测站的原始坐标时间序列;
对原始坐标时间序列进行小波分解和重构,得到低频部分和多层的高频部分;
对所述低频部分和每层高频部分分别进行奇异谱分解和重建;
将奇异谱分解和重建后的低频部分和每层高频部分进行合成得到原始坐标时间序列的拟合序列;
对所述拟合序列进行建模。
该方法的有益效果是:该方法先利用小波分解和重构将原始坐标时间序列分解成频率成分上相对单一、平滑的低频部分和各层高频部分,接着针对低频部分和各层高频部分分别进行奇异谱分解和重建,然后将奇异谱分解和重建后的各部分进行合成得到原始坐标时间序列的拟合序列,最后对合成后的拟合序列进行建模,在一定程度上降低了部分周期项如季节周期项、月周期项被当作噪声剔除的概率,从而提高了建模精度。
为了实现小波分解和重构,进一步地,采用haar小波、dbN小波、symN小波或bior小波模型对原始坐标时间序列进行小波分解和重构。
为了实现奇异谱分解和重建,进一步地,进行奇异谱分解和重建时选择的窗口长度为365。
本发明还提供了一种GNSS测站非线性运动建模装置,包括存储器和处理器,所述处理器用于运行存储在所述存储器中的程序指令,以实现如下方法:
获取GNSS测站的原始坐标时间序列;
对原始坐标时间序列进行小波分解和重构,得到低频部分和多层的高频部分;
对所述低频部分和每层高频部分分别进行奇异谱分解和重建;
将奇异谱分解和重建后的低频部分和每层高频部分进行合成得到原始坐标时间序列的拟合序列;
对所述拟合序列进行建模。
该装置的有益效果是:该装置先利用小波分解和重构将原始坐标时间序列分解成频率成分上相对单一、平滑的低频部分和各层高频部分,接着针对低频部分和各层高频部分分别进行奇异谱分解和重建,然后将奇异谱分解和重建后的各部分进行合成得到原始坐标时间序列的拟合序列,最后对合成后的拟合序列进行建模,在一定程度上降低了部分周期项如季节周期项、月周期项被当作噪声剔除的概率,从而提高了建模精度。
为了实现小波分解和重构,进一步地,采用haar小波、dbN小波、symN小波或bior小波模型对原始坐标时间序列进行小波分解和重构。
为了实现奇异谱分解和重建,进一步地,进行奇异谱分解和重建时选择的窗口长度为365。
附图说明
图1是本发明的GNSS测站非线性运动建模方法原理图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及具体实施例对本发明进行进一步详细说明。
方法实施例:
如图1所示,本实施例的GNSS测站非线性运动建模方法,包括以下步骤:
步骤1、获取GNSS测站的原始坐标时间序列S。具体过程如下:
选取IGS网站(该网站的网址为ftp://igs-rf.ensg.eu/pub/res/)上全球范围内某个纬度的GPS站坐标时间序列数据,作为GNSS测站的原始坐标时间序列S,设原始坐标时间序列S={s1,s2,s3…}。
步骤2、对原始坐标时间序列S进行小波分解和重构,得到低频部分和多层的高频部分。具体过程如下:
本实施例中,采用haar小波模型对原始坐标时间序列S={s1,s2,s3…}进行小波多尺度分解和重构,得到下式:
S=AB+D1+D2+D3+…+DN (1)
式中,AB={aB1,aB2,aB3…}为原始坐标时间序列S中低频信号的重构结果(即低频部分),D1={d11,d12,d13…},…,Di={di1,di2,di3…},…,DN={dN1,dN2,dN3…}分别为原始坐标时间序列S中第1层到第N层的高频信号的重构结果(即多层的高频部分)。作为其他实施方式,还可以采用dbN小波、symN小波或bior小波模型对原始坐标时间序列S进行小波多尺度分解和重构。
那么,ti时刻的原始坐标时间序列S可以用下式表示:
si=aBi+d1i+d2i+…+dNi (2)
式中,aBi为低频信号在ti时刻的值,d1i,d2i,…,dNi分别为各层高频信号在ti时刻的值,i表示变量。
步骤3、对低频部分和每层高频部分分别进行奇异谱分解和重建(即SSA建模)。
首先,利用SSA方法(即奇异谱分析方法)对低频部分AB={aB1,aB2,aB3…}进行奇异谱分解和重建,得到低频信号在ti时刻的拟合值为然后,利用SSA方法对每层高频部分Di={di1,di2,di3…},i=1,2,3,...,N分别进行奇异谱分解和重建,得到各层高频信号在ti时刻的拟合值分别为/>本实施例中,进行奇异谱分解和重建时,低频部分和每层高频部分对应的窗口长度、重构阶次是相同的。
下面以一维时间序列x=(x1,x2,x3,...,xN)为例,详细介绍利用SSA方法对其进行奇异谱分解和重建的过程,对低频部分和每层高频部分分别进行奇异谱分解和重建的过程与此类似,主要由以下4个步骤构成:
第一步,构建时滞矩阵。
首先,选择适当的窗口长度L,满足一般情况下,L的选取不宜超过整个时间序列数据长度N的1/3。其中,窗口长度L并没有统一的选择标准,而是需要根据具体应用而定,一般说来,L过大,会导致奇异值分解得到不同的部分产生混叠,而L太小,则无法使用信号从弱到强渐进划分,从而使得部分信号无法获取。如果能根据经验大致确定时间序列中数据的周期特征,则L一般取周期的公倍数。本实施例中,由于选取的原始坐标时间序列S中存在的已知周期有周年周期和半年周期,且本实施例所用数据的时间分辨率为一天,因此,可以选择已知周期的最小公倍数365作为窗口长度L。
其次,根据一维时间序列x=(x1,x2,x3,...,xN)构建时滞矩阵X,得到:
式(3)中,K=N-L+1,X为L×K阶时滞矩阵,且矩阵X的每一条副对角线上的元素相等,即若将矩阵X中第i行、第j列的元素用xi,j表示,则有xi,j=xi-1,j+1,因此矩阵X是一个Hankel矩阵(即汉克尔矩阵)。
第二步,SVD分解(即奇异值分解)。
对时滞矩阵X进行奇异值分解:
式(4)中,由矩阵XXT和矩阵XTX的前d(d=min{L,K})个最大非零特征值构成对角矩阵为矩阵X的奇异值,等价于矩阵XXT特征值的平方根;U为矩阵X的左奇异向量,等价于矩阵XXT的特征向量;V为矩阵X的右奇异向量,等价于XTX的特征向量。
由此推知:
式中,为矩阵X的奇异值,/>为奇异谱;Ui通常由经验正交函数(Empirical Orthogonal Functions,EOF)表示;Vi为主成分(Principle Components,PC);/>为矩阵X的第i个三重特征向量。最大特征值对应的特征向量,代表了信号的最大变化趋势,而较小的特征值对应的特征向量一般被当作噪声。
因此,时滞矩阵X的奇异值分解可以表示为:
X=X1+X2+…+Xd (7)
第三步,分组。
将矩阵Xi的下标{1,2,…,d}分成M个互不相交的集合I1,I2,…,IM,设I=[i1,i2,…,ip,那么与集合I相关的矩阵XI就可以表示为XI=Xi1+Xi2+…+Xip,时滞矩阵X就可以表示成:
X=XI1+XI2+…+XIM (8)
一般说来,在SVD分解中,前r(r<d)个特征向量对时滞矩阵X的贡献率ηr可以表示为:
第四步,对角平均化。
对角平均化的目的就是把第三步分解得到的矩阵XIM重新转换成为长度为N的新时间序列,称之为重建成分(Reconstruction component,RC),所有RC之和等于原始序列x。
假设z=z1,z2,…,zN为z经过对角平均化得到的时间序列,则对角平均化的公式可表示为:
这样,所有重建成分RC叠加后的和与原始序列x相同,即:
截取前P个贡献较大的成分近似表示原始序列x,则有:
其中,重构阶次P的选择需要根据具体应用而定,一般来说,重构阶次P主要根据奇异值的贡献率来确定,若P太小,则表示后面部分信号将被当作噪声剔除掉,若P太大,则会使得部分噪声被当作信号而提取出来。实际应用时,可以将贡献率较大的前几项之和与设定的某一阈值进行比较,来确定重构阶次P。
本实施例中,由于选取的原始坐标时间序列S中的周期特征主要集中在前几个特征值重构的序列中,所以重构阶次P一般取6~10即可。
步骤4、将奇异谱分解和重建后的低频部分和每层高频部分进行合成得到原始坐标时间序列S的拟合序列
根据步骤3中得到的原始坐标时间序列S中低频信号和各层高频信号在ti时刻的拟合值,得到原始坐标时间序列S在ti时刻的拟合序列为:
式(13)中,为低频信号在ti时刻的拟合值,/>分别为各层高频信号在ti时刻的拟合值。
步骤5、对拟合序列进行建模。
利用步骤4得到原始坐标时间序列的拟合序列后,可以根据实际需要建立多种模型,例如:(1)利用拟合序列建立时间序列趋势模型,用于分离和提取拟合序列中的趋势成分,进而得到GNSS测站三维方向存在的上升或移动趋势的量级;(2)利用拟合序列建立多路径效应改正模型,用于削弱原始坐标时间序列中多路径效应的影响;(3)利用拟合序列建立最优噪声模型,用于估计GNSS测站速度场的运动趋势;(4)利用拟合序列建立共模误差估计模型,用于提取原始坐标时间序列中的共模误差并剔除,进而提高原始坐标时间序列的精度。
基于奇异谱分析的非线性运动建模方法能有效地将半年及半年以上的周期项提取出来,而短周期项如季节周期项、月周期项则容易被忽略。本发明有效融合了小波多尺度分解和奇异谱分析两种方法的优点,先利用小波分解和重构将原始坐标时间序列分解成频率成分上相对单一、平滑的低频部分和各层高频部分,接着针对低频部分和各层高频部分分别进行奇异谱分解和重建,然后将奇异谱分解和重建后的各部分进行合成得到原始坐标时间序列的拟合序列,最后对合成后的拟合序列进行建模,在一定程度上降低了部分周期项如季节周期项、月周期项被当作噪声剔除的概率,从而提高了建模精度。
装置实施例:
本实施例的GNSS测站非线性运动建模装置,包括存储器和处理器,处理器用于运行存储在存储器中的程序指令,以实现如下方法:
获取GNSS测站的原始坐标时间序列;
对原始坐标时间序列进行小波分解和重构,得到低频部分和多层的高频部分;
对低频部分和每层高频部分分别进行奇异谱分解和重建;
将奇异谱分解和重建后的低频部分和每层高频部分进行合成得到原始坐标时间序列的拟合序列;
对拟合序列进行建模。
其中,该方法的具体实施过程与方法实施例中GNSS测站非线性运动建模方法的实施过程相同,此处不再赘述。
Claims (6)
1.一种GNSS测站非线性运动建模方法,其特征在于,包括以下步骤:
获取GNSS测站的原始坐标时间序列;
对原始坐标时间序列进行小波分解和重构,得到低频部分和多层的高频部分;
对所述低频部分和每层高频部分分别进行奇异谱分解和重建:首先,利用奇异谱分析方法对低频部分进行奇异谱分解和重建,得到低频信号在ti时刻的拟合值,然后,利用奇异谱分析方法对每层高频部分分别进行奇异谱分解和重建,得到各层高频信号在ti时刻的拟合值;进行奇异谱分解和重建时,低频部分和每层高频部分对应的窗口长度、重构阶次是相同的;
将奇异谱分解和重建后的低频部分和每层高频部分进行合成得到原始坐标时间序列的拟合序列;
对所述拟合序列进行建模。
2.根据权利要求1所述的GNSS测站非线性运动建模方法,其特征在于,采用haar小波、dbN小波、symN小波或bior小波模型对原始坐标时间序列进行小波分解和重构。
3.根据权利要求1所述的GNSS测站非线性运动建模方法,其特征在于,进行奇异谱分解和重建时选择的窗口长度为365。
4.一种GNSS测站非线性运动建模装置,其特征在于,包括存储器和处理器,所述处理器用于运行存储在所述存储器中的程序指令,以实现如下方法:
获取GNSS测站的原始坐标时间序列;
对原始坐标时间序列进行小波分解和重构,得到低频部分和多层的高频部分;
对所述低频部分和每层高频部分分别进行奇异谱分解和重建:首先,利用奇异谱分析方法对低频部分进行奇异谱分解和重建,得到低频信号在ti时刻的拟合值,然后,利用奇异谱分析方法对每层高频部分分别进行奇异谱分解和重建,得到各层高频信号在ti时刻的拟合值;进行奇异谱分解和重建时,低频部分和每层高频部分对应的窗口长度、重构阶次是相同的;
将奇异谱分解和重建后的低频部分和每层高频部分进行合成得到原始坐标时间序列的拟合序列;
对所述拟合序列进行建模。
5.根据权利要求4所述的GNSS测站非线性运动建模装置,其特征在于,采用haar小波、dbN小波、symN小波或bior小波模型对原始坐标时间序列进行小波分解和重构。
6.根据权利要求4所述的GNSS测站非线性运动建模装置,其特征在于,进行奇异谱分解和重建时选择的窗口长度为365。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910345177.8A CN110069868B (zh) | 2019-04-26 | 2019-04-26 | Gnss测站非线性运动建模方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910345177.8A CN110069868B (zh) | 2019-04-26 | 2019-04-26 | Gnss测站非线性运动建模方法与装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110069868A CN110069868A (zh) | 2019-07-30 |
CN110069868B true CN110069868B (zh) | 2024-04-02 |
Family
ID=67369141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910345177.8A Active CN110069868B (zh) | 2019-04-26 | 2019-04-26 | Gnss测站非线性运动建模方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110069868B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111160317B (zh) * | 2020-01-06 | 2023-03-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 微弱信号盲提取方法 |
CN112699608B (zh) * | 2020-12-31 | 2022-05-06 | 哈尔滨工业大学 | 适用于传感器掉电导致数据缺失的时间序列修复方法 |
CN113052053B (zh) * | 2021-03-18 | 2023-05-05 | 中国人民解放军战略支援部队信息工程大学 | 一种坐标时间序列周期项的提取方法 |
CN113642185A (zh) * | 2021-08-19 | 2021-11-12 | 东华大学 | 一种卫星高频颤振传递路径辨识方法 |
CN114332389B (zh) * | 2021-12-24 | 2022-11-08 | 中国测绘科学研究院 | 一种三维地壳形变模型的构建方法和装置 |
CN114844796B (zh) * | 2022-04-29 | 2024-07-09 | 济南浪潮数据技术有限公司 | 一种对时序kpi的异常检测的方法、装置及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106814378A (zh) * | 2017-01-17 | 2017-06-09 | 华东交通大学 | 一种gnss位置时间序列周期特性挖掘方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018090467A1 (zh) * | 2016-11-17 | 2018-05-24 | 丹阳华神电器有限公司 | 基于模糊熵的含噪信号处理方法及迭代奇异谱软阈值去噪方法 |
-
2019
- 2019-04-26 CN CN201910345177.8A patent/CN110069868B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106814378A (zh) * | 2017-01-17 | 2017-06-09 | 华东交通大学 | 一种gnss位置时间序列周期特性挖掘方法 |
Non-Patent Citations (2)
Title |
---|
基于奇异谱分析法的GPS时间序列周期项探测;汤文娟;;城市勘测(04);全文 * |
我国毫米级框架实现与维持发展现状和趋势;程鹏飞;成英燕;;测绘学报(10);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110069868A (zh) | 2019-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110069868B (zh) | Gnss测站非线性运动建模方法与装置 | |
CN110361778B (zh) | 一种基于生成对抗网络的地震数据重建方法 | |
Sun et al. | Automatic detection of volcanic surface deformation using deep learning | |
CN105572703B (zh) | 一种gps时间序列广义共模误差提取方法 | |
CN112348812B (zh) | 林分年龄信息测量方法及装置 | |
CN106202756B (zh) | 基于单层感知机的欠定盲源分离源信号恢复方法 | |
CN1656389A (zh) | 处理地震数据的高分辨率拉冬变换 | |
CN104599292A (zh) | 一种基于低秩矩阵分解的抗噪运动目标检测算法 | |
CN108809463B (zh) | 一种多基准时钟加权合成方法 | |
CN108830167A (zh) | 一种基于k-svd字典学习的海水温度场数据重构方法 | |
CN109725290B (zh) | 一种误差提取方法、装置、电子设备及可读存储介质 | |
CN108230375A (zh) | 基于结构相似性快速鲁棒的可见光图像与sar图像配准方法 | |
Chen et al. | ARU-net: Reduction of atmospheric phase screen in SAR interferometry using attention-based deep residual U-net | |
CN103020922A (zh) | 基于pca变换的sar图像相干斑抑制方法 | |
CN103870836A (zh) | 基于凸集投影超分辨率重构的合成孔径雷达图像目标识别方法 | |
CN116152206A (zh) | 一种光伏输出功率预测方法、终端设备及存储介质 | |
CN116992217A (zh) | 基于多尺度动态加权多维残差卷积的机械信号降噪方法 | |
CN106485074B (zh) | 一种基于自适应采样率的海洋温盐场采样方法 | |
CN104463245B (zh) | 一种目标识别方法 | |
CN111142134B (zh) | 一种坐标时间序列处理方法及装置 | |
CN103530857B (zh) | 基于多尺度的卡尔曼滤波图像去噪方法 | |
CN112946647A (zh) | 大气误差改正InSAR干涉图堆叠地质灾害普查方法和装置 | |
CN106971392A (zh) | 一种结合dt‑cwt和mrf的遥感图像变化检测方法与装置 | |
CN113487491B (zh) | 一种基于稀疏性与非局部均值自相似性的图像复原方法 | |
CN106934783B (zh) | 一种高频次遥感时间序列数据的时域重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |