谱模式格点怎么插值到经纬度格点上？

wheelar

某GRIB1数据，具体内容是：
wgrib -V A_HZJZ50RJTD300000_C_RJTD_20150130033932.bin
rec 1:43:date 2015013000 RELV kpds5=43 kpds6=100 kpds7=500 levels=(1,244) grid=38 500 mb 66hr fcst:
RELV=Relative vorticity
timerange 0 P1 66 P2 0 TimeU 1 nx -1 ny 73 GDS grid 0 num_in_ave 0 missing 0
center 34 subcenter 0 process 4 Table 3 scan: WE:SN winds(N/S)
thinned latlon: lat 0.000000 to 90.000000 by 1.250000 nxny 3447
long 60.000000 to 150.000000, 3447 grid pts (-1 x 73) scan 64 mode 128 bdsgrid 1
73 73 73 73 73 73 73 73 72 72 72 71 71 71 70
70 69 69 68 67 67 66 65 65 64 63 62 61 60 60
59 58 57 56 55 54 52 51 50 49 48 47 45 44 43
42 40 39 38 36 35 33 32 30 29 28 26 25 23 22
20 19 17 16 14 12 11 9 8 6 5 3 2
min/max data -0.000147048 0.000177202 num bits 12 BDS_Ref -0.000424567 DecScale 0 BinScale -21

据说是“谱模式格点”，每个纬度的数目不一样，赤道多一点，极地少，那串数字(73, 73, ... , 2) 共 73 个，加起来就 3447。第一个 73 就是 0 度的，最后一个 2 就是 90 度的。那请问用哪种插值方法插值到经纬度格点呢？是不是有专门的针对这种转换的方法？还是用普通的插值程序？插值好像有好多种的，比如线性插值，拉格朗日插值、三次样条等等，不知该用什么好？另外有Fortran程序吗？非常感谢！！

lzuatmos

是同一坐标系还是不同坐标系之间的转换？