はじめに

凹凸係数 $r_{i j}$ とは、各画素 $x_{i j}$ をその周辺の平均輝度 $R_{i j}$ で除することで得られる係数である。

カーネルを $(2 n + 1) \times (2 n + 1)$ としたとき、凹凸係数 $r_{i j}$ は次式で求まる。
$r_{i j} = \frac{x_{i j}}{R_{i j}}$ $R_{i j} = \frac{1}{(2 n + 1)^{2}} \sum_{p = - n}^{n} \sum_{q = - n}^{n} x_{i + p, j + q} = \underset{n \times n}{b l u r} (x_{i j})$
$b l u r$ は平均値ぼかしである。
（※論文で触れられてないが $R_{i j} = 0$ のゼロ除算はどう定義しよう？）

方法

画像をグレースケールで読み込む
平均輝度 $R_{i j}$ を求める
凹凸係数 $r_{i j}$ を求める
凹凸係数 $r_{i j}$ に二値化処理を行う

プログラム

using System;  
using OpenCvSharp;  

namespace 凹凸係数法 {  
    class Program {  
        static void Main(string[] args) {  
            //①画像をグレースケールで読み込み  
            Mat src = new Mat(@"D:\8f10.jpeg", ImreadModes.Grayscale);  
            //Cv2.ImShow("src", src);  
            src.ConvertTo(src, MatType.CV_32FC1);  

            //②平均輝度を算出  
            Mat R = new Mat();  
            Cv2.Blur(src, R, new Size(21, 21));  

            //③凹凸係数を算出  
            Mat r = new Mat();  
            Cv2.Divide(src, R, r);  
            r *= 127;  
            r.ConvertTo(r, MatType.CV_8UC1);  

            //④二値化する  
            Mat dst = new Mat();  
            Cv2.Threshold(r, dst, 0.98, 255, ThresholdTypes.Otsu);  
            dst.ConvertTo(dst, MatType.CV_8UC1);  

            dst.ImWrite(@"D:\dst.png");  
            Cv2.WaitKey();  
        }  
    }  
}