[R-sig-Geo] Re-projection with project raster results in rotated plot

2014-09-29 Thread Swagath Navin Manohar
Hello all, I have a problem concerning re-projecting a raster from EPSG:3035 to sr-org: 29, wrf-lambert-conformal-conic which results in a rotated plot. Here is the code I used: raster1=rasterFromXYZ(radon[, c("x", "y", "rn")]) http://spatialreference.org/ref/epsg/3035/proj4/ crs(raster

[R-sig-Geo] TIF with 3 bands to data frame

2013-05-03 Thread Swagath Navin
Hello all, I am new to rgdal. I downloaded a tif file from 'http://mrdata.usgs.gov/radiometric/data/NArad_U_geog83.tif' which i would like to convert into a regular data frame so that i can plot it using ggplot2. when i tried to read the file using readGDAL i get a warning message NArad_U_g

Re: [R-sig-Geo] reading .grd files

2013-04-30 Thread Swagath Navin
Thank you very much Barry. I converted the .grd file to a .flt file. Now i can work with the map. Also your blog at http://geospaced.blogspot.co.uk/2011/12/reading-afripop-data.html is very informative. Thank you very much for your help. Swagath On Mon, Apr 29, 2013 at 3:08 PM, Barry Rowlingson

[R-sig-Geo] gdalwarp gives different cellsize when re-projected

2013-03-07 Thread Swagath Navin
Dear all, I have a large data-frame with Lon, Lat , Value . Coordinates projected in EPSG:3035 with a resolution of 5 kmx5 km. I converted this into a SGDF. My goal was to project the file into EPSG:4326 with a resolution of 0.05. For re-projecting i converted the file to tif using writeGDA

[R-sig-Geo] creating a SGDF with transformed coordinates and spatial aggregation

2013-03-05 Thread Swagath Navin
Sorry..i am resending my mail again ..first mail i sent as HTML format Dear all, I have a large data frame in form of (Lon, Lat, Value) ((coordinates in EPSG:3035). I can convert the data frame into a object of class SGDF using the following commands: coordinates(df)=~x+y proj4string(df)=CRS(

[R-sig-Geo] creating a SGDF with transformed coordinates and spatial aggregation

2013-03-05 Thread Swagath Navin
Dear all, I have a large data frame in form of (Lon, Lat, Value) ((coordinates in EPSG:3035). I can convert the data frame into a object of class SGDF using the following commands: coordinates(df)=~x+y proj4string(df)=CRS("+init=epsg:3035") gridded(pts) = TRUE fullgrid(pts)=TRUE getGridTopolog

Re: [R-sig-Geo] create data frame after spDists (re-posted)

2012-08-27 Thread Swagath Navin
r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] on behalf of Rolf Turner [r.tur...@auckland.ac.nz] Sent: 24 August 2012 21:42 To: swagath navin Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] create data frame after spDists (re-posted) I believe that you would find the to

[R-sig-Geo] create data frame after spDists (re-posted)

2012-08-24 Thread swagath navin
Hello all, sorry for re-posting, but i got a scrub message for my first mail. I have two spatial points data frame df1 showing sample locations and df2 meteorological stations. I am interested in finding met stations closer to the sample stations and creating a separate data frame with sample loca

[R-sig-Geo] create data frame after spDists

2012-08-24 Thread swagath navin
Hello all, I have two spatial points data frame df1 showing sample locations and df2 meteorological stations. I am interested in finding met stations closer to the sample stations and creating a separate data frame with sample locations and the nearest met station. I calculated the distances betw