求助：风云2G卫星降水产品投影坐标

nlj917

MeteoInfo 发表于 2022-3-22 22:14
找一些插值方法试试吧

好的，老师。我读取出来的降水数据矩阵是2288*1144，经、纬度两个矩阵都是2288*2288，不知道怎么处理才能每个网格对应起来。

MeteoInfo

nlj917 发表于 2022-3-22 22:24
好的，老师。我读取出来的降水数据矩阵是2288*1144，经、纬度两个矩阵都是2288*2288，不知道怎么处理才能 ...

取北半球那一半坐标（纬度大于等于0）

nlj917

MeteoInfo 发表于 2022-3-22 22:27
取北半球那一半坐标（纬度大于等于0）

老师，请问我读取坐标的方法是否正确呢？我用以下方法读取了.dat文件，但是感觉评估的时候对应不上。
% clear;clc;
datainfo=hdf5info('C:\Users\Dell\Desktop\FY2G_PRE_024_NOM_20160101_TOAD.hdf');
CTT=hdf5read(datainfo.GroupHierarchy.Datasets(1));
filename='C:\Users\Dell\Desktop\NOM_ITG_2288_2288(0E0N)_LE.dat';
f1 = fopen(filename, 'r');
data = fread(f1,[2288,2288*2],'real*4'); %single 4bytes 读取
data(find(data==300))=nan;
lon=data(1:2288,1:2288)+104.5;
lat(1:1:2288,1:1:2288)=data(1:2288,2289:4576);
fclose(f1);

MeteoInfo

nlj917 发表于 2022-3-22 22:36
老师，请问我读取坐标的方法是否正确呢？我用以下方法读取了.dat文件，但是感觉评估的时候对应不上。
%  ...

matlab不怎么会用，参考MeteoInfoLab的代码：

fn = r'D:\Temp\FY\NOM_ITG_2288_2288(0E0N)_LE\NOM_ITG_2288_2288(0E0N)_LE.dat'
data = binread(fn, (2,2288,2288), datatype='float')
data[data==300] = nan
lat = data[0]
lon = data[1] + 104.5
lat = lat[:,1144:]
lon = lon[:,1144:]
lat = rot90(lat, k=3)
lon = rot90(lon, k=3)

nlj917

MeteoInfo 发表于 2022-3-22 23:10
matlab不怎么会用，参考MeteoInfoLab的代码：

fn = r'D:\Temp\FY\NOM_ITG_2288_2288(0E0N)_LE\NOM_ITG ...

好的，谢谢老师。还有最后一个疑问：请问老师怎样在大矩阵里挑选出某个小流域的经纬度矩阵呢？因为读取出来不是常规的等经纬度矩阵（小网格不是矩形），所以该如何挑选其中经纬度边界范围已知的小流域呢？谢谢老师。

MeteoInfo

nlj917 发表于 2022-3-23 09:49
好的，谢谢老师。还有最后一个疑问：请问老师怎样在大矩阵里挑选出某个小流域的经纬度矩阵呢？因为读取出 ...

判断每个数据点是否在流域多边形范围内，MatLab之类的软件应该都有相关的算法

nlj917

MeteoInfo 发表于 2022-3-23 09:56
判断每个数据点是否在流域多边形范围内，MatLab之类的软件应该都有相关的算法

好的，谢谢老师。