對稱場消矩測積法：不規則圖形的非成像總量式面積量測原理

版本：v0.2 重構稿 提出者：Neo.K 協作整理：Aletheia 審查參照：Theia 嚴格審查意見 核心名稱：SFM / DMR 中文名：對稱場消矩測積法；雙場消矩測積法 英文名：Symmetric Field Moment-Cancellation Method for Area Measurement 定位修正：本方法不主張作為一般數位面積演算法；其核心價值在於非成像、總量式、類比場量測條件下，利用互補代理場消去一階矩污染並恢復不規則區域面積。

v0.2 重構聲明

v0.1 將「雙環消矩測積法」放在不規則圖形數位計算方法論中心，並以旋轉體體積公式作為主要展示。經嚴格審查後，v0.2 對方法定位做出根本修正。

v0.1 中的代數關係是正確的：

\[ V_L=2\pi(RA+M_x) \]

\[ V_R=2\pi(RA-M_x) \]

\[ A=\frac{V_L+V_R}{4\pi R} \]

但若在數位環境中先由遮罩直接求出 \(A\) 與 \(M_x\)，再合成 \(V_L,V_R\)，再回推 \(A\)，此流程只能證明恆等式 \(A=A\)，不能證明一個獨立測量方法。即使改成逐像素加權求 \(V_L,V_R\)，它仍依賴對內容區像素的完整列舉；而像素列舉本身已足以直接求面積。因此，數位版 DCV/DMR 應被降級為「原理展示、sanity check、幾何矩解碼示範」，而非主方法。

v0.2 將核心翻轉為：本方法真正有生命的位置不是數位演算法，而是非成像、總量式、類比場量測。當儀器無法直接取得完整圖像、無法追蹤邊界、無法施加完美均勻場，卻能施加兩個互補權重場並取得總量訊號時，對稱場消矩測積法可透過兩次互補量測消除一階矩污染，恢復目標區域面積。

因此，本文從「雙環旋轉體求面積」改寫為「對稱代理場消矩測積」。旋轉體版本只是幾何例子；真正的 Kernel 是：

\[ Q_+=\int_D(a+b_xx+b_yy)\,dA \]

\[ Q_-=\int_D(a-b_xx-b_yy)\,dA \]

兩者相加：

\[ Q_++Q_-=2a\int_DdA=2aA \]

故：

\[ \boxed{A=\frac{Q_++Q_-}{2a}} \]

這是 v0.2 的主軸。

摘要

不規則圖形面積測量通常被理解為影像分割、像素計數、多邊形逼近、蒙地卡羅取樣、秤重或幾何積分問題。然而，許多實際量測場景並不允許取得完整圖像，也不允許追蹤圖形邊界；儀器可能只能取得一個總量訊號，例如總光通量、總吸收量、總熱反應、總電容耦合、總壓力響應或總質量變化。在這類條件下，面積測量的問題不再是「如何從完整形狀計算面積」，而是「如何從少量總量訊號中恢復目標面積」。

本文提出「對稱場消矩測積法」（Symmetric Field Moment-Cancellation Method, SFM；亦可稱 Dual Moment Removal, DMR）。其核心思想是：單一非均勻代理場對不規則區域的總量響應，通常不只包含面積，還會混入區域的一階矩，即質心偏移資訊。若場強沿空間呈線性梯度，則單次總量量測形式為：

\[ Q_+=aA+b_xM_x+b_yM_y \]

其中 \(A\) 是面積，\(M_x,M_y\) 是一階矩。若只量一次，面積與質心偏移無法分離。本文方法設計第二個互補場：

\[ Q_-=aA-b_xM_x-b_yM_y \]

兩次總量相加後，一階矩項抵消：

\[ Q_++Q_-=2aA \]

因此可由兩個互補場總量直接恢復面積：

\[ A=\frac{Q_++Q_-}{2a} \]

此方法的重點不是在數位遮罩上多做一次加權求和；那會退化為同義反覆。其真正價值在於非成像類比儀器：例如用兩個相反線性梯度光場照射不規則孔洞，只用總光量感測器接收通過光，而不拍攝影像、不追蹤邊界、不重建形狀。若兩個光場的均勻項相同、梯度項相反，則總光量之和消去質心污染，只保留孔洞面積。此原理可推廣到光學、熱場、電場、壓力場、化學反應場或其他可控代理場。

本文也重新定位旋轉體版本。將平面區域繞左右對稱外軸旋轉，得到：

\[ V_L=2\pi(RA+M_x),\qquad V_R=2\pi(RA-M_x) \]

這是對稱場消矩法的一個幾何實現，其中權重場為 \(2\pi(R\pm x)\)。但若體積由 \(A,M_x\) 反向合成，則只是恆等式驗證；若體積由獨立物理掃掠或獨立三維量測取得，才具有測量意義。

本文建立此方法的 Kernel、Runtime、Boundary、Error Model 與 Unit Test，並明確區分：數位版是教學/驗算工具；物理/類比場版才是方法本體。本文的主張不是「提出比像素法更有效的數位面積算法」，而是「提出一種在非成像、總量式量測條件下消除一階矩污染的對稱代理場原理」。

關鍵詞

不規則圖形、面積量測、非成像測量、總量式感測、代理場、對稱場、消矩、一階矩、質心污染、光通量、梯度場、旋轉體、帕普斯定理、積分幾何、SFM、DMR、雙場消矩測積法。

1. 問題重新定位：不是「算面積」，而是「從總量訊號恢復面積」

1.1 一般不規則圖形計算的預設

傳統面積計算方法通常預設至少能取得以下資訊之一：圖形完整邊界、完整內部遮罩、可列舉像素集合、可取樣判定函數、可剪下或稱重的實物區域、或可被掃描/成像的完整圖像。

若這些條件成立，最直接的方法通常已經足夠。多邊形有鞋帶公式，影像有像素計數，紙片有秤重法，解析區域有積分，複雜判定區域有蒙地卡羅。此時若額外引入旋轉體、三維體素或複雜代理場，往往只是把簡單問題複雜化。

因此，若把 SFM/DMR 放在「一般數位面積算法」的位置，它會被以下批評擊中：既然已經能列舉內容區，面積早就能直接得到；再用內容區去計算加權總量，最後回推面積，只是在做一個繞路的恆等式。

這個批評是正確的。

1.2 真正問題：只有總量，沒有圖像

本文真正要處理的是另一類問題：

給定一個不規則區域 \(D\)，儀器無法或不欲取得其完整圖像與邊界，只能取得某種總量響應 \(Q\)。如何從少數總量訊號中恢復面積 \(A\)？

這類問題與一般數位面積計算不同。其資料不是完整遮罩，而是類似：

\[ Q=\int_D w(x,y)\,dA \]

其中 \(w(x,y)\) 是儀器場、光場、熱場、電場、壓力場或反應場。儀器回傳的是一個總量，而非每一點的局部值。

若 \(w(x,y)\) 是常數：

\[ w(x,y)=a \]

則：

\[ Q=aA \]

面積可直接得到。但現實儀器很常無法提供完美均勻場。若場具有線性梯度：

\[ w(x,y)=a+b_xx+b_yy \]

則：

\[ Q=aA+b_xM_x+b_yM_y \]

此時單次量測不再只反映面積，而混入了一階矩，也就是圖形質心偏移。這就是本文所說的「一階矩污染」。

SFM/DMR 的目的不是從完整資料中重算面積，而是在總量式量測中設計第二個互補場，使污染項被消去。

2. 基本定義

令 \(D\subset\mathbb R^2\) 為可測不規則區域。面積定義為：

\[ A=M_{00}=\int_D1\,dA \]

一階矩定義為：

\[ M_x=M_{10}=\int_Dx\,dA \]

\[ M_y=M_{01}=\int_Dy\,dA \]

質心為：

\[ \bar x=\frac{M_x}{A},\qquad \bar y=\frac{M_y}{A} \]

代理場或權重場為：

\[ w:\mathbb R^2\to\mathbb R \]

總量式量測為：

\[ Q_w(D)=\int_Dw(x,y)\,dA \]

若 \(w\) 為仿射場：

\[ w(x,y)=a+b_xx+b_yy \]

則：

\[ Q_w(D)=aA+b_xM_x+b_yM_y \]

這裡的重點是：總量訊號 \(Q_w\) 同時包含面積與一階矩。若 \(M_x,M_y\) 未知，單一 \(Q_w\) 不能唯一決定 \(A\)。

3. 單通道代理場的失敗：面積與質心混疊

3.1 一維梯度情形

先考慮簡化情形：

\[ w_+(x)=R+x \]

則量測值：

\[ Q_+=\int_D(R+x)dA \]

展開：

\[ Q_+=R\int_DdA+\int_Dx\,dA \]

\[ Q_+=RA+M_x \]

若只知道 \(Q_+\)，則同時包含 \(A\) 與 \(M_x\)。不同形狀、不同面積、不同質心可能產生同一個 \(Q_+\)。因此單通道梯度場不是可靠面積量測器。

3.2 質心污染

若區域 \(D\) 整體偏向正 \(x\) 方向，\(M_x>0\)，則 \(Q_+\) 偏大。若偏向負 \(x\) 方向，\(M_x<0\)，則 \(Q_+\) 偏小。因此單一梯度場下，總量訊號不只反映「有多少面積」，也反映「面積分布在哪裡」。

這不是誤差，而是訊號本身的混合結構。若儀器場無法均勻，面積量測便會被位置分布污染。

3.3 為什麼平場不是總能解決？

若能施加完美均勻場，直接量測：

\[ Q_0=aA \]

