用MATLAB读取nc文件某一经纬度的数据

Lisahu

fname_in     = 'ec1.nc';
time         = ncread(fname_in,'time');
latitude          = ncread(fname_in,'latitude');  % 纬度
longitude          = ncread(fname_in,'longitude');  % 经度
r          = ncread(fname_in,'r');
longitude_index     = find(longitude ==117); % 选出经度在该区域的数据
latitude_index     = find(latitude == 34.5); % 选出纬度在该区域的数据
is_r_nan  = isnan(r);          % 返回与r一样大小的向量，判断r的值是否为NaN，是则为1，不是则为0
out_file_id   = fopen('data1.txt','w'); % 存放数据文件
for i =1 : length(time)
   
      
    if(is_r_nan(longitude_index ,latitude_index,i) == 0)        % 不为缺省值才做计算
        fprintf(out_file_id, '%.3f\t\n%.3f\t', longitude(longitude_index), latitude(latitude_index));
        
        fprintf(out_file_id, '%.3f\r\n', r(longitude_index,latitude_index,i));
    end
   
end
fclose(out_file_id);

Lisahu

所要读取的经纬度必须你的文件当中有 然后到txt文件当中查看你要的数值

五月雨

CMIP5的数据经度是0-360，请问0-180是东经吗？

zhouchen

我的总提示报错

xyyyy

谢谢楼主分享，学习了

易水

感谢楼主分享，学习了

阿童木

请问如果我是读取某一范围的经纬度怎么读取和编写程序呢？谢谢

李伟好气象

楼主，您能不能再发一份“.m”格式的文件，附带一下解释，您这个程序我就跑不动啊。

36N

五月雨 发表于 2019-2-27 17:48
CMIP5的数据经度是0-360，请问0-180是东经吗？

同学这个问题解决了嘛？0-180是东经度嘛？

深深

请问要提起某一范围的数据  怎么便是呢？ 比如北纬24-30，东经112-119？