Re: [R-sig-Geo] adapting spatial points and wrld_smpl to a reference system implicit in a .nc file

Folks: I went ahead and pushed a new version of gdalUtils to r-forge: install.packages("gdalUtils", repos="http://R-Forge.R-project.org;) gdal_translate has a new option to attempt to fix the issue above. Can you confirm if e.g.: gdal_translate(...,config='GDAL_NETCDF_BOTTOMUP="YES"') works

Re: [R-sig-Geo] Error running Mann-Kendall trend test on a raster stack

Dear Thiago, have a look at the gimms package, available on CRAN. The function of interest is significantTau(). https://github.com/environmentalinformatics-marburg/gimms/blob/master/R/significantTau.R Hope this helps, Tim On 27.02.2016 09:49, Thiago V. dos Santos wrote: Dear colleagues, I

[R-sig-Geo] Error running Mann-Kendall trend test on a raster stack

Dear colleagues, I have a raster stack with 89 layers, each layer representing yearly precipitation. I am trying to use the function rkt (from package "rkt") to detect a possible trend in my precipitation time series. Its usage is pretty simple, and this is how it runs on a data frame:

Re: [R-sig-Geo] adapting spatial points and wrld_smpl to a reference system implicit in a .nc file

Folks: Just noticed this thread -- I see I didn't include a "--config" option with any of the gdalUtils functions (it isn't one of the documented parameters on the individual utility website, but it seems it would have allowed you to run the GDAL_NETCDF_BOTTOMUP without setting a system

Re: [R-sig-Geo] resolution of openmap() raster layers

Agus, Mapview is using leaflet as engine. Due to this you will have the control icons on the map because first of all it is designed for interactive mapping within RStudio/R. I think there are two different approaches to save your maps: If you want to have a dump of the mapviewobject (but

Re: [R-sig-Geo] Problems with mapmisc::insetMap()

Hi. 1- `insetMap` returns a matrix of x and y coordinates, you needn't plot all of them. ``` nica <- getData("GADM", country="NIC", level=0) nicabg <- openmap(nica, path="landscape") map.new(nicabg,1,mar=c(2, 2, 2, 0) + 0.1) plot(nicabg,axes=TRUE) plot(nica,add=TRUE) loc <-

Re: [R-sig-Geo] odd behaviour of llgridlines() regarding axes labels

1) is ok with xpd clipping turned off op <- par(xpd = NA) plot(a,axes=FALSE) llgridlines(a, plotLines=TRUE, plotLabels=TRUE) par(op) 2) same op <- par(xpd = NA) plot(a,axes=FALSE,mar=c(4, 3, 2, 3)+0.1) llgridlines(a, plotLines=TRUE, plotLabels=TRUE) par(op) 3) Consider using graticule package

Re: [R-sig-Geo] resolution of openmap() raster layers

Stunning! Can I remove the buttons for saving to a bmp file? What attribution should be used for publishing? Agus On Thu, Feb 25, 2016 at 7:42 PM, Chris Reudenbach wrote: > Hi, > > if you just want to map the data, mapview could be an option that among > others avoid

[R-sig-Geo] odd behaviour of llgridlines() regarding axes labels

I see an odd behaviour of llgridlines() regarding axes labels: require(rgdal) require(raster) a <- extent(c(-90,-81.5625,8.407168,16.63619)) a <- as(a, 'SpatialPolygons') projection(a) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" No axis labels on x axis:

Re: [R-sig-Geo] Problems with mapmisc::insetMap()

Much better now! nica <- getData("GADM", country="NIC", level=0) nicabg <- openmap(nica, path="landscape") map.new(nicabg,1,mar=c(2, 2, 2, 0) + 0.1) plot(nicabg,axes=TRUE) plot(nica,add=TRUE) loc <- insetMap(nica,pos=c(-89.9,8.5), width=0.3, col=NA, lty=0,pch=".") points(loc) Just 2 questions:

Re: [R-sig-Geo] Multiplying a raster object

Thanks Loic, That helps a lot, actually. I used the wrong factor for transforming. Still adding the values of my raster don't add up to the points of the original ppp object. So I suppose there are better approaches than multiplying the raster by factor. If someone else wants to chip in on this