清晰的所有,关闭所有;
2 = imread (“BMB2.jpg”);
图;imshow (Ii);
标题(“对象模式”)
轴从
Ii0 =双(2);
为rgb = 1:3;
2 = Ii0 (:,:, rgb);
PH =兰德([256256]);
2 = Ii。* exp(2 *π* PH值);
M = 512;
我= 0 (512);
我(128:383,128:383)= 2;
z = 15;
w = 6500 * 10 ^ 8;
δ= 0.005;
r = 1: M;
c = 1: M;
[C R] = meshgrid (C, R);
p = exp(2我*π* z。* ((1 / w) ^ 2 - (1 / M /δ)^ 2 *(决定物价/ 2 - 1)。^ 2 - (1 / M /δ)^ 2 *(- / 2 - 1)。^ 2)。^ 0.5);
A0 = fftshift (ifft2 (fftshift(我)));
阿兹= A0。* p;
E = fftshift (fft2 (fftshift (Az)));
p = exp(2我*π* (- z)。* ((1 / w) ^ 2 - (1 / M /δ)^ 2 *(决定物价/ 2 - 1)。^ 2 - (1 / M /δ)^ 2 *(- / 2 - 1)。^ 2)。^ 0.5);
A1 = fftshift (ifft2 (fftshift (E)));
Az1 = A1。* p;
R1 = fftshift (fft2 (fftshift (Az1)));
R1 = (abs (R1)) ^ 2;
图;imshow (R1 / max (max (R1)));
标题(“重建图像(650海里)”)
轴从
dw = 50;
IMA = 0 (512512);
为g =今日;
w2 = (20000 - dw * g) * 10 ^ 8;
E2 = e * exp(我*π*信德(10)* (w-w2) / w / w2。* R *δ);
p = exp(2我*π* (- z)。* ((1 / w2) ^ 2 - (1 / M /δ)^ 2 *(决定物价/ 2 - 1)。^ 2 - (1 / M /δ)^ 2 *(- / 2 - 1)。^ 2)。^ 0.5);
Az2 = ifft2 (fftshift (E2)) * (fftshift (p));
R2 = fftshift (fft2 (Az2));
R = (abs (R2)) ^ 2;
IMA = IMA + R2;
结束
IMA = IMA / max (max (IMA));
图;imshow (IMA)
标题(“重建图像(白光)”)
轴从
结束