Re: [R-sig-Geo] creating a geotiff - and a new question

2009-08-13 Thread Monica Pisica
: roger.biv...@nhh.no > To: pisican...@hotmail.com > CC: b.rowling...@lancaster.ac.uk; r-sig-geo@stat.math.ethz.ch > Subject: Re: [R-sig-Geo] creating a geotiff - and a new question > > On Thu, 13 Aug 2009, Monica Pisica wrote: > >> >> WOW! thanks so much for the code. T

Re: [R-sig-Geo] creating a geotiff

2009-08-13 Thread Monica Pisica
Thanks! I will certainly try it. Monica > Date: Thu, 13 Aug 2009 18:14:21 +0200 > From: roger.biv...@nhh.no > To: pisican...@hotmail.com > CC: r-sig-geo@stat.math.ethz.ch > Subject: Re: [R-sig-Geo] creating a geotiff > > On Thu

Re: [R-sig-Geo] creating a geotiff - and a new question

2009-08-13 Thread Monica Pisica
o] creating a geotiff > From: b.rowling...@lancaster.ac.uk > To: pisican...@hotmail.com > CC: r-sig-geo@stat.math.ethz.ch > > On Thu, Aug 13, 2009 at 3:10 PM, Monica Pisica wrote: >> >> Thanks Barry! >> >> I will look into it - for now i have my clumsy code and bec

Re: [R-sig-Geo] creating a geotiff

2009-08-13 Thread Monica Pisica
ig-Geo] creating a geotiff > From: b.rowling...@lancaster.ac.uk > To: pisican...@hotmail.com > CC: r-sig-geo@stat.math.ethz.ch > > On Thu, Aug 13, 2009 at 1:35 PM, Monica Pisica wrote: >> >> Hi, >> >> I have some data as it is required for the function "image&

[R-sig-Geo] creating a geotiff

2009-08-13 Thread Monica Pisica
Hi, I have some data as it is required for the function "image" {graphics}, that means a vector x for x coord, a vector y for y coord and a matrix m with x rows and y cols. I want to save this as a geotiff. If i transform my data to take the form of a 3 column data.frame i can use gridded an

Re: [R-sig-Geo] read NITF save as geotif

2009-05-20 Thread Monica Pisica
eadGDAL('yourfile') > writeGDAL(imagery, 'destinationfile') > > GDALinfo('destinationfile') > > Hope this helps, > > Arien > > > On Tue, May 19, 2009 20:24, Monica Pisica wrote: >> >> Hi, >> >> I have some NITF (in lat

[R-sig-Geo] read NITF save as geotif

2009-05-19 Thread Monica Pisica
Hi, I have some NITF (in lat long) files and i would like to save them as geotif. Things would be almost straightforward but the pixels in the original image are rectangular, and not square, so i cannot grid. Do you have any idea how i can go around this? Thanks, i will really apprec

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

2009-01-28 Thread Monica Pisica
df(nc, nc$var[[6]])) > close.ncdf(nc) > > xm <- xr[1] - 0.5 * res[1] > xx <- xr[2] + 0.5 * res[1] > ym <- yr[1] - 0.5 * res[2] > yx <- yr[2] + 0.5 * res[2] > > r <- newRaster(xmn=xm, xmx=xx, ymn=ym, ymx=yx, ncols=dims[1], nrows=dims[2]) > r <- setValue

[R-sig-Geo] reading *.grd files - revisited

2009-01-27 Thread Monica Pisica
Hi, I've read already the thread about how to read *.grd files in R. What i want to do is to re-save the file in a geotiff format so i can load it in ArcGIS. I tried readRGDAL and i got the following error: x <- readRGDAL("test_grid.grd") Error in .local(.Object, ): GDAL Error 4: "test_g

Re: [R-sig-Geo] merging multiscale data

2008-12-04 Thread Monica Pisica
Hi Jeremy, Thank you very much for your thoughtful reply. The data i need to merge / interpolate is "topobathy" in the sense that i am trying to come up with a solution to do a seamless DEM for topography and bathymetry in a coastal area of Northern Gulf of Mexico. For this task i will have l

[R-sig-Geo] merging multiscale data

2008-12-03 Thread Monica Pisica
Hi everybody, I am very much interested in hearing what is your experience with merging different scale data to obtain a uniform surface at the finest scale possible, and if not at the most optimum scale. Now I know that "optimum" is in the eye of the beholder, but maybe we can agree on a def

Re: [R-sig-Geo] Displaying CIR images and overlay shape files

2008-05-21 Thread Monica Pisica
21:45:27 +0200> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]> CC: r-sig-geo@stat.math.ethz.ch> Subject: Re: [R-sig-Geo] Displaying CIR images and overlay shape files> > Monica Pisica wrote:> > I suppose that somehow I can use spplot, but I have to recognize that I don't k

[R-sig-Geo] Displaying CIR images and overlay shape files

2008-05-21 Thread Monica Pisica
Hi, I am still struggling to understand how to display images in R …. I have a CIR 3 band file (geotiff format with projection imbedded) and a polygon shape file with same projection. I want to display the CIR image and overlay on top only the contour of the poly shape file, and I want to dis

Re: [R-sig-Geo] displayDataset command

2008-05-20 Thread Monica Pisica
various vector layers.> > Tom Hengl> http://spatial-analyst.net > > > > -Original Message-> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of> Roger Bivand> Sent: dinsdag 20 mei 2008 16:38> To: Monica Pisica> Cc: r-sig-geo@stat.math.ethz.ch> S

Re: [R-sig-Geo] displayDataset command

2008-05-20 Thread Monica Pisica
Subject: Re: [R-sig-Geo] displayDataset command> > On Tue, 20 May 2008, Monica Pisica wrote:> > >> > Hi,> >> > I am trying to display a geotiff using displayDataset command from rgdal > > . which i can do pretty successfully but i cannot display axes fo

[R-sig-Geo] displayDataset command

2008-05-20 Thread Monica Pisica
Hi, I am trying to display a geotiff using displayDataset command from rgdal . which i can do pretty successfully but i cannot display axes for my map i tried to use axes = TRUE and i got the error that axes is matched by multiple arguments and if i try to set xaxs and yaxs as equal