Message-ID: <OFCD47AB0C.C6DFC7AC-ON85257865.0066E67C-85257865.006751D9@nps.gov>
Date: 2011-04-01T18:48:31Z
From: Steve_Friedman at nps.gov
Subject: filled contour plot with contour lines
I'm stumped, can anyone find my error in this sequence.
for(j in 1:(varsize[4]-1))
temp <- get.var.ncdf(nc=input,
varid="p_foraging",c(1,1,j),c(varsize[1],varsize[2],1))
filled.contour(x, y, temp, color = terrain.colors,
plot.title = title(main = paste("Everglades Wood Stork
Foraging Potential \nYear", (2000+j)),
xlab = "UTM East", ylab = "UTM North") ,
plot.axes = { contour(temp, add=T)
axis(1, seq(450000 , 580000, by = 10000))
axis(2, seq(2800000,4000000, by = 10000)) },
key.title = title(main="Probability") ,
key.axes = axis(4, seq(0 , 1 , by = 0.1))
The routine will work in a modified form without adding the coordinates
(the axis lines) but when I include these the routine produces various
errors, such as "dimension mismatch", or "unexpected end encountered."
I tried to follow the example on filled.contour help page.
Thanks in advance
Steve
Steve Friedman Ph. D.
Ecologist / Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
Steve_Friedman at nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147