matlab 計算灰度圖像的一階矩、二階矩、三階矩

​   一階矩，定義了每個顏色分量的平均強度

​  二階矩，反映待測區域顏色方差，即不均勻性

​  三階矩，定義了顏色分量的偏斜度，即顏色的不對稱性

 

close all;clear all;clc;
J = imread('lena.jpg');
K = imadjust(J,[70/255 160/255],[]);
figure;
subplot(121),imshow(J);
subplot(122),imshow(K);
[m,n] = size(J);
mm = round(m/2);
mn = round(n/2);
[p,q] = size(K);
pp = round(p/2);
qq = round(q/2);
J = double(J);
K = double(K);
colorsum = 0.0;
Javg = mean2(J)  %求原圖像一階矩
Kavg = mean2(K)  %求增強對比度后的圖像一階矩
Jstd = std(std(J))  %求原圖像的二階矩，因為一次std函數表示按列求標准差，兩次std表示求整個矩陣的標准差
Kstd = std(std(K))   %求增強對比度后的圖像二階矩
for i=1:mm
    for j=1:mn
        colorsum = colorsum+(J(i,j)-Javg)^3;
    end
end
Jske = (colorsum/(mm*mn))^(1/3)  %求原圖像的三階矩
colorsum = 0.0;
for i=1:pp
    for j=1:qq
        colorsum = colorsum + (J(i,j)-Kavg)^3;
    end
end
Kske = (colorsum/(pp*qq))^(1/3)  %求增強對比度后的圖像三階矩
        

部分函數說明：

mean2(A) : 求矩陣A的均值

std（x,flag,dim）: 求x的標准偏差

std(x,0,1) : 0表示求標准差時除n-1，1表示按列划分  

std(x,1,2) : 1表示求標准差時除n，2表示按行划分