站点资料插值成格点画降水图数据超出边界怎么办

subtropical

我把站点降水资料插值为格点资料画图，但是数据超过边界了，

请问要怎么解决？

然后搜了一些说用shapefile_mask_data ，下载了shapefile_utils.ncl这个脚本（在http://www.ncl.ucar.edu/Applications/shapefiles.shtml例子21里下载的），但是报错，请问这个怎么改？


我的绘图脚本：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
load "./shapefile_utils.ncl"
begin
min_lat = 15.0
max_lat =45.0
min_lon =105.0
max_lon =125.0
;-------------------------------------------------------------------------
f_data1=asciiread("D:/tf/rain/rain19610114.txt", (/2072,5/), "float")
f_data2=asciiread("D:/tf/rain/rain19610115.txt", (/2072,5/), "float")

sta_id=f_data1(:,0)
lon=f_data1(:,1)
lat=f_data1(:,2)
rian=f_data1(:,4)+f_data2(:,4)
;---------------------------------给变量赋给属性
lon@units="degrees_east"
lat@units="degrees_north"
;rian@units="mm"
rian@_FillValue=9999

;---------------------------------创建插值后的经纬度0.1*0.1
olon=fspan(105,145,410)   
olat=fspan(5,45,410)
olon@units="degrees_east"
olat@units="degrees_north"
;---------------------------------使用cressman插值函数
rscan=(/1,0.5,0.25/)
r_final=obj_anal_ic_Wrap(lon, lat, rian, olon, olat, rscan, False)
data_mask = shapefile_mask_data(r_final,"D:/cygwin/app/ncl/lib/ncarg/nclscripts/cnmap/cnmap.shp",True)    ;---shapefile_mask_data

。。。。。

subtropical

自己来顶一下{:eb303:}{:eb303:}

弓长铭

可以参考：

http://bbs.06climate.com/forum.p ... &extra=page%3D2file:///D:/hdpre/program/ex1.pdf

songyy

请问你这个问题解决了吗，我也用gc_inout 的时候报错了

songyy

说是纬度问题

duoci01

能否分享一下完整的脚本？期待您的回复

fly100516

请问楼主的问题解决了吗？谢谢！