[問題] 目前的問題已回收

看板MATLAB作者soiialan (雄哥)時間16年前 (2009/03/17 16:09)推噓0(0推 0噓 2→)

留言2則, 2人參與討論串1/1

要怎麼把這個錯改成可以跑的謝謝 ??? Error using ==> times Matrix dimensions must agree. Error in ==> fft_clocy at 53 cf=af1.*c; 以下是原程式碼 clc, clear, close all warning off clc; close all; % clear immean % img1=[zeros(60,30) ones(60,30)]; %fmain='5-sf=1-RPCScale=1.5'; b=5; %讀取影像並轉成灰階影像 [filename, pathname] = uigetfile({'*.tif';'*.jpg';'*.bmp'},'Select an image'); img = imread([pathname filename]); if ndims(img) == 3 % RGB image F_spatial = rgb2gray(img); elseif ndims(img) == 2 % gray image F_spatial = img; else disp('The image is not RGB or gray image!') end subplot(2,2,1), imshow(F_spatial), title('image in spatial domain') set(gcf, 'name', '計算影像的傅利葉頻譜並換算cpd','color','w') %img1=imread(strcat('E:\hdr\20090312\20090112-25張圖 \',fmain,'-gamma=0.33333.tif')); %img1=imread(strcat('E:\fft2spectra\20081229higher contrast exprement\1\1.tif')); %img1=srgb2lab(double(img1)); % figure,imshow(uint8(img1(:,:,1))),title(原圖LAB color) af1=fftshift(fft2(img(:,:,1))); figure,fftshow(af1,'log'),title(strcat('fftshow log')) F_size = size(F_spatial); %影像大小 aspect_ratio = F_size(1)/F_size(2); %影像高寬比 center_x = round(F_size(1)/2+1); % 頻率0的位置 center_y = round(F_size(2)/2); % 頻率0的位置 max_radius = ceil(sqrt((center_x/aspect_ratio)^2+ center_y^2))+1; %最大半徑(x 軸按高寬比縮放) meanP1d = zeros(max_radius,2); %一維log傅利葉頻譜平均值初始化 maxP1d = zeros(max_radius,1); %一維log傅利葉頻譜最大值初始化 a=1; for ii=1:30:301 % [x,y]=meshgrid(-300:299,-240:239); % [x,y]=meshgrid(-300:299,-225:224); [x,y]=meshgrid(-300:299,-300:299); z=sqrt(x.^2+y.^2); c1=(z<ii-30); c2=(z<ii); c=xor(c1,c2); figure,imshow(c); cf=af1.*c; % figure,fftshow(cf,'log'),title(strcat('fftshow log')) img2=img; img2(:,:,1)=ifft2(ifftshift(cf)); % immean(a,1)=mean(mean(img2(:,:,1))); % img2=(img2./max(img2(:)))*255; % img21=uint8((img2)); img21=lab2srgb(img2); immean(a,b)=mean(img21(:)); a=a+1; figure,imshow(img21) imwrite(img21,strcat('E:\fft2spectra\20081229higher contrast exprement\1\',' i=',num2str(ii),'.tif'),'tif') close all; end -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.137.51.66