查看“非局部平均”的源代码

'''非局部平均'''（Non-local means）是一種[[影像降噪]]的演算法，相較於局部（local）的演算法（如[[高斯模糊]]、[[非等向性擴散]]）只使用各個目標像素附近的點來將影像平滑化來去除雜訊，非局部平均演算法則對各個目標像素周圍定義一個區塊，並且對整個影像的所有像素依照該像素周圍區塊的區塊與目標像素區塊的相似度賦予權重、進行平均，如此可以使經過處理的影像更為清晰，並且損失較少的細節。<ref>{{Cite web|url=https://ieeexplore.ieee.org/document/1467423/|title=A non-local algorithm for image denoising - IEEE Conference Publication|accessdate=2018-06-22|work=ieeexplore.ieee.org|language=en-US}}</ref>

相較於其他的影像降噪演算法，非局部平均在的方法雜訊（method noise，定義為影像及其降噪後的結果的差值）與[[白雜訊]]較為相似，而通常認為方法雜訊應與白雜訊盡可能相似，因此這是一個正向的結果。<ref>{{Cite web|url=http://www.123seminarsonly.com/Seminar-Reports/029/42184313-10-1-1-100-81.pdf|title=On Image Denoise Method|accessdate=|author=Buades, Antoni|date=|publisher=}}</ref>非局部平均也被延伸至其他應用如[[去交錯]]及影像內插等。

== 雜訊模型 ==
如同其他影像降噪演算法，最簡單並且普遍的雜訊模型為在原始的圖片上加上白雜訊：

<math>v(i) = u(i) + n(i)</math>

其中<math>u(i)</math>為原始的圖片，<math>n(i)</math>為雜訊，在白雜訊的情況下為平均0標準差<math>\sigma</math>的i.i.d. [[高斯函数|高斯函數]]，<math>v(i)</math>為所觀測到包含雜訊的圖片。

== 演算法 ==
非局部平均的定義為：<ref name=":0">{{Cite journal|title=Non-Local Means Denoising|url=http://www.ipol.im/pub/art/2011/bcm_nlm/|last=Buades|first=Antoni|last2=Coll|first2=Bartomeu|date=2011-09-13|journal=Image Processing On Line|doi=10.5201/ipol.2011.bcm_nlm|volume=1|language=en|issn=2105-1232|last3=Morel|first3=Jean-Michel}}</ref>

<math>NL_u(x)=\frac{1}{C(x)}\int_{\Omega}^{} f(d(B(x), B(y)))\ u(y)\ dy</math>

其中<math>u</math>為要處理的圖片，<math>\Omega</math>為整張圖片的區域，<math>B(x)</math>為以<math>x</math>為中心的一個區塊

<math>d(x,y)</math>為<math>x</math>及<math>y</math>的[[欧几里得距离|歐幾里得距離]]：<math>d(x,y) = |x-y|_2^2</math>

<math>f</math>為一遞減函數，常見的如：<math>f(x)=exp(-x^2/h^2)</math>

<math>C(x)</math>為標準化的因子：<math>C(x)=\int_\Omega f(d(B(x),B(y)))\ dy</math>

此式可以解釋為圖上一點<math>x</math>經過降噪後的值為整張圖片上所有點的加權平均，其中每個點的權重為該點的附近區塊與<math>x</math>附近區塊的相似度（將兩個區塊各自的點以相同的排列視為一向量，計算歐幾里得距離），再經過一指數衰減的的函數（權重將落在(0,1]區間）。

相較於局部的演算法指考慮了每個點附近的點，非局部平均考慮了圖上所有的點，故稱為非局部。由於一張圖片中最相似的點不一定是距離近的點，反之亦然，故搜尋整張圖片上相似的點，利用週期性重複出現的部分如材質紋理或是伸長的邊緣等進行降噪可以得到更好的結果，但由於對每個點都要搜尋整張圖片上其他的點來比較相似度，故運算複雜度往往會比局部的演算法高。

== 實作 ==
上述的式子所定義的演算法為連續的，無法在實際的數位影像中使用，在實際應用中離散化的版本為：<ref name=":0" />

<math>NL_u[i]=\frac{1}{C(i)}\sum_{j\in\Omega}w(i,j)u(j)</math>

<math>C(i)=\sum_{j\in\Omega}w(i,j)</math>

其中<math>w(i,j)</math>為以<math>i</math>及<math>j</math>為中心的點的區塊的相似度，作為計算點<math>i</math>降噪後的值時點<math>j</math>的權重，常見的作法為計算以<math>i</math>及<math>j</math>為中心的正方形區塊的歐幾里得距離的平方，再經過一個指數遞減的函數，與前一節連續版本的例子相似，以彩色圖片為例，完整的式子如下：

<math>w(i,j)=e^{-\frac{max(d^2(B(i,f),B(j,f))-2\sigma^2, 0)}{h^2}}</math>

<math>d^2(B(i,f),B(j,f))=\frac{1}{3(2f+1)^2}\sum_{m=1}^3\sum_{n\in B(0,f)}(u_m(i+n)-u_m(j+n))^2</math>

<math>B(i,f)</math>為以點<math>i</math>為中心寬度為<math>2f+1</math>的正方形區域，即<math>[i_x-f,i_x+f]\times[i_y-f,i_y+f]</math>區間。

<math>u</math>的下標代表圖片每個像素值的三個分量（如RGB），黑白影像則無此項。

<math>\sigma</math>為雜訊的標準差，<math>h</math>為與<math>\sigma</math>相關的濾波器參數。

此權重函數可解釋為將區塊相似度大於某個程度（歐幾里得距離的平方小於等於<math>2\sigma^2</math>）的所有點都給予權重1，區塊相似度較小的點權重則由於指數遞減函數迅速的衰減。

== 區塊級實作 ==
除了上述以像素為單位的實作外（pixelwise），還有一種以區塊為單位的實作（patchwise implementation）：<ref name=":0" />

對於一個以<math>i</math>為中心的<math>(2f+1)\times(2f+1)</math>區塊<math>B(i,f)</math>，處理後的值

<math>\hat{B}_m(i,f)=\frac{1}{C}\sum_{Q(j,f)\in B(i,r)} u_m(Q)w(B,Q)</math>

其中<math>w</math>的算法與前一節的實作相同，<math>B(i,r)</math>以<math>r</math>的大小限制搜索區域以降低計算複雜度。

最後將這些值平均可得降噪後的圖片

<math>NL_u[i]=\frac{1}{(2f+1)^2}\sum_{j\in B(i,f)}\hat{B}(j)</math>

與前一節像素級的實作相比，由於最後一步集成的步驟降低了雜訊，PSNR較高，並且也減少了在邊緣附近震盪的雜訊，不過在保留圖片細節的方面則沒有明顯的進步。

== 參考來源 ==
{{reflist}}

[[Category:圖像處理]]