CN103871087B - 基于三维全变差稀疏先验的高光谱解混压缩感知方法 - Google Patents
基于三维全变差稀疏先验的高光谱解混压缩感知方法 Download PDFInfo
- Publication number
- CN103871087B CN103871087B CN201410102950.5A CN201410102950A CN103871087B CN 103871087 B CN103871087 B CN 103871087B CN 201410102950 A CN201410102950 A CN 201410102950A CN 103871087 B CN103871087 B CN 103871087B
- Authority
- CN
- China
- Prior art keywords
- matrix
- abundances
- follows
- total variation
- formula
- 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
Abstract
本发明公开了一种基于三维全变差稀疏先验的高光谱解混压缩感知方法,用于解决现有联合光谱解混的高光谱图像压缩感知算法精度低的技术问题。技术方案是采用随机观测矩阵从原始数据中抽取少量的样本作为压缩数据。重建过程,根据解混压缩感知模型,从光谱库中选择适当的光谱作为的模型中的端元矩阵,进而引入丰度值矩阵的三维全变差稀疏先验,通过求解受限的线性优化问题,精确地求解丰度值矩阵。最后使用线性混合模型重建原始数据。在HYDICE卫星拍摄的urban数据上当压缩比为1:20时,归一化的均方误差(normalized mean squared error,NMSE)小于0.09,当压缩比为1:10时,归一化均方误差同样小于0.08,相对于已有的压缩感知类算法精度提升10%以上。
Description
技术领域
本发明涉及一种高光谱解混压缩感知方法,特别涉及一种基于三维全变差稀疏先验的高光谱解混压缩感知方法。
背景技术
在高光谱图像中,丰富的光谱信息中包含着极大的数据冗余,严重地增加了高光谱图像采集、传输以及处理过程中的资源消耗。因此,设计一种高性能的高光谱图像压缩算法是非常必要的。已有的高光谱图像压缩算法主要分为两类,一类是基于信息编码的压缩方法,主要利用普通的图像压缩方法去除高光谱图像各个波段内部以及波段之间的冗余性来实现压缩,其中包括聚类的脉冲差分编码,三维小波变换,三维离散余弦变换等,然而这类压缩方法作用在图像获取之后,仍然需要消耗大量的资源来采集和存储数据,而且压缩率较低;另外一类是基于压缩感知(CompressiveSensing)的压缩方法,该类方法作用于数据采集端,通过采集原始稀疏信号的少量样本点实现数据的大幅度压缩,解压时利用压缩感知理论实现原始稀疏信号的精确重建,大大减少了采集和传输过程中的资源消耗。
文献“Acompressivesensingandunmixingschemeforhyperspectraldataprocessing,IEEETransactionsonImageProcessing,2012,21(3):1200–1210”公开了一种联合光谱解混的高光谱图像压缩感知算法。该方法首先使用随机观测矩阵对原始数据进行随机采样实现数据压缩;接着,在对应的光谱库中选择合适的端元;之后利用压缩感知理论对具有梯度稀疏性的丰度值矩阵进行精确重建;最后,结合选择的端元,利用线性混合模型实现原始高光谱图像的重建。然而,该方法仅仅考虑了丰度值矩阵空间维上的梯度稀疏性,忽略了丰度值矩阵光谱维上的潜在稀疏性,从而影响了该压缩算法的重建精度。
发明内容
为了克服现有联合光谱解混的高光谱图像压缩感知算法精度低的不足,本发明提供一种基于三维全变差稀疏先验的高光谱解混压缩感知方法。该方法采用随机观测矩阵从原始数据中抽取少量的样本作为压缩数据。重建过程,根据解混压缩感知模型,从光谱库中选择适当的光谱作为的模型中的端元矩阵,进而引入丰度值矩阵的三维全变差稀疏先验,通过求解受限的线性优化问题,精确地求解丰度值矩阵。最后使用线性混合模型重建原始数据。在HYDICE卫星拍摄的urban数据上当压缩比为1:20时,归一化的均方误差(normalizedmeansquarederror,NMSE)小于0.09,当压缩比为1:10时,归一化均方误差同样小于0.08,相对于已有的压缩感知类算法精度提升10%以上。
本发明解决其技术问题所采用的技术方案是:一种基于三维全变差稀疏先验的高光谱解混压缩感知方法,其特点是包括以下步骤:
步骤一、对于高光谱图像其中每一个像素的光谱xi表示成所有端元的线性组合如下:
xi=Whi(1)
其中,np表示空间上包含的像素数目,nb表示波段数量,为对应的丰度值向量。
整个数据X表示成丰度值矩阵和端元矩阵W的乘积:
X=WH(2)
在H中,行方向是光谱维,每一行代表不同像素的光谱在同一个端元上的投影;列方向是空间维,每一列代表一个像素的光谱在不同端元上的投影。
步骤二、采用满足高斯随机分布的归一化随机观测矩阵对原始数据进行随机采样,得到压缩数据如下:
F=AX=AWH(3)
其中,m表示对长度为nb的信号压缩后的长度,m<nb。
步骤三、对于有限的成像场景,根据场景信息从光谱库中抽取ne个光谱组成端元矩阵W。
步骤四、(1)在H的光谱维上应用一维的全变差稀疏先验,结合H空间维上的稀疏性,得到H的三维全变差稀疏先验,如下:
其中,ej和εj分别表示和空间中的第j个单位向量。TV(x)描述的是的全变差,Di(x)表示x梯度中的第i个分量。公式(4)中的第一部分表示H空间维上的二维全变差稀疏先验,其中对应的Di(·)为二维梯度;第二部分表示H光谱维上的一维全变差稀疏先验,其中对应的Di(·)为一维梯度。
(2)构建丰度值的其他先验。将线性解混模型中常用的丰度值先验引入,分别是混合光谱在不同端元上丰度值投影非负且全和为1的限制,如下:
其中,和是全部元素为1,长度分别为ne和np的向量。
(3)构建丰度值矩阵H的重建模型。结合公式(3)、(4)和(5)得到如下的重建模型:
为了方便后续求解,向公式(6)中引入分离变量vij=Di(Hej),得到:
(4)求解公式(7)得到丰度值矩阵H的估计具体求解过程如下:
①使用增广拉格朗日方法,根据公式(7)构建H的增广拉格朗日方程
其中,α=25,κ=25,β=213,γ=25为二次项惩罚系数,λij,πij,Π,υ为对应的拉格朗日乘子,初始化每个乘子的所有元素为0,||·||F表示Frobenius范数。
②固定拉格朗日乘子和H,更新分离变量vij,uij。形式如下:
③固定拉格朗日乘子和分离变量vij,uij,采用梯度下降法更新H。假定第k次更新,由Hk得到Hk+1,形式如下:
其中,为关于H一阶导数,形式如下:
式中,τ为梯度下降步长。其计算分为初始化和细化两步。在初始化过程中,当第一次更新H0时,τ采用最速下降法进行初始化,之后更新Hk,k=1,2,...时,对τ采用两点步长梯度法进行初始化。两点步长梯度法具体形式如下:
其中,tr(·)表示矩阵的迹。τ的细化过程具体如下:
(a)代入初始化的τ,根据公式(10)得到Hk+1,设置参数δ=3.2×10-4,η=0.6和计数器c=0;
(b)判断Hk+1是否满足如下的条件:
如果不满足,更新计数器c=c+1;
如果c<5,缩小步长τ=τ·η,继续循环判断是否满足(13);
否则τ由最速下降法确定,然后由公式(13)得到更新的Hk+1;
否则,得到更新的Hk+1。
④固定更新后的vij,uij和H,使用如下公式更新拉格朗日乘子:
⑤循环步骤②、③和④直至收敛,得到的最终估计的丰度值矩阵
步骤五、结合选择的端元矩阵W和线性混合模型公式(2)得到重建的高光谱数据
本发明的有益效果是:该方法采用随机观测矩阵从原始数据中抽取少量的样本作为压缩数据。重建过程,根据解混压缩感知模型,从光谱库中选择适当的光谱作为的模型中的端元矩阵,进而引入丰度值矩阵的三维全变差稀疏先验,通过求解受限的线性优化问题,精确地求解丰度值矩阵。最后使用线性混合模型重建原始数据。在HYDICE卫星拍摄的urban数据上当压缩比为1:20时,归一化的均方误差(normalizedmeansquarederror,NMSE)小于0.09,当压缩比为1:10时,归一化均方误差同样小于0.08,相对于已有的压缩感知类算法精度提升10%以上。
以下结合具体实施方式详细说明本发明。
具体实施方式
本发明基于三维全变差稀疏先验的高光谱解混压缩感知方法具体步骤如下:
在高光谱图像中,同一像素的不同波段的反射值构成了一个离散向量,称为该像素的光谱。通常,纯净物质在高光谱图像中具有独一无二的光谱,称作端元。由于存在地物混合和高光谱图像空间分辨率低等因素,像素的光谱往往是多种不同纯净地物光谱的混合。这种光谱混合现象通常可以使用线性混合模型来描述。该模型认为任何一个混合光谱是成像场景中所有端元的线性组合。端元在混合光谱中所占的比例称作丰度值。因此,对于高光谱图像(np表示空间上包含的像素数目,nb表示波段数量),其中每一个像素的光谱xi可以表示成所有端元的线性组合如下:
xi=Whi(1)
其中,为对应的丰度值向量。因此整个数据X可以表示成丰度值矩阵和端元矩阵W的乘积,如下:
X=WH(2)
在H中,行方向是光谱维,每一行代表不同像素的光谱在同一个端元上的投影;列方向是空间维,每一列代表一个像素的光谱在不同端元上的投影。本发明通过压缩感知技术重建丰度值矩阵,之后结合选择的端元和线性混合模型重建原始高光谱图像。
1、获取压缩数据。
采用满足高斯随机分布的归一化随机观测矩阵对原始数据进行随机采样,得到压缩数据如下:
F=AX=AWH(3)
其中,m表示对长度为nb的信号压缩后的长度,m<nb。
2、选择端元。
对于有限的成像场景,通常仅仅包含有限个端元。因此,可以根据场景信息从光谱库(如ASTER光谱库)中抽取ne个光谱组成端元矩阵W。
3、重建丰度值矩阵。
根据公式(3)求解丰度值矩阵是一个不适定问题,因此需要在公式(3)的基础上引入丰度值矩阵H的先验约束。本发明针对H设计了三维全变差稀疏先验,引入到公式(3)中,最后求解得到重建的具体过程如下:
(1)构建H的三维全变差稀疏先验。H的空间维上集中了不同像素在同一端元上的丰度值分量。因此,图像场景中的空间相似性在H的空间维上得以保持,而这种相似性可以使用空间上的二维梯度稀疏性来描述,即空间上的二维全变差稀疏先验。此外,H的光谱维描述了某个像素的光谱在不同端元上的丰度值投影。通常,在图像匀质区域内部的光谱,在该类物质对应的端元上具有较大的投影,其他端元上的投影近似为0,具有一定的稀疏性;而在匀质区域边界处的光谱,往往在两种或两种以上的端元产生相似的丰度值投影,具有梯度稀疏性。因此,为了满足了这两种情形下的稀疏性要求,本发明在H的光谱维上应用一维的全变差稀疏先验。结合H空间维上的稀疏性,得到H的三维全变差稀疏先验,如下:
其中,ej和εj分别表示和空间中的第j个单位向量。TV(x)描述的是的全变差,Di(x)表示x梯度中的第i个分量。公式(4)中的第一部分表示H空间维上的二维全变差稀疏先验,其中对应的Di(·)为二维梯度;第二部分表示H光谱维上的一维全变差稀疏先验,其中对应的Di(·)为一维梯度。
(2)构建丰度值的其他先验。由于本发明中引入了线性混合模型,因此需要将线性解混模型中常用的丰度值先验引入,分别是混合光谱在不同端元上丰度值投影非负且全和为1的限制,如下:
其中,和是全部元素为1,长度分别为ne和np的向量。
(3)构建丰度值矩阵H的重建模型。结合公式(3)、(4)和(5)得到如下的重建模型:
为了方便后续求解,向(6)中引入分离变量vij=Di(Hej),得到:
(4)求解公式(7)得到丰度值矩阵H的估计具体求解过程如下:
①使用增广拉格朗日方法,根据公式(7)构建H的增广拉格朗日方程
其中,α=25,κ=25,β=213,γ=25为二次项惩罚系数,λij,πij,Π,υ为对应的拉格朗日乘子,初始化每个乘子的所有元素为0,||·||F表示Frobenius范数。
②固定拉格朗日乘子和H,更新分离变量vij,uij。形式如下:
③固定拉格朗日乘子和分离变量vij,uij,采用梯度下降法更新H。假定第k次更新,由Hk得到Hk+1,形式如下:
其中,为关于H一阶导数,形式如下:
式中,τ为梯度下降步长。其计算分为初始化和细化两步。在初始化过程中,当第一次更新H0时,τ采用最速下降法进行初始化,之后更新Hk,k=1,2,...时,对τ采用两点步长梯度法进行初始化。两点步长梯度法具体形式如下:
其中,tr(·)表示矩阵的迹。τ的细化过程具体如下:
(c)代入初始化的τ,根据公式(10)得到Hk+1,设置参数δ=3.2×10-4,η=0.6和计数器c=0;
(d)判断Hk+1是否满足如下的条件:
如果不满足,更新计数器c=c+1;
如果c<5,缩小步长τ=τ·η,继续循环判断是否满足(13);
否则τ由最速下降法确定,然后由公式(13)得到更新的Hk+1;
否则,得到更新的Hk+1。
④固定更新后的vij,uij和H,使用如下公式更新拉格朗日乘子:
⑤循环步骤②、③和④直至收敛,得到的最终估计的丰度值矩阵
4、重建高光谱数据。
结合选择的端元矩阵W和线性混合模型公式(2)得到重建的高光谱数据
本发明采用随机观测矩阵从原始数据中抽取少量的样本作为压缩数据。重建过程,根据解混压缩感知模型,从光谱库中选择适当的光谱作为的模型中的端元矩阵,进而引入丰度值矩阵的三维全变差稀疏先验,通过求解受限的线性优化问题,精确地求解丰度值矩阵。最后使用线性混合模型重建原始数据。在HYDICE卫星拍摄的urban数据上当压缩比为1:20时,归一化的均方误差(normalizedmeansquarederror,NMSE)小于0.09,当压缩比为1:10时,归一化均方误差同样小于0.08,相对于已有的压缩感知类算法精度提升10%以上。
Claims (1)
1.一种基于三维全变差稀疏先验的高光谱解混压缩感知方法,其特征在于包括以下步骤:
步骤一、对于高光谱图像其中每一个像素的光谱xi表示成所有端元的线性组合,如下:
xi=Whi(1)
其中,np表示空间上包含的像素数目,nb表示波段数量,为对应的丰度值向量;整个数据X表示成丰度值矩阵和端元矩阵W的乘积:
X=WH(2)
在H中,行方向是光谱维,每一行代表不同像素的光谱在同一个端元上的投影;列方向是空间维,每一列代表一个像素的光谱在不同端元上的投影;
步骤二、采用满足高斯随机分布的归一化随机观测矩阵对原始数据进行随机采样,得到压缩数据如下:
F=AX=AWH(3)
其中,m表示对长度为nb的信号压缩后的长度,m<nb;
步骤三、对于有限的成像场景,根据场景信息从光谱库中抽取ne个光谱组成端元矩阵W;
步骤四、(1)在H的光谱维上应用一维的全变差稀疏先验,结合H空间维上的稀疏性,得到H的三维全变差稀疏先验,如下:
其中,ej和εj分别表示和空间中的第j个单位向量;TV(x)描述的是的全变差,Di(x)表示x梯度中的第i个分量;公式(4)中的第一部分表示H空间维上的二维全变差稀疏先验,其中对应的Di(·)为二维梯度;第二部分表示H光谱维上的一维全变差稀疏先验,其中对应的Di(·)为一维梯度;
(2)构建丰度值的其他先验;将线性解混模型中常用的丰度值先验引入,分别是混合光谱在不同端元上丰度值投影非负且全和为1的限制,如下:
其中,和是全部元素为1,长度分别为ne和np的向量;
(3)构建丰度值矩阵H的重建模型;结合公式(3)、(4)和(5)得到如下的重建模型:
为了方便后续求解,向公式(6)中引入分离变量vij=Di(Hej),得到:
(4)求解公式(7)得到丰度值矩阵H的估计具体求解过程如下:
①使用增广拉格朗日方法,根据公式(7)构建H的增广拉格朗日方程
其中,α=25,κ=25,β=213,γ=25为二次项惩罚系数,λij,πij,Π,υ为对应的拉格朗日乘子,初始化每个乘子的所有元素为0,||·||F表示Frobenius范数;
②固定拉格朗日乘子和H,更新分离变量vij,uij;形式如下:
③固定拉格朗日乘子和分离变量vij,uij,采用梯度下降法更新H;假定第k次更新,由Hk得到Hk+1,形式如下:
其中,为关于H一阶导数,形式如下:
式中,τ为梯度下降步长;其计算分为初始化和细化两步;在初始化过程中,当第一次更新H0时,τ采用最速下降法进行初始化,之后更新Hk,k=1,2,...时,对τ采用两点步长梯度法进行初始化;两点步长梯度法具体形式如下:
其中,tr(·)表示矩阵的迹;τ的细化过程具体如下:
(a)代入初始化的τ,根据公式(10)得到Hk+1,设置参数δ=3.2×10-4,η=0.6和计数器c=0;
(b)判断Hk+1是否满足如下的条件:
如果不满足,更新计数器c=c+1;
如果c<5,缩小步长τ=τ·η,继续循环判断是否满足(13);
否则τ由最速下降法确定,然后由公式(13)得到更新的Hk+1;
否则,得到更新的Hk+1;
④固定更新后的vij,uij和H,使用如下公式更新拉格朗日乘子:
⑤循环步骤②、③和④直至收敛,得到的最终估计的丰度值矩阵
步骤五、结合选择的端元矩阵W和线性混合模型公式(2)得到重建的高光谱数据
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410102950.5A CN103871087B (zh) | 2014-03-20 | 2014-03-20 | 基于三维全变差稀疏先验的高光谱解混压缩感知方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410102950.5A CN103871087B (zh) | 2014-03-20 | 2014-03-20 | 基于三维全变差稀疏先验的高光谱解混压缩感知方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103871087A CN103871087A (zh) | 2014-06-18 |
CN103871087B true CN103871087B (zh) | 2016-07-13 |
Family
ID=50909592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410102950.5A Active CN103871087B (zh) | 2014-03-20 | 2014-03-20 | 基于三维全变差稀疏先验的高光谱解混压缩感知方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103871087B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331880B (zh) * | 2014-10-20 | 2017-02-15 | 西安电子科技大学 | 基于几何空谱结构信息的高光谱混合像元分解方法 |
CN104734724B (zh) * | 2015-03-16 | 2017-11-24 | 西北工业大学 | 基于重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法 |
CN105513102B (zh) * | 2015-12-15 | 2018-11-16 | 西安电子科技大学 | 基于非局部全变差和低秩稀疏的高光谱压缩感知重建方法 |
CN105701845A (zh) * | 2016-01-14 | 2016-06-22 | 辽宁师范大学 | 协同稀疏测量和3d tv模型的高光谱影像压缩感知重构方法 |
CN111243043B (zh) * | 2019-11-21 | 2023-06-16 | 铜陵学院 | 基于修正线性混合模型的高光谱压缩感知方法、设备及系统 |
CN111260576B (zh) * | 2020-01-14 | 2022-07-05 | 哈尔滨工业大学 | 一种基于去噪三维卷积自编码网络的高光谱解混算法 |
CN111461087A (zh) * | 2020-06-16 | 2020-07-28 | 南京理工大学 | 基于光谱保持稀疏自编码器的高光谱异常检测方法 |
CN112364730B (zh) * | 2020-10-29 | 2023-01-17 | 济南大学 | 基于稀疏子空间聚类的高光谱地物自动分类方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102355268A (zh) * | 2011-08-31 | 2012-02-15 | 清华大学 | 基于稀疏化哈达玛矩阵的压缩感知观测矩阵构造方法 |
CN102661794A (zh) * | 2012-03-20 | 2012-09-12 | 清华大学 | 多光谱计算重构方法及系统 |
CN102722866A (zh) * | 2012-05-22 | 2012-10-10 | 西安电子科技大学 | 基于主成份分析的压缩感知方法 |
-
2014
- 2014-03-20 CN CN201410102950.5A patent/CN103871087B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102355268A (zh) * | 2011-08-31 | 2012-02-15 | 清华大学 | 基于稀疏化哈达玛矩阵的压缩感知观测矩阵构造方法 |
CN102661794A (zh) * | 2012-03-20 | 2012-09-12 | 清华大学 | 多光谱计算重构方法及系统 |
CN102722866A (zh) * | 2012-05-22 | 2012-10-10 | 西安电子科技大学 | 基于主成份分析的压缩感知方法 |
Non-Patent Citations (2)
Title |
---|
基于混合基稀疏图像表示的压缩传感图像重构;练秋生等;《自动化学报》;385-391;20100330;第36卷(第3期);第385-391页 * |
高光谱图像压缩感知投影与复合正则重构;冯燕等;《航空学报》;20120825;第33卷(第8期);第1466-1473页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103871087A (zh) | 2014-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103871087B (zh) | 基于三维全变差稀疏先验的高光谱解混压缩感知方法 | |
CN111260576B (zh) | 一种基于去噪三维卷积自编码网络的高光谱解混算法 | |
CN104978573B (zh) | 一种应用于高光谱图像处理的非负矩阵分解方法 | |
CN107525588B (zh) | 一种基于gpu的双相机光谱成像系统的快速重构方法 | |
CN104867124B (zh) | 基于对偶稀疏非负矩阵分解的多光谱与全色图像融合方法 | |
CN103020939B (zh) | 利用多时相数据去除光学遥感影像大面积厚云的方法 | |
CN103632385B (zh) | 基于空谱联合稀疏先验的卫星高光谱压缩感知重建方法 | |
CN105469360A (zh) | 基于非局部联合稀疏表示的高光谱图像超分辨率重建方法 | |
CN107657217A (zh) | 基于运动目标检测的红外与可见光视频的融合方法 | |
CN102708576A (zh) | 基于结构字典的分块图像压缩感知重建方法 | |
CN104050653A (zh) | 基于非负结构稀疏的高光谱图像超分辨率算法 | |
CN111598786B (zh) | 一种基于深度去噪自编码网络的高光谱图像解混方法 | |
CN103295197A (zh) | 基于字典学习和双边正则的图像超分辨率重建方法 | |
CN109087367B (zh) | 一种基于粒子群优化的高光谱图像快速压缩感知重构方法 | |
CN102682441A (zh) | 基于亚像元映射的高光谱图像超分辨重建方法 | |
CN108520495A (zh) | 基于聚类流形先验的高光谱图像超分辨重建方法 | |
CN103020912A (zh) | 一种结合波段聚类和稀疏表达的遥感影像复原方法 | |
CN112215267B (zh) | 一种面向高光谱图像的深度空谱子空间聚类方法 | |
CN109887050A (zh) | 一种基于自适应字典学习的编码孔径光谱成像方法 | |
CN108280486B (zh) | 一种基于端元簇的高光谱图像解混方法 | |
CN104091364A (zh) | 单幅图像超分辨重建方法 | |
CN104091368B (zh) | 基于空谱三维稀疏先验的高光谱解混压缩感知方法 | |
CN105701845A (zh) | 协同稀疏测量和3d tv模型的高光谱影像压缩感知重构方法 | |
CN102542547A (zh) | 一种基于光谱约束的高光谱图像融合方法 | |
CN102073875A (zh) | 基于稀疏表示的背景杂波量化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |