contour的一个问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

现在想做一图：
图形的输出范围固定，比如说xlim=c(0, 100), ylim=c(0, 100)
现在根据这幅图上一些数据，使用contour画等值线，现在想输入等值线的范围比如横坐标[20 80],纵坐标[30 60]
但是contour里边的xlim和ylim只是设定整幅图的输出范围，结果图形只能显示[20 80]和[30 60]的范围，而不是需要的[0 100], [0 100]的范围。

简单的说我的问题是，想在图形中只局部使用contour，可以做到吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：contour Tour con Our 我的问题

先画出你要的大范围，然后用contour的add = T，把contour加上去
比如：
x <- -6:16
plot(-1:2, -1:2, type = "n")
contour(outer(x, x), method = "edge", vfont = c("sans serif", "plain"), add = T)

qoiqpwqr 发表于 2012-2-21 12:36
先画出你要的大范围，然后用contour的add = T，把contour加上去
比如：
x

俺遇到的问题更为复杂，需要画好地图。请将site.xls先改写为site.csv存在R工作目录中才能运行下边的程序。

library("maps")
library("mapdata")

graphics.off()

data <- read.csv("site.csv",header=TRUE)
x    <- data$Lat
y    <- data$Long
NGEN <- data$NGEN

loess.NGEN    <- loess(NGEN~x*y, degree=2, span=0.3, normalize=F, family="symmetric")
NGEN.marginal <- list(x=seq(min(x),max(x),0.1),y=seq(min(y),max(y),0.1))
predict.NGEN  <- predict(loess.NGEN,expand.grid(NGEN.marginal))

tiff(file="NGEN.tif", width=85, height=72, units='mm',res=1200, compression='lzw', pointsize=5)
par(mar=c(6,7,1,1))
contour(NGEN.marginal$y, NGEN.marginal$x, t(predict.NGEN), levels=seq(4, 18, by=2), xlim=c(75, 134), ylim=c(19, 53), col="grey50", lwd=1, labcex=2, xlab="", ylab="", cex.lab=2, cex.axis=2)
mtext(expression(paste("Latitude  (",degree,"N)")), side = 2, line = 4, cex = 2)
mtext(expression(paste("Longitude  (",degree,"E)")), side = 1, line = 4, cex = 2)
map("china", add=TRUE, lwd=0.5, col="grey20")
contour(NGEN.marginal$y,NGEN.marginal$x,t(predict.NGEN),levels=seq(4, 18, by=2), xlim=c(75, 134), ylim=c(19, 53), col="grey5", lwd=1, labcex=2, xlab="", ylab="", cex.lab=2, cex.axis=2, add=TRUE)
mtext(expression(paste("Latitude  (",degree,"N)")), side = 2, line = 4, cex = 2)
mtext(expression(paste("Longitude  (",degree,"E)")), side = 1, line = 4, cex = 2)
dev.off