資源簡(jiǎn)介
NOAA18 AVHRR 白天和晚上云檢測(cè),之后求得日溫度差,同時(shí)求得改天有效海溫?cái)?shù)據(jù)

代碼片段和文件信息
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_;
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-11-16?11:37??多題實(shí)驗(yàn)7\
?????文件????????2377??2014-11-16?10:45??多題實(shí)驗(yàn)7\del_cloud.m
?????文件?????????910??2014-11-16?10:59??多題實(shí)驗(yàn)7\difference_maxarea.m
?????文件??????691608??2014-11-16?11:32??多題實(shí)驗(yàn)7\introduction.docx
?????目錄???????????0??2014-11-16?11:35??多題實(shí)驗(yàn)7\process_pic\
?????文件??????125327??2014-11-15?19:33??多題實(shí)驗(yàn)7\process_pic\T3_4.jpg
?????文件???????45343??2014-11-15?19:44??多題實(shí)驗(yàn)7\process_pic\T3_4_hist.jpg
?????文件??????155606??2014-11-15?19:33??多題實(shí)驗(yàn)7\process_pic\T3_5.jpg
?????文件???????46446??2014-11-15?19:44??多題實(shí)驗(yàn)7\process_pic\T3_5_hist.jpg
?????文件???????81321??2014-11-15?19:33??多題實(shí)驗(yàn)7\process_pic\T4_5.jpg
?????文件???????43819??2014-11-15?19:44??多題實(shí)驗(yàn)7\process_pic\T4_5_hist.jpg
?????文件?????1940491??2014-11-15?15:59??多題實(shí)驗(yàn)7\process_pic\ch3.tif
?????文件???????51598??2014-11-15?16:19??多題實(shí)驗(yàn)7\process_pic\ch3_hist.jpg
?????文件?????1940867??2014-11-15?15:59??多題實(shí)驗(yàn)7\process_pic\ch4.tif
?????文件???????51660??2014-11-15?16:19??多題實(shí)驗(yàn)7\process_pic\ch4_hist.jpg
?????文件?????1941195??2014-11-15?15:59??多題實(shí)驗(yàn)7\process_pic\ch5.tif
?????文件???????49486??2014-11-15?16:19??多題實(shí)驗(yàn)7\process_pic\ch5_hist.jpg
?????目錄???????????0??2014-11-16?11:31??多題實(shí)驗(yàn)7\result_pic\
?????文件??????121935??2014-11-16?10:47??多題實(shí)驗(yàn)7\result_pic\Diff.jpg
?????文件??????221013??2014-11-16?10:47??多題實(shí)驗(yàn)7\result_pic\Max_area.jpg
?????文件??????180220??2014-11-16?11:14??多題實(shí)驗(yàn)7\result_pic\N18D_sst.jpg
?????文件??????155733??2014-11-16?11:15??多題實(shí)驗(yàn)7\result_pic\N19N_sst.jpg
?????文件??????146334??2014-11-16?11:29??多題實(shí)驗(yàn)7\result_pic\eliminate.jpg
?????文件??????860812??2014-11-15?13:19??多題實(shí)驗(yàn)7\基于動(dòng)態(tài)閾值法的NOAA系列衛(wèi)星云檢測(cè)結(jié)果檢驗(yàn)_武艷?(1).caj
評(píng)論
共有 條評(píng)論