直方图特征的计算误差

Question

我只得到251块的峰度值，因为我把图像分割成块，我不知道出了什么问题，实际上它应该显示在252块上，请提前检查代码，谢谢。

clear all;
close all;

x =  imread('image009.png');
rgb = imresize(x,0.5);
green  = rgb(:,:,2);
[m,n]=size(green);
imnum=input('enter image number');
g=1;
q=1;
scaler=36;
scalec=37.5;
bno=1;
g=1;

for i=1:scaler:(m-scaler)
    if((i==1)||(i==(m-scaler))
        continue;    
    end

    for j=1:scalec:(n-scalec)

        if((j==1)||(j==(n-scalec)))
            continue;
        end
        tbpcolr = imcrop(green,[j,i,scalec,scaler])
        ip=input('press some key');
        mgreen = mean2(tbpcolr);
        stats = hisfeatures(tbpcolr);
        Mean(1,g) = stats(1)
        Variance(1,g) =stats(2)
        Skewness(1,g) =stats(3)
        Kurtosis(1,g)=stats(4)
        Energy(1,g) =stats(5)
        Entropy(1,g) =stats(6)
        Kurtosis = Kurtosis(~isnan(Kurtosis))
        Entropy = Entropy(~isnan( Entropy))
        Skewness = Skewness(~isnan(  Skewness))
        g=g+1;
        bno=bno+1; 
    end
end

blockno = [];
for  i = 1:252
    blockno(i)=i;
end

for  i = 1:252
    result_his=vertcat(Energy,Mean,Variance,Kurtosis,Skewness,Entropy,blockno); 
end

result_his=result_his.'

Answer 1

用NaN去除细胞会破坏基质结构。删除包含NaN的整行将丢弃实际数据。相反，统计工具箱有许多类似于其他MATLAB函数的函数，但它们将NaN值视为缺失，因此在计算中忽略它们。