[STATSGRASS] problems assesing mapsets

Roger Bivand Roger.Bivand at nhh.no
Thu Jan 13 06:50:31 EST 2005


On 13 Jan 2005, Marta Benito wrote:

> > Thank you Roger, yes, I am using GRASS 5.7 , and 
> > now it works:
> > 
> > 
> > > G <- gmeta(interp=TRUE)
> > > str(G)
> > > library(GRASS)
> > GRASS environment variables in: /tmp/grass6-marta-8901/gisrc
> > > G <- gmeta(interp=TRUE)
> > > str(G)
> > List of 17
> >  $ LOCATION: chr "espanya_modelo_devel"
> >  $ MAPSET  : chr "mapset1"
> >  $ proj    : chr "projection: 1 (UTM); zone:       30; datum:     
> > eur50; ellipsoid:  int
> > ernational"
> >  $ n       : num 4986657
> >  $ s       : num 3785668
> >  $ w       : num -176002
> >  $ e       : num 1259983
> >  $ nsres   : num 1000
> >  $ ewres   : num 1000
> >  $ Nrow    : num 1293
> >  $ Ncol    : num 1546
> >  $ Ncells  : num 2e+06
> >  $ xlim    : num [1:2] -176002 1259983
> >  $ ylim    : num [1:2] 3785668 4986657
> >  $ xseq    : num [1:1546] -175537 -174608 -173680 -172751 -171822 ...
> >  $ yseq    : num [1:1293] 3786132 3787061 3787990 3788919 3789848 ...
> >  $ ryseq   : num [1:1293] 4986193 4985264 4984335 4983406 4982477 ...
> >  - attr(*, "class")= chr "grassmeta"
> > 
> > but: 
> > 
> > >x <- rast.get(G, rlist="mapa")
> > >Error in rast.get(G, rlist = "mapa") : region for current mapset ??
> > run "g.region"
> > 
> > so, I have tried the same:
> > 
> > x <- rast.get(G, rlist="mapa", interp=T)
> > 
> > 
> > and now everything works, thanks a lot

Could other 5.7/6.0beta users report problems if possible - obviously, 
interp=TRUE is not the big performance hit it was 4-5 years ago, but it 
isn't a very robust solution. Does it also work for cygwin?

Roger

> > 
> > Marta
> 
> El jue, 13-01-2005 a las 12:15, Roger Bivand escribió:
> > On 13 Jan 2005, Marta Benito wrote:
> > 
> > > Hi, 
> > > I am having problems using the GRASS library in R. 
> > > I think I am doing something wrong when using gmeta but I don't know
> > > exactly, the problem began when I've done update of packages in R and
> > > GRASS. 
> > > 
> > > GRASS 5.7.cvs:~/grassdata/espanya_modelo_devel > R
> > > > library(GRASS)
> > > GRASS environment variables in: /tmp/grass6-marta-8901/gisrc
> > >  G<-gmeta()
> > > Error in gmeta() : region for current mapset MAPSET
> > > run "g.region"
> > 
> > You seem to be using GRASS 5.7 or 6.0beta1, is that correct?
> > 
> > Please try G <- gmeta(interp=TRUE), and post the result (str(G)) - this 
> > may be a WIND3 WINDOW not being recognised by the 5.0 version of the GRASS 
> > function used when run with interp=FALSE. 
> > 
> > Roger
> > 
> > > 
> > > I have already run g.region in GRASS, and also after the error. 
> > > Thank you in advance 
> > > 
> > > Marta
> > > 
> > > 
> > > _______________________________________________
> > > statsgrass mailing list
> > > statsgrass at grass.itc.it
> > > http://grass.itc.it/mailman/listinfo/statsgrass
> > > 
> > 
> > -- 
> > Roger Bivand
> > Economic Geography Section, Department of Economics, Norwegian School of
> > Economics and Business Administration, Breiviksveien 40, N-5045 Bergen,
> > Norway. voice: +47 55 95 93 55; fax +47 55 95 93 93
> > e-mail: Roger.Bivand at nhh.no
> > 
> > 
> > _______________________________________________
> > statsgrass mailing list
> > statsgrass at grass.itc.it
> > http://grass.itc.it/mailman/listinfo/statsgrass
> > 
> > 
> 
> 
> _______________________________________________
> statsgrass mailing list
> statsgrass at grass.itc.it
> http://grass.itc.it/mailman/listinfo/statsgrass
> 

-- 
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Breiviksveien 40, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 93 93
e-mail: Roger.Bivand at nhh.no





More information about the grass-stats mailing list