简介：这份资源围绕阿贝成像原理与空间滤波的MATLAB实现展开，面向学习傅里叶光学、信息光学课程的高校学生及需要完成相关实验仿真的研究者。它用空间频谱语言分析物光场结构，通过有意识地改变物频谱来获得期望的像，并以一维矩形光栅在4f光学系统中的频谱分析为核心案例，帮助读者理解从频域滤波到空域成像的完整链路。压缩包共6个文件，全部为m脚本，整体约4KB，涵盖主程序、滤波函数、卷积运算及正反傅里叶变换等模块，结构紧凑、便于逐行阅读与二次修改。资源重点复现两类典型情形：仅0级与0、±1级通过时的光栅像分布，以及去除0级空间频率后的像分布，读者可据此直观对比滤波前后图像变化，掌握频谱选择对成像结果的影响。目前已有1935人学习下载，适合作为光学实验仿真与MATLAB编程练习的参考素材。
1. 阿贝成像与空间滤波：从一次“像消失了”的实验说起
把一块透射光栅放在显微镜物镜下，理论上目镜里应该看到清晰的条纹。但如果你在物镜后焦面上放一个可调光阑，慢慢缩小它，会看到一件反直觉的事：条纹对比度先变好，然后突然整体消失，只剩一片均匀亮斑。这不是设备坏了，而是阿贝成像原理在起作用——成像本质上是衍射光被物镜收集后，在像平面干涉重建的结果。阿贝成像及空间滤波的 MATLAB 实现，就是用数值方法把“衍射—频谱—滤波—重建”这条链路完整跑一遍，让你能任意修改频谱面上的掩模，观察像怎么变。它适合光学、图像处理方向的学生和工程师，用来理解 4f 系统、傅里叶光学和频域滤波的物理对应关系。相比纯公式推导，代码能让你亲手把某个频率分量抠掉，看看到底哪部分信息在负责边缘、哪部分在负责轮廓。
2. 用 MATLAB 搭一个 4f 系统：从物面到频谱面再到像面
2.1 为什么选 FFT 而不是卷积来模拟衍射
阿贝成像的数学核心是菲涅耳衍射或夫琅禾费衍射。物面透射率函数经过一段距离传播后，在物镜后焦面上形成频谱分布。严格求解需要做卷积或傅里叶变换。MATLAB 里最直接的做法是用fft2把物面复振幅变换到频域，把频域坐标映射到后焦面的空间坐标，再乘上光瞳函数和滤波掩模，最后用ifft2回到像面。这样做的好处是：频谱面和像面都是矩阵，你可以像修图一样修改频谱面上的任意区域，物理意义清晰，计算量也可控。
常见做法是设物面采样数为 512×512 或 1024×1024，物面尺寸对应实际光栅周期。如果物面是周期结构，频谱面上会出现离散的衍射级次，这时滤波掩模只需要挡住或通过特定级次即可。如果物面是更复杂的图像，频谱连续分布，滤波掩模就相当于一个频域滤波器。
2.2 生成物面：光栅、字母和自定义图案
先构造一个透射光栅作为物面。下面代码生成一个周期为 16 像素、占空比 50% 的一维光栅，并扩展到二维。
% 参数设置 N = 512; % 采样点数 period = 16; % 光栅周期（像素） duty = 0.5; % 占空比 x = 1:N; grating1d = double(mod(x, period) < period * duty); % 扩展成二维光栅（只沿 x 方向变化） object = repmat(grating1d, N, 1); % 归一化到 0~1 object = object / max(object(:));这段代码里mod(x, period) < period * duty产生一个方波，repmat把它复制到每一行，形成沿 x 方向周期变化、y 方向均匀的透射光栅。如果你想要二维光栅，可以对 y 方向也做同样操作再相乘。物面矩阵object的每个元素代表该点的复振幅透射率，这里假设只有振幅调制，相位为 0。
如果想用字母或简单图案，可以直接读入二值图像或手动画一个。比如用insertText在空白矩阵上写一个字母，再二值化。物面越复杂，频谱越连续，滤波效果越像频域滤波。
2.3 传播到频谱面：fftshift 与坐标映射
物面到频谱面的传播用傅里叶变换实现。为了让零频分量在中心，需要做fftshift。同时要建立频率坐标与后焦面空间坐标的对应关系，方便后面画图。
% 傅里叶变换到频谱面 F = fft2(object); F = fftshift(F); % 零频移到中心 % 频率坐标（归一化） fx = (-N/2:N/2-1) / N; % 单位：1/像素 [FX, FY] = meshgrid(fx, fx); % 频谱强度（对数显示用） spectrum = log(1 + abs(F));fft2得到的是离散频谱，fftshift把零频搬到矩阵中心。fx是归一化频率，范围 -0.5 到 0.5。实际后焦面上的空间坐标与频率成正比：x_f = λ f fx，其中 λ 是波长，f 是物镜焦距。在仿真里我们只关心相对位置，所以直接用归一化频率作为坐标。
2.4 加载滤波掩模：低通、高通和方向滤波
滤波掩模就是一个与频谱矩阵同尺寸的矩阵，值在 0 到 1 之间。0 表示完全挡住，1 表示完全通过。下面生成三种典型掩模：低通、高通和方向滤波。
% 半径参数 r_low = 30; % 低通半径（像素） r_high = 20; % 高通半径（像素） % 计算每个点到中心的距离 R = sqrt(FX.^2 + FY.^2) * N; % 换算成像素距离 % 低通掩模：半径内通过 mask_low = double(R <= r_low); % 高通掩模：半径外通过 mask_high = double(R >= r_high); % 方向滤波：只通过水平方向频率（|fy| 小） mask_dir = double(abs(FY) < 0.02);R是频率面上每个点到零频的距离，乘以 N 后单位变成像素，方便设定半径。低通掩模保留低频，对应图像的大轮廓；高通掩模保留高频，对应边缘和细节；方向滤波只保留某一方向的频率分量，可以用来提取特定方向的条纹。
2.5 逆变换重建像面并对比
滤波后的频谱做逆变换回到像面，取模得到强度分布。
% 应用滤波 F_filtered = F .* mask_low; % 换成 mask_high 或 mask_dir 看不同效果 % 逆变换 image_filtered = ifft2(ifftshift(F_filtered)); image_filtered = abs(image_filtered); % 显示 figure; subplot(1,3,1); imshow(object, []); title('物面'); subplot(1,3,2); imshow(spectrum, []); title('频谱面（对数）'); subplot(1,3,3); imshow(image_filtered, []); title('滤波后像面');ifftshift把零频移回左上角，再ifft2得到复振幅，取模就是强度。注意ifft2的结果可能有微小虚部，abs后正常。显示时用[]自动缩放。
把mask_low换成mask_high，你会看到像面只剩下边缘；换成mask_dir，只剩水平方向的周期结构。这就是空间滤波的直观效果。
3. 参数怎么调：采样、周期、滤波半径与显示
3.1 采样数 N 和物面周期的关系
采样数 N 决定了频率分辨率。N 越大，频谱面上离散点越密，滤波掩模的边界越平滑。但 N 太大时计算变慢，而且如果物面周期不是整数像素，频谱会泄漏。建议让物面周期是整数像素，且 N 是周期的整数倍。例如周期 16，N 取 512，正好 32 个周期，频谱干净。如果周期是 15.3，频谱会出现大量旁瓣，滤波时容易误伤。
常见做法是先用小 N（如 256）快速试参数，确认效果后再用大 N 出图。如果物面是复杂图像，N 至少 512，否则高频细节被截断。
3.2 滤波半径的物理含义与选取
低通半径r_low对应允许通过的最高空间频率。半径越小，像越模糊；半径越大，细节越多。对于周期光栅，如果只保留零频和 ±1 级，像面会出现正弦条纹；如果只保留零频，像面变成均匀亮斑。你可以通过观察频谱面上衍射级次的位置来设定半径：先显示spectrum，找到各级次到中心的像素距离，然后让r_low刚好包含需要的级次。
高通半径r_high类似，但反过来。方向滤波的宽度参数0.02是归一化频率，对应实际角度。如果物面是水平条纹，水平方向频率集中在 fy=0 附近，用abs(FY) < 0.02就能提取。
3.3 显示动态范围：为什么频谱要取对数
频谱幅度动态范围很大，零频分量可能比高频高几个数量级。直接imshow(abs(F))只能看到中心一个亮点。取对数log(1+abs(F))压缩动态范围，才能看到衍射级次和连续分布。像面显示一般用线性，但如果像面整体偏暗，可以用imshow(image_filtered, [])自动拉伸。
注意log里加 1 是防止零值取对数出错。如果频谱中有零，log(0)会返回 -Inf，显示异常。
3.4 用循环批量测试不同掩模
手动改掩模效率低，可以写一个循环，把几种掩模依次应用并保存图像。
masks = {mask_low, mask_high, mask_dir}; names = {'低通', '高通', '方向'}; figure; for k = 1:3 F_f = F .* masks{k}; img = abs(ifft2(ifftshift(F_f))); subplot(1,3,k); imshow(img, []); title(names{k}); end这个循环把三个掩模的结果画在一张图上，方便对比。注意每次都要用ifftshift和ifft2，不能复用上一次的结果。
4. 避坑与排查：那些让像“消失”或“全黑”的细节
4.1 像面全黑或全白：fftshift 用反了
现象：滤波后像面要么全黑，要么全白，看不到任何结构。原因：fft2之后没有fftshift，或者逆变换时ifftshift和fftshift用混了。MATLAB 的fft2零频在左上角，fftshift把它移到中心。如果你在中心画了掩模，但没做fftshift，掩模就作用在错误的位置。解决：记住顺序——正变换fft2→fftshift→ 乘掩模 →ifftshift→ifft2。可以用一个简单测试：物面全 1，频谱只有零频一个点，滤波后应该还是全 1。
4.2 频谱显示只有中心一个亮点
现象：imshow(abs(F), [])只看到中心一个白点，周围全黑。原因：零频分量幅度远大于高频，线性显示时高频被压缩到看不见。解决：用log(1+abs(F))显示，或者用abs(F).^0.3做伽马校正。注意对数显示后零频仍然最亮，但衍射级次可见。
4.3 滤波后出现网格状伪影
现象：像面出现周期性网格或条纹，不是物面本身的结构。原因：物面周期与采样数不匹配，导致频谱泄漏，或者滤波掩模边缘太硬，产生振铃。解决：调整物面周期为整数像素，并让 N 是周期的整数倍；把硬边掩模换成高斯型软边掩模，例如exp(-(R/r_low).^2)。软边能减少吉布斯现象。
4.4 逆变换结果有虚部
现象：ifft2结果不是实数，取模后像面正常但数值有微小虚部。原因：浮点计算误差，或者滤波掩模不是共轭对称的。如果掩模在频域不是中心对称，逆变换结果会是复数。解决：对于实物面，滤波掩模应该满足中心对称（即mask(i,j) = mask(N-i+1, N-j+1)）。低通和高通掩模天然对称，方向滤波如果只取一侧，结果会是复数，取模即可，但物理上对应单边带滤波。
4.5 物面归一化问题导致像面亮度异常
现象：像面整体偏暗或偏亮，对比度不对。原因：物面矩阵没有归一化，或者ifft2后没有除以 N^2。MATLAB 的ifft2已经包含 1/N^2 因子，所以正反变换后能量守恒。但如果物面值在 0~255，频谱幅度会很大，滤波后像面也大。解决：把物面归一化到 0~1，或者显示时用[]自动缩放。不要手动再除 N^2，否则像面会暗 N^2 倍。
5. 进阶技巧：用软边掩模和相位滤波看更多细节
硬边掩模在频域相当于一个矩形窗，逆变换后会产生振铃。把低通掩模改成高斯型，像面过渡更自然。代码只需把mask_low = double(R <= r_low)换成mask_low = exp(-(R/r_low).^2)。高斯掩模没有严格截止频率，但能量集中在半径内，效果更平滑。
另一个技巧是相位滤波。前面只改了幅度，如果只改相位而保持幅度不变，可以看到相位对成像的影响。例如把频谱幅度设为 1，只保留原始相位，逆变换后像面会变成一种“边缘增强”的效果。代码：
F_phase = exp(1i * angle(F)); % 只保留相位 img_phase = abs(ifft2(ifftshift(F_phase)));这行代码把每个频谱点的幅度归一化，只留相位。逆变换后像面会突出边缘和轮廓，因为相位携带了位置信息。你可以对比原始像面和相位重建像面，理解幅度和相位各自的作用。
还有一个实用技巧：把滤波掩模做成可交互的。用imshow显示频谱，然后用drawcircle或roipoly让用户手动选区域，再生成掩模。这样不用每次改代码，适合探索性实验。例如：
figure; imshow(spectrum, []); h = drawcircle('Center', [N/2, N/2], 'Radius', 30); mask_custom = createMask(h);drawcircle返回一个 ROI 对象，createMask生成二值掩模。注意drawcircle的坐标是 (x,y)，而矩阵索引是 (行,列)，中心在 (N/2, N/2) 附近。生成掩模后同样做F .* mask_custom再逆变换。
我自己的习惯是：每次改完掩模，先看频谱面上哪些级次被挡住，再看像面变化，最后把物面、频谱、像面三张图并排保存。这样即使过了很久，也能一眼看出当时用的什么滤波。另外，物面周期不要设得太小，否则衍射级次间隔太远，滤波半径不好选。一般周期 8~32 像素比较合适。希望帮到你。
本文还有配套的精品资源，点击获取
