Coherent.rar_coherent_matlab

clc close all clear all [fn pn]=uigetfile('*.tif','!!! LLAMAR IMAGEN !!!'); Uo=imread([pn fn]); N=size(Uo,1) N=480; pantalla=get(0,'ScreenSize'); xw=pantalla(3);yw=pantalla(4); fig1=figure('Position',[150 150 0.5*xw 0.6*yw]); % figure imagesc(Uo) colormap(gray(256)) title(fn) lambda=0.0005; D=1; f=20; zo=-2000 m=f/(zo+f) zi=f-m*f rA=1.22*lambda*zi/D; n=4; du_mm=rA/n; du=du_mm/(lambda*zi); dx=(lambda*zi)/(N*du_mm); % D/dx rH=round((D/2)/du_mm); t=-N/2:1:N/2-1; x=dx*t; u=du*t; [tX,tY]=meshgrid(t,t); pH=double(sqrt(tX.*tX+tY.*tY)<=rH); % figure % imagesc(pH) tfUo=fft2(Uo); TfUo=tfUo/tfUo(1,1); MfUo=abs(fftshift(TfUo)); % figure % imagesc(MfUo) Ui=ifft2(TfUo.*fftshift(pH)); Ii=abs(Ui)/max(max(abs(Ui))); fig2=figure('Position',[1050 150 0.5*xw 0.6*yw]); imagesc(Ii) colormap(gray(256))