[STATSGRASS] v.variogram

Roger Bivand Roger.Bivand at nhh.no
Sun Oct 8 16:18:54 EDT 2006


On Sun, 8 Oct 2006, Tomas Lanczos wrote:

> > From: Roger Bivand [mailto:Roger.Bivand at nhh.no] 
> 
> > This looks as though the modification stops v.out.ogr -e from 
> > working, and results in no shapefile being written. Could you 
> > see what happens when you run v.out.ogr -e for your point 
> > vector data? My guess is that it will fail. It may also fail
> 
> Happens this, in fact the export performs:
> 
> GRASS_INFO_WARNING(15398,1): OGR can't parse PROJ.4-style parameter string:
> GRASS_INFO_WARNING(15398,1):  +proj=krovak +a=6377397.155 +rf=299.1528128
> +no_defs +to_meter=1.0000000000
> GRASS_INFO_WARNING(15398,1): (OGR Error code was 5)
> Exporting 89 points/lines...
> 89 features written 
> 

So the shapefile does get created? I do not think that we will get any 
further with v.variogram - please just use the R/GRASS interface 
interactively. Trying to encapsulate things that are changing on both 
sides probably wastes more time than it saves.

Roger

> 
> > with just v.out.ogr, because v.out.ogr tries to write a *.prj 
> > anyway for shapefiles, and that may just not work following 
> > the patch. Could you let me know if other OGR drivers do any better?
> 
> It's almost the same with the other drivers.
> 
> Tomas
> 
> > Roger
> > 
> > > Error in ogrInfo(dsn = dsn, layer = layer) : 
> > > 	Cannot open file
> > > Execution halted"
> > > 
> > > I suppose that this error occures when grass, resp. ogr is 
> > reading the 
> > > computed data back. Is it problem with the proj.4 style 
> > itself (how is 
> > > written back) or that I am using the patched proj4 (the patch was 
> > > prepared by Jachym Cepicky - I guess, as the proj4 does not 
> > handle the 
> > > krovak projection correctly). Am I right? Do You have any 
> > idea how to solve this?
> > > 
> > > > Apart from that, the script depends on using xpdf to view 
> > pdf files, 
> > > > an you don not have it installed or on your path.
> > > 
> > > Of course :-).
> > > 
> > > Thanks for any suggestion and/or advice in advance
> > > 
> > > Tomas
> > > 
> > > 
> > > 
> > > > Roger
> > > > 
> > > > > 
> > > > > Loading required package: sp
> > > > > Loading required package: maptools Loading required package: 
> > > > > foreign Loading required package: rgdal Geospatial Data 
> > > > > Abstraction Library extensions to R successfully loaded 
> > Error in 
> > > > > validObject(.Object) : invalid class
> > > > "GridTopology" object:
> > > > > cellsize has incorrect dimension
> > > > > Execution halted
> > > > > Done/usr/local/grass-6.3.cvs/scripts/v.variogram: line 
> > 135: xpdf: 
> > > > > command not found"
> > > > > 
> > > > > What can I do with it?
> > > > > 
> > > > > Note: If I tried to do kriging by the R I recognized that
> > > > there is no
> > > > > krige.G function for spgrass6.
> > > > > 
> > > > > Thank for any advice in advance.
> > > > > 
> > > > > Tomas
> > > > > 
> > > > > p.s: running Grass 6.3 cvs and R 2.3.0 on Ubuntu 6.0.6,
> > > > the databaze
> > > > > is PostgreSQL 8.1.3
> > > > > 
> > > > > _______________________________________________
> > > > > 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, Helleveien 30, 
> > > > N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
> > > > e-mail: Roger.Bivand at nhh.no
> > > > 
> > > 
> > 
> > --
> > Roger Bivand
> > Economic Geography Section, Department of Economics, 
> > Norwegian School of Economics and Business Administration, 
> > Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; 
> > fax +47 55 95 95 43
> > e-mail: Roger.Bivand at nhh.no
> > 
> > 
> 

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




More information about the grass-stats mailing list