1

エッジ指向のヒストグラムを実装しようとしています: しかし、コードに問題があり、 val=fix((angles(i,j)-a_min)/(BinSize))+1;Nanがどこにあるのかわかりません

私は8つのビンを取ります。1.最初に、ソーベル演算子を使用して勾配と角度を見つけています。2.後で、8つのビン（-pi/2からpi/2、atanはこれら2つの値の間の値を返すため）を持っており、角度に従って勾配の大きさをビニングしています。

    im=imread('cameraman.tif');
subplot(2,2,1);imshow(im);
im=im2double(im);

sob_x=[-1,0,1;-2,0,2;-1,0,1];
sob_y=[1,2,1;0,0,0;-1,-2,-1];

im_x=imfilter(im,sob_x);
im_y=imfilter(im,sob_y);
subplot(2,2,2);imshow(im_x);
subplot(2,2,3);imshow(im_y);

im_edge=zeros(size(im));
angles=zeros(size(im));

[r c]=size(im);
for i=1:r
    for j=1:c
        im_edge(i,j)=sqrt(im_x(i,j)^2+im_y(i,j)^2);
        angles(i,j)=atan(im_y(i,j)/im_x(i,j));%angles are -pi/2 to pi/2
    end
end

subplot(2,2,4);imshow(im_edge,[]);

im_temp=im2bw(im_edge,80/255);
figure();subplot(1,2,1);imshow(im_temp);
subplot(1,2,2);imshow(angles,[]);

%calculating the histogram



No_bins=8;
H=zeros(1,No_bins);
BinSize=pi/No_bins;


a_min=-pi/2;
for i=1:r
    for j=1:c
        val=fix((angles(i,j)-a_min)/(BinSize))+1;
        if (val>8)
            val=8;
        end
        if(val<1)
            val=1;
        end

        H(1,val)=H(1,val)+im_edge(i,j);


end
end

H
4

1 に答える 1

0

一部のピクセルでのこの angles(i,j)=atan(im_y(i,j)/im_x(i,j));%angles are -pi/2 to pi/2戻り値 Nan 値の実装。エッジ ピクセルでエッジの向きを抽出してみてください。

于 2015-01-13T13:48:26.100 に答える