即可。那 SFM/DMR 就沒有必要。但在以下情況中，均勻場可能不可得或不可靠：光源、熱源或電場天然有梯度；感測器靈敏度隨位置變化；樣品只能通過某個非均勻場區域；快速產線上無法做完整影像校正；場的均勻化成本高於做兩次互補量測；儀器幾何限制導致均勻場不可實現；量測目標很小、很快、很難成像；只能取得總量訊號，不能取得空間分布。

在這些場景中，方法的目標不是追求最完美的場，而是利用可控的非均勻性，設計互補場，使干擾項相互抵消。

4. 對稱場消矩原理

4.1 一維版本

設兩個互補場：

\[ w_+(x)=R+x \]

\[ w_-(x)=R-x \]

對同一區域 \(D\) 做兩次總量量測：

\[ Q_+=\int_D(R+x)dA=RA+M_x \]

\[ Q_-=\int_D(R-x)dA=RA-M_x \]

相加：

\[ Q_++Q_-=2RA \]

故：

\[ \boxed{A=\frac{Q_++Q_-}{2R}} \]

相減：

\[ Q_+-Q_-=2M_x \]

故：

\[ \boxed{M_x=\frac{Q_+-Q_-}{2}} \]

於是兩個量測通道同時給出面積與一階矩。體系的重點在於：未知的 \(M_x\) 不需要事先測出，它在相加時自動抵消，而在相減時可被讀出。

4.2 二維仿射場版本

設：

\[ w_+(x,y)=a+b_xx+b_yy \]

\[ w_-(x,y)=a-b_xx-b_yy \]

量測：

\[ Q_+=\int_D(a+b_xx+b_yy)dA \]

\[ Q_-= \int_D(a-b_xx-b_yy)dA \]

展開：

\[ Q_+=aA+b_xM_x+b_yM_y \]

\[ Q_-=aA-b_xM_x-b_yM_y \]

相加：

\[ Q_++Q_-=2aA \]

因此：

\[ \boxed{A=\frac{Q_++Q_-}{2a}} \]

相減：

\[ Q_+-Q_-=2(b_xM_x+b_yM_y) \]

這給出沿梯度方向的一階矩投影。若做多個方向的互補場，可恢復 \(M_x,M_y\)。

4.3 多方向版本

令有兩組互補場：

\[ w_{x+}=a+b x,\quad w_{x-}=a-b x \]

\[ w_{y+}=a+b y,\quad w_{y-}=a-b y \]

則：

\[ A_x=\frac{Q_{x+}+Q_{x-}}{2a} \]

\[ A_y=\frac{Q_{y+}+Q_{y-}}{2a} \]

理想情況：

\[ A_x=A_y=A \]

一階矩：

\[ M_x=\frac{Q_{x+}-Q_{x-}}{2b} \]

\[ M_y=\frac{Q_{y+}-Q_{y-}}{2b} \]

質心：

\[ \bar x=\frac{M_x}{A},\qquad \bar y=\frac{M_y}{A} \]

多方向版本具備自校驗意義，但只有在各通道為獨立物理量測時才有意義。若所有 \(Q\) 都由同一數位遮罩計算而來，則 \(A_x,A_y\) 的一致性只是恆等式，不是真正誤差檢查。

5. 旋轉體版本的正確定位

5.1 旋轉體是權重場的一個幾何實現

若將區域 \(D\) 繞左軸 \(x=-R\) 旋轉，點 \((x,y)\) 的旋轉半徑是 \(R+x\)。面元 \(dA\) 掃出的微小體積是：

\[ dV_L=2\pi(R+x)dA \]

故：

\[ V_L=\int_D2\pi(R+x)dA=2\pi(RA+M_x) \]

同理繞右軸 \(x=R\) 旋轉：

\[ V_R=\int_D2\pi(R-x)dA=2\pi(RA-M_x) \]

相加：

\[ V_L+V_R=4\pi RA \]

故：

\[ A=\frac{V_L+V_R}{4\pi R} \]

這個結果正確，但它只是 SFM/DMR 的一個特例，其中代理場為：

\[ w_L=2\pi(R+x),\qquad w_R=2\pi(R-x) \]

因此旋轉體不是方法的本體，而是幾何化的展示與可能的物理實作。

5.2 什麼時候旋轉體版本有測量意義？

只有當 \(V_L,V_R\) 是獨立於 \(A,M_x\) 被取得時，它才有測量意義。例如：真的把薄片掃掠成三維實體並測其體積；由 CAD 根據實物輪廓生成旋轉體，再用外部體積演算法量測；3D 掃描得到旋轉體，而非由已知面積合成；使用流體、沉積、材料成型等物理過程生成可量測體積。

若 \(V_L,V_R\) 是由以下方式得到：

\[ V_L=2\pi(RA+M_x) \]

\[ V_R=2\pi(RA-M_x) \]

其中 \(A,M_x\) 已由原圖直接算出，則該流程不是驗證方法，而是代數回路。

