Re: [R-sig-Geo] writeRaster to ascii file (.asc)

2013-07-11 Thread Mauricio Zambrano-Bigiarini
On 10/07/13 18:38, Manuel Spínola wrote: Dear list members, I am trying to export a raster file to a ascii file (.asc) but the min and max values are not the same. Any idea why? Hi Manuel, Could you provide your sessionInfo() ? bio1res class : RasterLayer dimensions : 382, 407,

Re: [R-sig-Geo] help for show.vgms

2013-07-11 Thread Edzer Pebesma
library(gstat) my.show.vgms = show.vgms my.show.vgms = edit(my.show.vgms) my.show.vgms() what exactly you have to change during the edit phase is beyond the context of this list, it is basic R programming. On 07/10/2013 07:51 PM, Saman Monfared wrote: Dear all, How can I add text and change

Re: [R-sig-Geo] Degrees of freedom in GWR

2013-07-11 Thread Roger Bivand
On Tue, 9 Jul 2013, Michael Blows wrote: Hi all, can any body clarify me why the degrees of freedom in R output of ?GWR comes in decimal? Mike Please read the references, where this is explained - why would you expect them to be integer when weighted proportions of observations are being

[R-sig-Geo] Bayesian Networks-GIS

2013-07-11 Thread Rosa Soledad Manrique Paredes
Dear all, Do exist a way to connect GIS data layers to a Bayesian Network? For example to use raster files as evidence variables? Thanks for your suggestions in advance Rosa. [[alternative HTML version deleted]] ___ R-sig-Geo mailing list

[R-sig-Geo] Data from intersecting Lines of shp files

2013-07-11 Thread ElSiebo
Hi Iam quite new to R and Iam trying to use it for GIS application. I used the maptools package to load a shapefile WVA-readShapeLines(WVA) and now I would like to fill gaps in the shape file (missing numbers) with the numbers from intersecting lines I have tried with gintersect but I can only

Re: [R-sig-Geo] Data from intersecting Lines of shp files

2013-07-11 Thread Barry Rowlingson
On Thu, Jul 11, 2013 at 12:55 PM, ElSiebo franz.tscheikner-gr...@uibk.ac.at wrote: Hi Iam quite new to R and Iam trying to use it for GIS application. I used the maptools package to load a shapefile WVA-readShapeLines(WVA) 1. Use readOGR from the rgdal package to read shapefiles - it will

Re: [R-sig-Geo] Data from intersecting Lines of shp files

2013-07-11 Thread Barry Rowlingson
On Thu, Jul 11, 2013 at 1:27 PM, Tscheikner-Gratl, Franz franz.tscheikner-gr...@uibk.ac.at wrote: Thanks for your reply. The problem is that I have data of a pipe network with missing attributes (construction year). So I would like to see what are the constructions years of the intersecting

[R-sig-Geo] writeRaster() writing float despite stating INT2U

2013-07-11 Thread Agustin Lobo
I'm having this strange problem: I write a raster as INT2U but get a geotif with Float64. Is perhaps a gdal problem? r - raster(matrix(1:100,nrow=10,ncol=10)) writeRaster(r, file=r.tif,format=GTiff,dataype=INT2U,overwrite=TRUE,NAflag=0) rgdal: version: 0.8-10, (SVN revision 478) Geospatial

Re: [R-sig-Geo] Data from intersecting Lines of shp files

2013-07-11 Thread Roger Bivand
On Thu, 11 Jul 2013, Barry Rowlingson wrote: On Thu, Jul 11, 2013 at 1:27 PM, Tscheikner-Gratl, Franz franz.tscheikner-gr...@uibk.ac.at wrote: Thanks for your reply. The problem is that I have data of a pipe network with missing attributes (construction year). So I would like to see what are

Re: [R-sig-Geo] writeRaster() writing float despite stating INT2U

2013-07-11 Thread Dutrieux, Loic
Hi Augustin, Could it be because of the typo in the code (dataype)? Cheers, Loïc -Original Message- From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-boun...@r-project.org] On Behalf Of Agustin Lobo Sent: jeudi 11 juillet 2013 15:19 To: r-sig-geo Subject: [R-sig-Geo]

