• 大小: 3.05MB
    文件类型: .zip
    金币: 1
    下载: 0 次
    发布日期: 2023-09-20
  • 语言: 其他
  • 标签: avhrr  NOAA  云检测  

资源简介

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

评论

共有 条评论