Hi Roger and colleagues, Everithing worked fine!
Remembering: I am working under VISTA 64bits, running grass65 that I compiled yesterday under MSYS. I installed the alpha version of R 2.9.0. But I had several issues regarding Security on Vista. I tryed several things, like install R 2.9.0 (also I tryed R 2.8.1), change the priority of security allowing R run As Administrator, and after installing the packages (sp, maptools, spatstat, adehabitat, rdgal and spgrass6), start MSYS>GRASS>R and require spgrass6. BUT I got error message: Loading required package: spgrass6 Warning message: In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, : there is no package called 'spgrass6' But if I start R from GRASS, and I install the packages after start R from GRASS/MSYS everything works fine. So I suggest (1) start MSYS (2) start GRASS (3) call R from GRASS/MSYS ( /c/Progra~2/R/R-2.9.0alpha/bin/Rgui.exe& ) (4) install the packages. Thanks Roger, issue solved! milton On Sun, Mar 29, 2009 at 11:04 AM, Roger Bivand <[email protected]> wrote: > On Sun, 29 Mar 2009, milton ruser wrote: > > Checking the result on the readRAST6, I noticed that >> the number of bytes are 2, and I need 4 bits to store >> values greather than 32565, I suppose. >> -------------- >> >>> mapa <- readRAST6(varint, cat=FALSE) >>> >> Creating BIL support files... >> Header File = >> >> C:/Data/Tibagi/vegPR_probio_proatlan_abr2009/tibagi_Probio_grass_sdb_apr2009/newLocation/PERMANENT/.tmp/uso_tibagi_fundo_soflor_ero30m_dila030m_mata_clump.hdr >> World File = >> >> C:/Data/Tibagi/vegPR_probio_proatlan_abr2009/tibagi_Probio_grass_sdb_apr2009/newLocation/PERMANENT/.tmp/uso_tibagi_fundo_soflor_ero30m_dila030m_mata_clump.wld >> Exporting raster as integer values (bytes=2) >> Using the current region settings... >> -------------------- >> But how can I chance the number of bytes, if it is not an >> argument on readRAST6? >> > > Thanks for including your sessionInfo() in the preceding mail. If you are > still using 0.5-10, you have no choice but to cast the data to floating > point using r.mapcalc (add 0.0 should do it). There have been many changes > in readRAST6() since that version, for example to use a GRASS GDAL plugin if > available, and to use GDAL to read and write GTiff temporary files. The > older route using r.out.bin only offers 2-byte integers, so was a reason for > making more formats available. > > So the answer in general is to update sp, spgrass6 and rgdal if you can. > Since you in addition on Windows, and probably using the Native Windows > version, until a new release you will hit a snag that more recent GRASS > r.out.gdal have a -c flag which is not in your release - here you can wait > for the next release, and use the trick of casting to floating point until > then. > > Roger > > Best wishes >> >> miltinho >> >> >> On Sun, Mar 29, 2009 at 12:43 AM, milton ruser <[email protected] >> >wrote: >> >> Hi there, >>> >>> I have a set of maps on GRASS 6.5 and I need to compute some basic stats >>> on R 2.7.1. For some maps I get the desired results, but for other I am >>> getting different values. My raster are integer values >= 1 (with some >>> NAs >>> where not mapped). >>> >>> I start R from GRASS session without problem, and can read the rasters >>> on R: >>> >>> require(spgrass6) >>> varint<-"uso_tibagi_fundo_soflor_ero30m_dila030m_mata_clump_AreaHA" >>> mapa <- readRAST6(varint, cat=FALSE) >>> mapa.valores<-m...@data[1] >>> mapa.valores<-mapa.valores[!is.na(mapa.valores)] >>> quantile(mapa.valores, probs=c(0.20,0.40,0.60,0.80)) >>> summary(mapa.valores) >>> rm(mapa, mapa.valores,varint) >>> But for one of my map I get different range of value between R and GRASS >>> ON GRASS - r.univar >>> n: 8336636 >>> minimum: 1 >>> maximum: 303396 >>> range: 303395 >>> mean: 125949 >>> >>> ON R >>> >>>> summary(m...@data[1]) >>>> >>> uso_tibagi_fundo_soflor_dila030m_clump_mata_limpa_AreaHA >>> Min. : -24284 >>> 1st Qu.: -24284 >>> Median : 31 >>> Mean : -6551 >>> 3rd Qu.: 2225 >>> Max. : 27846 >>> NA's :59217475 >>> >>> -----> sessionInfo() >>> R version 2.7.1 (2008-06-23) >>> i386-pc-mingw32 >>> locale: >>> >>> >>> LC_COLLATE=English_Canada.1252;LC_CTYPE=English_Canada.1252;LC_MONETARY=English_Canada.1252;LC_NUMERIC=C;LC_TIME=English_Canada.1252 >>> attached base packages: >>> [1] stats graphics grDevices utils datasets methods base >>> other attached packages: >>> [1] spgrass6_0.5-10 rgdal_0.5-25 sp_0.9-26 >>> loaded via a namespace (and not attached): >>> [1] grid_2.7.1 lattice_0.17-8 >>> Any help are welcome. >>> >>> Regards >>> >>> miltinho >>> Toronto, CA >>> >>> >>> >> [[alternative HTML version deleted]] >> >> _______________________________________________ >> R-sig-Geo mailing list >> [email protected] >> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >> >> > -- > 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: [email protected] > > [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