5.3 數位旋轉體驗證的降級

數位 voxel 旋轉體可保留，但其用途應明確限定為：視覺化、教學展示、幾何一致性檢查、軟體 pipeline 測試與 sanity check。

它不能被當作「比像素法更有效的面積算法」。因為若已有 2D 遮罩，像素法已經能直接給出面積；把它轉成 3D voxel 再回推只會增加計算成本與體素誤差。

6. 非成像光學版本：最小可行儀器模型

6.1 裝置構想

考慮一個不規則孔洞或遮罩 \(D\)。目標是測其開口面積，但不使用相機、不取得圖像、不追蹤邊界，只用總光量感測器。

裝置包含：可控制光場的照明模組、樣品平面、後方總光量感測器、兩次互補曝光控制、暗電流與背景校正、標準面積孔徑校準片。

第一次曝光施加：

\[ I_+(x,y)=I_0(a+b_xx+b_yy) \]

第二次曝光施加：

\[ I_-(x,y)=I_0(a-b_xx-b_yy) \]

透過孔洞的總光量：

\[ \Phi_+=\int_DI_+(x,y)dA \]

\[ \Phi_-=\int_DI_-(x,y)dA \]

展開：

\[ \Phi_+=I_0(aA+b_xM_x+b_yM_y) \]

\[ \Phi_-=I_0(aA-b_xM_x-b_yM_y) \]

相加：

\[ \Phi_++\Phi_-=2I_0aA \]

面積：

\[ \boxed{A=\frac{\Phi_++\Phi_-}{2I_0a}} \]

若 \(I_0a\) 不直接已知，可用已知面積 \(A_{\text{cal}}\) 的標準孔徑校準：

\[ K=2I_0a=\frac{\Phi_{+,\text{cal}}+\Phi_{-,\text{cal}}}{A_{\text{cal}}} \]

則未知面積：

\[ \boxed{A_{\text{unknown}}=\frac{\Phi_{+,\text{unknown}}+\Phi_{-,\text{unknown}}}{K}} \]

這個版本才真正避免循環：面積不是由像素數得到，而是由兩次獨立總光量量測恢復。

6.2 為什麼它不是普通光通量法？

普通光通量法要求均勻光場：

\[ \Phi=IA \]

若光場不均，需知道整個場分布或使用平場校正。SFM/DMR 的不同之處在於，它不要求每次光場均勻，而要求兩次光場互補，使一階不均項抵消。

普通光通量法要求：

\[ I(x,y)\approx I_0 \]

SFM/DMR 要求：

\[ I_+(x,y)+I_-(x,y)\approx 2I_0a \]

也就是說，單次不均可以接受；兩次總和必須均勻。這降低了單次場均勻性的要求，改為要求場的互補可控性。

6.3 真正優勢場景

此法可能有價值的場景包括：高速產線孔洞面積檢測、微孔/遮罩量測、高溫高壓或封閉環境中的非成像測量、低成本總量感測器、自然存在梯度的照明系統、只關心面積不關心形狀的快速檢測，以及樣品快速通過場區而只能取得時間序列總量訊號的情境。

7. 類比場的一般化

SFM/DMR 不限於光學。只要能建立總量式響應：

\[ Q=\int_Dw(x,y)dA \]

就可以使用互補場。

7.1 熱場版本

若區域 \(D\) 吸收熱量，且局部熱通量為：

\[ H_+(x,y)=a+b_xx+b_yy \]

\[ H_-(x,y)=a-b_xx-b_yy \]

總吸收熱量：

\[ Q_+=\int_DH_+dA \]

\[ Q_-=\int_DH_-dA \]

則：

\[ A=\frac{Q_++Q_-}{2a} \]

適合在無法成像、但可測總熱響應的場景中使用。

7.2 電場/電容版本

若一個不規則導電或介電區域在非均勻電場中產生總耦合訊號：

\[ C=\int_Dw(x,y)dA \]

則可用互補電極設計，讓一階位置偏差抵消。

7.3 壓力場版本

若薄片、孔洞或接觸區域在壓力敏感場中產生總力：

\[ F=\int_Dp(x,y)dA \]

且壓力場可做成互補梯度，則可由兩次總力恢復面積。

7.4 化學反應場版本

若某不規則區域暴露於濃度梯度場，總反應量為：

\[ Q=\int_Dc(x,y)dA \]

則兩個互補濃度場可消去位置偏差，恢復暴露面積。

7.5 數位濾波版本的正確用途

在數位影像中也能做：

\[ Q_+=\sum_{i\in D}w_+(x_i,y_i)\Delta A \]

\[ Q_-=\sum_{i\in D}w_-(x_i,y_i)\Delta A \]

但此版本不應宣稱是面積算法，因為遮罩已給出面積。它的合理用途是：模擬物理儀器、驗證場設計、測試誤差傳播、產生合成資料、教學展示幾何矩解碼。

