[STATSGRASS] v.variogram

[Markus Neteler]

On Sat, Oct 07, 2006 at 09:36:39PM +0200, Roger Bivand wrote:
> 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)

Hi,

probably "krovak" is the problem here. It was recently added to
PROJ (see that mailing list for details) and there were also suggestions
to rename it. I didn't follow that in detail. On the GRASS user
list should be some krovak users who may have an idea. Since OGR
fails, it seems to be a more general problem. It could also be that
GRASS needs krovak to be added, in GRASS 6.3 Glynn Clements recently
rewrote the projection definition part and simplified it a lot.


Markus

 
> 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
> 
> 
> _______________________________________________
> statsgrass mailing list
> statsgrass at grass.itc.it
> http://grass.itc.it/mailman/listinfo/statsgrass

-- 
Markus Neteler  <neteler itc it>  http://mpa.itc.it/markus/
ITC-irst -  Centro per la Ricerca Scientifica e Tecnologica
MPBA - Predictive Models for Biol. & Environ. Data Analysis
Via Sommarive, 18        -       38050 Povo (Trento), Italy