Re: [R-sig-Geo] some problems with krigeST

2013-10-08 Thread Edzer Pebesma
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On 10/09/2013 05:23 AM, 胡艺 wrote: > Dear all, > > I'm a newer to spatio-temporal analysis though I am quite familiar > with spatial one. When coming to spatio-temoral prediction using > krigeST, I have some queries as follows: > > 1. it seems that

[R-sig-Geo] some problems with krigeST

2013-10-08 Thread 胡艺
Dear all, I'm a newer to spatio-temporal analysis though I am quite familiar with spatial one. When coming to spatio-temoral prediction using krigeST, I have some queries as follows: 1. it seems that the krigeST can only deal with data with full layout format (e.g., STFDF class), however our data

Re: [R-sig-Geo] Creating a maximum raster layer from RasterStack

2013-10-08 Thread Agustin Lobo
Eddie, (I'm answering through the list, perhaps somebody else is interested) You need one element in "indices" for each input layer, so if you have 5 layers you should use mx <- stackApply(s,indices=c(1,1,1,1,1),max) The actual value of each index indicates to which output layer that input layer i

[R-sig-Geo] openmap on multiple plot

2013-10-08 Thread pierluigi de rosa
Dear all, I found that the code nm <- c("osm", "maptoolkit-topo", "mapquest", "mapquest-aerial", "bing", "stamen-toner", "stamen-watercolor", "esri", "esri-topo", "nps", "apple-iphoto", "skobbler") par(mfrow=c(3,4)) #Korea for(i in 1:length(nm)){ map <- openmap(c(43.468

[R-sig-Geo] change color palette to ggmap

2013-10-08 Thread pierluigi de rosa
Dear R-geo user, i wanto to change the default palette colour in map produced by ggmap. here the code I use: val_area <- openmap(c(43.275,12.500), c(43.227,12.588), minNumTiles=13, mergeTiles = TRUE, type='opencyclemap

Re: [R-sig-Geo] problem reading large ncdf on windows 7 64 bits

2013-10-08 Thread Ariel Ortiz-Bobea
I also ran into this problem. It was solved by installing the "ncdf4" package. Instructions can be found here: http://cirrus.ucsd.edu/~pierce/ncdf/ - Ariel Ortiz-Bobea Fellow at Resources for the Future -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/problem-readi

[R-sig-Geo] unable to load shared object

2013-10-08 Thread cpadwick padwick
Hi All, I'm trying to install rgdal on my centOS linux box and I'm running into a problem: Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared object '/home/cpadwick/R/x86_64-redhat-linux-gnu-library/2.15/rgdal/libs/rgdal.so': /home/user/R/x86_64-redhat-linux-gnu-library/

Re: [R-sig-Geo] problem reading hdf image

2013-10-08 Thread Agustin Lobo
Thanks! The following works: r <- raster("HDF4_EOS:EOS_GRID:/media/Iomega_HDD/VOLCAN/Modis_adls_MOD09A1/MOD09A1.A2011289.h11v11.005.2011298101214.hdf:MOD_Grid_500m_Surface_Reflectance:sur_refl_b01") And much beter, the following works as well: modislista <- getSds(file.path(modisdir,"MOD09A1.A20

Re: [R-sig-Geo] problem reading hdf image

2013-10-08 Thread Agustin Lobo
Yes, no problem from that side. I've used gdal_translate -sds and get all subdatasets saved as individual tif layers. But we should know what's wrong with the hdf flle as there are many thousands of files like this on the modis distribution data servers. Agus On Tue, Oct 8, 2013 at 1:58 PM, José

Re: [R-sig-Geo] ggmap error

2013-10-08 Thread Jari Oksanen
On 08/10/2013, at 13:40 PM, Barry Rowlingson wrote: > What version of R do you have? The paste0 function has been in R for a > while now... More exactly, paste0() has been in R release since April 2012. It was introduced in R 2.15.0 which was released on March 30, 2012. In my book that is such

Re: [R-sig-Geo] problem reading hdf image

2013-10-08 Thread Matteo Mattiuzzi
Section 'Subdatasets'ie: SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"MOD09A1.A2011289.h11v11.005.2011298101214.hdf":MOD_Grid_500m_Surface_Reflectance:sur_refl_b01 # this should work if R is in the HDF directory! if not redo gdalinfo with the file absolute path: gdalinfo /path/to/MOD09A1.A2011289.h11v

Re: [R-sig-Geo] problem reading hdf image

2013-10-08 Thread José Miguel Barrios
Have you tried converting the layer(s) of interest to something like ".tif" by using gdal_translate? 2013/10/8 Agustin Lobo > Here it goes, thaks a lot for your help: > > $ gdalinfo MOD09A1.A2011289.h11v11.005.2011298101214.hdf > Driver: HDF4/Hierarchical Data Format Release 4 > Files: MOD09A1.

[R-sig-Geo] Questions about creation of spsurvey.analysis object

2013-10-08 Thread Alessandro Samuel Rosa
Dear R users, I have a few questions about the use of the function spsurvey.analysis. First I'll describe my data and after I'll show you how I'm using the function. I hope someone has the time to help me out. My data is composed by n = 60 soil observations which were sampled in linear transec

Re: [R-sig-Geo] problem reading hdf image

2013-10-08 Thread Agustin Lobo
Here it goes, thaks a lot for your help: $ gdalinfo MOD09A1.A2011289.h11v11.005.2011298101214.hdf Driver: HDF4/Hierarchical Data Format Release 4 Files: MOD09A1.A2011289.h11v11.005.2011298101214.hdf Size is 512, 512 Coordinate System is `' Metadata: ASSOCIATEDINSTRUMENTSHORTNAME=MODIS ASSOCIAT

Re: [R-sig-Geo] problem reading hdf image

2013-10-08 Thread Matteo Mattiuzzi
Dear Agus, yes please send the output of gdalinfo (close to the end is the string you have to use). Matteo >>> Agustin Lobo 10/07/13 5:42 PM >>> Hi! I get the following error when try to read an hdf image (release 4): > modisima <- > stack(file.path(modisdir,"MOD09A1.A2011289.h11v11.0

Re: [R-sig-Geo] ggmap error

2013-10-08 Thread Barry Rowlingson
What version of R do you have? The paste0 function has been in R for a while now... On Tue, Oct 8, 2013 at 11:28 AM, pierluigi de rosa wrote: > Dear R-Geo user > I have a trouble when using any function from ggmap package... > It tells me > Error in message(paste0("Map from URL : ", url)) : > n

[R-sig-Geo] ggmap error

2013-10-08 Thread pierluigi de rosa
Dear R-Geo user I have a trouble when using any function from ggmap package... It tells me Error in message(paste0("Map from URL : ", url)) : not finding function "paste0" Some one experienced the same issue? thanks P -- Ing. Pierluigi De Rosa (PhD) tel: 075 791 0278 / fax: 075 7823038 cel: 3