[STATSGRASS] v.variogram

[Roger Bivand]

On Sat, 7 Oct 2006, Tomas Lanczos wrote:

> Hello,
> 
> > > v.variogram point=basement_surf layer=1 column=z resolution=70 
> > > maxdist=92330
> > > nlag=15
> > 
> > The script has not been updated to match changes in the 
> > spgrass6 package. 
> > west and south in line 117 of the script need to be changed 
> > to w and s from spgrass6 >= 0.3. The list in line 118 will 
> > need changing to data.frame for sp >= 0.9 and R >= 2.4.0. The 
> > script would be more robust if is used the mechanisms 
> > described in the GRASS book for argument passing. It could 
> > also condition on package releases.
> 
> Thank You very much for Your answer, I performed the changes You suggested.
> Now I got the following error message:
> 
> "GRASS_INFO_WARNING(12466,1): OGR can't parse PROJ.4-style parameter string:
> GRASS_INFO_WARNING(12466,1):  +proj=krovak +a=6377397.155 +rf=299.1528128
> +no_defs +to_meter=1.0000000000
> GRASS_INFO_WARNING(12466,1): (OGR Error code was 5)

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 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?

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