Ansys Zemax | 在 MATLAB 或 Python 中使用 ZOS-API 進(jìn)行光線追跡的批次處理
附件下載
聯(lián)系工作人員獲取附件
這篇文章會(huì)說明如何在 MATLAB 或 Python 中以 Zemax OpticStudio 應(yīng)用程式界面 (ZOS-API)處理光線數(shù)據(jù)庫(kù)(Ray Database, ZRD)檔案,過程中我們將使用ZRDLoader.dll。本文提供了在 Matlab 中批次處理序列光線追跡(一般、歸一化、偏振或非偏振),以及在 Matlab 和 Python 中使用方法 ReadNextSegmentFull() 批次處理非序列 ZRD 檔案的范例。
什么是 RayTrace.dll?
RayTrace.dll 結(jié)合了 ZRDLoader.dll 和 ReadNextSegmentFull,并在序列模式中與4個(gè) IBatchRayTrace 界面共同運(yùn)作。此外,RayTrace.dll 這個(gè)動(dòng)態(tài)連結(jié)函式庫(kù)涵蓋了以下的界面,并且繼承了 ZOSAPI.Tools.RayTrace 的命名空間:
1.IRayTraceDirectPolData (sequential):
在這個(gè)界面中,我們可以直接在 XYZ 坐標(biāo)中完成批次的偏振光線追跡,而這也和 DDE 光線追跡指令(模式3)十分相似。此外,該界面屬于 IBatchRayTrace 界面的一種。
2.IRayTraceDirectUnpolData (sequential):
在這個(gè)界面中,我們可以在 XYZ 坐標(biāo)中完成批次的非偏振光線追跡,而這與 DDE 光線追跡指令(模式1)相似。此外,該界面也同樣屬于 IBatchRayTrace 界面的一種。
3.IRayTraceNormPolData (sequential):
在這個(gè)界面中,我們可以使用歸一化光瞳坐標(biāo)(normalized pupil coordinate)進(jìn)行批次偏振光線追跡,這個(gè) DDE 光線追跡指令(模式2)相似。此外,這界面也同樣屬于 IBatchRayTrace 界面的一種。
4.IRayTraceNormUnpolData (sequential):
在這個(gè)界面中,我們可以使用歸一化光瞳坐標(biāo)進(jìn)行批次非偏振光線追跡,這與 DDE 光線追跡指令(模式0)相似。
5.IZRDReader (non-sequential):
我們使用這個(gè)界面在Matlab中設(shè)定RayDatabaseReader并使其運(yùn)作,同時(shí)讀取 C#dll 中每一個(gè)區(qū)段(segment)的值。
這個(gè) DLL 的設(shè)計(jì)是為了使光線或區(qū)段的原始值(raw value)可以作為宣告物件的屬性,而被傳遞到 MATLAB 或 Python,因此 DLL 運(yùn)作的過程中不會(huì)涉及任何的計(jì)算、合并和分析。這個(gè)DLL的主要目的是為了減少M(fèi)ATLAB或Python中十分耗時(shí)的for循環(huán),將原先效率極低的程式轉(zhuǎn)為原生且編譯完成的C#代碼。接著,當(dāng)原始值經(jīng)過所有循環(huán)后,C#代碼會(huì)將數(shù)值回傳到MATLAB或Python進(jìn)行最后的處理。
在序列模式的界面中有兩種情況會(huì)需要for循環(huán)的協(xié)助:AddRays()和ReadNextResults(),分別用來增加光線和讀取數(shù)據(jù)。為了提高效率,我們會(huì)以DLL達(dá)成這兩項(xiàng)任務(wù)。文章附件的其中兩個(gè)示例示范了如何在序列模式中高效率的產(chǎn)生方形網(wǎng)格或低量化誤差的光線(dithered rays)。根據(jù)預(yù)期達(dá)成的光線模式(若與預(yù)設(shè)模式不同的話),我們可以在MATLAB中使用原生數(shù)組索引(native array indexing)的方式增加DLL中的光線,這與原本的方法是十分相似的。無論是在ReadNextResult()或ReadNextResultFull(),最終回傳到MATLAB或Python中的結(jié)果均會(huì)與該界面最復(fù)雜輸出(verbose output)有著相同的物件屬性,包含了正確的相位和相對(duì)應(yīng)的案例。為了提高效率,我們會(huì)試著以 MATLAB 的陣列索引(array indexing)取代 for 或 while 循環(huán)的使用。
DLL 使用范例
我們使用三個(gè)基礎(chǔ)和三個(gè)進(jìn)階范例示范如何在一般序列模式(sequential direct)、歸一化序列模式(sequential normalized)和非序列 ZRD 光線追跡的情況下使用RayTrace.dll。這些示范附件都可以直接下載并使用,不需要對(duì)程式碼進(jìn)行任何更改。需要注意的是,RayTrace.dll 必須被存放在和其他 MATLAB 或Python 腳本(script)相同的資料夾中,才能確保使用時(shí)能被順利的呼叫。產(chǎn)生每個(gè) DLL 所使用的 C# 原始碼以及 Visual Studio 的解都可以在下方的文章中找到。
基本范例
· MATLAB_BatchRayTrace_Direct.m – (SEQ Direct).
這個(gè)范例使用了 IRayTraceDirectUnpolData 界面。此外,它還同時(shí)使用了RayTrace.dll (DirectUnpol.cs) 中的 ReadDirectUnpolData 函數(shù)。范例中,系統(tǒng)追跡了特定透鏡表面的光跡圖(footprint diagram),同時(shí)繪出系統(tǒng)的凈口徑 (Clear Semi-Diameter)和機(jī)械半直徑 (Mechanical Semi-Diameter)。

