【Matlab】繪製 MSTID 隨時間變化圖

tags: Matlab 高層大氣科學

Common 目標

1. 繪製2018年6月2日1200UT~1300UT
2. 每5分鐘之MSTID變化圖
3. 共13張圖，每張圖上標示日期與時間

Introduction 簡介

1. TID：Traveling Ionospheric Disturbances

2. TID 分類

   

3. MSTID

   

Matlab Code 程式碼

1. 繪製MSTID所需資料
   • 2018153.tar.gz
   • asian.mat
2. Code

clear all
close all

fndir='2018.153';
hmdir=pwd;
cd(fndir)
list=dir('v*.dat');
cd(hmdir)

% 重複繪製13次
hr=12;
for t=0:12    
    load asian.mat
    cd(fndir)
    mm=t*5;
    if t==12
        hr=13;
        mm=0;
    end
    % `~exist` means file not exist
    % `%02d` means two numbers, if not add 0
    if ~exist(['MSTID_' num2str(hr) num2str(mm,'%02d') 'UT.mat'])
        tindex1=(hr-1)*120+mm*2+1;
        tindex2=hr*120+mm*2+1;
        vtec=[];
        lons=[];
        lats=[];
        eles=[];
        for i=1:length(list)
            fn=list(i).name;
            tmp=load(['v' fn(2:end)]);
            tmp2=tmp(tindex1:tindex2,:);
            f=find(tmp2==0);
            tmp2(f)=nan;
            vtec=[vtec tmp2(end,:)-nanmean(tmp2)];
            tmp=load(['o' fn(2:end)]);
            lons=[lons tmp(tindex2,:)];
            tmp=load(['a' fn(2:end)]);
            lats=[lats tmp(tindex2,:)];
            tmp=load(['e' fn(2:end)]);
            eles=[eles tmp(tindex2,:)];
            %disp([num2str(i) '/' num2str(length(list))])
        end
        % eles<40 data delete
        f=find(lons==0 | lats==0 | eles<40);
        vtec(f)=[];
        lons(f)=[];
        lats(f)=[];
        eles(f)=[];
        eval(['save MSTID_' num2str(hr) num2str(mm,'%02d') 'UT.mat vtec lons lats eles'])
    else
        load(['MSTID_' num2str(hr) num2str(mm,'%02d') 'UT.mat'])
    end
    
    close all
    plot(as.lon,as.lat,'k');hold on
    axis equal
    axis([115 150 15 50])
    scatter(lons,lats,10,vtec,'filled');shading flat
    caxis([-2 2])
    title(['2018/06/02 ' num2str(hr) num2str(mm,'%02d') 'UT'])
    colorbar
    colormap(jet)
    
    cd(hmdir)
    print('-dpng', ['MSTID_' num2str(hr) num2str(mm,'%02d') 'UT.png'])
end

Result 結果

Reference 參考資料

1. 高層大氣科學 上課講義