资源简介
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
评论
共有 条评论