require(grDevices) # for coloursx <- y <- seq(-4*pi, 4*pi, len = 27)r <- sqrt(outer(x^2, y^2, "+"))image(z = z <- cos(r^2)*exp(-r/6), col = gray((0:32)/32))image(z, axes = FALSE, main = "Math can be beautiful ...", xlab = expression(cos(r^2) * e^{-r/6}))contour(z, add = TRUE, drawlabels = FALSE)
# Volcano data visualized as matrix. Need to transpose and flip# matrix horizontally.image(t(volcano)[ncol(volcano):1,])
# A prettier display of the volcanox <- 10*(1:nrow(volcano))y <- 10*(1:ncol(volcano))image(x, y, volcano, col = terrain.colors(100), axes = FALSE)contour(x, y, volcano, levels = seq(90, 200, by = 5), add = TRUE, col = "peru")axis(1, at = seq(100, 800, by = 100))axis(2, at = seq(100, 600, by = 100))box()title(main = "Maunga Whau Volcano", font.main = 4
转载于:https://www.cnblogs.com/arcserver/p/6437748.html
相关资源：GDAL进行shapefile数据栅格化