8. 誤差模型

8.1 背景與暗訊號

總量感測器通常有背景項 \(B\)：

\[ S_+=Q_++B_+ \]

\[ S_-=Q_-+B_- \]

需先做空場或遮蔽校正：

\[ Q_+=S_+-B_+ \]

\[ Q_-=S_--B_- \]

若背景估計誤差為 \(\Delta B_+,\Delta B_-\)，則面積誤差：

\[ \Delta A_B=\frac{\Delta B_++\Delta B_-}{2a} \]

8.2 均勻項校準誤差

面積公式：

\[ A=\frac{Q_++Q_-}{2a} \]

若 \(a\) 有相對誤差 \(\Delta a/a\)，則：

\[ \frac{\Delta A}{A}\approx-\frac{\Delta a}{a} \]

因此需用標準面積校準片測得有效係數：

\[ K=2a \]

或在光學版本中：

\[ K=2I_0a \]

8.3 互補不完全誤差

理想情況：

\[ w_++w_-=2a \]

若實際為：

\[ w_++w_-=2a+\eta(x,y) \]

則：

\[ Q_++Q_-=2aA+\int_D\eta(x,y)dA \]

面積估計偏差：

\[ \Delta A_\eta=\frac{1}{2a}\int_D\eta(x,y)dA \]

若 \(\eta\) 很小且平均接近零，誤差可控；若 \(\eta\) 有系統性空間結構，則需要校準或增加更多互補場。

8.4 二階場殘差

若場不是仿射，而含二階項：

\[ w_+(x,y)=a+b_xx+b_yy+c_{xx}x^2+c_{yy}y^2+c_{xy}xy \]

若 \(w_-\) 只反轉一階項，但二階項相同，則相加後：

\[ Q_++Q_-=2aA+2c_{xx}M_{20}+2c_{yy}M_{02}+2c_{xy}M_{11} \]

此時面積估計會被二階矩污染。解法包括：場設計時讓二階項也互補抵消；做多組場並解聯立方程；使用校準樣品估計二階殘差；限制樣品尺寸，使二階項可忽略；改用正交場基底。

8.5 樣品定位誤差

若樣品座標相對場中心偏移，在單通道量測中會造成一階偏差。SFM/DMR 的一階互補結構可抵消主要位置偏差，但前提是兩次曝光中樣品位置不變，或位置變化可忽略。

若兩次量測之間樣品移動，則：

\[ Q_+=aA+bM_x^{(1)} \]

\[ Q_-=aA-bM_x^{(2)} \]

相加：

\[ Q_++Q_-=2aA+b(M_x^{(1)}-M_x^{(2)}) \]

若 \(M_x^{(1)}\neq M_x^{(2)}\)，消矩不完全。因此高速或動態系統需同步曝光、快速切換或同時雙通道量測。

8.6 感測器非線性

若感測器輸出不是線性的：

\[ S=f(Q) \]

則不能直接相加。需在線性區操作，或先做反函數校正：

\[ Q=f^{-1}(S) \]

SFM/DMR 要求代理量在合成前可線性相加。

8.7 雜訊傳播

若兩次量測噪聲標準差為 \(\sigma_+,\sigma_-\)，且獨立，則：

\[ \sigma_A=\frac{\sqrt{\sigma_+^2+\sigma_-^2}}{2a} \]

若兩通道噪聲相關，需加入協方差：

\[ \sigma_A^2=\frac{\sigma_+^2+\sigma_-^2+2\operatorname{Cov}(Q_+,Q_-)}{(2a)^2} \]

9. 非循環驗證標準

v0.2 必須明確規定：什麼樣的驗證才算真正驗證方法，而不是驗證恆等式。

9.1 無效驗證

以下流程無效：

從遮罩直接算 \(A\)；
從遮罩直接算 \(M_x\)；
用公式合成 \(Q_+=aA+bM_x\)、\(Q_-=aA-bM_x\)；
再用 \(Q_++Q_-\) 回推 \(A\)。

這只是代數閉環。

9.2 弱驗證

以下流程可作為 sanity check，但不能證明方法有測量價值：

從遮罩逐像素計算 \(Q_+,Q_-\)；
用公式回推面積；
與像素面積比較。

原因是逐像素列舉已包含面積資訊。此流程可驗證程式實作與公式一致，但不能證明比像素法更有用。

9.3 有效數值模擬驗證

有效數值模擬應模擬儀器限制：

產生未知區域 \(D\)，但測量器不能直接讀出 \(A\)；
只允許測量器接收總量訊號 \(Q_+,Q_-\)；
加入場不均、暗電流、噪聲、二階殘差；
用校準片估計 \(K\)；
用兩次總量訊號回推 \(A\)；
最後才用真值 \(A\) 做誤差評估。

關鍵在於：演算法主流程不能直接使用 \(A\) 或完整像素面積，只能使用總量。

9.4 有效物理驗證

