CN103927725A - 基于分数阶微分的电影核磁共振图像序列运动场估计方法 - Google Patents
基于分数阶微分的电影核磁共振图像序列运动场估计方法 Download PDFInfo
- Publication number
- CN103927725A CN103927725A CN201410190540.0A CN201410190540A CN103927725A CN 103927725 A CN103927725 A CN 103927725A CN 201410190540 A CN201410190540 A CN 201410190540A CN 103927725 A CN103927725 A CN 103927725A
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- gamma
- centerdot
- resonance image
- 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.)
- Granted
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
基于分数阶微分的电影核磁共振图像序列运动场估计方法,属于核磁共振成像数据处理领域,本发明的目的是为了解决应用现有技术利用整数阶微分图像增强和直接建立光流方程的方法估计电影核磁共振图像序列运动场,存在丢失图像的纹理细节,估计结果受光照变化影响,不具有旋转不变性,抗噪声性能差,导致基于电影核磁共振图像运动估计的精确度低的问题。实现该方法的主要步骤为:一、对电影核磁共振图像利用分数阶微分进行纹理增强;二、通过Riesz变换提取图像的单演信号,即单演相位、单演方位、单演振幅;三、利用单演信号的相位向量建立光流方程;四、通过光流方程估计电影核磁共振图像序列的运动场。本发明用于电影核磁共振图像估计成像对象的运动。
Description
技术领域
本发明属于核磁共振成像数据处理领域。
背景技术
核磁共振成像(Magnetic Resonance Imaging,MRI)技术已经成为临床诊断的重要辅助手段。基于电影核磁共振图像(Cine-MRI)运动场估计是一个不可忽视的研究方向。由于Cine-MRI的灰度非常相近,难以找到密集的对应特征点,因而对于Cine-MRI的运动估计比对于加标记核磁共振图像(Tagged-MRI)的运动估计更具有难度,目前国内外研究较少。针对Cine-MRI的灰度相近,特征稀少的挑战,微分运算增强图像已得到广泛应用,应用整数阶微分图像增强方法(如Sobel、Prewitt、Laplacian算子)在增强图像的同时丢失了图像的纹理细节,尤其是平滑区域的纹理在整数阶微分后损失严重;而分数阶微分可以在提高图像信号高频信息的同时非线性地保留低频信息(纹理细节)。近年来,分数阶微积分逐渐应用于图像处理的各个方面。
单演信号是一种二维解析信号,是由Felsberg提出的,它是一维解析信号的一个扩展,原始的二维图像信息通过具有旋转不变性的Riesz变换,以一种非线性的方式被映射到了虚平面。与一维解析信号类似,该模型涵盖了原始图像在实平面和虚平面上的不同信息,可以据此提取出相应的单演幅度,单演相位和单演方向特征信息。与其他提取图像相位信息的方法如基于Gabor滤波器的相位计算方法相比,该模型的相位计算不需要对方向进行采样。相对于边缘轮廓信息,其中的单演相位特征包含了该图像的绝大部分重要特征信息,只需要这个单演相位向量信息就能建立光流方程并估计运动场。单演相位特征最大的优点是独立于亮度信息,不受光照变化的影响,而且该特征量具有旋转不变性同时也具有很好的抗噪声性能。
应用现有技术整数阶微分图像增强方法和直接建立光流方程的方法,估计电影核磁共振图像序列运动场,存在丢失图像的纹理细节,结果受光照变化影响,不具有旋转不变性,抗噪声性能差,基于电影核磁共振图像运动估计的精确度低。
发明内容
发明目的:为了解决应用现有技术利用整数阶微分图像增强和直接建立光流方程的方法估计电影核磁共振图像序列运动场,存在丢失图像的纹理细节,估计结果受光照变化影响,不具有旋转不变性,抗噪声性能差,导致基于电影核磁共振图像运动估计的精确度低的问题,本发明提出了一种基于分数阶微分的电影核磁共振图像序列运动场估计方法。
技术方案:本发明基于分数阶微分的电影核磁共振图像序列运动场估计方法,实现该方法的步骤如下:
一、利用分数阶微分对电影核磁共振图像序列进行纹理增强;
二、通过Riesz变换提取图像的单演信号,即单演相位、单演方位、单演振幅;
三、利用单演信号的相位向量建立光流方程;
四、通过光流方程估计电影核磁共振图像序列的运动场。
有益效果:
(1)利用分数阶微分对电影核磁共振图像进行纹理增强,克服了现有技术应用整数阶微分图像增强方法在增强图像的同时产生的丢失图像纹理细节的缺陷,可以在提高图像信号高频信息的同时非线性地保留低频信息,低频信息即图像纹理细节。
(2)利用单演信息构造光流方程,单演相位信息独立于亮度信息,不受光照变化的影响。
(3)本方法利用的单演信息特征量具有旋转不变性,具有很好的抗噪声性能,有效提高电影核磁共振图像序列的运动场的估计精度,比不利用单演信息构造光流方程估计精度高10%以上。
本发明用于电影核磁共振图像估计成像对象的运动。
附图说明
图1基于分数阶微分的电影核磁共振图像序列运动场估计方法的流程图;
图2任意一点(x,y)的8个对称方向,分别为x轴正方向,x轴负方向,y轴正方向,y轴负方向,左下方向,左上方向,右上方向,右下方向;其中水平向右为x轴正方向,水平向左为为x轴负方向,竖直向上为y轴正方向,竖直向下为y轴负方向;
图3分数阶微分增强掩模构造示意图,以兴趣点为中心,分别向8个对称方向构造掩模
图4单演信息提取示意图,p,q1,q2构成正交三维向量空间,Α=(p,q1,q2)T,q=(q1,q2)T;A与p的夹角为q与q1的夹角为θ。
具体实施方式
具体实施方式一:结合附图1,本实施方式基于分数阶微分的电影核磁共振图像序列运动场估计方法的具体实施步骤如下:
一、利用分数阶微分对电影核磁共振图像序列进行纹理增强;
二、通过Riesz变换提取图像的单演相位,单演方位,单演振幅,构造单演信号;
三、利用单演信号的相位向量建立光流方程;
四、通过光流方程估计电影核磁共振图像序列的运动场。
具体实施方式二:本实施方式是对具体实施方式一中的基于分数阶微分的电影核磁共振图像序列运动场估计方法的具体实施步骤一的进一步限定,结合附图2、附图3,步骤一所述利用分数阶微分对电影核磁共振图像序列进行纹理增强的过程为:
从连续函数整数阶导数的出发,将微分的阶数由整数扩展到分数,构建基本的v阶Grümwald–Letnikov分数阶微分方程,
这里连续函数s(u)表示一维图像信号,s(u)∈[a,u],a<u,a∈R,u∈R,u表示一维信号变化量,R表示实数集,[v]表示v的整数部分,表示整数集合;当v>0时,k不小于[v],表示信号的长度;是Gamma函数,表示Grümwald–Letnikov-分数阶微分操作符;当N值足够大时,推导出一维图像信号分数阶微分表达式,
针对二维图像s(x,y),有以下两个表达式,(x,y)表示图像中一点的坐标,
为了使分数阶微分掩模具有45°旋转不变性,针对二维图像中的任意一点(x,y)设计在8个对称方向下的分数阶微分,分别是X轴正方向,X轴负方向,Y轴正方向,Y轴负方向,左下方向,左上方向,右上方向,右下方向;
根据公式(2)和公式(3)构造分数阶微分增强掩模公式如下:
这里公式(5)中各式是兴趣点掩模公式,s0=s(x0,y0),(x0,y0)表示兴趣点坐标,用掩模依次卷积需增强图像的每个像素点。
具体实施方式三:本实施方式是对具体实施方式一中的基于分数阶微分的电影核磁共振图像序列运动场估计方法的具体实施步骤二的进一步限定,结合附图4,步骤二所述通过Riesz变换提取图像的单演相位,单演方位,单演振幅,构造单演信号的过程为:
a、通过3个二维空间正交滤波器将图像对应点的局部振幅、局部相位和局部方向分离出来;所采用的二维空间正交滤波器为差分泊松滤波器;空间正交滤波器由1个旋转不变的偶数带通滤波器be(Z)和2个奇数带通滤波器bo1(Z)、bo2(Z)组成;
b、求单演相位单演方向θ(Z)和单演振幅A(Z):
这里Z=(x,y),p(Z)=(I*be)(Z),q1(Z)=(I*bo1)(Z),q2(Z)=(I*bo2)(Z),q(Z)=[q1(Z),q2(Z)]T,“*”,符号表示二维卷积,I表示对应Z点灰度。
具体实施方式四:本实施方式是对具体实施方式一中的基于分数阶微分的电影核磁共振图像序列运动场估计方法的具体实施步骤三的进一步限定,步骤三所述利用单演信号的相位向量建立光流方程的方法为:
基于局部约束,假定以p点为中心的一个小区域Ω内各点的光流相同,对区域内不同的点给予不同的权重,离p点越近,权重越高;
列出光流的计算方程:
上式中,Ω代表2维图像中以p点为中心的一个小的区域,W(p)为窗函数,代表区域中各点的权重,离p点越近,权重越高,表示t时刻p点灰度的梯度,v[U V]T表示运动向量,U表示在X轴方向位移,V表示在Y轴方向位移,It(p,t)表示图像中p点在t时刻的灰度;
对于邻域Ω内的n点,其中,
W=diag(W(p1),W(p2),...,W(pn)),diag表示主对角矩阵是一个主对角线之外的元素皆为0的矩阵,
b=-(It(p1),It(p2),...,It(pn))T。
具体实施方式五:本实施方式是对具体实施方式四中的基于分数阶微分的电影核磁共振图像序列运动场估计方法的具体实施步骤四的进一步限定,步骤四通过光流方程估计电影核磁共振图像序列的运动场的方法为:
由光流计算方程:
解得运动向量v=(ATW2A)-1ATW2b;
其中ATW2A为2×2矩阵:
上式中所有的求和都是在邻域Ω内所有点上进行的;
Ix和Iy表示灰度在X轴和Y轴方向的变化率,本方法用单演信息相位沿方向n的变化率代替Ix和Iy;
代入单演信息求出单演相位信息下的运动向量v,对电影核磁共振图像进行区域划分,依次求出各个区域内单演相位信息下的运动向量v,根据各个运动向量构造出电影核磁共振图像的运动场。
Claims (5)
1.基于分数阶微分的电影核磁共振图像序列运动场估计方法,其特征在于,实现该方法的步骤如下:
一、利用分数阶微分对电影核磁共振图像序列进行纹理增强;
二、通过Riesz变换提取图像的单演信号,即单演相位、单演方位、单演振幅;
三、利用单演信号的相位向量建立光流方程;
四、通过光流方程估计电影核磁共振图像序列的运动场。
2.根据权利要求1所述的基于分数阶微分的电影核磁共振图像序列运动场估计方法,其特征在于,
步骤一所述利用分数阶微分对电影核磁共振图像进行纹理增强的过程为:
从连续函数整数阶导数的出发,将微积分的阶数由整数扩展到分数,构建基本的v阶Grümwald–Letnikov分数阶微分方程,
这里连续函数s(u)表示一维图像信号,s(u)∈[a,u],a<u,a∈R,u∈R,u表示一维信号变化量,R表示实数集,[v]表示v的整数部分,表示整数集合;当v>0时,k不小于[v],表示信号的长度;是Gamma函数,表示Grümwald–Letnikov-分数阶微分操作符;当N值足够大时,推导出一维图像信号分数阶微分表达式,
针对二维图像s(x,y),有以下两个表达式,(x,y)表示图像中一点的坐标,
针对二维图像中的任意一点(x,y)设计在8个对称方向下的分数阶微分,分别是X轴正方向,X轴负方向,Y轴正方向,Y轴负方向,左下方向,左上方向,右上方向,右下方向;
根据公式(2)和公式(3)构造分数阶微分增强掩模公式如下:
这里公式(5)中各式是兴趣点掩模公式,s0=s(x0,y0),(x0,y0)表示兴趣点坐标,用掩模依次卷积需增强图像的每个像素点。
3.根据权利要求1或2所述的基于分数阶微分的电影核磁共振图像序列运动场估计方法,其特征在于,步骤二通过Riesz变换提取图像的单演信号,即单演相位、单演方位、单演振幅的过程为:
a、通过3个二维空间正交滤波器将图像对应点的局部振幅、局部相位和局部方向分离出来;所采用的二维空间正交滤波器为差分泊松滤波器;空间正交滤波器由1个旋转不变的偶数带通滤波器be(Z)和2个奇数带通滤波器bo1(Z)、bo2(Z)组成;
b、求单演相位单演方向θ(Z)和单演振幅A(Z):
这里Z=(x,y),p(Z)=(I*be)(Z),q1(Z)=(I*bo1)(Z),q2(Z)=(I*bo2)(Z),q(Z)=[q1(Z),q2(Z)]T,“*”,符号表示二维卷积,I表示对应Z点灰度。
4.根据权利要求3所述的基于分数阶微分的电影核磁共振图像序列运动场估计方法,其特征在于,步骤三所述利用单演信号的相位向量建立光流方程的过程为:
基于局部约束,假定以p点为中心的一个小区域Ω内各点的光流相同,对区域内不同的点给予不同的权重,离p点越近,权重越高;
列出光流的计算方程:
上式中,Ω代表2维图像中以p点为中心的一个小的区域,W(p)为窗函数,代表区域中各点的权重,离p点越近,权重越高,表示t时刻p点灰度的梯度,v=[U V]T表示运动向量,U表示在X轴方向位移,V表示在Y轴方向位移,It(p,t)表示图像中p点在t时刻的灰度;
对于邻域Ω内的n点,其中,
W=diag(W(p1),W(p2),...,W(pn)),diag表示主对角矩阵是一个主对角线之外的元素皆为0的矩阵,
b=-(It(p1),It(p2),...,It(pn))T。
5.根据权利要求4所述的基于分数阶微分的电影核磁共振图像序列运动场估计方法,其特征在于,所述步骤四通过光流方程估计电影核磁共振图像序列的运动场的过程为:
由光流计算方程:
解得运动向量v=(ATW2A)-1ATW2b;
其中ATW2A为2×2矩阵:
上式中所有的求和都是在邻域Ω内所有点上进行的;
Ix和Iy表示灰度在X轴和Y轴方向的变化率,本方法用单演信息相位沿方向n的变化率代替Ix和Iy;
代入单演信息和求出单演相位信息下的运动向量v;
对电影核磁共振图像进行区域划分,依次求出各个区域内单演相位信息下的运动向量v,根据各个运动向量构造出电影核磁共振图像的运动场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410190540.0A CN103927725B (zh) | 2014-05-07 | 2014-05-07 | 基于分数阶微分的电影核磁共振图像序列运动场估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410190540.0A CN103927725B (zh) | 2014-05-07 | 2014-05-07 | 基于分数阶微分的电影核磁共振图像序列运动场估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103927725A true CN103927725A (zh) | 2014-07-16 |
CN103927725B CN103927725B (zh) | 2017-04-26 |
Family
ID=51145937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410190540.0A Expired - Fee Related CN103927725B (zh) | 2014-05-07 | 2014-05-07 | 基于分数阶微分的电影核磁共振图像序列运动场估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103927725B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104243976A (zh) * | 2014-09-23 | 2014-12-24 | 浙江科技学院 | 一种立体图像客观质量评价方法 |
CN105631897A (zh) * | 2015-12-22 | 2016-06-01 | 哈尔滨工业大学 | 基于单演信号特征距离和互相关变换光流算法的电影核磁共振图像序列运动估计方法 |
CN108805916A (zh) * | 2018-04-27 | 2018-11-13 | 沈阳理工大学 | 一种基于分数阶变分光流模型及对偶优化的图像配准方法 |
CN110133554A (zh) * | 2018-02-08 | 2019-08-16 | 深圳先进技术研究院 | 一种基于分数阶模型的磁共振指纹成像方法、装置及介质 |
CN114947951A (zh) * | 2022-08-01 | 2022-08-30 | 深圳华声医疗技术股份有限公司 | 超声成像控制方法、装置、超声设备及存储介质 |
CN117576334A (zh) * | 2024-01-16 | 2024-02-20 | 青岛市勘察测绘研究院 | 基于数字孪生技术的水下实景三维数据底板获取方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262778A (zh) * | 2011-08-24 | 2011-11-30 | 重庆大学 | 基于改进的分数阶微分掩模的图像增强方法 |
CN102693216A (zh) * | 2012-05-08 | 2012-09-26 | 苏州盛景空间信息技术有限公司 | 基于分数阶微分的点特征跟踪方法 |
EP2515133A1 (en) * | 2011-04-22 | 2012-10-24 | Eidgenössische Technische Hochschule (ETH) | Determining positions of a magnetic field probe in a magnetic resonance measurement |
CN102846326A (zh) * | 2011-06-03 | 2013-01-02 | 通用电气公司 | 用于处理选通图像数据的方法和系统 |
CN103247047A (zh) * | 2013-04-23 | 2013-08-14 | 重庆交通大学 | 一种基于分数阶偏微分的图像边缘检测方法 |
CN103390265A (zh) * | 2013-07-17 | 2013-11-13 | 蒲亦非 | 一种基于分数阶发展方程的纹理图像去噪滤波器 |
CN103400376A (zh) * | 2013-07-19 | 2013-11-20 | 南方医科大学 | 一种乳腺动态增强磁共振图像序列的配准方法 |
US20140002082A1 (en) * | 2007-07-11 | 2014-01-02 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
-
2014
- 2014-05-07 CN CN201410190540.0A patent/CN103927725B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140002082A1 (en) * | 2007-07-11 | 2014-01-02 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
EP2515133A1 (en) * | 2011-04-22 | 2012-10-24 | Eidgenössische Technische Hochschule (ETH) | Determining positions of a magnetic field probe in a magnetic resonance measurement |
CN102846326A (zh) * | 2011-06-03 | 2013-01-02 | 通用电气公司 | 用于处理选通图像数据的方法和系统 |
CN102262778A (zh) * | 2011-08-24 | 2011-11-30 | 重庆大学 | 基于改进的分数阶微分掩模的图像增强方法 |
CN102693216A (zh) * | 2012-05-08 | 2012-09-26 | 苏州盛景空间信息技术有限公司 | 基于分数阶微分的点特征跟踪方法 |
CN103247047A (zh) * | 2013-04-23 | 2013-08-14 | 重庆交通大学 | 一种基于分数阶偏微分的图像边缘检测方法 |
CN103390265A (zh) * | 2013-07-17 | 2013-11-13 | 蒲亦非 | 一种基于分数阶发展方程的纹理图像去噪滤波器 |
CN103400376A (zh) * | 2013-07-19 | 2013-11-20 | 南方医科大学 | 一种乳腺动态增强磁共振图像序列的配准方法 |
Non-Patent Citations (2)
Title |
---|
汪成亮 等: "自适应分数阶微分在图像纹理增强中的应用", 《重庆大学学报》 * |
黄赞: "基于计算机微视觉的微运动测量关键技术研究", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104243976A (zh) * | 2014-09-23 | 2014-12-24 | 浙江科技学院 | 一种立体图像客观质量评价方法 |
CN105631897A (zh) * | 2015-12-22 | 2016-06-01 | 哈尔滨工业大学 | 基于单演信号特征距离和互相关变换光流算法的电影核磁共振图像序列运动估计方法 |
CN105631897B (zh) * | 2015-12-22 | 2018-07-03 | 哈尔滨工业大学 | 基于单演信号特征距离和互相关变换光流算法的电影核磁共振图像序列运动估计方法 |
CN110133554A (zh) * | 2018-02-08 | 2019-08-16 | 深圳先进技术研究院 | 一种基于分数阶模型的磁共振指纹成像方法、装置及介质 |
CN110133554B (zh) * | 2018-02-08 | 2021-04-30 | 深圳先进技术研究院 | 一种基于分数阶模型的磁共振指纹成像方法、装置及介质 |
CN108805916A (zh) * | 2018-04-27 | 2018-11-13 | 沈阳理工大学 | 一种基于分数阶变分光流模型及对偶优化的图像配准方法 |
CN108805916B (zh) * | 2018-04-27 | 2021-06-08 | 沈阳理工大学 | 一种基于分数阶变分光流模型及对偶优化的图像配准方法 |
CN114947951A (zh) * | 2022-08-01 | 2022-08-30 | 深圳华声医疗技术股份有限公司 | 超声成像控制方法、装置、超声设备及存储介质 |
CN114947951B (zh) * | 2022-08-01 | 2022-10-25 | 深圳华声医疗技术股份有限公司 | 超声成像控制方法、装置、超声设备及存储介质 |
CN117576334A (zh) * | 2024-01-16 | 2024-02-20 | 青岛市勘察测绘研究院 | 基于数字孪生技术的水下实景三维数据底板获取方法 |
CN117576334B (zh) * | 2024-01-16 | 2024-04-16 | 青岛市勘察测绘研究院 | 基于数字孪生技术的水下实景三维数据底板获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103927725B (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103927725A (zh) | 基于分数阶微分的电影核磁共振图像序列运动场估计方法 | |
CN103455797B (zh) | 航拍视频中运动小目标的检测与跟踪方法 | |
CN102959584B (zh) | 功能磁共振图像配准方法 | |
CN103034982B (zh) | 一种基于变焦视频序列的图像超分辨率重建方法 | |
CN103106688B (zh) | 基于双层配准方法的室内三维场景重建方法 | |
WO2018024030A1 (zh) | 一种夜视红外图像中的基于显著性的道路目标提取方法 | |
CN103077505A (zh) | 基于字典学习和结构聚类的图像超分辨率重建方法 | |
CN107481315A (zh) | 一种基于Harris‑SIFT‑BRIEF算法的单目视觉三维环境重建方法 | |
CN103810755B (zh) | 基于结构聚类稀疏表示的压缩感知光谱图像重建方法 | |
CN103778636B (zh) | 一种用于无参考图像质量评价的特征构建方法 | |
CN107301664A (zh) | 基于相似性测度函数的改进局部立体匹配方法 | |
CN103218776A (zh) | 基于最小生成树的非局部的深度图超分辨率重建方法 | |
CN103654789A (zh) | 磁共振快速参数成像方法和系统 | |
CN107767425A (zh) | 一种基于单目vio的移动端AR方法 | |
CN103310453A (zh) | 一种基于子图像角点特征的快速图像配准方法 | |
CN104933678A (zh) | 一种基于像素强度的图像超分辨率重建方法 | |
CN103854262A (zh) | 基于结构聚类与稀疏字典学习的医学图像降噪方法 | |
CN102629374A (zh) | 基于子空间投影和邻域嵌入的图像超分辨率重建方法 | |
CN103700101A (zh) | 一种非刚性脑图像配准方法 | |
CN107330854A (zh) | 一种基于新型模板的图像超分辨率增强方法 | |
CN105427340B (zh) | 基于改进箱粒子滤波的视频目标跟踪方法 | |
CN104091364B (zh) | 单幅图像超分辨重建方法 | |
CN108510531A (zh) | 基于pcncc和邻域信息的sar图像配准方法 | |
CN107845107A (zh) | 一种图像透视变换的优化方法 | |
CN1948995A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170426 Termination date: 20180507 |