renyi熵是對通常的夏農熵的擴充套件，算是q階廣義熵。公式如下：
其中p和夏農熵公式中的p一樣，是概率。當q=1時公式退化為夏農熵公式。（如何證明？看wiki吧）
有用此公式尋找影象
最佳二值化閾值的。
首先定義前景區域a，背景區域b。
那麼前景與背景區域畫素相應的renyi熵就如下定義：
其中k是當前取的灰度級，p(a)是畫素在a區域的概率，p(b)類似。當然，這裡說區域，不是指空間區域，是畫素灰度級區域。
最後影象renyi熵求最佳閾值定義如下：
這裡得到的k就是分割閾值。
分割效果如下：
原圖：
分割後：
在程式設計時還需要確定階數q，我取的是2。
matlab**如下：
matlab練習程式（DBSCAN）
和kmeans相比，不需要事先知道資料的類數。以程式設計的角度來考慮，具體演算法流程如下 1.首先選擇乙個待處理資料。2.尋找和待處理資料距離在設定半徑內的資料。3.將找到的半徑內的資料放到乙個佇列中。4.拿佇列頭資料作為當前待處理資料並不斷執行第2步。5.直到遍歷完佇列中所有資料，將這些資料記為一...
matlab練習程式（c c 呼叫matlab）
就我目前了解的c 呼叫matlab有兩種方法。第一種是通過matlab引擎呼叫，也就是這裡用到的方法。第二種是用matlab將m檔案編譯為相應的h lib dll檔案再加以呼叫。使用engine所用到的h和lib檔案基本在d program files matlab r2010b extern裡面，...
matlab練習程式（SUSAN檢測）
susan運算元既可以檢測角點也可以檢測邊緣，不過角點似乎比不過harris，邊緣似乎比不過canny。不過思想還是有點意思的。主要思想就是 首先做乙個和原影象等大的目標影象。然後用乙個圓形的模板，用模板去遍歷原影象每個畫素，把模板內的每個畫素都和模板中心畫素比較，如果灰度小於乙個閾值，那麼就對目標...