NCL 站点作图插值

546243727

sunnygirlbo 发表于 2016-5-31 22:01
不好意思，打扰下，想请问一下，楼主那个插值半径是咋设置的，rscan那个，我对中国东部逐日降水插值，用的n ...

不好意思我也没理解透  (/10,7,4,1/)   我是试的 但没有插出过负值

广帅帅

这两天也写了点关于站点插值格点的程序，就楼主这个问题，发表一些粗略的见解。首先明白函数中的rscan( /10,7,4,1/)的中几个参数的意思。NCL 的插值方法采用的是Cressman 和  Barnes，楼主可以抽空看下Cressman插值原理，简单说在插值过程中，会确定插值过程中的影响因子和影响半径，然后先采用大的影响半径进行插值，再用小的影响半径进行继续插值和修订。那么这个rscan( /10,7,4,1/）的几个参数的意思就是前后插值的影响半径。NCL 官网没有对此作太多的说明，只说明括号内最多有四个参数，并且单调递减。根据我对脚本修改的影响来看，这个几个参数就是针对插值的背景场的”格距“，以你的脚本简单来说，你采用的分辨率是0.1*0.1,所以，你现在采用10、7、4、1就是代表第一次插值的影响半径是10除以0.1，也就是100个影响格点。所以，影响范围太大了，因此可能导致数据插值到省外的情况，你可以试试（0.5,0.1,0.05）试试效果会不会好一点~
（以上只是本人理解，如有不对还望见谅，勿责）

546243727

广帅帅 发表于 2016-6-30 20:05
这两天也写了点关于站点插值格点的程序，就楼主这个问题，发表一些粗略的见解。首先明白函数中的rscan( /10 ...

感谢耐心解答！

dannybear

546243727 发表于 2015-11-25 09:09
请教各位大神：
1.  olon = fspan(97,109,130） ，olat = fspan(26,35,100)  控制了分辨率吗 ， 是把分 ...

请问楼主问题解决了吗?我也遇到了和你类似的问题

dannybear

546243727 发表于 2015-11-25 09:09
请教各位大神：
1.  olon = fspan(97,109,130） ，olat = fspan(26,35,100)  控制了分辨率吗 ， 是把分 ...

rscan怎么选择

546243727

dannybear 发表于 2016-8-24 16:28
请问楼主问题解决了吗?我也遇到了和你类似的问题

楼上有回啊

千里草

546243727 发表于 2015-11-25 09:09
请教各位大神：
1.  olon = fspan(97,109,130） ，olat = fspan(26,35,100)  控制了分辨率吗 ， 是把分 ...

一、是不是控制分辨率你看一下你插值出来的数据是不是0.1X0.1的精度。
二、是因为你插值区域选择的就是这么大，肯定会插值到省外；
三、设置缺测值是没有影响的，或者说更好，因为有些区域（没有插值到）就设置成缺测了；
四、可以mask掉省外的
res@mpOutlineOn           = False        ; Turn on map outlines
  res@mpDataBaseVersion       = "MediumRes"
  res@mpDataSetName           = "Earth..4"
  res@mpAreaMaskingOn         = True
  res@mpMaskAreaSpecifiers    = (/"Si Chuan"/)   ；填写四川省区域，我不确定四川是不是这样表达的，你             可以去NCL官网搜搜
  res@mpLandFillColor         = "white"
  res@mpInlandWaterFillColor  = "white"
  res@mpOceanFillColor        = "white"
  res@mpOutlineBoundarySets   = "NoBoundaries"  
  res@mpFillDrawOrder  = "PostDraw"

546243727

千里草 发表于 2016-8-26 08:02
一、是不是控制分辨率你看一下你插值出来的数据是不是0.1X0.1的精度。
二、是因为你插值区域选择的就是 ...

谢谢耐心解答！

Eudora99

谢谢LZ分享！！！ps:四川省的底图是错的。

yeah...

数据缺测值应该ncl 自动识别的把