CN107154020A - 一种基于Curvelet变换的影像融合方法及系统 - Google Patents
一种基于Curvelet变换的影像融合方法及系统 Download PDFInfo
- Publication number
- CN107154020A CN107154020A CN201710301849.6A CN201710301849A CN107154020A CN 107154020 A CN107154020 A CN 107154020A CN 201710301849 A CN201710301849 A CN 201710301849A CN 107154020 A CN107154020 A CN 107154020A
- Authority
- CN
- China
- Prior art keywords
- image
- frequency coefficient
- component
- sar
- curvelet
- 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.)
- Pending
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000004927 fusion Effects 0.000 claims abstract description 140
- 238000011156 evaluation Methods 0.000 claims abstract description 42
- 230000009466 transformation Effects 0.000 claims abstract description 37
- 238000000844 transformation Methods 0.000 claims abstract description 36
- 238000000701 chemical imaging Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 11
- 230000000007 visual effect Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 abstract description 13
- 230000009286 beneficial effect Effects 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 4
- 238000007500 overflow downdraw method Methods 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
- G06T3/4061—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution by injecting details from different spectral ranges
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于Curvelet变换的影像融合方法及系统，该方法包括如下步骤：对原始的高光谱影像进行影像插值，新的高光谱影像；接着进行HIS变换，得到I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像；然后进行Curvelet变换，并根据融合规则进行融合，得到新的Curvelet变换系数；最后进行Curvelet逆变换，得到I*分量，同时将I*分量和S、H分量进行HIS逆变换，得到融合后的影像。还涉及一种系统，该系统包括：影像插值模块、影像处理模块、Curvelet变换模块、融合影像模块。通过本发明有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
Description
技术领域
本发明属于多源高分辨率卫星遥感影像融合领域，尤其涉及一种基于Curvelet变换的影像融合方法及系统。
背景技术
1982年，Haydn首次将HIS(Hue-Intensity-Saturation)方法应用到遥感影像融合中，HIS方法首先在多光谱影像中选择适当波段作为RGB三个通道数据，生成RGB彩色影像，然后将彩色影像变换到HIS变换空间，对待融合的高分辨率影像进行灰度拉伸替换I分量，最后进行HIS逆变换得到融合影像，HIS融合方法是一种实现简单、应用广泛的影像融合方法。
但是HIS融合方法要求对比度拉伸后的高分辨率影像不仅要与亮度分量高度相关，而且要与其光谱响应范围同多光谱影像的响应范围接近一致，在实际情况中，二者相关性很低，因此易发生光谱扭曲现象。对于高分辨率的SAR影像而言，由于SAR卫星成像系统基于相干成像原理，使得SAR影像出现不同程度的斑点噪声，高空间分辨率影像质量受其影响更严重，受探测机理影响，其影像无法探测地物目标的光谱信息，高光谱卫星影像数据含有丰富的光谱信息，但空间分辨率通常不会特别高，尤其是在拉伸后的高分辨率影像直接替换高光谱影像的亮度分量，易使SAR影像的噪声引入到融合影像中。
综上所述，由于SAR卫星成像系统基于相干成像原理，使得SAR影像出现不同程度的斑点噪声，高空间分辨率影像质量受其影响更严重，受探测机理影响，其影像无法探测地物目标的光谱信息。高光谱卫星影像数据含有丰富的光谱信息，但空间分辨率通常不会特别高。
因此，是否能够将HIS融合方法与SAR影像方法相融合，改善HIS方法存在的光谱失真和易引入噪声问题，从而得到高质量的融合影像是亟待需要解决的问题。
发明内容
本发明所要解决的技术问题是：现有技术中SAR影像高空间分辨率影像质量较差，其影像无法探测地物目标的光谱信息，并且其空间分辨率较低。
为解决上面的技术问题，本发明提供了一种基于Curvelet变换的影像融合方法，该方法包括如下步骤：
S1,对原始的高光谱影像进行影像插值，得到与SAR影像相同像元尺寸的新的高光谱影像；
S2，对新的高光谱影像进行HIS变换，得到I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像；
S3，对I分量和SAR*影像分别进行Curvelet变换，得到Curvelet变换系数，并根据融合规则对Curvelet变换系数进行融合，得到新的Curvelet变换系数；
S4，将新的Curvelet变换系数进行Curvelet逆变换，得到I*分量，同时将I*分量和S、H分量进行HIS逆变换，得到融合后的影像。
本发明的有益效果：本发明中对Curvelet变换并采用不同的融合规则进行融合，因此由于融合规则中对噪声的进行了抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
进一步地，所述S2中对SAR影像进行灰度拉伸，得到SAR*影像，其公式为：
其中，μ1，σ1分别为SAR影像的均值和标准差，μ0，σ0分别为I分量的均值和标准差。
进一步地，所述S3中包括：
S31，对I分量和SAR*影像分别进行Curvelet变换，分别得到与I分量对应的I分量低频系数、I分量中频系数、I分量高频系数和与新的SAR*影像对应的SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数；
S32，根据融合规则，将I分量低频系数、I分量中频系数、I分量高频系数分别对应地与SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数进行融合，对应地得到新的Curvelet变换的低频系数、中频系数、高频系数。
上述进一步地的有益效果：通过这些不同的融合规则进行的Curvelet变换，同时融合规则对噪声的进行了抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
进一步地，所述S4中包括：
S41，将新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到I*分量；
S42，将I*分量与S、H分量进行HIS反变换，得到融合后的影像。
进一步地，所述S2中包括：
在新的高光谱影像中选择三个波段作为R、G、B波段，并对该三个波段分别进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像。
上述进一步地有益效果：选择三个波段作为R、G、B波段，并非是所有的波段都进行，这样的目的减少非必要的影响的转换，同时也节省很多的时间，提高了精度。
进一步地，所述S32中的融合规则包括：区域标准差和局部方向信息熵的融合规则、特征量积的融合规则、高频系数取零的融合规则；
所述区域标准差和局部方向信息熵的融合规则，是用于I分量低频系数与SAR*影像低频系数进行融合，得到新的Curvelet变换的低频系数；
所述特征量积的融合规则，是用于I分量中频系数与SAR*影像中频系数进行融合，得到新的Curvelet变换的中频系数；
所述高频系数取零的融合规则，是用于I分量高频系数分与SAR*影像高频系数进行融合，得到新的Curvelet变换的高频系数。
上述进一步地的有益效果：通过这些不同的融合规则进行的Curvelet变换，低频采用区域标准差和局部方向信息熵的融合规则，中频采用特征量积的融合规则，高频采用高频系数取零的融合规则，同时融合规则对噪声的进行了抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
本发明还涉及一种基于Curvelet变换的影像融合系统，该系统包括：影像插值模块、影像处理模块、Curvelet变换模块、融合影像模块；
所述影像插值模块，其用于对原始的高光谱影像进行影像插值，得到与SAR影像相同像元尺寸的新的高光谱影像；
所述影像处理模块，其用于对新的高光谱影像进行HIS变换，得到I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像；
所述Curvelet变换模块，其用于对I分量和SAR*影像分别进行Curvelet变换，得到Curvelet变换系数，并根据融合规则对Curvelet变换系数进行融合，得到新的Curvelet变换系数；
所述融合影像模块，其用于将新的Curvelet变换系数进行Curvelet逆变换，得到I*分量，同时将I*和S、H分量进行HIS反变换，得到融合后的影像。
本发明的有益效果：通过这些不同的融合规则进行的Curvelet变换，低频采用区域标准差和局部方向信息熵的融合规则，中频采用特征量积的融合规则，高频采用高频系数取零的融合规则，同时融合规则对噪声的进行了抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
进一步地，所述Curvelet变换模块包括：第一Curvelet变换单元、第二Curvelet变换单元；
所述第一Curvelet变换单元，其用于对I分量和SAR*影像分别进行Curvelet变换，分别得到与I分量对应的I分量低频系数、I分量中频系数、I分量高频系数和与SAR*影像对应的SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数；
所述第二Curvelet变换单元，其用于根据融合规则，将I分量低频系数、I分量中频系数、I分量高频系数分别对应地与SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数进行融合，对应地得到新的Curvelet变换的低频系数、中频系数、高频系数。
上述进一步地有益效果：通过这些不同的融合规则进行的Curvelet变换，同时在融合规则中对噪声的进行抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
进一步地，所述融合影像模块中包括：第一融合影像单元、第二融合影像单元；
所述第一融合影像单元，其用于将新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到I*分量；
所述第二融合影像单元，其用于将I*分量与S、H分量进行HIS反变换，得到融合后影像。
进一步地，所述影像处理模块，其具体用于在新的高光谱影像中选择三个波段作为R、G、B波段，并对该三个波段分别进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像。
上述进一步地有益效果：选择三个波段作为R、G、B波段，并非是所有的波段都进行，这样的目的减少非必要的影响的转换，同时也节省很多的时间，提高了精度。
附图说明
图1为本发明的一种基于Curvelet变换的影像融合方法的流程图；
图2为本发明的一种基于Curvelet变换的影像融合方法的示意图；
图3为本发明的一种基于Curvelet变换的影像融合系统的示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
实施例1
如图1和图2所示，本实施例1中的一种基于Curvelet变换的影像融合方法，该方法包括如下步骤：
S1,对原始的高光谱影像进行影像插值，得到与SAR影像相同像元尺寸的新的高光谱影像；
在本实施例1中，我们是首先将原始的高光谱影像进行影像插值，这样通过影像插值可以得到与SAR影像相同像元尺寸的新的高光谱影像，也就是得到与SAR影像相同像元尺寸的HS影像。这样做的目的是将原始的高光谱影像的精度进一步提高。
S2，对新的高光谱影像进行HIS变换，得到I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像；
在本实施例1中我们是根据上述步骤S1得到的新的高光谱影像，对得到的新的高光谱影像进行HIS变换，将这新的高光谱影像进行分段，得到I、S和H分量，在进行HIS变换的同时，我们也对SAR影像进行灰度拉伸，得到SAR*影像。这样的目的是使得后续变换的时候能够更好地融合。
S3，对I分量和SAR*影像分别进行Curvelet变换，得到新的Curvelet变换系数；
在本实施例1中我们是对上述步骤S2中得到的SAR*影像和步骤S1中得到的I分量分别进行Curvelet变换，得到新的Curvelet变换系数，这样可以使得根据不同的融合规则，能够得到不同新的Curvelet变换系数。
S4，将新的Curvelet变换系数进行Curvelet逆变换，得到新的I*分量，同时将I*分量和S、H分量进行HIS逆变换，得到融合后的影像。
在本实施例1中我们是将上述步骤S3中得到的新的Curvelet变换系数进行Curvelet逆变换，通过Curvelet逆变换，我们可以得到I*分量，在我们进行Curvelet逆变换的同时，我们也将I*分量和S、H分量进行HIS逆变换，通过HIS逆变换，我们可以得到融合后的影像，从而实现本实施例1中的目的。通过上述实施例S1至S4步骤，我们对Curvelet变换的低频、中频和高频分量采用不同的融合规则，低频采用区域标准差和局部方向信息熵的融合规则，中频采用特征量积的融合规则，高频采用高频系数取零的融合规则，同时在融合规则中对噪声的进行抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
优选地，在本实施例1中我们对步骤S2进行了细化，其包括如下：所述S2中对SAR影像进行灰度拉伸，得到SAR*影像，其公式为：
其中，μ1，σ1分别为SAR影像的均值和标准差，μ0，σ0分别为I分量的均值和标准差。
优选地，我们还对上述步骤S3进行了细化，其包括如下：所述S3中包括：
S31，对I分量和SAR*影像分别进行Curvelet变换，分别得到与I分量对应的I分量低频系数、I分量中频系数、I分量高频系数和与SAR*影像对应的SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数；
S32，根据融合规则，将I分量低频系数、I分量中频系数、I分量高频系数分别对应地与SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数进行融合，对应地得到新的Curvelet变换的低频系数、中频系数、高频系数。
在本实施例1中我们是先对上述步骤S1中得到的I分量和上述步骤S2中得到的SAR*影像分别进行Curvelet变换，分别得到与I分量对应的I分量低频系数、I分量中频系数、I分量高频系数和与SAR*影像对应的SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数；在得到这些系数后，我们还根据不同的融合规则，将I分量低频系数、I分量中频系数、I分量高频系数分别对应地与SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数进行融合，对应地得到新的Curvelet变换的低频系数、中频系数、高频系数。通过这些不同的融合规则进行的Curvelet变换，低频采用区域标准差和局部方向信息熵的融合规则，中频采用特征量积的融合规则，高频采用高频系数取零的融合规则，同时融合规则对噪声的进行了抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
优选地，在本实施例1中我们还对上述步骤S4进行了细化，其包括如下：
S41，将新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到I*分量；
S42，将I*分量与S、H分量进行HIS反变换，得到融合后影像。
在本实施例1中我们是将上述步骤S32中得到的新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到新的I*分量；然后我们再将得到的I*分量与S、H分量进行HIS反变换，得到融合后的影像。
优选地，在本实施例1中我们还对上述步骤S2进行细化，其包括：在新的高光谱影像中选择三个波段作为R、G、B波段，并对该三个波段分别进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像。
在本实施例1中我们是在新的高光谱影像中选择三个波段作为R、G、B波段，并非是所有的波段都进行，这样的目的减少非必要的影响的转换，同时也节省很多的时间，提高了精度。我们再对选取的R、G、B波段进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到新的SAR*影像。
优选地，在本实施例1中我们还进一步地细化上述步骤S32，其中所述S32中的融合规则包括：区域标准差和局部方向信息熵的融合规则、特征量积的融合规则、高频系数取零的融合规则；
所述区域标准差和局部方向信息熵的融合规则，是用于I分量低频系数与SAR*影像低频系数进行融合，得到新的Curvelet变换的低频系数；
所述特征量积的融合规则，是用于I分量中频系数与SAR*影像中频系数进行融合，得到新的Curvelet变换的中频系数；
所述高频系数取零的融合规则，是用于I分量高频系数分与SAR*影像高频系数进行融合，得到新的Curvelet变换的高频系数。
在本实施例1中我们采用这样的不同的融合规则进行融合，对于低频系数，由于低频代表了影像的整体特征信息，为了尽可能多地保留地物的概貌信息，本实施例1中使用了以区域标准差和局部方向信息熵为主要判断依据的融合规则，更好地选择SAR影像与I分量的低频信息进行融合。
其中我们需要说明的是，所述区域标准差和局部方向信息熵的融合规则，局部方向信息熵可作为衡量区域显著性的准则，其值越大，表明该区域包含越丰富的信息，定义如下：
其中R表示以(k1,k2)为中心选取的局部区域，其大小为m×n。
记分别为PAN*和I的低频系数，SRPAN(k1,k2)、SRI(k1,k2)分别为的局部方向信息熵，σPAN(k1,k2)、σI(k1,k2)为的局部区域标准差。
另外，对于所述区域标准差和局部方向信息熵的融合规则，是用于I分量低频系数与SAR*影像低频系数进行融合，得到新的Curvelet变换的低频系数；其具体包括：
如果SRPAN(k1,k2)＞SRI(k1,k2)，
否则，
其中为融合后的低频系数，x1(k1,k2)、x2(k1,k2)分别为：
对于不同尺度及方向下的中频系数，由于绝对值取大的方法对噪声比较敏感，容易把噪声当作是有用信息注入到融合影像中，因此，融合中既需要考虑了局部窗口内整体特征，也需要考虑了每个像素的特征。特征量积使用了多个特征量表达细节信息，本实施例1中是采用基于特征量积的融合规则，尽可能多地保留有用的信息，减少噪声的影响。
其中，所述特征量积的融合规则，是用于I分量中频系数与SAR*影像中频系数进行融合，得到新的Curvelet变换的中频系数；其过程如下：
记分别为SAR*和I的中频系数，j为尺度参数，l为方向参数。分别为以(k1,k2)为中心，大小为n×n的局部区域的标准差，分别为(k1,k2)的梯度，则：
特征量积的定义如下：
其中分别为以(k1,k2)为中心，大小为n×n的局部区域的特征量积。
还有的是，得到的中频系数融合如下所示：
其中为融合后的中频系数。
最后我们还需要解析的是，所述高频系数取零的融合规则，是用于I分量高频系数分与SAR*影像高频系数进行融合，得到新的Curvelet变换的高频系数。其高频系数取零的融合规则是代表噪声的系数主要集中在高频系数中，本实施例1中采用高频系数取零的融合规则以尽可能地抑制SAR影像中的斑点噪声。
实施例2
如图3所示，本实施例2涉及一种基于Curvelet变换的影像融合系统，该系统包括：影像插值模块、影像处理模块、Curvelet变换模块、融合影像模块；
所述影像插值模块，其用于对原始的高光谱影像进行影像插值，得到与SAR影像相同像元尺寸的新的高光谱影像；
所述影像处理模块，其用于对新的高光谱影像进行HIS变换，得到I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像；
所述Curvelet变换模块，其用于对I分量和SAR*影像分别进行Curvelet变换，得到Curvelet变换系数，并根据融合规则对Curvelet变换系数进行融合，得到新的Curvelet变换系数；
所述融合影像模块，其用于将新的Curvelet变换系数进行Curvelet逆变换，得到I*分量，同时将I*和S、H分量进行HIS反变换，得到融合后的影像。
在本实施例2中对Curvelet变换的低频、中频和高频分量采用不同的融合规则，低频采用区域标准差和局部方向信息熵的融合规则，中频采用特征量积的融合规则，高频采用高频系数取零的融合规则，同时融合规则对噪声的进行了抑制，有效地改善了融合影像的光谱失真问题，极大地提高了空间分辨率。
优选地，所述Curvelet变换模块包括：第一Curvelet变换单元、第二Curvelet变换单元；
所述第一Curvelet变换单元，其用于对I分量和SAR*影像分别进行Curvelet变换，分别得到与I分量对应的I分量低频系数、I分量中频系数、I分量高频系数和与SAR*影像对应的SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数；
所述第二Curvelet变换单元，其用于根据融合规则，将I分量低频系数、I分量中频系数、I分量高频系数分别对应地与SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数进行融合，对应地得到新的Curvelet变换的低频系数、中频系数、高频系数。
优选地，所述融合影像模块中包括：第一融合影像单元、第二融合影像单元；
所述第一融合影像单元，其用于将新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到新的I*分量；
所述第二融合影像单元，其用于将I*分量与S、H分量进行HIS反变换，得到融合后的影像。
在本实施例2中我们需要说明的是，我们是先通过第一融合影像单元将将新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到I*分量；然后再通过第二融合影像单元将I*分量与S、H分量进行HIS反变换，得到融合后的影像。
优选地，所述影像处理模块，其具体用于在新的高光谱影像中选择三个波段作为R、G、B波段，并对该三个波段分别进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像。
在本实施例2中我们需要说明的是，我们是通过影像处理模块，在新的高光谱影像中选择三个波段作为R、G、B波段，并非是所有的波段都进行，这样的目的减少非必要的影响的转换，同时也节省很多的时间，提高了精度。我们再对选取的R、G、B波段进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像。
下面使我们根据上面的方法进行的一些数据试验，如下表所示：
由表1数据可以看出，与其他方法相比，研究的Curvelet_SAR方法得到的融合影像各波段的亮度均值和标准差数值均为最大，说明融合后的影像亮度较高，灰度级分布分散，融合影像可以反映出更多的特征信息；在信息熵指标上，Curvelet_SAR方法得到的融合影像各波段数值都是最大，说明与其他方法相比，该方法得到的融合影像中的信息量更为丰富。
在平均梯度指标上，由于Curvelet_SAR方法考虑了SAR影像中的斑点噪声对融合影像的影响，在高频系数融合规则中为了抑制斑点噪声而损失了部分细节信息，因此数值相对较低。
从表1中相关系数指标可以看出，文中方法得到的融合影像相关系数最高，其次是DWT_HIS方法，HIS方法得到的融合影像最低。因此说明Curvelet_SAR方法的光谱保持能力最好，HIS方法最差。在相对偏差指标上，DWT_HIS方法性能最优。
在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于Curvelet变换的影像融合方法，其特征在于,该方法包括如下步骤：
S1,对原始的高光谱影像进行影像插值，得到与SAR影像相同像元尺寸的新的高光谱影像；
S2，对新的高光谱影像进行HIS变换，得到I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像；
S3，对I分量和SAR*影像分别进行Curvelet变换，得到Curvelet变换系数，并根据融合规则对Curvelet变换系数进行融合，得到新的Curvelet变换系数；
S4，将新的Curvelet变换系数进行Curvelet逆变换，得到I*分量，同时将I*分量和S、H分量进行HIS逆变换，得到融合后的影像。
2.根据权利要求1所述的影像融合方法，其特征在于，所述S2中对SAR影像进行灰度拉伸，得到SAR*影像，其公式为：
<mrow>
<msup>
<mi>SAR</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mi>A</mi>
<mi>R</mi>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
</mrow>
其中，μ1，σ1分别为SAR影像的均值和标准差，μ0，σ0分别为I分量的均值和标准差。
3.根据权利要求1或2所述的影像融合方法，其特征在于，所述S3中包括：
S31，对I分量和SAR*影像分别进行Curvelet变换，分别得到与I分量对应的I分量低频系数、I分量中频系数、I分量高频系数和与SAR*影像对应的SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数；
S32，根据融合规则，将I分量低频系数、I分量中频系数、I分量高频系数分别对应地与SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数进行融合，对应地得到新的Curvelet变换的低频系数、中频系数、高频系数。
4.根据权利要求3所述的影像融合方法，其特征在于，所述S4中包括：
S41，将新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到I*分量；
S42，将I*分量与S、H分量进行HIS反变换，得到融合后的影像。
5.根据权利要求1或2所述的影像融合方法，其特征在于，所述S2中包括：
在新的高光谱影像中选择三个波段作为R、G、B波段，并对该三个波段分别进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像。
6.根据权利要求3所述的影像融合方法，其特征在于，所述S32中的融合规则包括：区域标准差和局部方向信息熵的融合规则、特征量积的融合规则、高频系数取零的融合规则；
所述区域标准差和局部方向信息熵的融合规则，用于I分量低频系数与SAR*影像低频系数进行融合，得到新的Curvelet变换的低频系数；
所述特征量积的融合规则，用于I分量中频系数与SAR*影像中频系数进行融合，得到新的Curvelet变换的中频系数；
所述高频系数取零的融合规则，用于I分量高频系数分与SAR*影像高频系数进行融合，得到新的Curvelet变换的高频系数。
7.一种利用权利要求1至6任一所述的影像融合方法的影像融合系统，其特征在于，该系统包括：影像插值模块、影像处理模块、Curvelet变换模块、融合影像模块；
所述影像插值模块，其用于对原始的高光谱影像进行影像插值，得到与SAR影像相同像元尺寸的新的高光谱影像；
所述影像处理模块，其用于对新的高光谱影像进行HIS变换，得到I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像；
所述Curvelet变换模块，其用于对I分量和SAR*影像分别进行Curvelet变换，得到Curvelet变换系数，并根据融合规则对Curvelet变换系数进行融合，得到新的Curvelet变换系数；
所述融合影像模块，其用于将新的Curvelet变换系数进行Curvelet逆变换，得到I*分量，同时将I*和S、H分量进行HIS反变换，得到融合后影像。
8.根据权利要求7所述的影像融合系统，其特征在于，所述Curvelet变换模块包括：第一Curvelet变换单元、第二Curvelet变换单元；
所述第一Curvelet变换单元，其用于对I分量和SAR*影像分别进行Curvelet变换，分别得到与I分量对应的I分量低频系数、I分量中频系数、I分量高频系数和与SAR*影像对应的SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数；
所述第二Curvelet变换单元，其用于根据融合规则，将I分量低频系数、I分量中频系数、I分量高频系数分别对应地与SAR*影像低频系数、SAR*影像中频系数、SAR*影像高频系数进行融合，对应地得到新的Curvelet变换的低频系数、中频系数、高频系数。
9.根据权利要求8所述的影像融合系统，其特征在于，所述融合影像模块中包括：第一融合影像单元、第二融合影像单元；
所述第一融合影像单元，其用于将新的Curvelet变换的低频系数、中频系数、高频系数进行Curvelet逆变换，得到I*分量；
所述第二融合影像单元，其用于将I*分量与S、H分量进行HIS反变换，得到融合后的影像。
10.根据权利要求7所述的影像融合系统，其特征在于，所述影像处理模块，其具体用于在新的高光谱影像中选择三个波段作为R、G、B波段，并对该三个波段分别进行HIS变换，得到对应的I、S和H分量，同时对SAR影像进行灰度拉伸，得到SAR*影像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710301849.6A CN107154020A (zh) | 2017-05-02 | 2017-05-02 | 一种基于Curvelet变换的影像融合方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710301849.6A CN107154020A (zh) | 2017-05-02 | 2017-05-02 | 一种基于Curvelet变换的影像融合方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107154020A true CN107154020A (zh) | 2017-09-12 |
Family
ID=59793512
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710301849.6A Pending CN107154020A (zh) | 2017-05-02 | 2017-05-02 | 一种基于Curvelet变换的影像融合方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107154020A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108764097A (zh) * | 2018-05-21 | 2018-11-06 | 北京理工大学 | 一种基于分段稀疏表示的高光谱遥感图像目标识别方法 |
| CN109118463A (zh) * | 2018-07-27 | 2019-01-01 | 中国科学院国家空间科学中心 | 一种基于hsl和图像熵的sar图像和光学图像融合方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110213556A1 (en) * | 2010-03-01 | 2011-09-01 | Bp Corporation North America Inc. | System and method for local attribute matching in seismic processing |
| CN102446351A (zh) * | 2010-10-15 | 2012-05-09 | 江南大学 | 一种多光谱与高分辨率全色图像融合方法研究 |
| CN102637297A (zh) * | 2012-03-21 | 2012-08-15 | 武汉大学 | 一种基于Curvelet变换的可见光与红外图像融合方法 |
| CN103198463A (zh) * | 2013-04-07 | 2013-07-10 | 北京航空航天大学 | 基于整体结构和空间细节信息融合的光谱图像全色锐化方法 |
| CN103198456A (zh) * | 2013-03-21 | 2013-07-10 | 西安电子科技大学 | 基于方向波域隐马尔可夫树模型的遥感图像融合方法 |
-
2017
- 2017-05-02 CN CN201710301849.6A patent/CN107154020A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110213556A1 (en) * | 2010-03-01 | 2011-09-01 | Bp Corporation North America Inc. | System and method for local attribute matching in seismic processing |
| CN102446351A (zh) * | 2010-10-15 | 2012-05-09 | 江南大学 | 一种多光谱与高分辨率全色图像融合方法研究 |
| CN102637297A (zh) * | 2012-03-21 | 2012-08-15 | 武汉大学 | 一种基于Curvelet变换的可见光与红外图像融合方法 |
| CN103198456A (zh) * | 2013-03-21 | 2013-07-10 | 西安电子科技大学 | 基于方向波域隐马尔可夫树模型的遥感图像融合方法 |
| CN103198463A (zh) * | 2013-04-07 | 2013-07-10 | 北京航空航天大学 | 基于整体结构和空间细节信息融合的光谱图像全色锐化方法 |
Non-Patent Citations (7)
| Title |
|---|
| MOHAMED ELHABIBY: "：Second Generation Curvelet Transforms Vs Wavelet transforms and Canny Edge Detector for Edge Detection from WorldView-2 data，Mohamed Elhabiby", 《INTERNATIONAL JOURNAL OF COMPUTER SCIENCE & ENGINEERING SURVEY》 * |
| 何贵青: "一种基于特征量积的遥感图像融合方法", 《计算机工程与应用》 * |
| 利用快速离散CURVELET变换的遥感影像融合: "利用快速离散Curvelet变换的遥感影像融合", 《武汉大学学报》 * |
| 李晖晖: "基于Curvelet变换的SAR与可见光图像融合研究", 《光电子•激光》 * |
| 杨旸: "多传感器图像融合方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 肖化超: "基于 IHS 变换和 Curvelet 变换的卫星遥感图像融合方法", 《华南理工大学学报》 * |
| 马苗: "基于Curvelet变换的图像质量多尺度灰评价", 《计算机工程与应用》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108764097A (zh) * | 2018-05-21 | 2018-11-06 | 北京理工大学 | 一种基于分段稀疏表示的高光谱遥感图像目标识别方法 |
| CN108764097B (zh) * | 2018-05-21 | 2021-07-20 | 北京理工大学 | 一种基于分段稀疏表示的高光谱遥感图像目标识别方法 |
| CN109118463A (zh) * | 2018-07-27 | 2019-01-01 | 中国科学院国家空间科学中心 | 一种基于hsl和图像熵的sar图像和光学图像融合方法 |
| CN109118463B (zh) * | 2018-07-27 | 2021-10-19 | 中国科学院国家空间科学中心 | 一种基于hsl和图像熵的sar图像和光学图像融合方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109118431B (zh) | 一种基于多记忆及混合损失的视频超分辨率重建方法 | |
| CN108304755B (zh) | 用于图像处理的神经网络模型的训练方法和装置 | |
| Battiato et al. | A locally adaptive zooming algorithm for digital images | |
| CN104272739B (zh) | 用于不同颜色取样格式的基于模式的系数扫描和对角变换 | |
| US8724894B1 (en) | Colorization of digital imagery | |
| WO2016045242A1 (zh) | 一种图像放大方法、图像放大装置及显示设备 | |
| JPWO2005067294A1 (ja) | 画像処理方法、画像処理装置、サーバクライアントシステム、サーバ装置、クライアント装置および画像処理システム | |
| DE102018119625A1 (de) | Reduzieren von strukturierten IR-Mustern bei der Stereoskopischen Tiefensensorbildgebung | |
| CN101295399A (zh) | 分辨率增强设备和方法 | |
| CN102893609B (zh) | 图像处理设备和图像处理设备的控制方法 | |
| CN111260580A (zh) | 一种基于图像金字塔的图像去噪方法、计算机装置及计算机可读存储介质 | |
| CN105340268A (zh) | 图像处理装置、图像处理方法以及图像处理程序 | |
| CN116051428B (zh) | 一种基于深度学习的联合去噪与超分的低光照图像增强方法 | |
| CN113781320A (zh) | 一种图像处理方法、装置、终端设备及存储介质 | |
| CN104156911A (zh) | 一种图像融合处理方法及系统 | |
| CN102053804A (zh) | 图像处理装置及控制方法 | |
| CN111800629A (zh) | 视频解码方法、编码方法以及视频解码器和编码器 | |
| CN106056565B (zh) | 一种基于多尺度形态学双边滤波分解和对比度压缩的mri与pet图像融合方法 | |
| CN107154020A (zh) | 一种基于Curvelet变换的影像融合方法及系统 | |
| CN106875370B (zh) | 一种全色图像和多光谱图像的融合方法及装置 | |
| CN109859153B (zh) | 一种基于自适应光谱-空间梯度稀疏正则化的多光谱图像融合方法 | |
| CN103517008B (zh) | 保持锐度不变的同时增强视频数据的对比度的方法和系统 | |
| CN106485664B (zh) | 一种基于小波变换和Wallis变换的卫星图像色彩平衡方法 | |
| Ghimire et al. | Color image enhancement in HSV space using nonlinear transfer function and neighborhood dependent approach with preserving details | |
| CN101917624A (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170912 |