TWI776848B - 用於獲得電磁場的波前的層析分佈的方法和光學系統 - Google Patents
用於獲得電磁場的波前的層析分佈的方法和光學系統 Download PDFInfo
- Publication number
- TWI776848B TWI776848B TW107103790A TW107103790A TWI776848B TW I776848 B TWI776848 B TW I776848B TW 107103790 A TW107103790 A TW 107103790A TW 107103790 A TW107103790 A TW 107103790A TW I776848 B TWI776848 B TW I776848B
- Authority
- TW
- Taiwan
- Prior art keywords
- wavefront
- optical
- dimensional
- images
- plane
- Prior art date
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 213
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000009826 distribution Methods 0.000 title claims abstract description 32
- 230000005672 electromagnetic field Effects 0.000 title abstract description 6
- 238000005315 distribution function Methods 0.000 claims abstract description 37
- 210000001747 pupil Anatomy 0.000 claims description 23
- 230000001186 cumulative effect Effects 0.000 claims description 19
- 238000003325 tomography Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 description 16
- 238000006073 displacement reaction Methods 0.000 description 12
- 239000011159 matrix material Substances 0.000 description 7
- 230000004075 alteration Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000001902 propagating effect Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
- G01J2009/002—Wavefront phase distribution
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本發明涉及一種用於光學系統(100)中光的波前(104)的二維重建的方法,包含:測量不同光學平面(101、102)具有光程差。特別地,這種方法適用於用圖像探測器探測電磁場波前的層析分佈。任何標準的二維相機。測量在不同光學平面(101、102)具有一光程差的至少二圖像中光強度的分佈函數。特別地,這種方法適用於用圖像偵測器(例如任何標準的二維相機)探測電磁場波前的層析分佈。
Description
本發明涉及如申請專利範圍第1項的前言部分所述類型的方法、申請專利範圍第10項的前言部分中所述的電腦系統、申請專利範圍第12項的前言部分中所述的存儲媒體、以及如申請專利範圍第13項的前言部分所述的光學系統。
當電磁波通過不均勻介質時,其波前相對於原來的形狀會產生扭曲或變形。波前扭曲會導致光學系統中出現像差,從而降低光學系統的表現。
舉例來說,在地面用望遠鏡進行天文成像時,由於與動蕩的地球大氣的相互作用以及與光學元件(例如望遠鏡的鏡頭或鏡子)的相互作用,來自遠處的天文觀測目標的光波前會變形,從而導致圖像質量下降。
然而,光學系統在其他技術領域的表現,例如攝影、醫學影像(如組織或顯微鏡的斷層掃描),會因為波前與非均勻介質的相互作用而受到沿光路徑發生的波前扭曲的不利影響。
為了校正這種波前變形或像差並重建波前的初始形狀,尤其是使用自適應光學系統,其旨在通過所謂的波前感測器來測量波前變形或像差。
這種波前感測器的典型例子包含Shack-Hartmann感測器、金字塔形感測器(pyramidal sensors)、剪切干涉測量(shearing interferometry)和曲率感測器(curvature sensors)。
已知的波前感測器的缺點在於它們在技術上是複雜的,例如,它們包含多個如鏡頭的光學元件。此外,藉由這些波前感測器產生的數據進行解析以重建初始的波前是需要高密集計算和具挑戰性的。
因此,本發明的目的是提供用於光學系統中的波前重建的改進的裝置。具體而言,例如,本發明的目的是簡化和加速波前扭曲的測量,特別是加速對波前斜率的估計,以及促進和加速初始波前的重建,即重建初始的波前形狀。
根據本發明,此目的可通過根據申請專利範圍第1項的方法、申請專利範圍第10項的電腦系統、申請專利範圍第12項的存儲介質以及申請專利範圍第13項的光學系統來實現。
有利的實施例和進一步的發展是附屬項的標的。
用於光學系統的波前、光波前的二維重建的示例性方法可以包含以下步驟中的一者、部分或全部:
在具有光程差的不同光學平面上的至少二個圖像中,測量光強度的分佈函數(distribution function),例如光強度的二維分佈函數;
該測量的步驟包含在各平面內的不同角度範圍內確定各平面中的光強度的多個一維累積分佈函數(cumulative distribution function);
將所確定的不同光學平面的一維累積分佈函數進行匹配,以便於在一平面(例如位於所述不同的光學平面之間的中間平面)中導出二維波前斜率估計值;以及
將該些二維波前斜率估計積分,以重建位於該些不同光學平面之間的該平面中的波前的二維形狀。
參考下面的結合圖式詳細描述的實施例,本發明的優點和特徵以及獲得它們的方式將變得顯而易見。然而,本發明不限於下文揭示的實施例,並且可以以許多不同的形式實施。另外,提供這些示例性實施例,使得本發明將徹底和完整的揭示,並且將向本領域技術人員充分傳達該範圍。本發明的範圍應由申請專利範圍界定。
示例性的光學系統或光學採集系統由此可以包含數位相機,例如常見的二維相機、望遠鏡、顯微鏡、積分顯示器等成像系統。
尤其是,在光學系統的不同光學平面上拍攝的圖像可以被例如感光耦合元件(charge-coupled device,CCD)擷取。
於此測量圖像中光強度的分佈函數的術語,例如,圖像中光強度的二維分佈函數,根據圖像的光學平面的已知特徵,還可以理解為確定光強度的分佈函數,例如,光強度的二維分佈函數。例如,如果光學系統的瞳孔具有足夠的特徵,例如通過光瞳函數,光瞳平面上光強度的分佈函數可由光學系統的參數直接確定。
此外,由於在這樣的示例性情況下,光瞳平面圖像與通過的波前面無關地保持相同,所以瞳孔圖像(即光瞳平面中的光強度的分佈函數)只需要確定或評估一次即可。
此外,所確定的不同光學平面的一維累積分佈函數的匹配,例如可通過使用沿著直線或沿著軸(如,橫座標,距離函數的值相同)測量確定的一維累積分佈函數之間的距離的度量或距離函數對不同光學平面的一維累積分佈函數的成對匹配來執行。
應注意的是,在此假定波前(即電磁波前沿)或光子的傳播遵循幾何光學的法則,即假設傳播方向垂直於波前。為了完整性,注意到波前可以由具有相同相位的一組波點來定義,即波前或波前形狀可以由相位圖(phase map)來描述,例如二維相圖。
此外,假定傳播波前的光強度的分佈函數可以由概率密度函數(probability density function,PDF)表示光子到達的概率。進一步假定,波前傳播節省了通量(flux),即指一維PDF曲線下的總面積保持不變。
需要進一步指出的是,波前重建的光學平面可以位於沿著光學系統的光路上的任何位置。換句話說,其中波前將被重建的光學平面不需要與任何特定的光學平面重合,例如,光學系統的孔徑平面或光瞳平面。
具體地說,例如可以設想,在具有光程差的不同光學平面上拍攝的圖像,即可測量光強度的分佈函數的圖像可以位於光學系統的孔徑或光瞳平面之前或之後,使得波前將被重建的光學平面也可以位於光學系統的孔徑或光瞳平面之前或之後。
也可能的是,在不同的光學平面上拍攝的具有光程差的圖像,即其中光強度的分佈函數被測量的圖像也可以位於相對於可能的光圈面或光瞳平面的不同距離處。
上面描述的方法提供了更有效和更快的方法以估計二維波前斜率以及從來自波前感測器或例如具有光程差的圖像對的數據導出的所述二維波前斜率估計來恢復波前的二維初始形狀。特別地,所提出的方法可以更有效地恢復波前的初始形狀或波前的初始相位,其中波前變化或光強度變化發生在二維(即沿著不同的方向)中。
另外,根據本發明,例如可以實現在xy平面中拍攝的具有光程差的光路(例如,沿著z軸)上的波前解析度,可以達到幾皮米pm(picometer)。
與已知技術相比,除了其他之外,需要涉及多個光學元件的複雜的波前感測器設計,並且還需要困難和繁瑣的集成技術和波前參數化,例如,使用Zernike多項式(Zernike polynomials)或Radon轉換(Radon transformations),相對於可以計算要恢復的波前的速度,這裡介紹的方法省去了這樣的計算密集型技術和參數化,從而顯著地促進了波前的重建。換句話說,該方法克服了現有技術的問題,可精確地二維波前重建需要計算資源密集和復雜的波前參數化和轉換。
由於其計算效率的提高,該方法適於實際應用,例如在自適應光學系統中控制反饋迴路。
此外,減少計算負擔使得該方法可在普通的中央處理器單元(CPU)或圖形處理器單元(GPU)等可以在智慧型手機、數位相機、平板電腦、筆記型電腦、桌上型電腦和其他普通電子設備中找到的普通電子元件上實現。
這裡描述的方法的另一個有利效果在於,因為不需要額外的光學元件來測量或感測波前,例如在用波前感測器(如Shack-Hartmann感測器或金字塔式感測器)進行波前感測的情況下,不需要額外的空間或時間調製被添加到作為測量的波前。
除其他外,所述確定前述確定各平面內不同角度範圍內各平面內的光強度的多個一維累積分佈函數,可以理解為確定各平面內的光強度在一定範圍的不同方向上的光強度的多個一維累積分佈函數,即在每個平面內的不同直線的範圍或至少一條直線上,其中所述直線各自具有可以由角度表示的斜率。
換句話說,每個平面內的可能角度可以定義每個平面中的直線的斜率,沿著該平面可以確定每個平面中的光強度的一維累積分佈函數。
用於確定每個平面中的光強度的多個一維累積分佈函數,限定每個平面中所述直線的斜率的每個平面內的不同角度的可能範圍可以例如在每個平面中包含至少兩個不同的角度,例如兩個正交角度(如 0°和90°)。
也可能確定每個平面中的光強度的一維累積分佈函數所沿的角度或直線與波前或波前變化的主要方向或軸線相匹配。
但是,可以想到的是,當假定波前變化僅在特定角度或特定方向上發生時,僅用於所述特定角度或特定方向的一維累積分佈函數被確定,以恢復初始波前形狀或波前相位。
而且,上述示例性步驟可以提供執行電磁場的波前的三維分佈(即,體積分佈)的層析成像(tomography)的基礎,電磁場的波前的三維分佈例如波前沿光學系統的光軸的三維分佈,如,在一堆圖像上。換句話說,本文中呈現的方法允許例如允許從一疊或一系列重建的二維波前導出沿著光學系統的光軸的波前的三維分佈。具體而言,本文描述的方法與當前的波前層析成像方法相比,提供了計算更快且更有效的方式來執行所述電磁場波前的三維分佈的層析成像。
在上述示例性方法步驟中,可以例如在光瞳平面(如出瞳平面、或光學系統的孔徑平面)中拍攝至少兩個圖像中的一個。
在光瞳平面中拍攝圖像的可能的示例性步驟也可以被理解為基於光學系統的瞳孔的明確定義的表徵模型來表徵光瞳平面中的圖像。換句話說,至少兩個圖像中的一者可以是可通過理論模型和/或經驗模型充分表徵的圖像。換句話說,這樣的特徵化圖像可以被計算,並且不需要由照相機擷取或測量。
如上所述,這具有以下優點:確定光瞳平面或孔徑平面中的多個一維累積分佈函數可以基於或來源於瞳孔或孔徑的充分定義的模型,即明確定義和/或充分表徵的光瞳函數或孔徑函數。
進一步可能的是,至少一個圖像可能失焦,即在光學系統的散焦平面上。
具體而言,可以想到的是一圖像被拍攝於焦內或焦前和一圖像被拍攝焦外或焦後,即一圖像在光學系統的焦平面之前及一圖像超出光學系統的焦平面。換句話說,至少兩個圖像都可被散焦(defocused)。
所有拍攝的圖像都可能失焦,也就是說所有的圖像都可能會散焦。
但是,也可能會聚焦來自不同光學平面的部分或全部圖像。
此外,注意到,這裡描述的方法也可以應用於無焦光學系統(afocal system),即沒有焦點的光學系統。只需要在測量或確定的圖像之間存在光程差。
可以使用測量或確定的圖像之間的任何光程差。例如,光學系統的分數(fraction)或總長度可以用作光程差。藉此光程差的大小或長度用可以被確定的波前斜率估計或波前形狀或波前相位,可以確定解析度。特別是增加光程差,即增加測量或確定的圖像之間的光學平面之間的距離,可以增加要恢復的波前的波前相位解析度或波前形狀解析度。
電磁場波前的三維分佈(即體積分佈),例如沿光學系統的光軸的波前的三維分佈。通過在具有光程差的不同光學平面上拍攝的圖像,可以提供導出這種相位或波前層析成像的基礎,並且其中例如所述圖像是由具有感光耦合元件(CCD)的圖像檢測器(例如, 一個普通的二維數位相機)拍攝的。因此,所述圖像可以提供用於提供三維波前相位分佈的層析測量的基礎。
換言之,例如可能的是,在不同的光學平面上拍攝多個圖像(例如兩個以上的圖像),以在多個光學平面上重建出多個二維波前形狀。
這有利地允許執行上述的相位或波前層析成像,以在保持圖像的二維解析度的同時,從不同的角度觀看成像的物體。
可替換地或另外地,波前相位層析成像可以基於估計波前相位、或者估計波前斜率,不在所測量或確定的圖像的整個尺寸上,而是從估計波前相位或僅估計部分波前斜率或測量或確定的圖像的部分來執行。
換句話說,可以從多個圖像(例如在不同的光學平面(如在光路位置上)上拍攝兩個以上的圖像)中導出示例性的三維波前相位分佈或層析成像,在探測三維波前相位分佈的多個光學平面上重建多維二維波前形狀,和/或可以從針對給定的光學平面從測量或確定的圖像的不同部分或部分重建的多個二維波前,獲得示例性的三維波前相位分佈或層析成像。
因此,圖像可以被分成多個部分,其中每個部分重建波前的二維形狀。更精確地說,二維波前或波前相位可以在一個平面(例如,位於測量或確定圖像的不同光學平面之間的部分之間的中間平面)的部分重建。
由於從要觀察的物體的不同部分和/或不同深度發出或散射光或光線(或從不同物體發出或散射的光或光線)將以不同的角度或位置射中測量或確定的圖像的圖像平面,即可以在測量或確定的圖像的光學平面中的不同位置處形成圖像,測量或確定的圖像的不同區域可被映射回到待觀察的物體(或不同物體的部分)的不同部分或不同深度。
換言之,通過計算或恢復用於測量或確定的圖像的不同區域或部分的波前或波前相位,可以獲得待觀察物體的不同視角。
換句話說,本文描述的方法可以應用於測量或確定的圖像的多個不同的部分或子部分,從不同部位或不同深度的待觀察物體或來自不同物體的部分恢復初始的波前。
因此,這裡描述的方法可以包含以下步驟中的其中一者、部分或全部:
在具有光程差的不同光學平面上的至少二個圖像中,測量光強度的分佈函數;
其中該測量的步驟可以包含在每個截面內的不同角度範圍內確定光強度的多個一維累積分佈函數;
將所確定的不同光學平面的不同截面的一維累積分佈函數進行匹配,以便於在一平面(例如位於所述不同的光學平面之間的中間平面)的截面上導出二維波前斜率估計值;以及
將該些二維波前斜率估計積分,以重建位於該些不同光學平面之間的該平面的截面中的波前的二維形狀。
對於假定賦予要觀察的物體的波前的多個測量或確定的圖像和可能的多個選擇的相幕(phase screen)的可能的示例性分佈,相當於該區域或其所代表的物體的部分內的相位變化,所述可能的相幕可以從所述多個部分或區域中恢復或確定的波前相位來計算,所述測量或確定的圖像可以被劃分成多個部分或區域。
假設,例如要被觀察的物體或目標介質至少部分透明,則所述物體或目標介質可被建模為沿著光軸(例如,z軸)分佈的一組離散相幕,並且其中所述一組階段螢幕可以基於測量或確定的圖像(可以被分割成多個部分或區域),來計算或恢復一組相幕。
從計算的觀點來看,可以用矩陣模擬相幕,其中不同的矩陣元素代表不同的相變值。如前所述,至少部分透明的物體或目標體積可接著建模為一組或一堆的相幕。
通過將測量或確定的圖像分成多個部分或區域,所述多個部分或區域可以例如被理解為以特定角度捕捉至少部分透明的物體或目標體積的投影(或線積分)。
因此,測量或確定的圖像的特定部分或區域,可以對應於在特定角度相幕的投影(或線積分)。
換句話說,測量或確定的圖像可以分成多個部分或區域,這些部分或區域可以反映特定角度的相幕的假設分佈的投影或線積分。
從這些投影到所測量或確定的圖像的多個部分或區域上,可以定義至少一個方程式系統,由此可以確定相幕的所述矩陣元素的未知值。
例如並且為了簡單起見,相幕由具有代表不同相位變化值的元素x1、x2、x3、x4的以下示例性矩陣表示:
該矩陣的以下兩個投影,例如0°和90°;
對於0度的投影:
對於90度的投影:
因此,從這四個方程可以確定未知值x1、x2、x3、x4。這個例子可以擴展到更多的矩陣大小或尺寸。
換句話說,通過恢復相幕來解決所述方程式系統,可以執行要被觀察的物體或目標體積的相位層析成像。
還值得注意的是,一旦波前(即位於所述不同光學平面之間的波前)被恢復,所述恢復的波前可以根據瑞利-索末菲繞射(Rayleigh-Sommerfeld)的原理傳播。
由此可以模擬未被擷取的光學平面中的光強度,光學系統只需拍攝兩個圖像即可作為光場相機和相位層析感測器。
如上所述,該方法可以用電腦執行,即電腦系統可以包含至少一個處理器單元,該處理器單元被配置成執行如本文所述的用於波前的二維重建的方法。
所述至少一個處理器單元可以是例如常見的中央處理器單元(CPU)或圖形處理器單元(GPU)。
由於該方法涉及圖像(例如可以擷取、存儲和處理作為基於像素的圖像),圖形處理器單元(GPU)可以特別適合於執行該方法。然而,最後但並非最不重要的是,由於這裡描述的方法比當前已知的技術計算更有效率,所以它也可以容易地在普通的中央處理器單元(CPU)上執行。
例如這裡介紹的方法的計算時間可以是ms級或更小。具體而言,例如,對於256×256像素的示例性圖像尺寸,已經獲得具有1.14毫秒的GPU的計算時間,以及對於140×140的示例性圖像尺寸,具有0.73毫秒的GPU的計算時間。
這裡描述的方法可以在電腦可讀指令中執行,可讀指令可以存儲在一個或多個電腦可讀存儲介質上或者在另一個電腦可讀數據文件中。
圖像構成了本文所述的光學平面中波前的二維重建的基礎,並且圖像可選地還可以提供波前的三維分佈(即體積分佈)的層析成像,可以例如由圖像檢測器(如包含光學系統的感光耦合元件(CCD)的照相機和/或包含照相機的專用波前感測器)擷取,照相機例如一個普通的二維數位相機。
換句話說,示例性的光學系統可以包含至少一個圖像檢測器和/或至少一個波前感測器,並且光學系統可以被配置為在波前重建和/或層析成像中,基於來自至少一個圖像檢測器的數據和/或基於來自至少一個波前感測器的數據,執行如本文所述的波前的三維分佈(即體積分佈)。
特別地,可以設想的是,可能的波前感測器可以是曲率感測器,和/或其中包含光學採集系統(例如包含感光耦合元件(CCD)的照相機)的波前感測器,感光耦合元件(CCD)的照相機例如一個普通的二維數位相機。
以下,使用示例性框架給出光學系統的光學平面中拍攝的圖像的光強度或光強度分佈的分佈函數(例如二維分佈函數)的測量或估計的示例。
其中,CDF表示f
沿著斜率p
和截距t
的直線的累積分佈函數。
該轉換V
可以表示為:
另外,對於相同的x值,我們可以用表示距離度量或者距離函數D來表示兩條曲線a和b在x軸上的距離。這裡,兩條曲線a和b可以被理解為沿著如上定義的用於V轉換的圖像的線(即具有斜率或角度和截距的線)的累積分佈函數CDF。
函數可用於確定在不同光學平面上的兩個圖像之間的測量的光子到達位置、或光子射線位置、或測量的局部光線強度位置的空間位移,其中空間位移由波前像差引起並取決於波前像差。換句話說,通過將具有光程差的光學系統的兩個不同圖像中確定的光子射線位置、或測量的局部光強度進行比較,可以限制波前像差的出現,並且允許重建初始的波前形狀。
此外,我們假設在示例性光學系統中,我們已經在具有光程差的不同光學平面上擷取了兩個圖像I1
、I2
,即兩個不同的光強度分佈。
因此,位於圖像I1
、I2
的不同光學平面之間的波前重建平面中的像素(x, y)的水平一階導數可以由下式給出:
且位於圖像I1
、I2
的不同光學平面之間的波前重建平面中的像素(x, y)的垂直一階導數可以由下式給出:
利用波前的斜率與光子的位移之間的幾何光學的線性關係,並利用所確定的一階導數與局部波前斜率(即在波前重建平面中的像素(x,y)處)的等價性,位於圖像I1
、I2
的不同光學平面之間的平面中重建的二維波前,可以通過沿著x和y對進行積分來獲得。
所述導出的波前斜率的積分通過估計來重建位於所述不同光學平面之間的光學平面中的波前的二維形狀可以使用,例如快速傅立葉轉換(fast Fourier transform,FFT)演算法來執行,其例如包含以下步驟中的其中一者、部分或全部:
計算傅立葉逆轉換以獲得重建波前的估計值
但是,也可以應用其他集成方案。
下圖舉例說明:
圖1示出簡化的光學系統100的例子,其示出了簡化的一維情況下示例性光子106a、106b、107a、107b、108a、108b的位移與波前104的斜率(即局部斜率)之間的關係,其中所有的波前和強度變化僅在一個方向上發生,並且對於具有三個探測到的光子的示例性情況,例如,在具有光程差(如沿z軸有一個光程差)的每個光學平面101、102中擷取圖像。
示例性光子106a、106b、107a、107b、108a、108b在它們各自的光學平面101、102中的分佈由此可以被解釋為表示光強度分佈,並且光學平面101、102可以被解釋為具有光程差的不同圖像。
附圖標記109表示光學系統的示例性方向,其中z軸與光學系統100的光軸(未示出)相同或平行。
假定光子106a、106b、107a、107b、108a、108b以圖像平面或光學平面102和101之間的直線105a、105b、105c行進,並且假設光子的傳播方向垂直於它們相應的局部波前,則光子沿x軸的位移由(局部)波前斜率乘以兩個光學平面101、102之間的距離給出。
因此波前104的局部波前斜率104a、104b、104c,通過分別將一個光學平面101的光子106b、107b、108b匹配到另一個光學平面102的光子106a、107a、108a上,可在光子位置之間的光學平面103處或在光學平面102、101之間的中途被估計或重建。
其中測量光子分佈(即光強度分佈)的光學平面102、101可以位於沿光學系統100的光路的任何地方。因此,波前將被重建的光學平面103也可以位於沿著光學系統100的光路的任何地方。換句話說,將被重建的波前的光學平面103不需要與光學系統100的任何特定的光學平面(例如孔徑平面或光瞳平面)重合。
如前所述,可以想到,在具有光程差的不同光學平面上拍攝的圖像,即光強度的分佈函數被測量的圖像,可以位於光學系統100的孔徑或光瞳平面之前或之後,使得波前重建的光學平面103也可位於光學系統100的孔徑或光瞳平面之前或之後。
因此,有可能在不同的光學平面上拍攝具有光程差的圖像,即其中光強度的分佈函數被測量的圖像可以位於相對於可能的光圈面或光瞳平面的不同距離處。
上面描述的根據本發明的方法,現在也允許在更複雜的情況下恢復波前的形狀,其中波前和強度變化發生在二維和不同方向上。
圖2示例性示出了測量射線位置偏移202的誤差203的誤差曲線201的曲線圖200,即測量光線位置或空間位移的空間位移或偏移的誤差203、或測量局部光強度位置的空間位移或偏移202。
如上所述,可以假定波前的斜率與垂直於波前傳播的光子或光子射線的空間位移之間的關係可以遵循線性關係。
光線位置或局部光強位置的位移可以測量越好,初始的波前形狀或波前相位可以恢復越好。
可以實現的波前解析度的估計可以由下式給出:
其中d
表示在其上執行在此呈現的方法的兩個圖像之間的距離,即光程差,以m為單位,p
是物體空間中的像素尺寸。
從上面的示例性估計可以進一步看出,可實現的波前解析度可以隨著光程差的增加而增加,這是因為較長的距離可以放大給定的光子射線角度(例如光子射線相對於光軸的傳播角度)的偏移或位移。
在所示出的示例性情況下,小於0.5像素的位移或偏移的誤差相當小,例如小於10%,這意味著,例如幾公分(cm)的光程差,可以獲得小至皮米(picometer)級的波前解析度。
為了完整性,我們注意到,可以通過atan (p
/d
)來估計最小可測量角度或波前斜率或波前相位斜率,其中d表示在其上執行在此呈現的方法的兩個圖像之間的距離,即光程差,p
是圖像平面或圖像感測器的像素大小。
圖3示出了具有光軸315的簡化的光學系統300的示例,示出了基於在具有光程差的不同光學平面或測量平面316、317處,對測量或確定的圖像301、302進行分割或劃分的用於波前相位層析成像的光學系統300的可能的示例性配置。
此外,示例性的光學系統300包含選用的光學元件304,例如能夠將從待觀察的示例性物體307、308傳播的光線或光線束309、310聚焦到示例性焦平面303上的透鏡。所述示例性物體307、308可以是位於距光軸315的不同距離處的兩個不同的物體(如圖所示),或者它們可以是單個物體的不同部分。
如示例性所示,光線或光線束309、310可以在每個平面309、310上的不同分離位置處,即在測量或確定的圖像301、303的不同位置處,射中示例性光學平面316、317。
如上所述,測量或確定的圖像301、303可被分割或劃分或分隔成不同的部分或區域,其中區域或部分可以重疊或可以分開。例如圖像301可以被分割成兩個區域311、314,並且圖像303可以被分割成兩個區域312、313。
其他分割方案也是可以想像的。這裡顯示的簡單分割只是說明性的。如示例所示,物體307的光線在光學平面316的區域311處射中圖像301,並在光學平面317的區域312處射中圖像302,而物體308的光線在光學平面316的區域314處射中圖像302,並在光學平面317的區域313處射中圖像302。
在圖像301、302的整個尺寸或圖像感測器(例如感光耦合元件(CCD))的整個尺寸上,無法應用這裡描述的用於波前或波前相位二維重建的方法,重建波前或波前相位的方法只能應用於其中每個圖像301、302或測量平面316、317或圖像感測器(未示出)被分成部分或區域上。
換句話說,波前相位不是在整個圖像301、302上被恢復,而是每個圖像301、302的每個部分或區域311、314、312、313的波前相位被恢復。
更確切地說,位於相應圖像301、302(即光學平面316、317之間)的部分或區域311、314、312、313之間的平面中的部分的波前形狀或波前相位可以被恢復。
假設,例如至少部分透明的物體或目標介質體積318,所述至少部分透明的物體或目標介質體積318可以被建模為一組不同的離散相幕305、306,其中如上所述的相幕可以通過矩陣建模,其中不同的矩陣元素表示不同相位變化值,用於藉由不同地區的物體或目標介質體積,通過所述至少部分透明的物體或目標介質體積318的波前傳播的相位變化。
通過將測量或確定的圖像301、302分成多個部分或區域311、314、312、313,所述多個部分或區域可以例如被理解為擷取至少部分透明的物體或目標介質體積318的在特定角度的投影(或線積分)。
因此,測量或確定的圖像301、302的特定部分或區域311、314、312、313可對應於特定角度的相幕305、306的投影(或線積分)。
然後,如上所述,將測量或確定的圖像301、302劃分成多個部分或區域311、314、312、313可形成定義方程式系統的基礎,從中可以計算或恢復多個相幕305、306。
計算或恢復的相幕305、306的集合,然後可以尤其允許執行例如至少部分透明的物體或要觀察的目標介質體積318的波前相位層析成像,其用於在不同視角和/或不同深度下的物體的不同部分或不同的物體。
以下是包含圖1、圖2和圖3的三張表格,並且其中附圖標記標識以下部件:
100:示例性光學系統
101:在具有(第一)光強度分佈的(第一)光路位置處的示例性(第一)像平面或(第一)光學平面
102:在具有(第二)光強度分佈的(第二)光路位置處的示例性(第二)圖像平面或(第二)光學平面
103:在所述第一和第二光學平面之間的示例性光學平面,其中波前將被重建,例如光學系統的光圈平面
104:示例性的波前被重建
104a:示例性局部波前段具有(第一)局部斜率
104b:示例性局部波前段具有(第二)局部斜率
104c:示例性局部波前段具有(第三)局部斜率
105a:示例性光子傳播軌跡/光子傳播方向/局部波前傳播指南
105b:示例性光子傳播軌跡/光子傳播方向/局部波前傳播指南
105c:示例性光子傳播軌跡/光子傳播方向/局部波前傳播指南
106a:示例性光子表示光學平面102中的局部光強度
106b:示例性光子表示光學平面101中的局部光強度
107a:示例性光子表示光學平面102中的局部光強度
107b:示例性光子表示光學平面101中的局部光強度
108a:示例性光子表示光學平面102中的局部光強度
108b:示例性 光子表示光學平面101中的局部光強度
109:示例性光學系統的方向,z軸是相同的或平行於光學系統的光軸(未示出)
200:示例性 射線偏移測量的誤差圖
201:示例性誤差曲線
202:示例性橫座標軸,例如 y軸,例如 像素中的光線偏移
203:示例性座標軸,例如 x軸,例如 標準化尺度從0到1的射線偏移測量的誤差
300:示例性 光學系統
301:示例性(第一)光路位置處的(第一)圖像
302:示例性 (第二)光路位置處的(第二)圖像
303:示例性可能的焦平面
304:示例性光學系統的光學元件,例如透鏡
305:示例性(第一)相幕
306:示例性 (第二) 相幕
307:示例性 (第一) 要觀察的物體
308:示例性 (第二) 要觀察的物體
309:示例性(第一)物體發出的光線(光線束)
310:示例性 (第二)物體發出的光線(光線束)
311:示例性(第一)圖像301的(第一)部分或區域
312:示例性 (第二) 圖像 302的(第一)部分或區域
313:示例性 (第二) 圖像 302的(第二) 部分或區域
314:示例性 (第一)圖像301的(第二) 部分或區域
315:示例性光軸
316:示例性(第一)光學平面或(第一)測量平面
317:示例性(第二)光學平面或(第二)測量平面
318:示例性至少部分透明的目標(介質)量
100‧‧‧光學系統101、102、103‧‧‧光學平面104‧‧‧波前104a、104b、104c‧‧‧波前斜率105a、105b、105c‧‧‧直線106a、106b、107a、107b、108a、108b‧‧‧光子109‧‧‧方向200‧‧‧曲線圖201‧‧‧誤差曲線202‧‧‧偏移203‧‧‧誤差300‧‧‧光學系統301、302、303‧‧‧圖像303‧‧‧焦平面304‧‧‧光學元件305、306‧‧‧相幕307‧‧‧物體308‧‧‧物體309、310‧‧‧光線或光線束311、314、312、313‧‧‧部分或區域315‧‧‧光軸316、317‧‧‧光學平面或測量平面318‧‧‧目標介質體積
圖1是示例性的光學系統。 圖2是射線偏移測量精度的示例性圖。 圖3是示例性的波前相位層析成像配置圖。
100‧‧‧光學系統
101、102、103‧‧‧光學平面
104‧‧‧波前
104a、104b、104c‧‧‧波前斜率
105a、105b、105c‧‧‧直線
106a、106b、107a、107b、108a、108b‧‧‧光子
109‧‧‧方向
Claims (13)
- 一種用於一光學系統中的光波前的二維重建方法,包含:在具有一光程差的多個不同光學平面上的至少二圖像中測量一光強度的一分佈函數,包含:確定各該不同光學平面內該光強度的多個一維累積分佈函數;匹配確定的該些不同光學平面的該些一維累積分佈函數以得到一光學平面中的多個二維波前斜率估計,該光學平面位於所述不同光學平面之間的一中間平面,其中根據位於不同光學平面內的該些一維累積分佈函數其中二者在正交於該光學系統的光軸之方向上的距離得到該些二維波前斜率估計其中一者;以及將該些二維波前斜率估計積分,以重建位於該些不同光學平面之間的該光學平面中的一波前的一二維形狀;其中,該些不同光學平面沿該光學系統的光軸排列,且該至少二圖像之間的距離定義出該光程差。
- 如請求項1所述之方法,其中該至少二圖像中的其中一者是通過理論模型和/或經驗模型充分表徵的被計算圖像。
- 如請求項1所述之方法,其中該至少二圖像中的其中一者被拍攝在該光學系統的一光瞳平面中。
- 如請求項1所述之方法,其中一該圖像位於焦點內和另一該圖像位於焦點外。
- 如請求項1所述之方法,其中在該些不同的光學平面上拍攝多個圖像以在該些不同的光學平面上重建多個二維波前形狀。
- 如請求項1其中任一項所述之方法,其中該些圖像被分成多個部分,且其中各該部分重建該波前的該二維形狀。
- 如請求項1其中任一項所述之方法,其中恢復的該波前根據瑞利-索末菲(Rayleigh-Sommerfeld)繞射進行傳播。
- 一種包含至少一處理器單元的電腦系統,該處理器單元用以執行根據如請求項1-7中任一項所述的用於光波前的二維重建方法。
- 如請求項8所述之電腦系統,其中該至少一該處理器單元是圖形處理器單元(graphical processing unit,GPU)。
- 如請求項8所述之電腦系統,包含一或多個其中存儲有多個指令的電腦可讀存儲介質,該些指令在由一個或多個處理器執行時指導該至少一處理器執行根據請求項1所述的方法。
- 一種光學系統,包含至少一圖像檢測器,其中該光學系統用以基於來自至少一該圖像檢測器且包含該些不同光學平面上拍攝之圖像的一數據執行根據前述請求項1中所述的用於光波前的二維重建方法和/或用以 執行基於來自至少一圖像檢測器的該數據的多個波前的一三維分佈的一層析成像。
- 如請求項11所述之光學系統,更包含至少一波前感測器,其中該光學系統更用以基於來自至少一該波前感測器且包含波前曲率的一數據在一波前重建處執行根據前述請求項1的方法和/或用以基於來自至少一該波前感測器的該數據,執行多個波前的一三維分佈的一層析成像。
- 如請求項12所述之光學系統,其中該波前感測器是一曲率感測器和/或其中該波前感測器包含一光學採集系統。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP17382047.3 | 2017-02-03 | ||
??17382047.3 | 2017-02-03 | ||
EP17382047.3A EP3358321B1 (en) | 2017-02-03 | 2017-02-03 | Method and optical system for acquiring the tomographical distribution of wave fronts of electromagnetic fields |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201830338A TW201830338A (zh) | 2018-08-16 |
TWI776848B true TWI776848B (zh) | 2022-09-11 |
Family
ID=58264456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW107103790A TWI776848B (zh) | 2017-02-03 | 2018-02-02 | 用於獲得電磁場的波前的層析分佈的方法和光學系統 |
Country Status (10)
Country | Link |
---|---|
US (1) | US11015981B2 (zh) |
EP (1) | EP3358321B1 (zh) |
JP (1) | JP6968895B2 (zh) |
KR (1) | KR102343320B1 (zh) |
CN (1) | CN110476043B (zh) |
CA (1) | CA3051969C (zh) |
IL (1) | IL268304B2 (zh) |
RU (1) | RU2733822C1 (zh) |
TW (1) | TWI776848B (zh) |
WO (1) | WO2018141853A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3736796A1 (en) * | 2019-05-07 | 2020-11-11 | Wooptix S.L. | Method and optical system for characterizing displays |
CN113203485B (zh) * | 2021-04-27 | 2022-08-05 | 浙江大学 | 一种通过单次曝光实现轴向相位差波前重建的装置及方法 |
EP4134644A1 (en) | 2021-08-11 | 2023-02-15 | Wooptix S.L. | System and method for extracting information on the spatial distribution of wavefronts |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5454004A (en) * | 1994-05-06 | 1995-09-26 | Regents Of The University Of Minnesota | Phase grating and mode-selecting mirror for a laser |
GB9820664D0 (en) * | 1998-09-23 | 1998-11-18 | Isis Innovation | Wavefront sensing device |
US7609388B2 (en) | 2002-01-24 | 2009-10-27 | Icos Vision Systems Nv | Spatial wavefront analysis and 3D measurement |
JP2007533977A (ja) * | 2004-03-11 | 2007-11-22 | アイコス・ビジョン・システムズ・ナムローゼ・フェンノートシャップ | 波面操作および改良3d測定方法および装置 |
ES2372515B2 (es) * | 2008-01-15 | 2012-10-16 | Universidad De La Laguna | Cámara para la adquisición en tiempo real de la información visual de escenas tridimensionales. |
GB0813907D0 (en) * | 2008-07-30 | 2008-09-03 | Univ Durham | Sub-micron 3d holographic lithpgraphy |
EP2553531B1 (en) * | 2010-04-01 | 2019-10-16 | SeeReal Technologies S.A. | Method and device for encoding three-dimensional scenes which include transparent objects in a holographic system |
US8907260B2 (en) * | 2011-01-14 | 2014-12-09 | The United States Of America, As Represented By The Secretary Of The Navy | Extended source wavefront sensor through optical correlation with a change in centroid position of light corresponding to a magnitude of tip/tilt aberration of optical jitter |
US9247874B2 (en) * | 2013-02-01 | 2016-02-02 | Carl Zeiss Meditec, Inc. | Systems and methods for sub-aperture based aberration measurement and correction in interferometric imaging |
CN103226059A (zh) * | 2013-04-09 | 2013-07-31 | 中国科学院西安光学精密机械研究所 | 光学***波前测量装置及测量方法 |
CN105283118B (zh) * | 2013-06-10 | 2018-01-02 | 埃西勒国际通用光学公司 | 用于确定待测试光学***的波前像差数据的方法 |
-
2017
- 2017-02-03 EP EP17382047.3A patent/EP3358321B1/en active Active
-
2018
- 2018-02-01 US US16/483,385 patent/US11015981B2/en active Active
- 2018-02-01 RU RU2019126068A patent/RU2733822C1/ru active
- 2018-02-01 JP JP2019542229A patent/JP6968895B2/ja active Active
- 2018-02-01 WO PCT/EP2018/052533 patent/WO2018141853A1/en active Application Filing
- 2018-02-01 CN CN201880023516.4A patent/CN110476043B/zh active Active
- 2018-02-01 KR KR1020197025611A patent/KR102343320B1/ko active IP Right Grant
- 2018-02-01 CA CA3051969A patent/CA3051969C/en active Active
- 2018-02-02 TW TW107103790A patent/TWI776848B/zh active
-
2019
- 2019-07-28 IL IL268304A patent/IL268304B2/en unknown
Non-Patent Citations (2)
Title |
---|
期刊 Marcos A. Van Dam et al., "Wave-front sensing from defocused images by use of wave-front slopes", APPLIED OPTICS VOL.41 NO.26 September 2002 pages 5497~5502; * |
期刊 Saloni Pal et al., "Practical application of the geometric wavefrontsensor for adaptive optics" 2016 International Conference on Image and Vision Computing New Zealand (IVCNZ) 1 November 2016 pages 1-5 * |
Also Published As
Publication number | Publication date |
---|---|
IL268304B1 (en) | 2023-01-01 |
EP3358321B1 (en) | 2024-07-03 |
CA3051969A1 (en) | 2018-08-09 |
KR20190138632A (ko) | 2019-12-13 |
IL268304B2 (en) | 2023-05-01 |
CN110476043B (zh) | 2021-09-07 |
JP6968895B2 (ja) | 2021-11-17 |
CN110476043A (zh) | 2019-11-19 |
KR102343320B1 (ko) | 2021-12-23 |
US11015981B2 (en) | 2021-05-25 |
IL268304A (en) | 2019-09-26 |
RU2733822C1 (ru) | 2020-10-07 |
WO2018141853A1 (en) | 2018-08-09 |
US20200011738A1 (en) | 2020-01-09 |
TW201830338A (zh) | 2018-08-16 |
CA3051969C (en) | 2021-11-16 |
EP3358321A1 (en) | 2018-08-08 |
JP2020506390A (ja) | 2020-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8305485B2 (en) | Digital camera with coded aperture rangefinder | |
US8432479B2 (en) | Range measurement using a zoom camera | |
DK2239706T3 (en) | A method for real-time camera and obtaining visual information of three-dimensional scenes | |
US8773550B2 (en) | Range measurement using multiple coded apertures | |
EP2597421B1 (en) | Stereo distance measurement apparatus and stereo distance measurement method | |
US20110267485A1 (en) | Range measurement using a coded aperture | |
US8330852B2 (en) | Range measurement using symmetric coded apertures | |
TWI776848B (zh) | 用於獲得電磁場的波前的層析分佈的方法和光學系統 | |
TWI687661B (zh) | 用於確定場景相關電磁場的複振幅的方法及裝置 | |
JP2018081378A (ja) | 画像処理装置、撮像装置、画像処理方法および画像処理プログラム | |
TW202307526A (zh) | 光學系統、用於擷取波前的空間分佈上的資訊的方法及系統以及電腦可讀取儲存媒介 | |
Raj | Accurate depth from defocus estimation with video-rate implementation | |
Zhang et al. | Exemplar-Based Elemental Images Restoration in Integral Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
GD4A | Issue of patent certificate for granted invention patent |