clear;
%數(shù)據(jù)處理部分
[FileNamePathName]?=?uigetfile（‘*.*‘‘電離層網(wǎng)格產(chǎn)品‘）;?%選擇電離層產(chǎn)品文件
f?=?fopen（fullfile（PathNameFileName）?‘r‘）;?%讀取文件數(shù)據(jù)
????Data?=?fscanf（f‘%c‘）;?%以字符形式儲存于Data中
fclose（f）;
TECstart?=?strfind（Data?‘START?OF?TEC?MAP‘）;?%查找TEC?MAP開始、結(jié)束標(biāo)識
TECend?=?strfind（Data?‘END?OF?TEC?MAP‘）;?
Dateflag?=?strfind（Data?‘?EPOCH?OF?FIRST?MAP?‘）;??%提取日期
Date?=?str2num（strtrim（Data（Dateflag-60:?Dateflag）））;
Year?=?Date（1）;?Month?=?Date（2）;?Day?=?Date（3）;?%月年日
TEC（111）?=?0;
flag?=?0;
PaintData=zeros（14）;
for?i?=?1:13??%0-24時(shí)每兩小時(shí)一個(gè)TEC?MAP
????Time=2*（i-1）;
????CurrentMap?=?Data（TECstart（i）:?TECend（i））;?%當(dāng)前時(shí)間電離數(shù)據(jù)圖
????TECflag?=?strfind（CurrentMap?‘LAT/LON1/LON2/DLON/H‘）;?%電離層數(shù)據(jù)開始位置
????for?j?=?1:71??%緯度87.5至-87.5間隔2.5
????????LAT=87.5-2.5*（j-1）;
????????TECdata?=?str2num（char（strsplit（strtrim（CurrentMap（TECflag（j）+20:?TECflag（j）+390）））））;
????????for?k?=?1:73??%經(jīng)度-180.0至180.0間隔5
????????????LON=-180+5*（k-1）;
????????????%TEC（ijk）?=?TECdata（k1）;
????????????flag?=?flag?+?1;
????????????PaintData（flag1:4）?=?[Time?LAT?LON?TECdata（k1）];
????????end
????end
end

%繪圖部分
ScreenSize?=?get（0?‘ScreenSize‘）;??%獲得屏幕大小
set（gcf?‘Position‘?[ScreenSize（3）/2-650?ScreenSize（4）/2-200?1300?400]）;??%繪圖窗口居中
FileName?=?strcat（fullfile（PathName?strtok（FileName?‘.‘））?‘.gif‘）;
for?i=1:13
????Time?=?（i-1）?*?2;??%時(shí)間
????x?=?PaintData（1+5183*（i-1）:5183*i3）;??%讀取數(shù)據(jù)
????y?=?PaintData（1+5183*（i-1）:5183*i2）;
????z?=?PaintData（1+5183*（i-1）:5183*i4）;
????[XYZ]?=?griddata（xyzlinspace（-180180100）‘linspace（-87.587.5100）‘v4‘）;?%插值處理

????subplot（121）;
????mesh（XYZ）;??%mesh格網(wǎng)繪圖
????shading?interp;??%對圖形對象的顏色著色進(jìn)行色彩的插值處理使色彩平滑過渡
????colorbar;??%色標(biāo)
????axis（[-200?200?-100?100?0?500]）;?%坐標(biāo)軸大小
????caxis（[0500]）;??%調(diào)節(jié)顏色標(biāo)尺范圍
????title（sprintf（‘%d年%d月%d日%d時(shí)VTEC分布圖（三維顯示）‘?Year?Month?Day?Time‘））;?%標(biāo)題
????xlabel（‘經(jīng)度‘）ylabel（‘緯度‘）zlabel（‘VTEC值‘）;%各個(gè)坐標(biāo)軸標(biāo)簽
????
????subplot（122）;
????pcolor（XYZ）;??%二維繪圖
????shading?interp;??%對圖形對象的顏色著色進(jìn)行色彩的插值處理使色彩平滑過渡
????colorbar;??%色標(biāo)
????axis?tight;??%坐標(biāo)軸的最大值和最小值和的數(shù)據(jù)范圍一致
????caxis（[0500]）;??%調(diào)節(jié)顏色標(biāo)尺范圍
????set（gca‘nextplot‘‘replacechildren‘）;??%設(shè)定axis覆蓋重畫模式
????title（sprintf（‘%d年%d月%d日%d時(shí)VTEC分布圖（平面顯示）‘?Year?Month?Day?Time‘））;?%標(biāo)題
????xlabel（‘經(jīng)度‘）ylabel（‘緯度‘）zlabel（‘VTEC值‘）;%各個(gè)坐標(biāo)軸標(biāo)簽
????
????M（i）=getframe（gcf）;??%保存各個(gè)時(shí)刻繪圖窗口到每一幀
????%保存至GIF文件
????im?=?frame2im（M（i））;
????[imindcm]?=?rgb2ind（im256）;
????if?i?==?1
????????imwrite（imindcmFileName?‘gif‘?‘Loopcount‘?inf）;
????else
????????imwrite（imindcmFileName?‘gif‘?‘WriteMode‘?‘a(chǎn)ppend‘）;
????end
end
close?all;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????809376??2017-10-31?23:47??VTEC模型建立\VTEC.exe

