Robert,
I'm using
packageDescription(raster)
Package: raster
Type: Package
Title: Geographic analysis and modeling with raster data
Version: 0.9.9-14
Date: 8-Febrary-2010
(Maybe I should wait few more hours to get the fixed version?
The date above seems earlier than your message, but just to
Dear all,
I've plotted my data using the following
image.plot(data, zlim=c(0,500), add=FALSE, nlevel=64, horizontal=T,
legend.shrink=0.5, legend.width=.8, col=tim.colors(64))
However, I see white spaces in my image which I attribute them to values
that is 500 than that defined in zlim().
Hi!
After
r3 - readGDAL(test2.tif)
test2.tif has GDAL driver GTiff
and has 256 rows and 256 columns
image(r3,col=rainbow(128))
I get the image N-S reversed.
I've put an screenshot and the actual geotif image here:
http://sites.google.com/site/eospansite/dummy/test2_screenshot2.jpeg
On Wed, 10 Feb 2010, Agustin Lobo wrote:
Hi!
After
r3 - readGDAL(test2.tif)
test2.tif has GDAL driver GTiff
and has 256 rows and 256 columns
image(r3,col=rainbow(128))
I get the image N-S reversed.
I've put an screenshot and the actual geotif image here:
Thanks Sarah. But I think that package sp has an specific method of
image() for objects SpatialGridDataFrame, which is the class returned
by readGDAL(), and this specific method takes care of this problem. Check
image.SpatialGridDataFrame in package {sp}
Agus
2010/2/10 Sarah Goslee
Hi,
What if you make an spplot, still the same problem?
cheers,
Paul
Agustin Lobo wrote:
Thanks Sarah. But I think that package sp has an specific method of
image() for objects SpatialGridDataFrame, which is the class returned
by readGDAL(), and this specific method takes care of this
Yes, same problem
Agus
2010/2/10 Paul Hiemstra p.hiems...@geo.uu.nl:
Hi,
What if you make an spplot, still the same problem?
cheers,
Paul
Agustin Lobo wrote:
Thanks Sarah. But I think that package sp has an specific method of
image() for objects SpatialGridDataFrame, which is the class
FYI
This problem was actually identified because I was getting wrong
results at calculating
averaged raster values for intersecting line segments.
I include the last message sent by Robert. His trick is a good
workaround by now.
Agus
-- Forwarded message --
From: Robert J.
On Wed, 10 Feb 2010, Agustin Lobo wrote:
I think I tripped on this same stone, but cannot find anything in my
records, perhaps it's on the lost computer.
Sorry, can't find such a posting - maybe my search keys were wrong - you
could try the nabble archive.
The file test2.tif was made by
On Wed, 10 Feb 2010, Agustin Lobo wrote:
FYI
This problem was actually identified because I was getting wrong
results at calculating
averaged raster values for intersecting line segments.
I include the last message sent by Robert. His trick is a good
workaround by now.
And there is a similar
I forward to qgis devel and Marco Hugentobler.
I really cannot understand, Qgis uses gdal for writing and reading geotif files
as far as I know.
I think the point is that QGIS reports the origin in the NW corner and then uses
negative resolution for Y. But do not understand the discrepancy in
the
Dear researchers
I download the new version of R 2.10.1 and I have a problem to read maptools
Help as the older version
thanks for attention
Gianni
library(maptools)
Carico il pacchetto richiesto: foreign
Carico il pacchetto richiesto: sp
Carico il pacchetto richiesto: lattice
?maptools
No
Dear Reseachers
I have a problem to read the Help file of maptools and now I see I have a
problem to read.shape fucntion. I download the new version of R (2.10.1)
library(maptools)
x - read.shape(sample.shp, package=maptools)
Error: I don't find the function read.shape
it's the first time I
On Wed, 10 Feb 2010, gianni lavaredo wrote:
Dear researchers
I download the new version of R 2.10.1 and I have a problem to read maptools
Help as the older version
thanks for attention
Gianni
library(maptools)
Carico il pacchetto richiesto: foreign
Carico il pacchetto richiesto: sp
Carico
On Wed, 10 Feb 2010, gianni lavaredo wrote:
Dear Reseachers
I have a problem to read the Help file of maptools and now I see I have a
problem to read.shape fucntion. I download the new version of R (2.10.1)
library(maptools)
x - read.shape(sample.shp, package=maptools)
Error: I don't find the
dear Rearchers
I know this is easy but after 5 hours on R, I am honest and I don't
rememeber to get X, Y coordinate from point shapefile
sample.shp - readShapeSpatial(sample.shp)
proj4string(sample.shp) - CRS(+init=epsg:26911)
str(sample.shp)
str(sample.shp)
Formal class
xy - coordinates(sample.shp)
will give you a 2-column matrix xy of the X and Y coordinates.
For more complete reporting for this list, please provide the packages
you are using and the version. The easiest way is to use sessionInfo()
See here for more information, and the posting guide:
Hello,
On Thu, Feb 11, 2010 at 9:17 AM, gianni lavaredo
gianni.lavar...@gmail.com wrote:
Thanks Mike,
i have this question. Is It possible to get X,Y and have two new columns
inside the dbf of shape? as GET X and GET Y of ArcMap
I have no idea about any of ArcMap's functions, and I'm sure
dear R users,
there is same possibility to convert a points shape file in a polygon shape
file? where the corner are the last points
thanks
gianni
[[alternative HTML version deleted]]
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
Hello, this example came up in discussion today - this will retrieve
all the raw coordinates from a SpatialPolygonsDataFrame.
I cannot find any helper functions for this in maptools or sp, but if
there are some I'd like to know. This doesn't do anything to tag each
vertex in terms of its origin
On Wed, 10 Feb 2010, Wouter Buytaert wrote:
Hi list,
I found some strange behaviour in the latest version of spDistsN1:
spDistsN1(cbind(c(1,2,3,4),c(0,0,0,0)), c(0,0), longlat = TRUE)
[1] 0 0 0 0
spDistsN1(cbind(c(1,2,3,4),c(0,0,0,0)), c(0,0), longlat = FALSE)
[1] 1 2 3 4
It seems that
21 matches
Mail list logo