name='d:/';
imgname=[name,'01_test'];
cdata=imread(imgname,'bmp');
cdata = double(cdata);
Size = 17;
step = [1 1.5 2 2.5 3 3.5 4];
%step = [0.1 0.3 0.5 1 1.5 2 3];
for is = 1:6
sigma = step(is);
period = step(is)*2.0*sqrt(2);
orient = 0;
[height, width] = size(cdata);
bcAcu = -Inf*ones(height, width);
Acu = zeros(height, width);
bsAcu = -Inf*ones(height, width);
for i = 0:20:360
[cmask, smask] = gabormask(Size, sigma, period, (orient+i)*pi/180);
bc = imfilter(cdata,cmask);
%bcAcu = bcAcu + bc;
bcAcu = max(bcAcu, bc);
% figure,imshow(bc,[]);
bs = imfilter(cdata,smask);
bsAcu = max(bsAcu, bs);
%figure,imshow(nnormal(bs),[]);
end
% figure,imshow(nnormal(bcAcu),[]);
% figure,imshow(nnormal(bsAcu),[]);
figure,imshow(bcAcu,[]);
% figure,imshow(bsAcu,[]);
fr = uint8(nnormal(bcAcu,10,50));%uint8(nnormal(1*(bcAcu)+1.5*bsAcu/step(is)));
frs = uint8(nnormal(bsAcu,min(min(bsAcu)),max(max(bsAcu))));
str = [name,'.IM0.dcm_0013.gabor',num2str(is),'.bmp'];
%imwrite(fr,str,'bmp');
%Acu = max(bcAcu,Acu);
%Acu = max(1.1*Acu,nnormal(bcAcu,10,50))
??? Attempt to reference field of non-structure array.
Error in ==> imagesci\private\readbmp at 17
map = info.Colormap;
Error in ==> imread at 306
[X, map] = feval(fmt_s.read, filename, extraArgs{:});

up

真心求教

能问下你的程序还在吗？？