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]> 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<-open.ncdf("gridded_obs.monthly.Tmin.1999.nc") > test.1.1999<-get.var.ncdf(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] > > _______________________________________________ > R-sig-Geo mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
