Project

General

Profile

« Previous | Next » 

Revision 96629c37

Added by Jim Regetz about 13 years ago

  • ID 96629c3791b63f3b8c5d4b40e11f947be2f2da02

added histogram plot of ASTER-SRTM differences

View differences:

terrain/dem/dem-assessment.R
260 260
    xlab="SRTM elevation (m)", ylab="ASTER elevation (m)", cex=0.5)
261 261
abline(median(delta.vals), 1, col="red", cex=0.5)
262 262
abline(0, 1, col="blue", lty=2, cex=0.5)
263
# add inset histogram of differences
264
opar <- par(fig=c(0.55, 0.95, 0.1, 0.6), new=TRUE)
265
h <- hist(delta.vals[abs(delta.vals)<60], breaks=48, xlab=NA, main=NULL,
266
    col=grey(0.8), border=grey(0.3), yaxt="n", ylab=NA, cex.axis=0.5,
267
    cex.lab=0.5, tcl=-0.25, mgp=c(3,0,0))
268
text(10, 0.4*max(h$counts), labels=paste("Entire range:\n(",
269
    min(delta.vals), ", ", max(delta.vals), ")", sep=""), cex=0.6,
270
    adj=c(0,0))
271
mtext("SRTM-ASTER (m)", side=1, cex=0.5, line=0.6)
263 272
dev.off()
264 273

  
265 274

  

Also available in: Unified diff