MATLAB的答案

有效的多重网格插值的方法

1视图(30天)
FeAde
FeAde 2020年11月28日
回答: Sindar 2020年11月30日
我有分散的点雨量计站随机缺失值。我也有NDVI值在这些站(完整的记录)。从2016年1月-Decemeber,数据。
我想插入这两个变量具有独立性在日常时间范围不包括任何南车站为这两个变量(即南排除在降雨和相应的归一化植被指数)。他们将被保存为Geotiff。这是我的代码一个网格。
我想做一个循环与南在每个时间步排除在外。
% % %加载数据
负载data_now.dat
(经度、纬度)= ndgrid(经度、纬度);
lonobs = data_now (: 1);
latobs = data_now (:, 2);
J = data_now (: 3);
G = data_now (: 4)
%经度和纬度netcdf文件读取
朗= ncread (“important_gh.nc”,“朗”);
lat = ncread (“important_gh.nc”,“纬度”);
%插入
Rainq = griddata (latobs lonobs, J,纬度、经度);
NDVIq = griddata (latobs lonobs, G,纬度、经度);
谢谢

在回答这个问题。

答案(1)

Sindar
Sindar 2020年11月30日
我不知道它会有多快,但我认为这将工作
%经度和纬度netcdf文件读取
朗= ncread (“important_gh.nc”,“朗”);
lat = ncread (“important_gh.nc”,“纬度”);
(经度、纬度)= ndgrid (dl.lon dl.lat);
data_file =…
为印第安纳州= 1:长度(data_file)
%加载到一个变量
data_now =负载(data_file{印第安纳});
lonobs = data_now (: 1);
latobs = data_now (:, 2);
J = data_now (: 3);
G = data_now (: 4);
%找到有效的位置
idx = (~ isnan (J)) & (~ isnan (G));
%插入,删除南的位置
Rainq = griddata (latobs (idx) lonobs (idx), J (idx)、纬度、经度);
NDVIq = griddata (latobs (idx) lonobs (idx), G (idx)、纬度、经度);
%保存
结束

类别

找到更多的在插值在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!