Thanks Roger, you are right, I had missed that part (don't ask me how).
projectRaster() is indeed the function in raster to use (but it is a bit slow).
Something like
r<-raster(paste(asc.in.dir, "\\", asc.in.files[1], sep=""))
projection(r) <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0
Thank you Roger. Using
do.call("spRbind", myListOfSpatialPointsObjects)
do not work, but I can work around.
Renaud
2009/5/8 Roger Bivand
> On Fri, 8 May 2009, Renaud Lancelot wrote:
>
> Thank you, here is the full process trying to follow Edzer's suggestion:
>>
>> library(RODBC)
>> library(r
On Fri, 8 May 2009, Robert Hijmans wrote:
Tim,
"proj" is not a valid argument in raster(x, ... ) when x (the first
argument) is of class character (i.e. interpreted as a filename); but
I can add that argument. For now, what should work is:
r<-raster(paste(asc.in.dir, "\\", asc.in.files[i], sep
The raster package depends on rgdal and sp. It is under development in
R-forge (that is, it is not on CRAN yet). You can download the source
or install it directly on Win and Linux with this command:
install.packages("raster", repos="http://R-Forge.R-project.org";)
Robert
On Fri, May 8, 2009 at 1
On Fri, 8 May 2009, Renaud Lancelot wrote:
Thank you, here is the full process trying to follow Edzer's suggestion:
library(RODBC)
library(rgdal)
channel <- odbcConnectAccess("farms.mdb")
sigal <- sqlQuery(channel, "select * from T_SIGAL")
Liste <- by(sigal,
list(zone = sigal$zone),
Hello.
Raster package seems to be very interesting since it
doesn't require external applications.
However currently, no binaries prepared for Win32.
Will this Win32 binary available?
Regards.
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.et
Tim,
"proj" is not a valid argument in raster(x, ... ) when x (the first
argument) is of class character (i.e. interpreted as a filename); but
I can add that argument. For now, what should work is:
r<-raster(paste(asc.in.dir, "\\", asc.in.files[i], sep=""))
projection(r) <- "+proj=tmerc +lat_0=0
I need to take an ESRI grid (.asc) which is in geographic coordinates
(lat/lon) and project it to UTM coordinates. The following command is
giving me an error saying I'm giving it an undefined argument. I can read
the raster without the projs arguement, but I'm stuck when adding that on.
> r<-ra
Thank you, here is the full process trying to follow Edzer's suggestion:
library(RODBC)
library(rgdal)
channel <- odbcConnectAccess("farms.mdb")
sigal <- sqlQuery(channel, "select * from T_SIGAL")
Liste <- by(sigal,
list(zone = sigal$zone),
function(x){
strg
Its all quite simple I think (i.e. "hope"). Have
you used spTransform(rdgal)?
## here I have a file from a local group delivered
## in UTM format that I transformed OK to long-lat.
library(rgdal) # for map projection support; automatically loads sp
## transform internal measurements from UTM to
Thank you Edzer. I was hoping there was a simpler solution, but it's still
fairly simple.
Renaud
2009/5/7 Edzer Pebesma
> Yes, but you will not like that call.
>
> In R (package sp), each spatial data set can only have one single
> projection. I'd suggest, in pseudo-code
>
> - import the data i
Thank you Milton. I have already explored this source, without success.
Renaud
2009/5/7 milton ruser
> Hi Renaud,
>
> may be you find some solution reading last posts.
>
>
> http://search.r-project.org/cgi-bin/namazu.cgi?query=utm+to+latlong&max=20&result=normal&sort=score&idxname=Rhelp02a&idxn
Yes, but you will not like that call.
In R (package sp), each spatial data set can only have one single
projection. I'd suggest, in pseudo-code
- import the data into a data.frame
- for each UTM zone do:
select the locations
set their CRS
reproject to long/lat
- cbind the locations to a
Hi Renaud,
may be you find some solution reading last posts.
http://search.r-project.org/cgi-bin/namazu.cgi?query=utm+to+latlong&max=20&result=normal&sort=score&idxname=Rhelp02a&idxname=functions&idxname=docs
cheers
milton
On Thu, May 7, 2009 at 5:11 PM, Renaud Lancelot wrote:
> Dear all,
>
>
Dear all,
I have a big dataset (> 200,000 lines) of georeferenced locations where
the coordinates are given in UTM, spanning over several zones. Is is
possible to convert all the UTM coordinates into longlat WGS84 with a
single call ?
Renaud
--
Renaud Lancelot
EDEN Project, coordinator
http:
I would like to change the color of a binary mask when I plot it. The
default is red and beige. I've tried setting col parameter, but it is
only working with polygonal windows.
data(demopat)
plot(demopat$window, col=3)
#Not OK :
plot(as.mask(demopat$window), col=3)
Thanks,
Etienne
[[al
Dear R-mappers and analyzers of spatial data,
just wanted to bring your attention to a new package on CRAN, RgoogleMaps,
which allows you to plot on Google map tiles in R.
The package serves two purposes:
(i) Provide a comfortable R interface to query the Google server for static
maps
(ii) Use the
Dear R-sig-geo users,
I have the output of a watershed segmentation in vector format (shapefile)
which has it's attribute table populated with statistics regarding spectral
reflectance of each polygon object. The attribute data was sourced from a
geographically co-incident aerial photograph. I
18 matches
Mail list logo