· MATLAB_BatchRayTrace_Normalized.m – (SEQ Norm).
這個(gè)范例使用了 IRayTraceNormUnpolData 界面,而該界面還包含了RayTrace.dll (NormUnpol.cs) 中的 ReadNormUnpolData 函數(shù)。范例中,系統(tǒng)以超過10,000條光線追跡一個(gè)方形圖案投射至像面的結(jié)果。設(shè)定上只使用了單一視場(chǎng),且整個(gè)過程耗時(shí)不到一秒。

MATLAB_ZRDLoaderFull.m or PythonNET_ZRDLoaderFull.py – (NSC ZRD)
這個(gè)范例使用了 IZRDReader 界面,且同時(shí)使用了RayTrace.dll (ZRDLoaderFull.cs)中的 ReadZRDData 函數(shù)。范例中,系統(tǒng)解析了 ZRD 檔案并決定了特定探測(cè)器上接收到的總能量。

進(jìn)階范例
· MATLAB_BatchRayTrace_Surface_AOI.m – (SEQ Direct).
這個(gè)范例使用了 IRayTraceDirectUnpolData 界面,且同時(shí)使用了RayTrace.dll (DirectUnpol.cs)的ReadDirectUnpolData 函數(shù)。范例中,系統(tǒng)追跡了光線網(wǎng)格投射至一個(gè)特定表面的結(jié)果。過程中使用了 LMN 方向的 cosine 向量和 Nxyz 法向量計(jì)算光線在表面上每個(gè) XY 點(diǎn)的入射角。

· MATLAB_BatchRayTrace_Ex22_Performance_Comparison.m – (SEQ Norm).
這個(gè)范例使用了 IRayTraceNormUnpolData 界面,同時(shí)還使用了RayTrace.dll (NormUnpol.cs)的ReadNormUnpolData 函數(shù)。范例中,系統(tǒng)使用低量化誤差光線圖型(dithered ray pattern) 繪制了整個(gè)視場(chǎng)的點(diǎn)列圖(spot diagram)。模擬過程中,系統(tǒng)可以在一秒內(nèi)追跡10,000條光線。相同情況下,若使用原生 MATLAB 的 for 循環(huán),繪制961條光線將耗時(shí)超過4秒。根據(jù)電腦的性能差異,使用 DLL 最高可將運(yùn)算過程加速為原先的40倍。

MATLAB_ZRD_Pixelated_Detector_xybin.m – (NSC ZRD).
這個(gè)范例使用了IZRDReader界面,同時(shí)也使用了RayTrace.dll (ZRDLoaderFull.cs)的 ReadZRDData 函數(shù)。范例中的系統(tǒng)有矩形光源、矩形探測(cè)器,以及一個(gè)由 “LETTERF.BMP”產(chǎn)生的幻燈片(Slide) 物件,最后在探測(cè)器上顯現(xiàn)出圖型 “F”。接著,我們使用 ZRD 檔案中的“xybin”將每個(gè)像素的光通量加總,重新在 ZRD 檔案中產(chǎn)生一個(gè)新的探測(cè)器結(jié)果。這樣的做法可使我們免于使用探測(cè)器查看器(Detector Viewer)的分析功能,直接在 MATLAB 中看到結(jié)果。


疑難排解
若是使用了范例提供的 Visual Studio 解,但卻在編譯的過程中出現(xiàn)錯(cuò)誤的話,請(qǐng)先照著以下步驟操作以確保 ZOSAPI 和 ZOSAPI_Interfaces 的來源(Reference)是正確的。我們可以打開 Solution Explorer,并點(diǎn)選 References。假如這時(shí) ZOSAPI 和 ZOSAPI_Interfaces 旁邊出現(xiàn)黃色的警示標(biāo)語,你將會(huì)需要進(jìn)行以下操作:
1.移除現(xiàn)有的 references (鼠標(biāo)右鍵>移除(Remove))
2.右鍵選取 Reference 并選擇 Add References
3.點(diǎn)選 Browse 并找到 OpticStudio 的安裝資料夾
4.以快捷鍵 Shift+Click 選擇 ZOSAPI 和 ZOSAPI_Interfaces,接著點(diǎn)選 Add
5.點(diǎn)選 OK 加入新的解
6.選取所有的 references,打開 Properties 接著將 Copy Local 由True 改為 False