?????文件???????2907??2017-10-31?23:50??VTEC模型建立\VTEC.m

?????文件????2021370??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\1-格式轉(zhuǎn)換結(jié)果\2020-result.txt

?????文件????2021429??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\1-格式轉(zhuǎn)換結(jié)果\2020-說明結(jié)果.txt

?????文件????2021370??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\1-格式轉(zhuǎn)換結(jié)果\3460-result.txt

?????文件????2021429??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\1-格式轉(zhuǎn)換結(jié)果\3460-說明結(jié)果.txt

?????文件???????1482??2017-10-28?14:34??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\1.f90

?????文件????2021370??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\2020-result.txt

?????文件????2021429??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\2020-說明結(jié)果.txt

?????文件????2021370??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\3460-result.txt

?????文件????2021429??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\3460-說明結(jié)果.txt

?????文件??????11649??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\Debug\1.obj

?????文件???????3934??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\Debug\BuildLog.htm

?????文件??????53248??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\Debug\vc100.pdb

?????文件?????913920??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\Debug\VTEC分布圖-數(shù)據(jù)提取.exe

?????文件????????381??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\Debug\VTEC分布圖-數(shù)據(jù)提取.exe.intermediate.manifest

?????文件????2165760??2017-10-28?11:44??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\Debug\VTEC分布圖-數(shù)據(jù)提取.pdb

?????文件?????869995??2017-10-28?02:58??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\igsg2020.16i

?????文件?????869428??2017-10-28?00:34??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\igsg3460.16i

????..A..H.???????250??2017-10-30?20:43??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\VTEC分布圖-數(shù)據(jù)提取.u2d

?????文件???????1967??2017-10-26?20:36??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取\VTEC分布圖-數(shù)據(jù)提取.vfproj

?????文件????????944??2017-10-26?20:25??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取.sln

????..A..H.?????13312??2017-10-30?20:43??VTEC模型建立\程序-格式轉(zhuǎn)換\VTEC分布圖-數(shù)據(jù)提取.suo

?????文件???????1373??2017-10-28?03:09??VTEC模型建立\程序-靜態(tài)動態(tài)繪圖\2020三維gif.m

?????文件???????1377??2017-10-28?03:09??VTEC模型建立\程序-靜態(tài)動態(tài)繪圖\2020平面gif.m

?????文件???????1373??2017-10-28?02:47??VTEC模型建立\程序-靜態(tài)動態(tài)繪圖\3460三維gif.m

?????文件???????1377??2017-10-28?02:48??VTEC模型建立\程序-靜態(tài)動態(tài)繪圖\3460平面gif.m

?????文件?????507267??2017-10-28?16:56??VTEC模型建立\結(jié)果-VTEC分布圖\2016年12月12日\3460三維-gif.gif

?????文件?????295630??2017-10-28?16:53??VTEC模型建立\結(jié)果-VTEC分布圖\2016年12月12日\3460平面-gif.gif

?????文件?????123618??2017-10-28?17:03??VTEC模型建立\結(jié)果-VTEC分布圖\2016年12月12日\靜態(tài)\3460三維-1.jpg

............此處省略67個(gè)文件信息