AVHRR数据云检测

NOAA18 AVHRR 白天和晚上云检测，之后求得日温度差，同时求得改天有效海温数据

clear;clc;
ch1=load（‘sst1.txt‘）;
ch2=load（‘sst2.txt‘）;
ch3=load（‘sst3.txt‘）;
ch4=load（‘sst4.txt‘）;
ch5=load（‘sst5.txt‘）;
sst=load（‘sst.txt‘）;
save ch1;save ch2;save ch3;save ch4;save ch5;save sst;

ch1_=reshape（ch140004000）;
ch2_=reshape（ch240004000）;
ch3_=reshape（ch340004000）;
ch4_=reshape（ch440004000）;
ch5_=reshape（ch540004000）;
sst_=reshape（sst40004000）;
sst_（sst_>30 | sst_<0）=NaN;
save ch1_;save ch2_;save ch3_; save ch4_;save ch5_;save sst_;

figure（1）;imagesc（ch1_‘）;colorbar;caxis（[0 30]）; 
figure（2）;imagesc（ch2_‘）;colorbar;caxis（[0 30]）;
figure（3）;imagesc（ch3_‘）;colorbar;caxis（[0 30]）;
figure（4）;imagesc（ch4_‘）;colorbar;caxis（[0 30]）;
figure（5）;imagesc（ch5_‘）;colorbar;caxis（[0 30]）;
figure（15）;imagesc（sst_‘）;colorbar;caxis（[0 30]）;

figure（6）;hist（ch3_）;
figure（7）;hist（ch4_）;
figure（8）;hist（ch5_）;
figure（16）;hist（sst_）;

saveas（figure（1）‘ch1‘‘tif‘）;
saveas（figure（2）‘ch2‘‘tif‘）;
saveas（figure（3）‘ch3‘‘tif‘）;
saveas（figure（4）‘ch4‘‘tif‘）;
saveas（figure（5）‘ch5‘‘tif‘）;
saveas（figure（15）‘sst‘‘tif‘）;

saveas（figure（6）‘ch3_hist‘‘jpg‘）;
saveas（figure（7）‘ch4_hist‘‘jpg‘）;
saveas（figure（8）‘ch5_hist‘‘jpg‘）;
saveas（figure（16）‘sst_hist‘‘jpg‘）;


%两通道比较
T3_4=ch3_-ch4_;
T3_5=ch3_-ch5_;
T4_5=ch4_-ch5_;

figure（9）;imagesc（T3_4‘）;colorbar;
figure（10）;imagesc（T3_5‘）;colorbar;
figure（11）;imagesc（T4_5‘）;colorbar;
figure（12）;hist（T3_4）;%直方图
figure（13）;hist（T3_5）;
figure（14）;hist（T4_5）;
saveas（figure（9）‘T3_4‘‘jpg‘）;
saveas（figure（10）‘T3_5‘‘jpg‘）;
saveas（figure（11）‘T4_5‘‘jpg‘）;
saveas（figure（12）‘T3_4_hist‘‘jpg‘）;
saveas（figure（13）‘T3_5_hist‘‘jpg‘）;
saveas（figure（14）‘T4_5_hist‘‘jpg‘）;

cloud2=zeros（40004000）;%局部方差
for i=2:3999
    for j=2:3999
        temp=[ch4_（i-1j-1）ch4_（i-1j）ch4_（i-1j+1）...
            ch4_（ij-1）ch4_（ij）ch4_（ij+1）...
            ch4_（i+1j-1）ch4_（i+1j）ch4_（i+1j+1）];
        cloud2（ij）=nanvar（temp）;
    end
end
figure（17）;imagesc（cloud2‘）;colorbar;caxis（[0 0.5]）;
figure（18）;hist（cloud2）;
saveas（figure（18）‘cloud2_hist‘‘jpg‘）;
save cloud2;

cloud=zeros（40004000）;%云
cloud（T3_4>10）=1;
cloud（cloud2>0.5）=1;
figure（20）;imagesc（cloud‘）;colorbar;


%去云
index1=find（cloud2>0.3）;
index2=find（T3_4>10）;
sst_（index1）=NaN;
sst_（index2）=NaN;
figure（21）;imagesc（sst_‘）;colorbar;
save sst_;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-11-16 11:37  多题实验7\
     文件        2377  2014-11-16 10:45  多题实验7\del_cloud.m
     文件         910  2014-11-16 10:59  多题实验7\difference_maxarea.m
     文件      691608  2014-11-16 11:32  多题实验7\introduction.docx
     目录           0  2014-11-16 11:35  多题实验7\process_pic\
     文件      125327  2014-11-15 19:33  多题实验7\process_pic\T3_4.jpg
     文件       45343  2014-11-15 19:44  多题实验7\process_pic\T3_4_hist.jpg
     文件      155606  2014-11-15 19:33  多题实验7\process_pic\T3_5.jpg
     文件       46446  2014-11-15 19:44  多题实验7\process_pic\T3_5_hist.jpg
     文件       81321  2014-11-15 19:33  多题实验7\process_pic\T4_5.jpg
     文件       43819  2014-11-15 19:44  多题实验7\process_pic\T4_5_hist.jpg
     文件     1940491  2014-11-15 15:59  多题实验7\process_pic\ch3.tif
     文件       51598  2014-11-15 16:19  多题实验7\process_pic\ch3_hist.jpg
     文件     1940867  2014-11-15 15:59  多题实验7\process_pic\ch4.tif
     文件       51660  2014-11-15 16:19  多题实验7\process_pic\ch4_hist.jpg
     文件     1941195  2014-11-15 15:59  多题实验7\process_pic\ch5.tif
     文件       49486  2014-11-15 16:19  多题实验7\process_pic\ch5_hist.jpg
     目录           0  2014-11-16 11:31  多题实验7\result_pic\
     文件      121935  2014-11-16 10:47  多题实验7\result_pic\Diff.jpg
     文件      221013  2014-11-16 10:47  多题实验7\result_pic\Max_area.jpg
     文件      180220  2014-11-16 11:14  多题实验7\result_pic\N18D_sst.jpg
     文件      155733  2014-11-16 11:15  多题实验7\result_pic\N19N_sst.jpg
     文件      146334  2014-11-16 11:29  多题实验7\result_pic\eliminate.jpg
     文件      860812  2014-11-15 13:19  多题实验7\基于动态阈值法的NOAA系列卫星云检测结果检验_武艳 （1）.caj