有效物理驗證應包含：標準面積孔徑校準、未知不規則孔洞樣品、兩個互補場、總量感測器、背景扣除、場互補性測試、與獨立方法比較如顯微影像像素法或秤重法，以及誤差與適用範圍報告。

若 \(Q_+,Q_-\) 來自真實儀器，而非由已知面積合成，則方法才真正成立為量測方法。

10. 實驗設計草案

10.1 光學最小實驗

器材

可調亮度螢幕或投影面板；
不規則孔洞遮罩；
光電二極體或照度感測器；
暗箱；
已知面積校準孔；
控制程式。

步驟

建立第一張光場圖案 \(I_+(x)=I_0(a+bx)\)；
建立第二張光場圖案 \(I_-(x)=I_0(a-bx)\)；
用已知面積孔校準 \(K\)；
放入未知不規則孔洞；
量測通過總光量 \(\Phi_+,\Phi_-\)；
扣除暗電流與背景；
計算：

\[ A=\frac{\Phi_++\Phi_-}{K} \]

用相機拍攝後像素法作為獨立真值比較。

成功標準

單次梯度場估計會隨孔洞位置偏移；
雙場相加後對位置偏移不敏感；
面積估計接近像素法或標準孔徑結果；
誤差隨互補場校準改善而下降。

10.2 位置偏移測試

將同一孔洞遮罩在 \(x\) 方向移動。理論上：

單通道：

\[ Q_+=aA+bM_x \]

會隨位置改變。

雙通道：

\[ Q_++Q_-=2aA \]

應近似不變。

這是 SFM/DMR 最關鍵的實驗展示：它不是單純測面積，而是抵抗一階位置污染。

10.3 非均勻場測試

設計故意不均勻的光場，讓普通單通量法失效。比較三種方法：單一梯度場、假設均勻場的普通光通量法、SFM/DMR 雙場法。

若雙場法在一階不均下恢復面積，則其獨特性得到支持。

10.4 二階殘差測試

加入二階場：

\[ I(x)=a+bx+cx^2 \]

若只做一階互補，殘差會出現。測試目標是估計方法邊界，而不是假裝它萬能。

11. 方法比較

| 方法 | 是否需成像 | 是否需邊界 | 是否需均勻場 | 是否抗一階位置污染 | 核心限制 | |---|---:|---:|---:|---:|---| | 像素法 | 是 | 間接需要 | 否 | 是 | 需完整影像與分割 | | 多邊形法 | 是/需輪廓 | 是 | 否 | 是 | 需乾淨邊界 | | 秤重法 | 否 | 否 | 不適用 | 是 | 材料面密度需均勻 | | 普通光通量法 | 否 | 否 | 是 | 否 | 光場需均勻 | | 單梯度場量測 | 否 | 否 | 否 | 否 | 混入質心 | | SFM/DMR 雙場法 | 否 | 否 | 不要求單次均勻，只要求互補 | 是，一階 | 需互補場與線性總量感測 |

這張表說明：SFM/DMR 的價值不在取代像素法，而在取代「不均勻場下的單通道總量量測」。

12. Kernel / Runtime / Boundary / Unit Test

12.1 Kernel Layer

核心形式：

\[ Q_+=\int_D(a+\ell(x,y))dA \]

\[ Q_-=\int_D(a-\ell(x,y))dA \]

其中 \(\ell(x,y)\) 是一階線性項：

\[ \ell(x,y)=b_xx+b_yy \]

則：

\[ Q_++Q_-=2aA \]

\[ A=\frac{Q_++Q_-}{2a} \]

若 \(\ell\) 可被方向化設計，差分可取得一階矩投影：

\[ Q_+-Q_-=2\int_D\ell(x,y)dA \]

12.2 Runtime Layer

類比量測 Runtime

Input:
    unknown aperture/region D
    complementary fields w+, w-
    total-response sensor
    calibration aperture

Process:
    1. measure background B+, B-
    2. measure calibration responses C+, C-
    3. compute K = (C+ + C-) / A_cal
    4. place unknown D
    5. measure S+, S-
    6. subtract background: Q+ = S+ - B+, Q- = S- - B-
    7. estimate A = (Q+ + Q-) / K
    8. optionally estimate moment projection from Q+ - Q-
    9. report uncertainty

數位模擬 Runtime

Input:
    hidden test region D
    simulated fields w+, w-
    sensor model with noise and background

Restriction:
    measurement module may not directly output pixel count area

Process:
    1. generate Q+ = sensor_integral(D, w+)
    2. generate Q- = sensor_integral(D, w-)
    3. apply calibration
    4. estimate A
    5. compare to ground truth only after estimate is produced

12.3 Boundary Layer

下界必要條件

區域 \(D\) 的總量響應可被定義；
場響應近似線性可加；
兩個場具有可校準的互補關係；
均勻項 \(a\) 或總校準係數 \(K\) 可取得；
樣品在兩次量測中不發生不可忽略變形或位移；
背景訊號可扣除；
二階以上場殘差可忽略、可校準或可建模。

