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

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

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

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

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

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

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

版本：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 可補：可執行的非循環感測模擬程式、光學實驗圖示、場互補誤差的數值測試、二階殘差補償、與平場光通量法的實驗比較，以及更嚴格的命題/證明格式。