# 對稱場消矩測積法：不規則圖形的非成像總量式面積量測原理 **版本**：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 無效驗證 以下流程無效： 1. 從遮罩直接算 \(A\)； 2. 從遮罩直接算 \(M_x\)； 3. 用公式合成 \(Q_+=aA+bM_x\)、\(Q_-=aA-bM_x\)； 4. 再用 \(Q_++Q_-\) 回推 \(A\)。 這只是代數閉環。 ## 9.2 弱驗證 以下流程可作為 sanity check，但不能證明方法有測量價值： 1. 從遮罩逐像素計算 \(Q_+,Q_-\)； 2. 用公式回推面積； 3. 與像素面積比較。 原因是逐像素列舉已包含面積資訊。此流程可驗證程式實作與公式一致，但不能證明比像素法更有用。 ## 9.3 有效數值模擬驗證 有效數值模擬應模擬儀器限制： 1. 產生未知區域 \(D\)，但測量器不能直接讀出 \(A\)； 2. 只允許測量器接收總量訊號 \(Q_+,Q_-\)； 3. 加入場不均、暗電流、噪聲、二階殘差； 4. 用校準片估計 \(K\)； 5. 用兩次總量訊號回推 \(A\)； 6. 最後才用真值 \(A\) 做誤差評估。 關鍵在於：演算法主流程不能直接使用 \(A\) 或完整像素面積，只能使用總量。 ## 9.4 有效物理驗證 有效物理驗證應包含：標準面積孔徑校準、未知不規則孔洞樣品、兩個互補場、總量感測器、背景扣除、場互補性測試、與獨立方法比較如顯微影像像素法或秤重法，以及誤差與適用範圍報告。 若 \(Q_+,Q_-\) 來自真實儀器，而非由已知面積合成，則方法才真正成立為量測方法。 --- # 10. 實驗設計草案 ## 10.1 光學最小實驗 ### 器材 1. 可調亮度螢幕或投影面板； 2. 不規則孔洞遮罩； 3. 光電二極體或照度感測器； 4. 暗箱； 5. 已知面積校準孔； 6. 控制程式。 ### 步驟 1. 建立第一張光場圖案 \(I_+(x)=I_0(a+bx)\)； 2. 建立第二張光場圖案 \(I_-(x)=I_0(a-bx)\)； 3. 用已知面積孔校準 \(K\)； 4. 放入未知不規則孔洞； 5. 量測通過總光量 \(\Phi_+,\Phi_-\)； 6. 扣除暗電流與背景； 7. 計算： \[ A=\frac{\Phi_++\Phi_-}{K} \] 8. 用相機拍攝後像素法作為獨立真值比較。 ### 成功標準 1. 單次梯度場估計會隨孔洞位置偏移； 2. 雙場相加後對位置偏移不敏感； 3. 面積估計接近像素法或標準孔徑結果； 4. 誤差隨互補場校準改善而下降。 ## 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 ```text 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 ```text 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 ### 下界必要條件 1. 區域 \(D\) 的總量響應可被定義； 2. 場響應近似線性可加； 3. 兩個場具有可校準的互補關係； 4. 均勻項 \(a\) 或總校準係數 \(K\) 可取得； 5. 樣品在兩次量測中不發生不可忽略變形或位移； 6. 背景訊號可扣除； 7. 二階以上場殘差可忽略、可校準或可建模。 ### 上界排除條件 以下情況不應宣稱 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：非循環驗證偽代碼 ```text 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 的定位修正版。其主要變更為： 1. 將方法本體從旋轉體積改為互補代理場； 2. 將數位版降級為 sanity check； 3. 將類比場版本升為核心； 4. 新增非循環驗證標準； 5. 新增光學最小實驗設計； 6. 修正命名，降低「Dual Curvature」用語； 7. 明確指出適用邊界與不應宣稱的場景。 後續 v0.3 可補：可執行的非循環感測模擬程式、光學實驗圖示、場互補誤差的數值測試、二階殘差補償、與平場光通量法的實驗比較，以及更嚴格的命題/證明格式。