簡介傅里葉光學(xué)4f系統(tǒng)邊緣提取資源包面向光學(xué)信息處理與圖像處理方向的課程設(shè)計(jì)、畢業(yè)設(shè)計(jì)及學(xué)習(xí)參考人員。資源從理論推導(dǎo)出發(fā)研究邊緣提取算子尺寸與空間復(fù)濾波器空間分布之間的嚴(yán)格關(guān)系并給出濾波器振幅與位相分布的設(shè)定依據(jù)同時(shí)利用多張圖片驗(yàn)證邊緣提取效果幫助學(xué)習(xí)者完整掌握從理論建模到系統(tǒng)實(shí)現(xiàn)的鏈路。壓縮包共19個(gè)文件含Matlab源碼、14張實(shí)驗(yàn)圖像、PDF報(bào)告及說明文檔源碼覆蓋邊緣提取核心實(shí)現(xiàn)圖像展示輸入、頻譜、相位及邊緣結(jié)果報(bào)告詳述推導(dǎo)過程與參數(shù)設(shè)計(jì)包體僅1.2MB結(jié)構(gòu)清晰、便于下載。該資源已有679人學(xué)習(xí)適合具備基礎(chǔ)傅里葉光學(xué)知識(shí)、希望動(dòng)手實(shí)現(xiàn)光學(xué)圖像邊緣提取的本科生和研究生深入?yún)⒖寂c復(fù)用。1. 傅里葉光學(xué)中的4f系統(tǒng)邊緣提取一次透鏡變換取代逐像素卷積邊緣提取在數(shù)字圖像處理里是卷積運(yùn)算Sobel 核一滑至少 O(N2) 的乘加換成 Canny還得先高斯導(dǎo)數(shù)、非極大值抑制。而傅里葉光學(xué)給出的做法是把圖像放在 4f 系統(tǒng)的物面第一個(gè)透鏡把圖像變換到頻譜面那里正好是光學(xué)圖像的空間頻率坐標(biāo)。邊緣對(duì)應(yīng)高頻成分頻譜中心是零頻和低頻于是只需要在頻譜面擋掉中心區(qū)域第二個(gè)透鏡做逆變換輸出面上直接就是邊緣增強(qiáng)圖像。這個(gè)過程在光速量級(jí)并行完成沒有逐像素循環(huán)?；诟道锶~光學(xué)中的4f系統(tǒng)實(shí)現(xiàn)光學(xué)圖像的邊緣提取關(guān)鍵是濾波器傳遞函數(shù)與系統(tǒng)參數(shù)匹配適用于光計(jì)算加速、顯微成像、機(jī)器視覺預(yù)檢。從透鏡相位公式到 Python 數(shù)值仿真再到實(shí)驗(yàn)調(diào)參可以一步步落地。2. 4f系統(tǒng)成像鏈路的前提透鏡傅里葉變換、頻譜面坐標(biāo)與數(shù)值模擬2.1 薄透鏡的傅里葉變換性質(zhì)從光程差到焦平面空間頻譜薄透鏡對(duì)光場的相位調(diào)制近似為t(x, y) exp[ -j (π / (λ f)) (x2 y2) ]當(dāng)一個(gè)復(fù)振幅分布 U?(x, y) 放在透鏡前焦面在傍軸條件下透鏡后焦面上的光場 U_f(u, v) 可以寫成對(duì) U? 的二維傅里葉變換U_f(u, v) ∝ ? U?(x, y) exp[ -j 2π (u x v y) / (λ f) ] dx dy把這里的 u / (λ f) 和 v / (λ f) 解釋成空間頻率 ξ、η就得到結(jié)論后焦面各點(diǎn)位置與圖像的空間頻率一一對(duì)應(yīng)光是圖像頻譜的“物理顯示器”。4f 系統(tǒng)就是在這個(gè)事實(shí)之上再接一個(gè)透鏡把頻譜面上的場做逆變換恢復(fù)成一幅圖像。整個(gè)系統(tǒng)長度嚴(yán)格來說是 4 個(gè)焦距物面在透鏡 L1 前焦面L1 后焦面是頻譜面該平面同時(shí)作為透鏡 L2 的前焦面L2 后焦面輸出處理后的圖像。中間這段“頻譜面”就是放濾波器的位置。2.2 4f系統(tǒng)的坐標(biāo)、波長與孔徑約束頻譜面上的物理坐標(biāo)與空間頻率的換算關(guān)系是x_f ξ λ f y_f η λ f如果拍某個(gè)方向上的周期條紋周期 d則空間頻率 ν 1/d。在頻譜面上它對(duì)應(yīng)的半徑是r λ f ν這里非常容易把單位搞錯(cuò)。常見做法是把 ν 用 cycles/mm 表示λ 用 mm 表示f 用 mm 表示算出來的 r 仍然是 mm。例如 λ 532 nm 0.000532 mmf 100 mmν 20 cycles/mm則 r 0.000532 × 100 × 20 1.064 mm。也就是要擋住零頻、保留 20 cycles/mm 以上頻率時(shí)遮光圓斑半徑大約是 1 毫米。透鏡口徑和頻譜面濾波器的物理尺寸都要按這個(gè)式子預(yù)估否則濾出來的邊緣可能只是光斑邊緣而不是圖像邊緣。2.3 用 Python 數(shù)值模擬 4f 系統(tǒng)的邊緣提取最小流程在沒有光路時(shí)先用numpy.fft把 4f 系統(tǒng)的兩次變換完整模擬一遍。下面這段代碼可以直接跑通輸入一個(gè)簡單二值圖形在頻譜面做高通濾波輸出邊緣增強(qiáng)圖。import numpy as np import matplotlib.pyplot as plt N 512 # 采樣點(diǎn)數(shù) pitch 4.8e-6 # 模擬像元尺寸單位 m lam 532e-9 # 波長單位 m f 150e-3 # 透鏡焦距單位 m # 構(gòu)造模擬物面兩個(gè)灰色矩形作為待檢測物體 obj np.zeros((N, N)) obj[96:192, 128:384] 1.0 obj[256:384, 96:256] 0.7 obj[32:64, 32:64] 0.0 # 頻域坐標(biāo)零頻在數(shù)組中心 fx np.fft.fftshift(np.fft.fftfreq(N, dpitch)) fy np.fft.fftshift(np.fft.fftfreq(N, dpitch)) FX, FY np.meshgrid(fx, fy) FU np.sqrt(FX**2 FY**2) # 空間頻率幅度單位 cycles/m # 物體到頻譜面一次傅里葉變換 U_f np.fft.fftshift(np.fft.fft2(obj)) # 高通濾波器擋掉 20 cycles/mm 以下分量 # 20 cycles/mm 20000 cycles/m r0 20e3 H (FU r0).astype(np.float32) # 頻譜面濾波后再還原到像面 U_filt U_f * H out np.abs(np.fft.ifft2(np.fft.ifftshift(U_filt))) # 歸一化輸出方便觀察 out (out - out.min()) / (out.ptp() 1e-12) # 打印有效高頻分量對(duì)應(yīng)的頻譜面半徑 radius_mm lam * f * r0 * 1e3 print(f頻譜面遮擋半徑: {radius_mm:.2f} mm)邏輯說明np.fft.fftfreq(N, dpitch)生成一維空間頻率坐標(biāo)fftshift讓零頻落在數(shù)組中央FU可以理解成頻譜面上各點(diǎn)的空間頻率。U_f是物面的二維傅里葉變換等效于經(jīng)過透鏡 L1 到達(dá)頻譜面。H是一個(gè)以零頻為中心的圓形高通閾值掩模把低頻乘成 0第二個(gè)透鏡的逆變換就把剩下頻譜還原成邊緣圖像。參數(shù)說明r0的單位是 cycles/m寫成20e3是 20 cycles/mm。如果改變波長或焦距要注意同步調(diào)整r0否則可能出現(xiàn)“低通濾波”或“全黑”的結(jié)果。pitch只影響頻域坐標(biāo)刻度卻決定FU與真實(shí)空間頻率的對(duì)應(yīng)是調(diào)參時(shí)最容易錯(cuò)的一條。3. 邊緣提取的空間濾波核高通、拉普拉斯、螺旋相位的傳遞函數(shù)選擇3.1 振幅高通擋零頻的最直接做法最常見的 4f 邊緣提取實(shí)現(xiàn)是在頻譜面中心放一個(gè)不透光圓點(diǎn)。傳遞函數(shù)寫成H(ξ, η) 0, sqrt(ξ2 η2) ν0 H(ξ, η) 1, sqrt(ξ2 η2) ν0這個(gè)濾波器去除直流和低頻輸出的強(qiáng)度圖邊緣非常亮但需要注意邊緣顯示成“雙線”。原因是圖像中臺(tái)階邊緣經(jīng)過高通后兩側(cè)都會(huì)出現(xiàn)沖激輪廓人眼看到的是邊緣兩側(cè)各一圈亮線。這在工業(yè)缺陷檢測里問題不大但如果要做亞像素定位后續(xù)還需要數(shù)字細(xì)化。優(yōu)點(diǎn)是光學(xué)實(shí)現(xiàn)簡單一塊膠合玻璃上鍍不透明圓斑就能用缺點(diǎn)是低頻信息全部丟失均勻區(qū)域輸出接近零邊緣強(qiáng)度也受照明背景影響。3.2 拉普拉斯濾波器頻域傳遞函數(shù)與邊緣位置的對(duì)應(yīng)數(shù)字圖像處理中拉普拉斯算子對(duì)應(yīng)頻域乘 -4π2(ξ2 η2)。在 4f 系統(tǒng)里如果能制作一塊振幅透過率隨半徑平方變化的濾波片那么輸出可以直接等效于光學(xué)拉普拉斯算子H_lap(ξ, η) -4π2(ξ2 η2)實(shí)際工程里很少做這種漸變玻璃因?yàn)楣に囌`差會(huì)直接影響頻響。更常見的做法是在數(shù)字仿真里用它作為“標(biāo)準(zhǔn)答案”用來評(píng)估實(shí)際濾波器的性能。仿真代碼只需把上一節(jié)的H替換成H_lap -4.0 * np.pi**2 * (FX**2 FY**2) # 零頻附近直流噪聲會(huì)被放大先擋掉 H_lap[FU 20e3] 0這個(gè)濾波器放大高頻的同時(shí)也放大探測器噪聲。頻譜面高頻區(qū)如果有塵埃或灰塵點(diǎn)會(huì)在輸出面形成亮斑因此實(shí)驗(yàn)時(shí)后面一定要接低通或中值濾波。3.3 螺旋相位濾波各向同性邊緣增強(qiáng)的光學(xué)實(shí)現(xiàn)螺旋相位濾波的傳遞函數(shù)是H_spiral(ξ, η) exp(jθ), θ arctan2(η, ξ)它不遮擋任何頻率而是給頻譜引入一個(gè)隨方向旋轉(zhuǎn)的相位延遲。物面上的邊緣經(jīng)過逆變換后不同方向的邊緣由于頻譜相位差會(huì)被同等增強(qiáng)并且輸出保留直流分量看到的邊緣是“浮雕”而不是純黑背景。這種濾波器近年經(jīng)常出現(xiàn)在顯微成像和定量相位成像中光學(xué)實(shí)現(xiàn)是螺旋相位板。THETA np.arctan2(FY, FX) H_spiral np.exp(1j * THETA) U_spiral np.fft.fftshift(np.fft.fft2(obj)) * H_spiral out_spiral np.abs(np.fft.ifft2(np.fft.ifftshift(U_spiral)))注意螺旋相位濾波后輸出是復(fù)數(shù)場取模后可以看到邊緣但相位信息里還包含更多邊緣方向線索。如果只取強(qiáng)度會(huì)丟掉方向信息后續(xù)要做高精度定位建議同時(shí)保留實(shí)部和虛部。3.4 三種濾波核的選型表濾波器傳遞函數(shù)輸出邊緣特征光學(xué)實(shí)現(xiàn)難度對(duì)低頻噪聲敏感度振幅高通H0 當(dāng) νν0雙線亮邊背景暗低光闌/鍍膜點(diǎn)低拉普拉斯H∝ν2二階微分邊緣兩側(cè)正負(fù)峰高漸變膜難均勻高螺旋相位Hexp(jθ)浮雕式各向同性邊緣中一體化相位板中選型時(shí)優(yōu)先考慮物面照明穩(wěn)定性。照明均勻性不夠好時(shí)振幅高通的輸出背景波動(dòng)會(huì)遠(yuǎn)大于邊緣信號(hào)這時(shí)螺旋相位濾波的優(yōu)勢更明顯它保留直流分量不會(huì)把背景照明梯度變成假邊緣。4. 4f系統(tǒng)中的3個(gè)必調(diào)參數(shù)截止半徑、焦距像元匹配、相干噪聲抑制4.1 截止半徑的物理換算與微調(diào)依據(jù)實(shí)驗(yàn)臺(tái)上調(diào)整高通濾波圓點(diǎn)首先要知道“擋多大”對(duì)應(yīng)“濾掉多低頻”。換算公式為r0_mm λ_mm × f_mm × ν_cutoff_cycles_per_mm以 532 nm 波長為例不同焦距下截止頻率和圓點(diǎn)半徑的關(guān)系如下截止頻率f100mmf150mmf200mm10 cycles/mm0.532 mm0.798 mm1.064 mm20 cycles/mm1.064 mm1.596 mm2.128 mm50 cycles/mm2.660 mm3.990 mm5.320 mm注意這個(gè)半徑是頻譜面上不透明圓點(diǎn)的半徑。激光光斑如果不到 1 mm直接用直徑 2 mm 的點(diǎn)就能擋掉幾乎所有低頻。調(diào)試時(shí)不要只看邊緣亮度還要觀察均勻區(qū)域的背景波動(dòng)。背景波動(dòng)大通常是圓點(diǎn)太小或照明中心沒有對(duì)準(zhǔn)頻譜面中心。順帶提醒4f 系統(tǒng)頻譜面上的圓點(diǎn)不能理解成“擋住零頻”就完事它對(duì)頻譜的衰減是硬截?cái)唷＝財(cái)鄷?huì)導(dǎo)致輸出圖像產(chǎn)生衍射振鈴邊緣周圍出現(xiàn)明暗交替條紋。如果實(shí)驗(yàn)里看到邊緣附近一圈一圈的紋路優(yōu)先懷疑這個(gè)原因而不是鏡頭污點(diǎn)。4.2 透鏡焦距與像元尺寸匹配讓邊緣方向落在奈奎斯特頻率內(nèi)光學(xué)邊緣提取最終要送到 CCD 或 CMOS 上量化。探測器對(duì)空間頻率有硬上限ν_nyquist 1 / (2 × pixel_size)以 4.8 μm 像元為例奈奎斯特頻率約為 104 cycles/mm。4f 系統(tǒng)輸出面上可用的最大空間頻率由透鏡口徑?jīng)Q定近似為ν_max D / (2 λ f)如果 ν_max 高于探測器奈奎斯特頻率邊緣細(xì)節(jié)會(huì)產(chǎn)生混疊出現(xiàn)摩爾紋或假邊緣。常規(guī)做法是讓 ν_max 略低于 ν_nyquist或者直接在輸出面加一個(gè)光學(xué)低通濾波片。具體到系統(tǒng)設(shè)計(jì)先量出物面需要檢測的最小特征尺寸 d_min用它的倒數(shù)作為目標(biāo)截止頻率然后倒推f r0 / (λ ν_cutoff)如果算出的焦距太長頻譜面濾波器物理尺寸會(huì)很大系統(tǒng)占據(jù)桌面空間迅速增加焦距太短頻譜面上高頻和低頻擠在一起難以用普通光闌分離。我一般會(huì)把 r0 控制在 13 mm 之間以此倒推焦距。4.3 激光散斑與照明孔徑的折中用激光照明時(shí)4f 系統(tǒng)輸出的邊緣圖像上常有一層顆粒狀散斑這種散斑在頻譜上表現(xiàn)為寬帶噪聲和高頻邊緣混在一起。高通濾波后散斑顆粒被放大邊緣圖像看起來像蒙了一層砂紙。解決思路有三個(gè)使用旋轉(zhuǎn)毛玻璃或漫射器降低空間相干性但會(huì)降低系統(tǒng)截止頻率需要重新調(diào)截止半徑。在光源后加一個(gè)小孔光闌把照明光束的數(shù)值孔徑控制在頻譜面濾波器尺寸的 1/3 以內(nèi)。改用可調(diào)諧半導(dǎo)體激光或相干長度較短的激光二極管讓散斑對(duì)比度下降。這三點(diǎn)沒有絕對(duì)優(yōu)先級(jí)。需要做邊緣精確測量的場景優(yōu)先降低空間相干性只做大面積缺陷檢測的保留激光照明并加大濾波器圓點(diǎn)直徑即可抑制部分散斑。模擬調(diào)參時(shí)可以做一個(gè)快速掃描用 2.3 節(jié)的對(duì)象對(duì)多個(gè)r0計(jì)算輸出邊緣強(qiáng)度for nu_cut in [10e3, 20e3, 40e3, 80e3]: H (FU nu_cut).astype(np.float32) out_tmp np.abs(np.fft.ifft2(np.fft.ifftshift( np.fft.fftshift(np.fft.fft2(obj)) * H))) edge out_tmp[obj 0.5].mean() bg out_tmp[obj 0.2].mean() print(nu_cut, edge - bg)這里edge - bg是邊緣區(qū)域與背景區(qū)域的輸出強(qiáng)度差。掃描結(jié)果一般會(huì)在某個(gè)nu_cut處出現(xiàn)峰值繼續(xù)增大截止半徑邊緣內(nèi)部細(xì)節(jié)被濾掉差值反而下降?，F(xiàn)場調(diào)光路時(shí)就按這個(gè)趨勢找靈敏區(qū)。5. 光學(xué)邊緣提取的驗(yàn)證方法階躍響應(yīng)、閾值對(duì)接與工程檢查5.1 用一維階躍目標(biāo)驗(yàn)證邊緣保持率在做整機(jī)標(biāo)定時(shí)先在物面放置一張高反差分標(biāo)定板拍一條與邊緣方向垂直的線。理想情況下4f 系統(tǒng)輸出的是階躍信號(hào)的導(dǎo)數(shù)即一個(gè)尖峰。測量這個(gè)尖峰的半高寬就能評(píng)估系統(tǒng)對(duì)邊緣的保持能力。step np.ones(256) step[64:] 0.0 f_step np.fft.fftshift(np.fft.fft(step)) freq np.fft.fftshift(np.fft.fftfreq(step.size, dpitch)) flt (np.abs(freq) 20e3).astype(float) out_step np.abs(np.fft.ifft(np.fft.ifftshift(f_step * flt))) fwhm np.sum(out_step 0.5 * out_step.max()) * pitch * 1e6 print(f邊緣響應(yīng)FWHM: {fwhm:.2f} um)這段代碼把一維階躍當(dāng)成物面透過率做高通濾波然后通過輸出峰值的半高寬評(píng)估邊緣被“拖寬”了多少。FWHM 越小邊緣提取越銳但在實(shí)際光路中FWHM 不會(huì)小于照明波長和數(shù)值孔徑?jīng)Q定的愛里斑尺寸。5.2 從光學(xué)邊緣到數(shù)字閾值與精定位算法對(duì)接光學(xué)邊緣提取的輸出是強(qiáng)度值不是二值化邊緣。后續(xù)接數(shù)字閾值時(shí)需要先做歸一化否則不同照明亮度下閾值漂移。常見做法是在輸出圖像上取一個(gè)不含目標(biāo)的背景窗口計(jì)算均值 μ 和標(biāo)準(zhǔn)差 σ然后設(shè)閾值threshold μ k σk 通常取 3 到 5。如果 k 太小散斑和灰塵會(huì)被當(dāng)成邊緣k 太大低對(duì)比度邊緣丟失。這個(gè)方法比直接取固定灰度值可靠因?yàn)樗詣?dòng)適應(yīng)照明波動(dòng)。4f 系統(tǒng)輸出的邊緣是多像素寬亞像素定位建議用灰度重心法而非梯度法。因?yàn)楣鈱W(xué)邊緣由于衍射振鈴并非單峰直接取梯度最大值會(huì)鎖定在振鈴環(huán)上重心計(jì)算則更穩(wěn)定。5.3 實(shí)驗(yàn)調(diào)試檢查清單先確認(rèn)物面在透鏡 L1 前焦面并用刀口法找頻譜面位置再放濾波器。遮擋中心時(shí)觀察像面是否出現(xiàn)對(duì)稱的衍射條紋不對(duì)稱說明光軸沒有對(duì)準(zhǔn)圓點(diǎn)中心。用已知分辨率板代替測試圖拍攝邊緣輸出確認(rèn) 20 cycles/mm 附近有穩(wěn)定響應(yīng)。檢查輸出圖像四個(gè)角的亮度是否一致不一致可能是透鏡 L2 孔徑切割導(dǎo)致漸暈。記錄光源波長、焦距、截止半徑、探測器曝光時(shí)間四個(gè)參數(shù)保證同一產(chǎn)線可復(fù)現(xiàn)。最后一點(diǎn)經(jīng)驗(yàn)頻譜面上阻擋低頻的小圓點(diǎn)盡量貼在透明玻璃片上而不是用裸線或大頭針。裸線邊緣衍射會(huì)形成明顯的十字或圓環(huán)衍射條紋在邊緣圖像上表現(xiàn)為固定圖案用圖像處理很難完全去除。本文還有配套的精品資源點(diǎn)擊獲取