100行で作る定数時間バイラテラルフィルタ - Qiita

J_p = \frac{\sum_{q \in \mathcal{N}_p} \exp\left(\frac{\|p-q\|_2^2}{-2\sigma_{s}^{2}}\right)\exp\left(\frac{\|I_p-I_q\|_2^2}{-2\sigma_{r}^{2}}\right)I_q}{{\sum_{q \in \mathcal{N}_p} \exp\left(\frac{\|p-q\|_2^2}{-2\sigma_{s}^{2}}\right)\exp(\frac{\|I_p-I_q\|_2^2}{-2\sigma_{r}^{2}})}} 畳み込みカーネルの範囲は集合で$\mathcal{N}_p$で表しています（2重ループで書くのめんどくさい）． $p,q$に対する位置変化に対する重みをスペースカーネル（一個目のexp），$I_p,I_q$に対する輝度変化に対する重