%形狀特征，Hu的七個不變矩
%The?Author?:?LongJie?Shi
%Data?:?2018.9.4
function?Hu?=?Hu_Moment（gray_image）%函數命名，輸入灰度圖像

format?short?e%有效數字5位+指數3位
image?=?double（gray_image）;%圖像矩陣運算時的數據類型是double類型
m00?=?sum（sum（image））;%0階矩圖像灰度值的和，第一個sum先計算列向量的和第二個sum計算總和
%1階中心矩賦值為0
m10?=?0;
m01?=?0;
[rowcol]?=?size（image）;%計算圖像尺寸儲存在[rowcol]
for?i?=?1:row?????????????????%計算1階中心矩的值
????for?j?=?1:col
????????m10?=?m10?+?i*image（ij）;
????????m01?=?m01?+?j*image（ij）;
????end
end
%歸一化的1階矩，代表圖像重心
n10?=?m10/m00;????
n01?=?m01/m00;

%2階3階中心矩賦初值為0
n20?=?0;n02?=?0;n11?=?0;n30?=?0;n12?=?0;n21?=?0;n03?=?0;
for?i?=?1:row??????????????????????????%計算2階3階中心矩
????for?j?=?1:col
????????n20=n20+i^2*image（ij）;
????????n02=n02+j^2*image（ij）;
????????n11=n11+i*j*image（ij）;
????????n30=n30+i^3*image（ij）;
????????n03=n03+j^3*image（ij）;
????????n12=n12+i*j^2*image（ij）;
????????n21=n21+i^2*j*image（ij）;
????end
end
%歸一化的2階3階中心矩
n20?=?n20/m00^2;
n02=n02/m00^2;
n11=n11/m00^2;
n30=n30/m00^2.5;
n03=n03/m00^2.5;
n12=n12/m00^2.5;
n21=n21/m00^2.5;

%Hu的七個不變矩由二階三階矩的線性組合構成
h1?=?n20?+?n02;???????????????????
h2?=?（n20-n02）^2?+?4*（n11）^2;
h3?=?（n30-3*n12）^2?+?（3*n21-n03）^2;??
h4?=?（n30+n12）^2?+?（n21+n03）^2;
h5?=?（n30-3*n12）*（n30+n12）*（（n30+n12）^2-3*（n21+n03）^2）+（3*n21-n03）*（n21+n03）*（3*（n30+n12）^2-（n21+n03）^2）;
h6?=?（n20-n02）*（（n30+n12）^2-（n21+n03）^2）+4*n11*（n30+n12）*（n21+n03）;
h7?=?（3*n21-n03）*（n30+n12）*（（n30+n12）^2-3*（n21+n03）^2）+（3*n12-n30）*（n21+n03）*（3*（n30+n12）^2-（n21+n03）^2）;
%利用Hu的不變矩作為特征向量
Hu?=?[h1?h2?h3?h4?h5?h6?h7];???

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1719??2018-09-04?10:06??matlab源碼計算Hu的七個不變矩\Hu_Moment.m
?????文件??????118742??2018-09-04?10:12??matlab源碼計算Hu的七個不變矩\Hu_Moment.mat