Hi Steve,
(I'm replying to your e-mail at the mailing list to get it included
there also)
I've written some stuff on graphs with sp at the R wiki:
http://wiki.r-project.org/rwiki/doku.php?id=tips:spatial-data:spatial_data_visualization
or have a look at the gallery containing plots with the code that
created them:
http://r-spatial.sourceforge.net/gallery/
cheers!
Paul
[EMAIL PROTECTED] wrote:
Hello
I'm working with a large hydrological data set stored in a netCDF format.
The file stores x and y coordinates in the UTM projected coordinate system,
yet when I use image to graphically display the z variable, the image is
distorted in the sense that it does not plot the map in the correct spatial
organization.
I'm wondering if I need to define the projection of the netCDF file with
rgdal or proj4 routines first before I send it to the graphics device.
My code is as follows:
q1_2001 <- open.ncdf("H:\\SKF_DESKTOP FILES\\My
Documents\\EDEN\\EDEN\\Surfaces\\2000_q1.nc", readunlimi=FALSE) #opens ncdf
file for reading
wat.data2000q1 <- get.var.ncdf(q1_2001, verbose=FALSE ) # gets the real
information
# GENERAL EXAMINATION OF HEADER DATA in the wat.data file
day <- get.var.ncdf(q1_2001, "time") # length(day) 91 days in quarter
UTMx <- get.var.ncdf(q1_2001, "x") # columns (eastings) # should
return 405
UTMy <- get.var.ncdf(q1_2001, "y") # rows (northings) #
should return 287
# plot first 91 days (3 months of the year)
for(i in 1:91) {
!is.na( image(UTMx, UTMy, z = wat.data2001q1[,,i], col=brewer.pal(8,
"YlGnBu"),
axes=T, pty="s", ylab="UTM Northing", xlab="UTM Easting",
main = "First Quater 2001") )
}
As I indicated above the map is displayed on the graphics device. However
the orientation is distorted pulling the x axis to wide and the y axis too
tall. How can I set the graphics device to know the orientation and
scaling (if these are the correct terms) in order to display this map
correctly?
All insights will be greatly appreciated.
Thanks
Steve
Steve Friedman Ph. D.
Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
Office (305) 224 - 4282
[EMAIL PROTECTED]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
--
Drs. Paul Hiemstra
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone: +31302535773
Fax: +31302531145
http://intamap.geo.uu.nl/~paul
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
