关于用站点数据画图出的问题，希望大神们帮忙

羊羊羊

我手里有已经转化好的二进制站点数据，和map的映射文件，写了站点数据的ctl，以上部分都没有错，接下来用fortran编了画格点的程序，又为格点写了ctl，之后写gs运行就报错，我把过程贴出来，希望大家能帮忙解决~~谢谢~~
这是站点数据的ctl
rainfall.ctl
DSET    f:\data\prepitation\1.dat
UNDEF   -999.9
TITLE   rainfall
DTYPE   station
STNMAP  f:\data\prepitation\1.map
TDEF    11993 linear 01Jan1980 1dy
VARS 1
p 0 99 rainfall
ENDVARS

这是画格点的fortran
Program  stn2grd
      parameter(nx=71, ny=41)
      real lat(ny), lon(nx)
      real s(nx,ny)
      open(1, file= 'f:/data/prepitation/2.dat', form='binary')
      lat(1)=15.0
      lon(1)=70.0
           do j=1,ny-1
                lat(j+1)=lat(j)+1.0
           enddo
           do i=1,nx-1
                lon(i+1)=lon(i)+1.0
           enddo

            do i=1,nx
                 do j=1,ny
                       s(i,j)= 1
                 enddo
            enddo
            write(1)    s
            end

这是格点对应的grid.ctl
dset f:/data/prepitation/2.dat
undef  -999.0
title Grid data
xdef  71     linear   70  1
ydef  41     linear   15   1
zdef  1      linear   1000    1
tdef  11993  linear   01jan1980  1dy
vars
g 1   99   grid data
endvars

这是运行的rain1.gs
'open  d:/grid.ctl'
'open  d:/rainfall.ctl'
'enable print d:/index.gmf'
'set lon 90 112'
'set lat 15 36'
'set time 3jun1992 14jun1992'
'set mpdset hires  cnworld'
'define   a=oacres(g,p.2)'
'set gxout shaded'
'd a'
'set gxout contour'
'd a'
'print'
'disable print'
'reinit'

                               
登录/注册后可看大图

f117hqh

格点文件2.dat和ctl对应不上，ctl中描述有11993个时间层，2.dat中显然只有一个时间层

羊羊羊

f117hqh 发表于 2015-3-15 09:29
格点文件2.dat和ctl对应不上，ctl中描述有11993个时间层，2.dat中显然只有一个时间层

可是我改完了grid.ctl的时间层，改成了1，还是提错
open error:missing or invalid variable count
the invalid description file record is:
vars
the data file was not opened

f117hqh

羊羊羊 发表于 2015-3-15 09:52
可是我改完了grid.ctl的时间层，改成了1，还是提错
open error:missing or invalid variable count
the ...

提示的很清楚了
vars 1
g 1   99   grid data
endvars

羊羊羊

f117hqh 发表于 2015-3-15 09:29
格点文件2.dat和ctl对应不上，ctl中描述有11993个时间层，2.dat中显然只有一个时间层

我重新找了一个云南底图差值，发现能出图了，但是我的数据有西藏，云南，广西一片，效果与实际有差距，有没有办法找到三个省，或者这一块的底图差值？

羊羊羊

f117hqh 发表于 2015-3-15 10:17
提示的很清楚了
vars 1
g 1   99   grid data

恩恩，这个问题解决了~谢谢~~我重新用了个底图文件，现在能出了，只是效果有点差，我数据是云南广西西藏一片，只用云南底图做会有差距，请问要怎么做呀？

f117hqh

羊羊羊 发表于 2015-3-15 10:32
我重新找了一个云南底图差值，发现能出图了，但是我的数据有西藏，云南，广西一片，效果与实际有差距，有 ...

可以同时用三个省的底图 记得家园里就有全国分省的底图

羊羊羊

f117hqh 发表于 2015-3-15 10:44
可以同时用三个省的底图 记得家园里就有全国分省的底图

谢谢~各省的底图我都有，只是我刚刚发现，即使我用了云南底图的格点，不改ctl，出图的效果与用全国底图一样，改了ctl的经纬度范围后，发现出的图打不开了，我应该怎么写ctl呢？
这是原来中国底图的ctl
dset d:\grid.grd
undef -999.0
title Sample GRIB Data
xdef 71 linear 70 1
ydef 41 linear 15 1
zdef 1 linear 500 1
tdef 1 linear DEC1979 1yr
vars  1
g   0  99  grid data prepared for oacres function
endvars
这是我后来改的，图打不开的ctl
dset d:\yunn.grd
undef -999.0
title Sample GRIB Data
xdef 11 linear 97 1
ydef 31 linear 21 1
zdef 1 linear 500 1
tdef 1 linear JAN1980 1dy
vars  1
g   0  99  grid data prepared for oacres function
endvars

qxtlyf

羊羊羊 发表于 2015-3-15 11:10
谢谢~各省的底图我都有，只是我刚刚发现，即使我用了云南底图的格点，不改ctl，出图的效果与用全国底图一 ...

vars  1
g   0  99  grid data prepared for oacres function
中g   0  99还是g   1  99？

羊羊羊

qxtlyf 发表于 2015-3-15 11:17
vars  1
g   0  99  grid data prepared for oacres function
中g   0  99还是g   1  99？

是g   0  99  ，开始那个写错了，后来改对了，现在问题是不知道怎么写底图对应的ctl，才能把数据成功插值到云南，云南底图文件我已经有了，谢谢~~