Question

您好我是Matlab的新手......我正在尝试不使用histeq进行直方图均衡...

但由于某种原因，我总是得到一个错误：???指数超过矩阵维度。

这是我的代码........................................... .......感谢您的帮助

  clc

  I = imread ('Machine-Edge.PNG');
  I2 = rgb2gray(I);
  colormap gray;

  y = imhist(I2);


  %using hist eq. built in fn
  I3= histeq(I2);
  z= imhist(I3);

  %my equalization
  r = size(I2,1);
  c = size(I2,2);
  A= zeros(1,256);

  %counting number of pixels of the image and putting the count in Array A
  for j=1:r
    for x=1:c
        v=I2(j,x);
        A(v+1)=A(v+1)+1;
    end
  end

  %pi=n/size
    for y=1;256
      pi(y)= ((A(y))/(r*c));
    end

   %calculate CI (cumulated pi )
    ci(1)=pi(1);
    for yy=2;256
         ci(yy) = ci(yy-1)+ pi(yy);
    end

    %calculate T=range *Ci
      for b=1;256
        T(b)=ci(b)*255;
       end

    %equilization..replacing each pixel with T value
     for j=1:r
         for x=1:c
             I4(j,x) =T(I2(j,x));

          end
     end




    vv= imhist(I4);



   figure
   subplot(3,2,1)
    imagesc(I2)
   subplot(3,2,2)
   plot(y)

  subplot(3,2,3)
  imagesc(I3)
   subplot(3,2,4)
   plot(z)

   subplot(3,2,5)
   imagesc(I4)
   subplot(3,2,6)
    plot(vv)

Answer 1

这是一个老帖子，但OP使用;而不是：在它们的for循环中（即对于y = 1; 256应该读取y = 1：256）。更正后的代码如下：

 clc

 I = imread ('Machine-Edge.PNG');
 I2 = rgb2gray(I);
 colormap gray;

 y = imhist(I2);


 %using hist eq. built in fn
 I3= histeq(I2);
 z= imhist(I3);

 %my equalization
 r = size(I2,1);
 c = size(I2,2);
 A= zeros(1,256);

 %counting number of pixels of the image and putting the count in Array A
 for j=1:r
   for x=1:c
       v=I2(j,x);
       A(v+1)=A(v+1)+1;
   end
 end

 %pi=n/size
   for y=1:256
     pi(y)= ((A(y))/(r*c));
   end

  %calculate CI (cumulated pi )
   ci(1)=pi(1);
   for yy=2:256
        ci(yy) = ci(yy-1)+ pi(yy);
   end

   %calculate T=range *Ci
     for b=1:256
       T(b)=ci(b)*255;
      end

   %equilization..replacing each pixel with T value
    for j=1:r
        for x=1:c
            I4(j,x) =T(I2(j,x));

         end
    end




   vv= imhist(I4);



  figure
  subplot(3,2,1)
   imagesc(I2)
  subplot(3,2,2)
  plot(y)

 subplot(3,2,3)
 imagesc(I3)
  subplot(3,2,4)
  plot(z)

  subplot(3,2,5)
  imagesc(I4)
  subplot(3,2,6)
   plot(vv)

Matlab直方图均衡，无需使用内置函数

1 个答案: