clc
clear
fid1=fopen（‘‘）;???????%讀取相干系數
Coh?=?fread（fid1[1000700]‘uint8‘）;

Coh=Coh‘;%轉置
fclose（fid1）;
fid1=fopen（‘‘）;%讀取幅度值
Amp=fread（fid1[1000700]‘uint32‘）;
Amp=Amp‘;
fclose（fid1）;
fid1=fopen（‘‘）;%讀取相位
Pha=?fread（fid1[1000700]‘uint8‘）;
Pha=Pha‘;
fclose（fid1）;
imagesc（Pha）;
%row=2640;col=2450*2;
%fid1=fopen（‘G:\畢業設計\導師資料\INSAR入門資料\濾波程序及數據\程序及數據\master_slave.flt‘）;
%[fltcount]=fread（fid1[49002640]‘float32‘‘ieee-be‘）;
%fclose（fid1）;
%flt=flt‘;
%flt_real=flt（:1:2:end-1）;
%flt_image=flt（:2:2:end）;
%IFG_shanghai=flt_real+i*flt_image;
%Interferogram=IFG_shanghai;%干涉圖
%A1=angle（IFG_shanghai）;%求相位角
%Pha=A1
%enta?coherence?data?from?gamma

%fid1=fopen（‘G:\畢業設計\導師資料\INSAR入門資料\濾波程序及數據\程序及數據\master_slave.cc‘）;
%[Cohcount]=fread（fid1[24502640]‘float32‘‘ieee-be‘）;
%fclose（fid1）;
%Coh=Coh‘;



%b=max（max（Pha））;
%a=min（min（Pha））;
%c=Pha（1:7001:1000）
%imagesc（Pha）
Pha=Pha*pi*2/255-pi;

%imagesc（Pha）
Coh=Coh/255;
Interferogram=Amp.*exp（i*Pha）;?
[LP]=size（Interferogram）;
%REAL=Amp.*cos（Pha）
%imagesc（REAL）
%IMAGE=Amp.*sin（Pha）
%figure（2）
%imagesc（IMAGE）
%figure（3）
%sumf=REAL.*IMAGE
%imagesc（Pha）
%figure（4）
%imagesc（sumf）
Block=32;
Overlap=28;
increment=Block-Overlap;
Start_point=Overlap/2+1;%從15開始
Start_point_row=2;%開始行（從第二行開始）
Start_point_col=2;%開始列
End_point_row=L-?（Block-Overlap/2）;%終止行（倒數第二行）要加括號嗎
End_point_col=P-?（Block-Overlap/2）;%終止列（倒數第二列）
grid=zeros（size（Interferogram））;%格網值賦0

for?ii=Start_point:increment:End_point_row%ii從15開始
?????for?jj=Start_point:increment:End_point_col
?????????fetch_int=Interferogram（（ii-14）:（ii-14）+31（jj-14）:（jj-14）+31）;%取整個圖塊的干涉值
?????????fetch_coh=Coh（（ii-14）:（ii-14）+31（jj-14）:（jj-14）+31）;
?????????namda1=mean2（fetch_coh（15:1815:18））;%對中間位置相干系數兩次取平均（結果為一個數）
?????????namda2=namda1/100;%化為百分數
?????????H?=?fft2（fetch_int）;????	%與下一行一起組成FFT變換		
?????????H?=?fftshift（H）;
?????????S=abs（H）;
?????????S=S/max（max（S））;%將化為百分比（比上最大值）
?????????S=S.^（1-namda1）;%對應位置相乘得平滑值
?????????H=H.*S;
?????????H?=?ifftshift（H）;???%與下一行一起組成反FFT變換		???
?????????fetch_int?=?ifft2（H）;???
????????