颜色距求法代码

close all;% 关闭所有图形窗口
clear all;%清除工作空间所有变量
clc;%清空命令行
J=imread('C:\F\MATLAB\meinv.jpg');%读入要处理的清晰图像，并赋值给J
K=imadjust(J,[70/255,160/255],[]);%灰度级调整将【70 160】扩展到【0 255】，增强对比度；
figure;
subplot(121),imshow(J);%显示原图像
subplot(122),imshow(K);%显示对比度增强的图像
[m,n]=size(J);%求图像J数据矩阵的大小赋值给[m,n]表示维矩阵；
mm=round(m/2);%对m/2取整赋值给mm
mn=round(n/2);
[p,q]=size(k);
pp=round(p/2);
qq=round(q/2);
J=double(J);%将数据变为double类型
K=double(K);
colorsum=0.0;%给灰度值总和赋0值
Javg=mean2(J);%求原图像一阶矩
Kavg=mean2(K);%求增强对比度后的图像的一阶矩
Jstd=std(std(J));%求原图像二阶矩
Kstd=std(std(K));%求增强对比度后的图像的二阶矩
for i=1:mn%循环求解灰度值总和
    for j=1:mn
        colorsum=colorsum+(J(i,j)-Javg)^3;
    end
end
Jske=(colorsum/mm*mn)^(1/3)%求原图像三阶矩
colorsum=0.0;%给灰度值总和赋0值
for i=1:pp%循环求解灰度值总和
    for j=1:qq
        colorsum=colorsum+(J(i,j)-Kavg)^3;
    end
end