上界排除條件

以下情況不應宣稱 SFM/DMR 有效：兩個場不互補；感測器嚴重非線性且未校正；樣品在兩次量測間移動；場含大量二階以上殘差且未處理；用已知 \(A\) 合成 \(Q\) 再回推 \(A\)；已有完整數位遮罩卻宣稱 SFM 是更有效面積算法；把旋轉體視覺化當作獨立實驗驗證；把一階抵消誇大成任意誤差抵消。

12.4 Unit Test

Test 1：標準孔徑

已知面積 \(A_{\text{cal}}\)，檢查：

\[ K=\frac{Q_++Q_-}{A_{\text{cal}}} \]

是否穩定。

Test 2：同面積不同位置

同一圖形平移。單通道 \(Q_+\) 應變化，雙通道和 \(Q_++Q_-\) 應近似不變。

Test 3：不同形狀同面積

測試不同形狀但同面積樣品。雙通道估計應接近相同面積；差分則反映不同質心。

Test 4：二階殘差

故意加入二階場，觀察估計偏差。用於界定方法上界。

Test 5：與獨立方法比較

使用像素法、秤重法或標準孔徑作為外部參照，不得用 SFM 自己生成真值。

13. 與帕普斯定理的關係

帕普斯定理說明，平面區域繞不相交外軸旋轉，體積等於面積乘以質心路徑長度：

\[ V=2\pi A\bar r \]

若旋轉軸為 \(x=-R\)，則：

\[ \bar r=R+\bar x \]

\[ V_L=2\pi A(R+\bar x) \]

因為：

\[ M_x=A\bar x \]

故：

\[ V_L=2\pi(RA+M_x) \]

單軸版本需要質心，因此不能單靠 \(V_L\) 求 \(A\)。雙軸版本引入右軸：

\[ V_R=2\pi(RA-M_x) \]

相加消去 \(M_x\)。因此，旋轉體形式可理解為：

帕普斯關係的雙通道消矩化。

但 v0.2 更進一步指出，帕普斯只是其中一個幾何例子。SFM/DMR 的本體是互補代理場，不必真的旋轉，也不必真的生成三維體。

14. 與積分幾何與矩理論的關係

SFM/DMR 可被放在幾何矩框架中理解。面積是零階矩：

\[ M_{00}=\int_D1\,dA \]

質心資訊來自一階矩：

\[ M_{10}=\int_Dx\,dA \]

\[ M_{01}=\int_Dy\,dA \]

單一仿射場量測得到的是零階矩與一階矩的線性組合：

\[ Q=aM_{00}+b_xM_{10}+b_yM_{01} \]

兩個互補場相當於建立一組線性方程，使 \(M_{00}\) 可被分離。若增加更多場，可進一步求解更高階矩。

因此，SFM/DMR 的更一般表述是：

透過可控代理場取得幾何矩的線性觀測，再透過對稱設計或聯立解碼分離目標矩。

這可連接到更廣義的形狀量測、矩特徵、非成像感測與壓縮感知式幾何測量。但本文 v0.2 僅主張一階消矩測積，不主張已完成完整形狀重建理論。

15. 命名修正

v0.1 使用 DCV（Dual Curvature Volume Method）。經審查後，此名稱不夠準確，原因是方法真正雙的不是曲率，而是代理場、量測通道與一階矩符號。

建議名稱：

中文正式名：對稱場消矩測積法
中文簡稱：雙場消矩法
英文正式名：Symmetric Field Moment-Cancellation Method for Area Measurement
英文縮寫：SFM 或 DMR

本文保留 DMR 作為方法族縮寫，SFM 作為具體場量測版本縮寫。

16. 方法論判決：它到底是不是一個方法？

答案取決於它所處的資料環境。

16.1 在完整數位遮罩環境中

若已知每個像素是否屬於 \(D\)，則面積可直接計算：

\[ A=N\Delta A \]

此時 SFM/DMR 不是必要方法，只是展示：

\[ \int_D1\,dA \]

可被寫成兩個互補權重積分的平均。

判決：教學工具，不是主算法。

16.2 在旋轉體由公式合成時

若 \(V_L,V_R\) 由已知 \(A,M_x\) 合成，則是恆等式。

判決：sanity check，不是驗證。

16.3 在獨立物理旋轉體量測時

若 \(V_L,V_R\) 由真實掃掠體積量得，則可回推面積。

判決：可作為物理測量方法，但工程成本需評估。

16.4 在非成像總量式互補場量測時

若 \(Q_+,Q_-\) 由兩個互補物理場獨立量得，且儀器無法直接取得完整圖像，則 SFM/DMR 能消去一階矩污染並恢復面積。

判決：方法本體成立。

17. 與 v0.1 的差異摘要

