CN103530876B - 一种基于傅里叶变换的红外图像分析方法 - Google Patents
一种基于傅里叶变换的红外图像分析方法 Download PDFInfo
- Publication number
- CN103530876B CN103530876B CN201310469916.7A CN201310469916A CN103530876B CN 103530876 B CN103530876 B CN 103530876B CN 201310469916 A CN201310469916 A CN 201310469916A CN 103530876 B CN103530876 B CN 103530876B
- Authority
- CN
- China
- Prior art keywords
- image
- infrared
- fourier transform
- frequency spectrum
- spectrogram
- 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
Landscapes
- Spectrometry And Color Measurement (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明公开了一种基于傅里叶变换的红外图像分析方法,它主要应用于红外辐射强度图像和偏振度图像分析。它首先对待分析的红外辐射强度图像和红外偏振度图像进行傅里叶变换,得到各自的频谱。然后对得到的频谱进行中心化,得到中心化的频谱。然后根据所创建的分区模型,根据不同的n值,对辐射强度图像和偏振度图像的频谱特点进行分析。实验结果表明该方法能够较好说明红外图像的频谱分布特征,可以为红外图像进一步处理提供理论依据。
Description
技术领域:
本发明属于图像探测与处理领域,具体涉及一种基于傅里叶变换的红外图像分析方法。
背景技术:
目前,红外图像融合已经成为一种重要的红外图像处理技术,涉及到可见光与红外图像融合、中波红外与长波红外图像融合、红外辐射强度图像与偏振度图像融合等领域。其中数据层融合通常用于多源图像组合,多源图像组合是将不同传感器获得的同一场景的图像经过匹配、合成处理,最终获得一幅克服单一传感器图像在几何、光谱和空间分辨率上存在局限性的优质图像,以利于进一步对其有效应用。
对不同的红外图像进行融合,首先要对图像的特点进行分析,因此设计一种红外图像分析方法非常有必要。采用傅里叶分析手段,对红外图像进行分析,可以得出红外图像频谱分布特点,为红外图像融合提供理论依据。
发明内容:
针对红外偏振图像融合的实际需要,本发明提出一种基于傅里叶变换的红外图像分析方法,事实表明用此方法可以对红外图像频谱分布特点进行分析。
实现本发明的技术解决方案为:一种基于傅里叶变换的红外图像分析方法,步骤如下:
步骤1:对待分析红外图像X进行二维离散傅里叶变换,得到图像的频率域F(μ,ν),它为复数阵列,可以表示为F(μ,ν)=R(μ,ν)+jI(μ,ν),然后根据公式|F(μ,ν)|=[R2(μ,ν)+I2(μ,ν)]1/2,得到其频谱图|F(μ,ν)|。
步骤2:中心化该频谱图,将|F(μ,ν)|的原点移到频谱矩形的中心,得到|FC(μ,ν)|。
步骤3:如下图所示,对得到的|FC(μ,ν)|,我们定义了频谱能量集中度p。以频谱图中心像素点为中心,以n个像素为半径的圆内所有像素灰度值之和比上整幅图像的像素灰度值之和。其中n为整数,对一幅M×N的图像,其取值范围为:1~min(M/2,N/2)。
具体方法如下:
步骤4:对红外辐射强度图像和红外偏振度图像,我们分别求得它们在n取不同值时的p值,从而估计它们的频谱分布情况,得出红外辐射强度图像和红外偏振度图像的频谱特点。
附图说明
图1:基于傅里叶变换红外图像分析方法框图。
图2:一幅256×256红外辐射强度图像。
图3:一幅256×256红外偏振度图像。
图4:频谱图分区方式。
具体实施方式:
下面根据附图对本发明的具体实施方式作进一步的详细说明。
图1是基于傅里叶变换的红外图像分析方法框图。
图2是一幅256×256红外辐射强度图像X。
图3是一幅256×256红外偏振度图像Y。
图4是频谱图分区方式。
步骤1:对图2红外辐射强度图像X进行二维离散傅里叶变换,得到其频率域FX(μ,ν),它为复数阵列,可以表示为FX(μ,ν)=R(μ,ν)+jI(μ,ν),然后根据公式|FX(μ,ν)|=[R2(μ,ν)+I2(μ,ν)]1/2,得到其频谱图|FX(μ,ν)|。对图3红外偏振度图像Y,按照上述步骤得到FY(μ,ν)和|FY(μ,ν)|。其中图像X和图像Y大小均为256×256个像素。
步骤2:中心化该频谱图,将|FX(μ,ν)|和|FY(μ,ν)|的原点移到频谱矩形的中心,得到|FCX(μ,ν)|和|FCY(μ,ν)|。
步骤3:如图4所示,对得到的|FCX(μ,ν)|和|FCY(μ,ν)|,我们定义了频谱能量集中度p。以频谱图中心像素点为中心,以n个像素为半径的圆内所有像素灰度值之和比上整幅图像的像素灰度值之和,即以A中所有像素灰度值之和比上整幅图像灰度值之和。其中n为整数,对图2和图3两幅256×256的图像,其取值范围为:1~128。
具体方法步骤如下:
步骤4:对红外辐射强度图像和红外偏振度图像,我们分别求得它们在n取不同值时的p值,从而估计它们的频谱分布情况,得出红外辐射强度图像和红外偏振度图像的频谱特点。下表给出了n取10、20、30、40、50时的p值。从中可以看出红外辐射强度图像和红外偏振度图像频谱分布情况,如下表所示。
P | 红外辐射强度频谱图 | 红外偏振度频谱图 |
P,n=10 | 0.6263 | 0.0740 |
P,n=20 | 0.8092 | 0.1490 |
P,n=30 | 0.8780 | 0.2273 |
P,n=40 | 0.9156 | 0.3162 |
P,n=50 | 0.9390 | 0.4021 |
Claims (1)
1.一种基于傅里叶变换的红外图像分析方法,其特征在于包括以下步骤:
步骤1:对待分析红外图像X进行二维离散傅里叶变换,得到图像的频率域F(μ,ν),它为复数阵列,可以表示为F(μ,ν)=R(μ,ν)+jI(μ,ν),然后根据公式:
|F(μ,ν)|=[R2(μ,ν)+I2(μ,ν)]1/2
得到其频谱图|F(μ,ν)|;
步骤2:中心化该频谱图,将|F(μ,ν)|的原点移到频谱矩形的中心,得到|FC(μ,ν)|;
步骤3:如下图所示,对得到的|FC(μ,ν)|,定义频谱能量集中度p,以频谱图中心像素点为中心,以n个像素为半径的圆内所有像素灰度值之和比上整幅图像的像素灰度值之和,其中n为整数,对一幅M×N的图像,其取值范围为:1~min(M/2,N/2);
具体方法如下:
步骤4:对红外辐射强度图像和红外偏振度图像,分别求得它们在n取不同值时的p值,从而估计它们的频谱分布情况,得出红外辐射强度图像和红外偏振度图像的频谱特点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310469916.7A CN103530876B (zh) | 2013-10-10 | 2013-10-10 | 一种基于傅里叶变换的红外图像分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310469916.7A CN103530876B (zh) | 2013-10-10 | 2013-10-10 | 一种基于傅里叶变换的红外图像分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103530876A CN103530876A (zh) | 2014-01-22 |
CN103530876B true CN103530876B (zh) | 2016-03-30 |
Family
ID=49932855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310469916.7A Active CN103530876B (zh) | 2013-10-10 | 2013-10-10 | 一种基于傅里叶变换的红外图像分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103530876B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104103073A (zh) * | 2014-07-14 | 2014-10-15 | 中国人民解放军国防科学技术大学 | 一种红外偏振图像边缘检测方法 |
CN105373809B (zh) * | 2015-11-06 | 2018-10-09 | 重庆大学 | 基于非负最小二乘稀疏表示的sar目标识别方法 |
CN109472805B (zh) * | 2018-10-09 | 2020-10-02 | 哈工大机器人(山东)智能装备研究院 | 一种基于傅里叶变换的红外图像自适应分层方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7620265B1 (en) * | 2004-04-12 | 2009-11-17 | Equinox Corporation | Color invariant image fusion of visible and thermal infrared video |
CN102222339A (zh) * | 2011-06-17 | 2011-10-19 | 中国科学院自动化研究所 | 一种基于纹理和强度特征融合的多尺度背景建模方法 |
CN102881010A (zh) * | 2012-08-28 | 2013-01-16 | 北京理工大学 | 基于人眼视觉特性的融合图像感知清晰度评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6843176B2 (en) * | 2002-04-26 | 2005-01-18 | Kodak Polychrome Graphics, Llc | Method to remove unwanted, unexposed, radiation-sensitive layer in a lithographic printing plate |
-
2013
- 2013-10-10 CN CN201310469916.7A patent/CN103530876B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7620265B1 (en) * | 2004-04-12 | 2009-11-17 | Equinox Corporation | Color invariant image fusion of visible and thermal infrared video |
CN102222339A (zh) * | 2011-06-17 | 2011-10-19 | 中国科学院自动化研究所 | 一种基于纹理和强度特征融合的多尺度背景建模方法 |
CN102881010A (zh) * | 2012-08-28 | 2013-01-16 | 北京理工大学 | 基于人眼视觉特性的融合图像感知清晰度评价方法 |
Non-Patent Citations (6)
Title |
---|
A wavelet-based adaptive fusion algorithm of infrared polarization imaging;Yang weil等;《International Symposium on Photoelectronic Detection and Imaging 2011: Advances in Infrared Imaging and Applications》;20110908;第8193卷;第819324-1页-第819324-9页 * |
Fusion of infrared polarization and intensity images using support value transform and fuzzy combination rules;Yang Fengbao等;《Infrared Physics & Technology》;20130525;第60卷;第235-243页 * |
Modulation transfer function measurement of an infrared focal plane array by use of the self-imaging property of a canted periodic target;Nicolas Guerineau等;《APPLIED OPTICS》;19990201;第38卷(第4期);第631-637页 * |
不同视场的多源图像融合技术研究;田思等;《兵工学报》;20120131;第33卷(第1期);第26-31页 * |
基于边缘特征和能量的图像融合方法;罗南超等;《计算机工程》;20100831;第36卷(第15期);第202-203页 * |
基于频谱分析和图像融合的红外弱小目标增强;王晓文等;《中国电子科学研究院学报》;20130630;第8卷(第3期);第258-265页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103530876A (zh) | 2014-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yue et al. | Deep recursive super resolution network with Laplacian Pyramid for better agricultural pest surveillance and detection | |
CN102436652B (zh) | 一种多源遥感图像自动配准方法 | |
CN101510309B (zh) | 基于复小波纹理区域合并的改进分水岭sar图像分割方法 | |
CN104240244B (zh) | 一种基于传播模式和流形排序的显著性物体检测方法 | |
CN104200452A (zh) | 基于谱图小波变换的红外与可见光图像融合方法及其装置 | |
CN103455813A (zh) | 一种ccd图像测量系统光斑中心定位的方法 | |
CN103530876B (zh) | 一种基于傅里叶变换的红外图像分析方法 | |
CN102184557A (zh) | 一种复杂场景的显著区域检测方法 | |
CN103578091A (zh) | 一种基于拉普拉斯金字塔的红外偏振图像融合方法 | |
CN105335972A (zh) | 基于小波轮廓波变换与视觉显著性的经编织物疵点检测方法 | |
CN103325105A (zh) | 一种高精度合成孔径雷达图像自动配准方法及设备 | |
CN104574301A (zh) | 一种太赫兹图像重构方法及系统 | |
CN103473749B (zh) | 一种基于全变分图像融合的方法及装置 | |
CN104252703A (zh) | 基于小波预处理和稀疏表示的卫星遥感图像超分辨率重构方法 | |
CN103345742B (zh) | 基于一种改进马尔可夫随机场模型的遥感图像变化的检测方法 | |
CN104102900A (zh) | 一种车辆识别系统 | |
CN103310461A (zh) | 基于块卡尔曼滤波的图像边缘提取方法 | |
CN103530860B (zh) | 基于自适应的自回归模型的高光谱图像超分辨率方法 | |
CN102945550A (zh) | 一种基于高斯尺度空间构建遥感影像语义的方法 | |
CN102903087B (zh) | 基于gpu编程的sar图像去噪方法 | |
CN102542547A (zh) | 一种基于光谱约束的高光谱图像融合方法 | |
CN110473145B (zh) | 被动式太赫兹图像融合超高分辨率重建的检测方法及装置 | |
CN104574304A (zh) | 一种毫米波图像重构方法及系统 | |
CN102855612A (zh) | 基于灰度线阵ccd图像的自适应增强算法 | |
CN103065310A (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 |