Re: [R-sig-Geo] writeRaster() writing float despite stating INT2U

2013-07-11 Thread Agustin Lobo
LLoic, Thanks and I apologize! Agus On Thu, Jul 11, 2013 at 4:46 PM, Dutrieux, Loic loic.dutri...@wur.nl wrote: Hi Augustin, Could it be because of the typo in the code (dataype)? Cheers, Loïc -Original Message- From: r-sig-geo-boun...@r-project.org

Re: [R-sig-Geo] R or R studio 3d manipulate

2013-07-11 Thread Greg Snow
Here is an example that works with the scatterplot3d function (from the package of the same name). I only found 2 arguments that control the projection, angle and scale.y: library(scatterplot3d) library(TeachingDemos) # from example 1 of scatterplot3d zz - seq(-10, 10, 0.01) xx -

[R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Michael Treglia
Hi All, I am using he raster package, and in particular the raster function. I am dealing with a number of layers, and it would be easiest to set the CRS upon import. Looking through the documentation, it seems that it is doable, and I have been using the following code to get started, though it

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Hodgess, Erin
Hi Michael: Did you try this; test-raster(us_ppt_1981_2010.02.asc, crs=+ellps=GRS80) with the plus sign, please? That might do it. Thanks, Erin From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] on behalf of Michael Treglia

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Michael Treglia
Hi Erin, Thanks for the quick response - I did try it with the plus sign, and same result. (I also tried it with 'crs=+proj=longlat +ellps=GRS80' with the same result too.) test-raster(us_ppt_1981_2010.02.asc, crs=+ellps=GRS80) test class : RasterLayer dimensions : 3105, 7025,

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Michael Treglia
I just tried running it in an older version of R, 2.13.0 (64-bit) and it ran with no problem... the exact same code run from Notepad++ Anybody know what has changed/if I need to adjust my for R 3.x.x? t-raster(us_ppt_1981_2010.02.asc, crs=+proj=longlat +ellps=GRS80) test class :

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Michael Sumner
Please report on the actual code you ran, you've just sent two messages, one with test-raster(us_ppt_1981_2010.02.asc, crs=+ellps=GRS80) this correctly cannot set the CRS and so it is NA and the other with t-raster(us_ppt_1981_2010.02.asc, crs=+proj=longlat +ellps=GRS80) test We have no idea

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Michael Treglia
Hi, I apologize - I must have mistakenly copied/pasted only the last t Here is how it was run, and the results, in R 3.0.0 (64-bit): test-raster(us_ppt_1981_2010.02.asc, crs=+proj=longlat +ellps=GRS80) test class : RasterLayer dimensions : 3105, 7025, 21812625 (nrow, ncol, ncell)

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Hodgess, Erin
Hi again. It looks like the crs argument no longer works when creating a raster from a file. I've tried all kinds of variations on the theme also, but no luck. Weird. Thanks Erin From: Michael Treglia [mtreg...@gmail.com] Sent: Thursday, July 11, 2013 8:21 PM

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

2013-07-11 Thread Michael Treglia
Thanks for trying Erin! I appreciate knowing its not just me :-) -Mike On Thu, Jul 11, 2013 at 6:40 PM, Hodgess, Erin hodge...@uhd.edu wrote: Hi again. It looks like the crs argument no longer works when creating a raster from a file. I've tried all kinds of variations on the theme also,

[R-sig-Geo] map topology problem and solution (off-topic)

2013-07-11 Thread Michael Sumner
Hello, This explanation and solution to a common problem is really excellent: http://bost.ocks.org/mike/example/ While it's not R, I thought it would be of general interest to many on this list, apologies for the noise otherwise. Cheers, Mike. -- Michael Sumner Hobart, Australia e-mail:

[R-sig-Geo] NMDS R2 goodness of fit statistic

2013-07-11 Thread Federico Tomasetto
Dear members, I hope the question make sense to you. Using envfit function in vegan package, is the goodness of fit statistic R2 based on the fit across all NMDS axes or just one? Many thanks, Federico [[alternative HTML version