A = imread('fig1c.jpg'); Green = A(:,:,2); s = size(Green); FF = zeros(s); for i=1:s(1) FF(i,:) = fft(Green(i,:)); end unit = 50000/1086; y = 0:unit/1000:unit*(s(1)-1)/1000; x = 0:1/unit:1/unit*(s(2)-1); surf(x,y,abs((FF))); axis([0.15 0.62 0 2]); view(0,90); shading flat; colorbar; caxis([0 5000]);