clear;
%數據處理部分
[FileNamePathName]?=?uigetfile（‘*.*‘‘電離層網格產品‘）;?%選擇電離層產品文件
f?=?fopen（fullfile（PathNameFileName）?‘r‘）;?%讀取文件數據
????Data?=?fscanf（f‘%c‘）;?%以字符形式儲存于Data中
fclose（f）;
TECstart?=?strfind（Data?‘START?OF?TEC?MAP‘）;?%查找TEC?MAP開始、結束標識
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時每兩小時一個TEC?MAP
????Time=2*（i-1）;
????CurrentMap?=?Data（TECstart（i）:?TECend（i））;?%當前時間電離數據圖
????TECflag?=?strfind（CurrentMap?‘LAT/LON1/LON2/DLON/H‘）;?%電離層數據開始位置
????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??%經度-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;??%時間
????x?=?PaintData（1+5183*（i-1）:5183*i3）;??%讀取數據
????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格網繪圖
????shading?interp;??%對圖形對象的顏色著色進行色彩的插值處理使色彩平滑過渡
????colorbar;??%色標
????axis（[-200?200?-100?100?0?500]）;?%坐標軸大小
????caxis（[0500]）;??%調節顏色標尺范圍
????title（sprintf（‘%d年%d月%d日%d時VTEC分布圖（三維顯示）‘?Year?Month?Day?Time‘））;?%標題
????xlabel（‘經度‘）ylabel（‘緯度‘）zlabel（‘VTEC值‘）;%各個坐標軸標簽
????
????subplot（122）;
????pcolor（XYZ）;??%二維繪圖
????shading?interp;??%對圖形對象的顏色著色進行色彩的插值處理使色彩平滑過渡
????colorbar;??%色標
????axis?tight;??%坐標軸的最大值和最小值和的數據范圍一致
????caxis（[0500]）;??%調節顏色標尺范圍
????set（gca‘nextplot‘‘replacechildren‘）;??%設定axis覆蓋重畫模式
????title（sprintf（‘%d年%d月%d日%d時VTEC分布圖（平面顯示）‘?Year?Month?Day?Time‘））;?%標題
????xlabel（‘經度‘）ylabel（‘緯度‘）zlabel（‘VTEC值‘）;%各個坐標軸標簽
????
????M（i）=getframe（gcf）;??%保存各個時刻繪圖窗口到每一幀
????%保存至GIF文件
????im?=?frame2im（M（i））;
????[imindcm]?=?rgb2ind（im256）;
????if?i?==?1
????????imwrite（imindcmFileName?‘gif‘?‘Loopcount‘?inf）;
????else
????????imwrite（imindcmFileName?‘gif‘?‘WriteMode‘?‘append‘）;
????end
end
close?all;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????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模型建立\程序-格式轉換\1-格式轉換結果\2020-result.txt

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

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

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

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

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

?????文件????2021429??2017-10-28?11:44??VTEC模型建立\程序-格式轉換\VTEC分布圖-數據提取\2020-說明結果.txt

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

?????文件????2021429??2017-10-28?11:44??VTEC模型建立\程序-格式轉換\VTEC分布圖-數據提取\3460-說明結果.txt

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

............此處省略67個文件信息