On Tue, Mar 6, 2012 at 1:49 AM, Zia Ahmed z...@cornell.edu wrote:
Is there anyway to calculate population weighted centroids of a spatial
polygon in R? Any idea will be appreciated. Thanks
Some more clues would be appreciated! What's your data? A set of
points representing cities with
Hello everyone,
This is going to sound quite silly, but I'm an amateur R user and after a
while I haven't been able to find the solution for my problem...
I've got a shapefile coming from ArcMap 10, with a Projected Coordinate
System: World_Cylindrical_Equal_Area (I chose that one when I was
Hi!
Is there some file with same name of the shapefile, but with the prj
extension?
On Tue, Mar 6, 2012 at 10:08 AM, Ester Polaina Lacambra e.pola...@gmail.com
wrote:
Hello everyone,
This is going to sound quite silly, but I'm an amateur R user and after a
while I haven't been able to find
Hi,
I have several environmental rasters (loaded as stack/brick) in R which I want
to use as predictors for Species Distribution Models (preferably Maxent) with
the dismo-package. Some of the predictors are habitat specific variables and
are restricted to a river (like a line of raster cells,
Hi Johannes,
On 06.03.2012, at 14:45, Johannes Radinger wrote:
Hi,
I have several environmental rasters (loaded as stack/brick) in R which I
want to use as predictors for Species Distribution Models (preferably Maxent)
with the dismo-package. Some of the predictors are habitat specific
Hi everybody,
My problem is this : I have a raster layer of world temperature and I have
also a shapefile representing a region of the world. I want obtain a new raster
layer of world temperature without the values at region shapefile.
Thank for your help,
Camille Leclerc
Hi,
I need to load some data from the Worldclim Database and reproject them
to some other, already existing rasters. Because I don't want to use the
reprojected rasters from the memory I want to save the brick as a GeoTiff (like
the other rasters) using writeRaster. But somehow writeRaster()
On 05/03/12 14:40, Edzer Pebesma wrote:
This has now been changed/repared in sp on r-forge (svn), and will
appear from sp 0.9-97 on. A test run would be:
library(sp)
xyz = data.frame(expand.grid(x=1:10,y=1:10),rnorm(100))
coordinates(xyz)=~x+y
spplot(xyz, cex=10)
which used to plot points in
Hi Johannes,
did you transform your SpatialPointsDataFrame to a regular data.frame?
point.df=data.frame(point.spdf)
subset(points.df, points.df[,value] 0 , select = c(1, 2))
coordinates(points.df) = ~X + Y ##You should keep those columns (X,Y) in
the subset command in order to backtransform the
Yes, look for argument alpha e.g. in ?rgb or ?bpy.colors
On 03/06/2012 04:49 PM, Mauricio Zambrano-Bigiarini wrote:
On 05/03/12 14:40, Edzer Pebesma wrote:
This has now been changed/repared in sp on r-forge (svn), and will
appear from sp 0.9-97 on. A test run would be:
library(sp)
xyz =
Hi Erin,
I don't know if you fix your problem yet but I had the same error yesterday
and the problem is that ArcGIS do not accept the caracter . in a field
name..
hope this helps
Sam
On Thu, Mar 1, 2012 at 9:38 PM, Hodgess, Erin hodge...@uhd.edu wrote:
Hi R Sig Geo People:
I'm having some
Hi,
Original-Nachricht
Datum: Tue, 6 Mar 2012 10:53:11 -0500
Von: Samuel Turgeon sturgeo...@gmail.com
An: Johannes Radinger jradin...@gmx.at
CC: r-sig-geo@r-project.org
Betreff: Re: [R-sig-Geo] subset() for spatialPointsDataFrame
Hi Johannes,
did you transform your
On 06/03/12 16:55, Edzer Pebesma wrote:
Yes, look for argument alpha e.g. in ?rgb or ?bpy.colors
Thank you very much Edzer. It works perfectly:
library(sp)
xyz = data.frame(expand.grid(x=1:10,y=1:10),rnorm(100))
coordinates(xyz)=~x+y
spplot(xyz, cex=10, alpha=0.7)
Cheers,
Mauricio
--
Am 06.03.2012 16:53, schrieb Samuel Turgeon:
Hi Johannes,
did you transform your SpatialPointsDataFrame to a regular data.frame?
point.df=data.frame(point.spdf)
subset(points.df, points.df[,value] 0 , select = c(1, 2))
Indexing with the correct number of dimensions as done here (in contrast
Johannes, it works for me:
library(raster)
crs - +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=350 +y_0=0
+datum=potsdam +units=m +no_defs +ellps=bessel +towgs84=606.0,23.0,413.0
r - raster(nc=10, nr=10, crs=crs)
r[] - 1:ncell(r)
b - brick(r,r,r)
x - writeRaster(b, 'datum.tif', overwrite=TRUE)
x
Hi
I am getting the following error when trying to read and stack a .nc file.
PJ - stack(2030-2060-Tmax.nc,varname=TASMAX)
Warning message:In .stackCDF(x, varname = varname, bands = bands) :
TASMAXhas 4 dimensions, I do am using the last one
PJ
[1] TASMAXhas 4 dimensions, I do am using
Sajid,
You are not getting an error (it is a warning), and you are not saying what
you want to achieve, that makes it hard to say how you can handle this?
For efficiency, you should probably use 'brick' rather than 'stack' when
using a single ncdf file. You can use the can use lvar argument to
Yes perfect thanks.
-Original Message-
From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-boun...@r-project.org]
On Behalf Of Edzer Pebesma
Sent: Tuesday, 6 March 2012 7:56 p.m.
To: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] adjusting and moving scale
Kanika, maybe something
Hi all R gurus,
I'm trying to produce simple returns on a series of market prices but am
receiving zeros, even when using the formula directly. The timeSeries
function returns is the only method, of the methods that I've tried,
that works for me. My preference is to use PerformanceAnalytics
I feel like I have done this before, so I am sorry if I have made a stupid
mistake, however, I am out of ideas.
I am trying to make a shapefile out of a list of polygons. I want to
include the a data frame as well, but I can't seem to create a
SpatialPolygons-Class longer than one object. The one
Dear Robert
Thanks for the reply.. Yes I should have explained more. Actually I want to
subset certain bands representing seasonal months to create rasters and
followed by summary stats.
I was not able to subset with those warnings.
But your tip on using brick instead of stack worked and now
21 matches
Mail list logo
