csi_signal_simulate — INTERFEROMETRY `simulate` op

• 資料種類:無 → sweep(僅由參數決定的運算子 —— 不接受影像或資料輸入)

• 呼叫:import interferometry; interferometry.csi_signal_simulate(surface_um=6.0, z_start_um=0.0, z_step_um=0.05, n_planes=241, wavelength_um=0.6, envelope_fwhm_um=2.8, envelope_sigma_um=None, bias=0.5, amplitude=0.4, reflectivity=1.0, noise=0.0, seed=0)(或 opsinterferometry.get("csi_signal_simulate"))

用法

合成已知表面高度處一個像素的 z 掃描干涉圖。

> 以下的詳細說明為原文 —— 摘要與標題已翻譯。

The coherence-scanning forward model, and the reason every other operator here

has an exact answer to be checked against::

I(z) = bias + amplitude*reflectivity * exp(-(z-z0)^2 / 2 sigma^2)

• cos(4*pi*(z-z0)/wavelength)

with `z0 = surface_um`. The 4 is the double pass — light goes down to the

surface and back, so one fringe is `wavelength/2` of *height*, not a whole

wavelength. Getting that factor wrong is a clean 2x in every height this

module produces, which is why it is written out here rather than hidden in a

constant.

surface_um: the true surface height `z0`, in the scan's own

coordinate. Need not land on a scan plane — the

sub-step case is the interesting one and the tests use

it deliberately.

z_start_um/z_step_um/n_planes: the scan grid,

`z_k = z_start_um + k*z_step_um`.

wavelength_um: mean wavelength of the source.

envelope_fwhm_um: the FWHM of the envelope along the scan axis. Give

this *or* `envelope_sigma_um`, never both. It is

half the source coherence length, because the

double pass makes OPD = 2z;

:func:csi_design returns both under separate names

for exactly that reason.

bias/amplitude: the intensity pedestal `a and fringe amplitude b`.

reflectivity: per-pixel scale on the fringe amplitude (>= 0). It

scales the envelope and therefore

:func:csi_contrast_map, and — this is the honest part

— it does not move the envelope peak, so it does not

bias :func:csi_peak_position. A *spatially varying*

reflectivity biases nothing either; what does bias the

centroid is where the peak sits in the window, and that

is documented on :func:csi_peak_position.

noise: additive Gaussian sigma (0 = the exact model).

seed: integer seed for that noise (no `None`).

Returns a 1-D float64 array of `n_planes` intensities.

Ground truth: with `noise=0` and the surface centred in the scan, the

`"gaussian" estimator of :func:csi_peak_position` returns *surface_um* to

1.43e-07 um over sub-step offsets, and to 2.9e-14 um when the envelope is

given analytically instead of through the Hilbert transform (both pinned in

the tests).

Raises `ValueError`: a non-real / non-finite / string / bool parameter,

a non-positive `z_step_um / wavelength_um` / envelope width, a negative

`bias / amplitude / reflectivity / noise, n_planes` outside

`[3, MAX_SCAN_POINTS], a z_step_um at or past the wavelength_um/4`

Nyquist ceiling, and a *surface_um* outside the scan range (which is the case

a real instrument reports as "no surface found", not as a height).

詳細使用指南

• coherence_scanning 族使用指南

參考(範例資料・文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者・年份)與用途見上面的族使用指南。

可執行的範例(實際呼叫該運算子並已驗證的樣例)

• coherence_scanning — py -3.11 examples/coherence_scanning.py

型別可銜接的下一個運算子(可接受 `sweep` 作為輸入)

csi_envelope · csi_peak_position · chromatic_confocal_height

同類別(`simulate`)

csi_stack_simulate · chromatic_confocal_simulate

*Provenance: interferometry.py — INTERFEROMETRY 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

csi_signal_simulate — INTERFEROMETRY simulate op

用法