全球电离层图GIM（IONEX)数据读取matlab代码

克里普crirp

已于 2024-05-25 17:49:59 修改

阅读量1.2k

点赞数 3

文章标签：人工智能

于 2024-05-25 17:48:18 首次发布

本文链接：https://blog.csdn.net/weixin_44392003/article/details/139200887

版权

全球电离层图目前应用非常广泛，目前提供该数据产品的单位包括CODE、CAS、UPC、WHU等等。国内推荐中科院空天院李子申、王宁波团队的产品：ftp://ftp.gipp.org.cn/product/，包括IONEX、BDGIM、ION文件、DCB文件等。

下面提供IONEX数据读取的matlab代码，采用正则式进行读取，非常简洁高效。

function [ION,TimeTk] = readionexfile(ionfile)
str = fileread(ionfile);
ca1 = regexp(str, '(?<=START OF TEC MAP).+?(?=END OF TEC MAP)', 'match');
ca1 = ca1(1:end-1);
ION = nan( 73, 71, length(ca1) );
lat2ix  = @(lat) round((lat+87.5)/2.5)+1;
% lon2ix  = @(lon) round((lon+180)/5.0)+1;
for jj = 1 : length( ca1 )
    buf = regexp( ca1{jj}, '\n', 'split', 'once' );
    buf = regexp( buf{2} , '\n', 'split', 'once' );
    ut  = textscan( buf{1}, '%*f%*f%*f%f%*[^\n]' );
    ut  = ut{1};
    ca2 = regexp( buf{2}, 'LAT/