Re: [R-sig-Geo] reading .grd file in R?

2010-10-08 Thread Roman Luštrik
Have you tried raster package? On Sat, Oct 9, 2010 at 1:32 AM, wrote: > > > Can someone let me know what package and function should be used to read > *.grd files in R? I tried readGDAL, open.ncdf functions already which threw > errors for file format mismatch. > > Thank you all, > Maha >

[R-sig-Geo] reading .grd file in R?

2010-10-08 Thread govindas
Can someone let me know what package and function should be used to read *.grd files in R? I tried readGDAL, open.ncdf functions already which threw errors for file format mismatch.  Thank you all, Maha [[alternative HTML version deleted]] ___

Re: [R-sig-Geo] World map in magnetic coordinates

2010-10-08 Thread Clint Bowman
Why not plot data and map in original coordinates and then overlay with grid lines showing magnetic coordinates? -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE:

Re: [R-sig-Geo] World map in magnetic coordinates

2010-10-08 Thread Barry Rowlingson
On Fri, Oct 8, 2010 at 7:21 PM, ogbos okike wrote: > library(maps) >  data(worldMapEnv) >  map("world") > The three steps above in R gives you the world map I mean. if you do: w = map("world") then w will be a simple list with $x and $y components, as well as $range and $names. If you write

[R-sig-Geo] transparent png images for map overlays

2010-10-08 Thread eick
I have created a SpatialPoints object, converted it to a ppp object, and then run spatialstat's density function create a density object. I can use the png function to make a density plot. I would like the density plot to be transparent so that I can include it as a layer on a map. Is there any

Re: [R-sig-Geo] World map in magnetic coordinates

2010-10-08 Thread ogbos okike
On 8 October 2010 19:52, Barry Rowlingson wrote: > On Fri, Oct 8, 2010 at 5:32 PM, ogbos okike wrote: > > > My data points are given in latitude and longitude (geographic). Using > AACGM > > geomagnetic coordinate software, I converted them to magnetic > coordinates. > > World map (such as implem

Re: [R-sig-Geo] Question about back-transformation(geoR)

2010-10-08 Thread Paulo Justiniano Ribeiro Jr
Samuel No, it is not exp(prediction on the normal scale) by exp(porediction in the normal scale + 0.5*variance or the prediction in the normal scale) geoiR folows lcosely what is describe in Diggle and Ribeiro (2007) The Section 3.8 (transformed Gaussian models) it describes the basis for th

Re: [R-sig-Geo] World map in magnetic coordinates

2010-10-08 Thread Barry Rowlingson
On Fri, Oct 8, 2010 at 5:32 PM, ogbos okike wrote: > My data points are given in latitude and longitude (geographic). Using AACGM > geomagnetic coordinate software, I converted them to magnetic coordinates. > World map (such as implemented in R, which can be generated using maps has > no data poi

[R-sig-Geo] Question about back-transformation(geoR)

2010-10-08 Thread Samuel Turgeon
Hi, I'm working with the function krige.conv of the library geoR. I'm wondering how the predicted result are back-transformed when the lambda parameter is set to 0 (log-transformation)? Is it a simply transformation like exp(result$prediction)? I saw some code example on the web where the exp() fu

Re: [R-sig-Geo] World map in magnetic coordinates

2010-10-08 Thread ogbos okike
Hi Barry, On 8 October 2010 18:05, Barry Rowlingson wrote: > On Fri, Oct 8, 2010 at 4:42 PM, ogbos okike wrote: > > > However, I am interested in organizing my data in magnetic coordinates. I > > have converted my data points from geographic coordinates to magnetic > > coordinates. > > > My pro

Re: [R-sig-Geo] World map in magnetic coordinates

2010-10-08 Thread Barry Rowlingson
On Fri, Oct 8, 2010 at 4:42 PM, ogbos okike wrote: > However, I am interested in organizing my data in magnetic coordinates. I > have converted my data points from geographic coordinates to magnetic > coordinates. > My problem now is to generate world map in magnetic coordinates. > Can anybody p

[R-sig-Geo] World map in magnetic coordinates

2010-10-08 Thread ogbos okike
Dear All, I am using raster to do my geographic plots. I also use maps package to generate world map in geographic coordinates. However, I am interested in organizing my data in magnetic coordinates. I have converted my data points from geographic coordinates to magnetic coordinates. My problem

Re: [R-sig-Geo] spplot: 1) filling background polygon without overwriting foreground polygon; 2) left justifying sp.layout text

2010-10-08 Thread Lyndon Estes
Dear Edzer, Very many thanks for the rapid response. What you posted is exactly what I was hoping for! I will download the updated sp package and reproduce my maps. Thanks again! Best, Lyndon On Fri, Oct 8, 2010 at 10:51 AM, Edzer Pebesma wrote: > Thanks for the reproducable example. I can s

[R-sig-Geo] problem to read ASCII file in Raster or readAsciiGrid (different stretched value and classified value)

2010-10-08 Thread gianni lavaredo
dear Researchers, I have a ASCII file with classified value (ex: 0-28) and stretched value (0-29871) in Arcmap. When I try to read read with Raster or readAsciiGrid the final result is with stretched value. Some people know some solution? thanks in advance Gianni [[alternative HTML vers

[R-sig-Geo] spplot: 1) filling background polygon without overwriting foreground polygon; 2) left justifying sp.layout text

2010-10-08 Thread Lyndon Estes
Greetings, I have two questions related to spplot: 1. I am trying to plot one SpatialPolygonsDataFrame over another, and fill the background polygon with grey. I can't seem to fill the background polygon without having it cover over my supposed foreground polygon, even if I specify "first = T" (w

Re: [R-sig-Geo] points, raster

2010-10-08 Thread steven mosher
see pointsToRaster() On Fri, Oct 8, 2010 at 2:12 AM, Arnald Marcer wrote: > Hi, > > Can anybody help me a little on how to transfer values from a data frame > (coord.x, coord.y, value) to > a raster file. I have a raster initialized to 0 values and I want to > substitute zeros for the specifi

[R-sig-Geo] points, raster

2010-10-08 Thread Arnald Marcer
Hi, Can anybody help me a little on how to transfer values from a data frame (coord.x, coord.y, value) to a raster file. I have a raster initialized to 0 values and I want to substitute zeros for the specific values at coord.x and coord.y Thank you very much, Arnald Marcer CREAF __