1

誰かがこれを解決する方法を教えてもらえますか？

grayImage=dicomread('028tm.R.dcm');
sx=32;
sy=32;
theta=[0 pi/4 2*pi/4 3*pi/4 4*pi/4 5*pi/4 6*pi/4 7*pi/4];

gamma=1;
psi=0;
sigma=6; % ????12
lambda=[5 6 7 8 9 7];

# Creating 40 Gabor Filters
G = cell(5,8);
for i = 1:5
  for j = 1:8
    G{i,j}=zeros(65,65);
  end
end
for i = 1:5
  for j = 1:8
    f=1/lambda(i);      

    [T,gabout] = gaborfilter1(grayImage,sx,sy,f,theta(j));
    u{i,j}=mean2(gabout);
    del{i,j}=std2(gabout);  
    skew{i,j} =skewness(gabout);
    sk=[skew{:,:}];
    std=[u{:,:}];
    mn=[del{:,:}];
    disp(sk)             
  end
end

これが私のワークスペースを表示します。

mn <1x40 double>
sk<1x10280 double>
std<1x40 double>

なぜsk（歪度の値）が1X102802倍になるのですか????? それは1x40標準偏差と平均のようでなければなりません。sk（歪度）をどのように作成します1x40か？

誰かが私を助けてくれることを願っています。

4

1 に答える 1

2

あなたの問題はそれでmean2ありstd2、画像全体を操作しますが、画像の列skewnessを操作します。新しい関数を定義する必要があります

function sk = skewness2(X)
    sk = skewness(X(:));
end

そして、の代わりにそれを使用してskewnessください。

于 2013-03-15T11:35:30.337 に答える