simcal_fit — SIMCAL `fit` op

• 数据种类:simmodel → calresult

• 调用: import fullseye as fs; fs.ledger.simcal_fit(model=None, method='lm', n_starts=16, fix=None, threshold=10.0, allow_undetermined=False, seed=0) -> 'dict'(要直接调用实现,import simcal; simcal.simcal_fit(model=None, method='lm', n_starts=16, fix=None, threshold=10.0, allow_undetermined=False, seed=0) -> 'dict';从台账取用则 opssimcal.get("simcal_fit"))

用法

将带观测的 simmodel 的自由旋钮拟合到其图像（加权最小二乘）-> calresult。

> 以下的详细说明为原文 —— 摘要与标题已翻译。

Minimises chi^2 = sum over valid pixels of (observed − simulated)^2 / variance, the variance coming from the noise model at the current prediction (re-weighted each

iteration: the quasi-likelihood estimating equation, unbiased for a signal-dependent variance). Knobs live in u in [0, 1] (log knobs on a log axis).

Starting point: the initial knobs and n_starts Latin-hypercube points over the ranges are scored by chi^2 and the best is used (the design's "start from the best point of the plan"). method: "lm" (Levenberg-Marquardt, Marquardt-scaled damping; a raw step that leaves the range is not taken — its projection onto the range is tried, a knob on a range end whose descent points outward is held there, and a point the forward model refuses rejects the step and raises the damping, as in lensopt), "nm" (Nelder-Mead; a point outside the range scores the metric's worst value via metriccontract.value_or_worst), "multistart" (LM from the best n_starts // 4 plan points, at least 2, ranked with metriccontract.best_of — for multi-modal problems such as thin-film fringe orders). fix = {knob: value} pins knobs for this fit.

Refusal (metriccontract.MetricContractError, a ValueError) when, at the end point, a knob has no influence (a full-range move changes chi^2 by less than 3.84), is confounded with others (collinearity index gamma_K > threshold, sensitivity.mv_confounding), or sits within 1 % of a range end — the message names the knobs, the reason and the combinations the data do fix. With allow_undetermined=True those knobs get no value (undetermined holds the reason, status = partial) and determined_combos carries the combinations with their values (a product of powers for log knobs). No noise model / no observed images:

refused.

详细使用指南

• sim_to_real_calibration 族使用指南

参考(示例数据、文献)

• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。

• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。

• 算法的正典(作者、年份)与用途见上面的族使用指南。

可运行的示例(实际调用该算子并已验证的样例)

• poc_pocket_interreflection — py -3.11 examples/poc_pocket_interreflection.py

• poc_simcal_thin_film — py -3.11 examples/poc_simcal_thin_film.py

• simcal_checks_tour — py -3.11 examples/simcal_checks_tour.py

• simcal_tour — py -3.11 examples/simcal_tour.py

类型可衔接的下一个算子(可接受 `calresult` 作为输入)

simcal_identifiability · simcal_uncertainty · simcal_validate · simcal_discrepancy · simcal_report

同类别(`fit`)

simcal_screen

*Provenance: simcal.py — SIMCAL 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

simcal_fit — SIMCAL fit op

用法