怎么得到一副sar复数图像的幅度图像?

clc;
clear;
fid=fopen('C:\','r','ieee-be') %文件路径名字 读写方式 字节序
for i=1:512
for j=1:512
I(i,j)=fread(fid,1,'int16'); %实数
Q(i,j)=fread(fid,1,'int16'); %虚数
A(i,j)=(Q(i,j)^2+I(i,j)^2)^0.5; %计算振幅
end
end
fclose(fid);
figure;imshow(A,(0,1000)); %显示振幅强度图像