如何使用样条插值为栅格数据与一个特定的条件?

2视图(30天)

显示旧的评论

• 2018年3月9日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/387387-how-to-use-spline-interpolation-for-raster-data-with-a-specific-condition

评论道: Unai San Miguel 2018年7月9日

新folder.zip

嗨,我有一个从卫星降雨测量NetCDF文件(请看附件)。每个网格位于细胞的中心。我想做样条插值与一个特定的条件。

样条插值在默认情况下实施以下两个条件,1)表面经过完全通过数据点和2)表面必须有最小曲率。

但在这种情况下,我也想使用一个额外的约束。插值的每个像素值的约束条件是必须符合相应的卫星像素值。换句话说,每个卫星像素的值必须等于所有内插值像素的平均值。

                         %的降水地图原始卫星数据
                        
                         文件=' 3 b-day.ms.mrg.3imerg.20150903-s000000-e235959.v05.nc4.nc ';
                        
                         pre = ncread(文件,“precipitationCal”);
                        
                         朗= ncread(文件,“朗”);
                        
                         lat = ncread(文件,“纬度”);
                        
                         图
                        
                         地图= pcolor(经度、纬度、pre);
                        
                         colorbar

每个网格的卫星之间的距离是0.1度,插值后的距离是0.001度,LatitudeLimits = (46.6 - 49.1), LongitudeLimits = (13.9 - 17.3)。

样条插值的代码到目前为止我所做的可以看到下图:

                         %的降水插值降水数据的地图
                        
                         [xq, yq] = ndgrid(1:0.01:长度(lat), 1:0.01:长度(朗));
                        
                         F = griddedInterpolant(前样条的);
                        
                         Vq1 = F (xq, yq);
                        
                         图
                        
                         xq, map_2 = pcolor (yq Vq1);
                        
                         集(map_2,“EdgeAlpha”,0);
                        
                         colorbar