傅里叶单像素成像
文章目录
张子邦, 陆天傲, 彭军政,等.傅里叶单像素成像技术与应用[J]. 红外与激光工程, 2019, 48(6).、
数学原理
通过获取图像的傅里叶谱来重建图像本身。图像的傅里叶谱通过对目标物体投影一些列不同的傅里叶基地图案,并通过单像素探测器测量所得的光强值来获得。图像的重建过程仅仅是依次二维离散逆傅里叶变换。利用自然图像在傅里叶域的稀疏性以实现高效率的成像。
对一个目标物体进行二维成像,本质上就是一个对该目标物体的光强分布的空间信息获取和再现过程。二维单像素成像,就是通过仅适用一个没有空间分辨能力的探测器——单像素探测器——来获取物体的空间信息,并对所获取的信息施以响应的算法来重建物体图像。由此可见,二维单像素成像本身包含了一对关于空间分辨能力的矛盾,即:成像过程对空间分辨能力的需求而单像素探测器缺乏空间分辨能力的矛盾。
基于傅里叶分析,1807年,J.Fourier提出傅里叶分析理论,指出任意一个来纳许的信号都可以通过傅里叶基得线性组合构成(傅里叶正变换的意义)。
换言之,任意一个连续信号都可以适用无穷多个空间频率不同的正弦或余弦波进行线性叠加(傅里叶逆变换的物理意义)
二维傅里叶正变换和逆变换:
C ( f x , f y ) = ∫ − ∞ + ∞ ∫ − ∞ + ∞ I ( x , y ) e x p [ − j ⋅ 2 π ( f x x + f y y ) ] d f x d f y C(f_x,f_y)=\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}I(x,y)exp\left[-j\cdot2\pi(f_xx+f_yy)\right]df_xdf_y C(fx,fy)=∫−∞+∞∫−∞+∞I(x,y)exp[−j⋅2π(fxx+fyy)]dfxdfy
I ( x , y ) = ∫ − ∞ + ∞ ∫ − ∞ + ∞ C ( f x , f y ) e x p [ − j ⋅ 2 π ( f x x + f y y ) ] d x d y I(x,y)=\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}C(f_x,f_y)exp\left[-j\cdot2\pi(f_xx+f_yy)\right]dxdy I(x,y)=∫−∞+∞∫−∞+∞C(fx,fy)exp[−j⋅2π(fxx+fyy)]dxdy
- x,y是空间域中的直角坐标;
- f x , f y f_x,\;\;f_y fx,fy是傅里叶域中的直角坐标,对应x,y方向的空间频率;
- I ( x , y ) I(x,y) I(x,y)表示二维图像
- C ( f x , f y ) C(f_x,f_y) C(fx,fy)表示二维图像的傅里叶谱
- j是虚数单位
扩展到离散信号:
C ( u , v ) = ∑ x = 0 M − 1 ∑ y = 0 N − 1 I ( x , y ) e x p [ − j 2 π ( u x M + v y N ) ] C(u,v)=\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}I(x,y)exp\left[-j2\pi\left(\frac{ux}{M}+\frac{vy}{N}\right)\right] C(u,v)=x=0∑M−1y=0∑N−1I(x,