Great that works.
Thanks,
Sam
Robert J. Hijmans wrote:
Hi Sam,
This should work
m <- t(test.1.1999)
m <- m[nrow(m):1,]
r <- raster(m)
plot(r)
Robert
On Fri, Sep 4, 2009 at 11:25 AM, Sam Veloz <[email protected]
<mailto:[email protected]>> wrote:
I am working with some climate data in Netcdf format that are
available here:
http://hydro.engr.scu.edu/files/gridded_obs/monthly/ncfiles/
I can successfully read in the files and plot individual bands:
library(ncdf)
test.nc
<http://test.nc><-open.ncdf("gridded_obs.monthly.Tmin.1999.nc
<http://gridded_obs.monthly.Tmin.1999.nc>")
test.1.1999<-get.var.ncdf(test.nc
<http://test.nc>,varid="Tmin",start=c(1,1,1),count=c(462,222,1))
library(raster)
image(test.1.1999) #looks good
test.raster<-raster(test.1.1999)
image(test.raster) #now North America is on its side, east coast
of NA is now on the south side of the image
Is it possible to convert this file to a raster file without
having the axes switched? Any guesses what I am doing wrong?
I had a similar problem using gdal_translate where the output was
a mirror image of what it should have been.
Thanks,
Sam
--
****************************************************
Sam Veloz
Postdoctoral Researcher
Department of Environmental Science and Policy
University of California, Davis
[email protected] <mailto:[email protected]>
_______________________________________________
R-sig-Geo mailing list
[email protected] <mailto:[email protected]>
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
****************************************************
Sam Veloz
Postdoctoral Researcher
Department of Environmental Science and Policy
University of California, Davis
[email protected]
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
