Michael Sumner wrote:
Here's one way. I'm not sure about how to control the aspect ratio
(something I've been meaning to check for a while now).
somewhat of a non-newby answer to Michael: asp is set in plot.Spatial,
and if the default value of NA is passed, set to:
if (is.na(asp)) asp
Here's one way. I'm not sure about how to control the aspect ratio
(something I've been meaning to check for a while now).
library(sp)
data(volcano)
x <- 10*(1:nrow(volcano))
y <- 10*(1:ncol(volcano))
## image xyz list
imlist <- list(x = x, y = y, z = volcano)
## SpatialGridDataFrame (image2Grid
At the suggestion of Walmes Zeviani I downloaded and am taking a look at the sp
package.
In an attempt to get going, can someone offer a suggestion about how to convert
the below over to a version that uses the sp package?
data(volcano)
x <- 10*(1:nrow(volcano))
y <- 10*(1:ncol(volcano))
imag
On Mon, 4 Jan 2010, Peter Adler wrote:
Dear R-sig-geo,
I am new to RGDAL and PROJ4, and am having trouble projecting a raster
dataset.
The dataset is a 721 x 721 cell raster of snow cover from the NSIDC. The
projection information is here:
http://nsidc.org/data/atlas/epsg_3408.html
My com
Dear R-sig-geo,
I am new to RGDAL and PROJ4, and am having trouble projecting a raster
dataset.
The dataset is a 721 x 721 cell raster of snow cover from the NSIDC. The
projection information is here:
http://nsidc.org/data/atlas/epsg_3408.html
My command:
D=readGDAL("annual_snow_trend_ease
Hi,
Check out the spatial task view [1] and the R wiki [2].
cheers,
Paul
[1] http://cran.r-project.org/web/views/Spatial.html
[2] http://wiki.r-project.org/rwiki/doku.php?id=tips:spatial-data
GRIGNOLA, FERNANDO E [AG/6042] wrote:
Hello,
I'm working on a project using weather data mostly on t
Hello,
I'm working on a project using weather data mostly on the statistical
aspects of it, i.e. environmental classification. I have a bunch of
shape files (*.dbf, *.prj, *.sbn, *.sbx, *.shp, *.shx) from different
world regions and I would like to start using R to do the mapping of
different clust
Ashton Shortridge wrote:
On Monday 16 November 2009 08:53:03 Carson Farmer wrote:
... In addition, I use linux (Debian) and the desktop GIS tools that
I tried (svSig, Qgis, and some more) did not provide publication quality
maps.
I would strongly suggest you take another look at some
On Monday 16 November 2009 08:53:03 Carson Farmer wrote:
> >... In addition, I use linux (Debian) and the desktop GIS tools that
> > I tried (svSig, Qgis, and some more) did not provide publication quality
> > maps.
>
> I would strongly suggest you take another look at some of these
> packages (es
Roger Bivand wrote:
...
> Including a code snippet demonstrating your problem and posting a link
> to a simple data set (shape file plus election data) would be a more
> constructive way to proceed, no? The code in Ch 5 in our book on
> www.asdar.org covers a very similar case,
...
The book lin
On Sun, 15 Nov 2009, tom sgouros wrote:
Roger Bivand wrote:
I'm trying to make some simple maps of election results and I am having
some trouble. I'm hoping to use Grass (or QGis) to make the maps, but
apparently this requires getting my data into the shape files, and I
can't figure out ho
On Mon, 16 Nov 2009, Karl Ove Hufthammer wrote:
On Sun, 15 Nov 2009 19:52:48 +0100 (CET) Roger Bivand
wrote:
Do use readShapeSpatial() (or better readOGR() in rgdal)
Why (or how) is readOGR better?
Because it uses OGR, with an updated rather than frozen shapelib, and
reads the *.prj file
>... In addition, I use linux (Debian) and the desktop GIS tools that
> I tried (svSig, Qgis, and some more) did not provide publication quality
> maps.
I would strongly suggest you take another look at some of these
packages (especially QGIS).
These are relatively young projects, and are growing
Hello all,
An alternative to create (possibly simple) shapefiles is the
"shapefiles" package.
http://cran.r-project.org/web/packages/shapefiles/index.html.
Also the "PBSmapping" package is a very good alternative to make maps
within R, and to perform basic GIS operations (i.e. overlap of poly
Raphael Saldanha wrote:
Hi Tom!
I'm a geographer and R enthusiast too, but what I can say is: R is a great
software to statistics. The final production of maps, although it's
possible, it's simple and faster made with GIS softwares, like GRASS,
Quantum GIS or ArcGIS.
Hi,
As you already stat
On Sun, 15 Nov 2009 19:52:48 +0100 (CET) Roger Bivand
wrote:
> Do use readShapeSpatial() (or better readOGR() in rgdal)
Why (or how) is readOGR better?
--
Karl Ove Hufthammer
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.c
tom sgouros wrote:
> Alex Mandel wrote:
>
>> This could be done easily in either QGIS or R. I would leave GRASS out
>> of the discussion for now if we're just talking about making a map and
>> bring it back in if you start wanting to analyze the spatial
>> relationship between your data and other
tom sgouros wrote:
> tom sgouros wrote:
>
>> Dan Putler wrote:
>
>>> Er, what exactly is your data? Do you have data that is linked to zip
>>> codes and you want to do a "thematic" (choropleth) map of the data? Do
>>> you have point data (say a lat/lon values), and want to plot the points?
>>>
I'd go with the shapefiles, the E00 files are ArcInfo export format
files, and it is very likely you would convert to shapefile sets along
the way.
Dan
On Sun, 2009-11-15 at 23:27 -0500, tom sgouros wrote:
> tom sgouros wrote:
>
> >
> > Dan Putler wrote:
>
> > > Er, what exactly is your data
tom sgouros wrote:
> Dan Putler wrote:
>
>> Tom,
>>
>> Er, what exactly is your data? Do you have data that is linked to zip
>> codes and you want to do a "thematic" (choropleth) map of the data? Do
>> you have point data (say a lat/lon values), and want to plot the points?
>> Given what you've s
tom sgouros wrote:
>
> Dan Putler wrote:
> > Er, what exactly is your data? Do you have data that is linked to zip
> > codes and you want to do a "thematic" (choropleth) map of the data? Do
> > you have point data (say a lat/lon values), and want to plot the points?
> > Given what you've said
Dan Putler wrote:
> Tom,
>
> Er, what exactly is your data? Do you have data that is linked to zip
> codes and you want to do a "thematic" (choropleth) map of the data? Do
> you have point data (say a lat/lon values), and want to plot the points?
> Given what you've said, hard to figure out how
Tom,
Er, what exactly is your data? Do you have data that is linked to zip
codes and you want to do a "thematic" (choropleth) map of the data? Do
you have point data (say a lat/lon values), and want to plot the points?
Given what you've said, hard to figure out how to get you going in the
right di
Raphael Saldanha wrote:
> Hi Tom!
>
> I'm a geographer and R enthusiast too, but what I can say is: R is a great
> software to statistics. The final production of maps, although it's
> possible, it's simple and faster made with GIS softwares, like GRASS,
> Quantum GIS or ArcGIS.
This is exactl
Hi Tom!
I'm a geographer and R enthusiast too, but what I can say is: R is a great
software to statistics. The final production of maps, although it's
possible, it's simple and faster made with GIS softwares, like GRASS,
Quantum GIS or ArcGIS.
My tip is: organize data with spreadsheet and GIS sof
Roger Bivand wrote:
> > I'm trying to make some simple maps of election results and I am having
> > some trouble. I'm hoping to use Grass (or QGis) to make the maps, but
> > apparently this requires getting my data into the shape files, and I
> > can't figure out how to do that. I've found re
On Sun, 15 Nov 2009, Tom Sgouros wrote:
Hello all:
I'm trying to make some simple maps of election results and I am having
some trouble. I'm hoping to use Grass (or QGis) to make the maps, but
apparently this requires getting my data into the shape files, and I
can't figure out how to do that
Hello all:
I'm trying to make some simple maps of election results and I am having
some trouble. I'm hoping to use Grass (or QGis) to make the maps, but
apparently this requires getting my data into the shape files, and I
can't figure out how to do that. I've found read.shape and that's
great,
AIL PROTECTED]>
Subject: [R-sig-Geo] Newbie question: Reclassifying raster (ndvi)
using fisher method
To: [EMAIL PROTECTED], "r-sig-geo@stat.math.ethz.ch"
Message-ID:
<[EMAIL PROTECTED]>
Content-Type: text/plain
Hi all,
I am not sure as to which list
Hi all,
I am not sure as to which list this best fits (Grass or R-sig-geo).
I am trying to reclass a raster map (ndvi) into 5 different zones using the
fisher method, ultimatley i am trying to reclass as natural breaks.
The point of this exercise is for precision agriculture. I want to determine
Oehler, Friderike (AGPP) wrote:
> Thanks Edzer, I shall try again:
>
> 1) I would like to use spplot to map my factor "TYP" as dots of different
> colours, however the resulting plot uses the same colours for the first and
> last value (10,40). I guess that my use of the "cuts" argument is wrong, b
Thanks Edzer, I shall try again:
1) I would like to use spplot to map my factor "TYP" as dots of different
colours, however the resulting plot uses the same colours for the first and
last value (10,40). I guess that my use of the "cuts" argument is wrong, but
I can't find any better:
LAT <-c(-6.
Oehler, Friderike (AGPP) wrote:
> Dear list members,
>
> I am still struggeling with spplot. Sorry for being that unprofessional, but
> honestly I find the spplot help page quite hard to understand. (Many thanks
> Edzer for the very helpful hint a few days ago.)
>
> 1) My first plot actually works
Dear list members,
I am still struggeling with spplot. Sorry for being that unprofessional, but
honestly I find the spplot help page quite hard to understand. (Many thanks
Edzer for the very helpful hint a few days ago.)
1) My first plot actually works all fine, but the colours seem to loop, i.e.
Roger Bivand wrote:
> It may also be worth asking on the PostGIS list whether you could upload
> your geometries to PostGIS and have GEOS do the measurement - if you have
> access to the software and experience with it. It may be that it can
> handle geographical coordinates.
>
I have to be
On Sat, 16 Jun 2007, David Hugh-Jones wrote:
> Thank you for this help. I realised after writing this that indeed,
> spatstat won't work because I have world-scale data in unprojected
> latitude/longitude form... so my distance measure is no use. This
> leads to more questions
>
> - is there a can
Thank you for this help. I realised after writing this that indeed,
spatstat won't work because I have world-scale data in unprojected
latitude/longitude form... so my distance measure is no use. This
leads to more questions
- is there a canonical way to convert projections from within "sp"? My
cu
On Fri, 15 Jun 2007, David Hugh-Jones wrote:
> Hi R-spatialists,
>
> I'm working with a big list of polygons (imported into R with package
> sp as a SpatialPolygonDataFrame, and I am trying to calculate minimum
> distances between each pair of polygons. For this, I need to get the
> lines for each
Hi R-spatialists,
I'm working with a big list of polygons (imported into R with package
sp as a SpatialPolygonDataFrame, and I am trying to calculate minimum
distances between each pair of polygons. For this, I need to get the
lines for each polygons. At the moment I can't seem to find a way to
do
Hallo,
just a couple hints:
> In other words, I would like to know if there is a relationship between yield
> and phosphorus, taking into account the spatial dependence. But how can I
> test this?
I think you should look at the cross-variograms. Or did you already
reject this option?
> Anothe
Montclair State University
Montclair, NJ, 07043
Tel: 973-655-4313
Fax: 973-655-4072
email: [EMAIL PROTECTED]
- Original Message -
From: "GRETA G. GRAMIG" <[EMAIL PROTECTED]>
Date: Wednesday, December 20, 2006 11:57 am
Subject: [R-sig-Geo] Newbie Questions
> I have lim
I have limited experience with spatial analysis, but I am trying to learn.
Perhaps someone with patience for a newbie will field my questions? If this is
not the appropriate forum for such questions perhaps someone could suggest a
more appropriate forum.
I have a small dataset that consists of
On Fri, 3 Nov 2006, Jeffrey Stratford wrote:
> Hi.
>
> Does anyone have an example of estimating Moran's I from a file with
> point (not polygon) data? I'm getting close(r) but I still have a few
> problems.
>
> Here's my code so far:
>
> library(spdep)
> ctpoints <- read.csv("c:\\Georgia\\xy.
Hi.
Does anyone have an example of estimating Moran's I from a file with
point (not polygon) data? I'm getting close(r) but I still have a few
problems.
Here's my code so far:
library(spdep)
ctpoints <- read.csv("c:\\Georgia\\xy.csv", header=TRUE)
# only two columns containing the long-lat data
Hi!
My name is Raphael Saldanha, I'm a student in Brazil.
I'm using R by the first time now (and maybe my english too! hehehe), it's
really a very good program for spatial data analysis. But I have some
questions:
1. How can I discard some variables using princomp with a dataframe?
2. Watch is t
Jeff:
SPDEP might be the package where you can find numerous functions and tests for
your data.
Hope this helps.
Danlin
Jeffrey Stratford <[EMAIL PROTECTED]>:
> Greetings! I would appreciate some feedback as to where you think a
> newbie should start when learning spatial stats with R.
>
>
Greetings! I would appreciate some feedback as to where you think a
newbie should start when learning spatial stats with R.
I have avian point count data and landscape data (Landsat) over a 70 by
40 km area that I would like to start analyzing. I have already,
developed some relationships betw
Just an addition to the question I have asked before:
one easy way would be to fit:
- a global constant (to account for the average of the points of the
matrix);
- a plane to account for linear terms (in the example below, it should be a
plane with both negative partial derivatives with respect "
Hello,
I am new in this list. I have a problem and I think that the solution to
this problem can be found using spatial statistics, but I do not know how.
My problem is the following: I have a matrix of data, of around 6 rows and 5
columns. In reality, I have a time series, where at each point i
49 matches
Mail list logo