| 項目 | v0.1 | v0.2 | |---|---|---| | 主名稱 | 雙環消矩測積法 / DCV | 對稱場消矩測積法 / SFM-DMR | | 主體 | 旋轉體積 | 互補代理場總量量測 | | 數位版定位 | 主要驗證之一 | 降級為 sanity check | | 類比場定位 | 未來工作 | 核心章節 | | 創新點 | 雙旋轉體消質心 | 雙代理場消一階矩污染 | | 驗證要求 | 公式與 voxel 測試 | 非循環、總量式、獨立量測 | | 應用場景 | 不規則圖形計算 | 非成像面積量測、抗場梯度污染 | | 命名問題 | Dual Curvature 不準 | Symmetric Field / Dual Moment Removal |

18. 結論

v0.2 對方法做出根本修正。本文不再把雙環旋轉體版本宣稱為一般不規則圖形數位面積算法，而是將其重新定位為「對稱場消矩測積法」的一個幾何例子。方法的真正核心是：當總量式量測訊號混入一階矩，即質心偏移污染時，設計兩個互補代理場，使一階項以相反符號出現，並在相加時抵消，只保留零階面積項。

核心公式為：

\[ Q_+=aA+b_xM_x+b_yM_y \]

\[ Q_-=aA-b_xM_x-b_yM_y \]

\[ A=\frac{Q_++Q_-}{2a} \]

這使方法從「漂亮的面積恆等式」轉為「抗一階場不均/位置偏移污染的非成像總量式量測原理」。

它不應與像素法在完整數位圖像場景中競爭；在那裡它必然多餘。它應該進入的戰場，是無法或不欲成像、無法追邊界、只能取得總量訊號、且單通道總量會被場梯度污染的類比測量場景。若在這樣的場景中，\(Q_+,Q_-\) 能由獨立物理量測取得，並透過校準消除背景與互補誤差，則 SFM/DMR 成立為一種真正的測量方法。

因此，本文最終主張是：

不規則圖形面積可以被視為零階矩。單一非均勻代理場會把零階矩與一階矩混合；對稱互補場可以使一階矩湮滅，將面積從總量訊號中解碼出來。

這是從「計算面積」走向「設計代理場與量測結構」的升級。

附錄 A：最小公式集

\[ A=\int_DdA \]

\[ M_x=\int_Dx\,dA \]

\[ M_y=\int_Dy\,dA \]

\[ Q_+=\int_D(a+b_xx+b_yy)dA \]

\[ Q_-= \int_D(a-b_xx-b_yy)dA \]

\[ Q_+=aA+b_xM_x+b_yM_y \]

\[ Q_-=aA-b_xM_x-b_yM_y \]

\[ A=\frac{Q_++Q_-}{2a} \]

\[ b_xM_x+b_yM_y=\frac{Q_+-Q_-}{2} \]

一維版本：

\[ Q_+=RA+M_x \]

\[ Q_-=RA-M_x \]

\[ A=\frac{Q_++Q_-}{2R} \]

旋轉體版本：

\[ V_L=2\pi(RA+M_x) \]

\[ V_R=2\pi(RA-M_x) \]

\[ A=\frac{V_L+V_R}{4\pi R} \]

附錄 B：非循環驗證偽代碼

Given:
    hidden region D
    calibration region C with known area A_C
    complementary fields w_plus, w_minus
    sensor model Sensor(D, w) -> scalar

Calibration:
    B_plus  = Sensor(empty, w_plus)
    B_minus = Sensor(empty, w_minus)

    C_plus  = Sensor(C, w_plus)  - B_plus
    C_minus = Sensor(C, w_minus) - B_minus

    K = (C_plus + C_minus) / A_C

Unknown:
    S_plus  = Sensor(D, w_plus)  - B_plus
    S_minus = Sensor(D, w_minus) - B_minus

Estimate:
    A_hat = (S_plus + S_minus) / K

Validation:
    Compare A_hat with external ground truth only after estimate is produced.

附錄 C：一句話版本

對稱場消矩測積法不是用數位遮罩繞路算面積，而是在只能取得總量訊號的非成像量測中，用兩個互補場讓質心偏移項互相湮滅，從而把不規則區域的面積解碼出來。

附錄 D：版本註記

v0.2 是對 v0.1 的定位修正版。其主要變更為：

將方法本體從旋轉體積改為互補代理場；
將數位版降級為 sanity check；
將類比場版本升為核心；
新增非循環驗證標準；
新增光學最小實驗設計；
修正命名，降低「Dual Curvature」用語；
明確指出適用邊界與不應宣稱的場景。

後續 v0.3 可補：可執行的非循環感測模擬程式、光學實驗圖示、場互補誤差的數值測試、二階殘差補償、與平場光通量法的實驗比較，以及更嚴格的命題/證明格式。

原始檔（供 RAG/下載）：/raw/lm-000451.md [md] · id: lm-000451