natgrid插值的问题

玉面飞猪

总是报错，提示fatal:Dimension sizes of left hand side and right hand side of assignment do not match
百思不得其解。资料用的micaps资料

脚本如下：
filepath="./06071508.000"

   argu    = asciiread(filepath,-1,"string") ;以字符串形式读取参数文件入数组argu
   station = stringtofloat(str_get_field(argu(11::),1," ")) ;从数组lines中获取站号
   lon     = stringtofloat(str_get_field(argu(11::),2," ")) ;从数组lines中获取经度值lon
   lat     = stringtofloat(str_get_field(argu(11::),3," ")) ;从数组lines中获取纬度值lat
   height  = stringtofloat(str_get_field(argu(11::),4," ")) ;从数组lines中获取海拔高度
   rain    = stringtofloat(str_get_field(argu(11::),5," ")) ;从数组lines中获取站点数据值
;插值后的范围*********************************************************
     data1 = new((/40,66/),"float")
     data1!0="lat"
     data1!1="lon"
     data1&lat=fspan(20,32,40)
     data1&lon=fspan(110,122,66)
     data1&lat@units="degrees_north"
     data1&lon@units="degrees_east"
     rain@_FillValue = 999999.000000
    data1 = natgrid(lon,lat,rain,data1&lon,data1&lat)

兰溪之水

在哪一行开始报错了？把具体的报错信息贴出来吧~

玉面飞猪

问题已解决，是因为natgrid函数中lon，和lat标反了

上帝也厕所

兰溪之水 发表于 2012-3-18 20:19

                               
登录/注册后可看大图

在哪一行开始报错了？把具体的报错信息贴出来吧~

麻烦问一下格点插值到站点是哪个函数

上帝也厕所

玉面飞猪 发表于 2012-3-20 21:06

                               
登录/注册后可看大图

问题已解决，是因为natgrid函数中lon，和lat标反了

麻烦问一下格点插值到站点是哪个函数

孤独卡齐莫多

上帝也厕所 发表于 2012-9-24 16:22
麻烦问一下格点插值到站点是哪个函数

同问,请问您的问题解决了么?

四月雨之音

兰溪之水 发表于 2012-3-18 20:19
在哪一行开始报错了？把具体的报错信息贴出来吧~

您好，我用natgrid的时候把缺测的踢掉了，但是程序走到插值这里就停了，报错内容Natgrid - two input triples have the same x/y coordinates
            but different data values:

                First triple:  109.000000 25.000000 251.000000
                Second triple: 109.000000 25.000000 268.000000
不知道是什么原因造成？确定经纬度信息没有重叠的站点

Clover

四月雨之音 发表于 2020-8-1 22:01
您好，我用natgrid的时候把缺测的踢掉了，但是程序走到插值这里就停了，报错内容Natgrid - two input tri ...

遇到同样的问题，请问你解决了吗？

lfy222

Clover 发表于 2020-11-2 17:29
遇到同样的问题，请问你解决了吗？

您好，我现在也遇到同样问题，您解决了吗