提示:本站收集StackOverFlow近2千万问答，支持中英文搜索，鼠标放在语句上弹窗显示对应的参考中文或英文， 本站还提供 中文繁体 英文版本 中英对照 版本，有任何建议请联系yoyou2525@163.com。
我有几个空间多边形数据帧，我要依次添加到R 4.0.2中的绘图中。 多边形代表美国各州，国家森林，国家公园和野火边界。 虽然我显然不能提供一个工作示例的所有数据，但是以下代码是我尝试制作的绘图的简化版本：
# equal area projection
proj_coords <- "+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0"
# load libraries
library(rgdal)
library(rgeos)
# Make a bbox of the total extent
e <- as(raster::extent(-1400000,-850000,2300000,2750000), "SpatialPolygons")
proj4string(e) <- proj_coords
# Make a smaller box to add to the plot for Stack overflow example
e2 <- as(raster::extent(-1300000,-900000,2400000,2650000), "SpatialPolygons")
proj4string(e2) <- proj_coords
# Plot
lt.green <- adjustcolor("lightgreen", alpha.f = 0.5)
plot(e,axes=T,xlim=c(-1150000,-1100000),ylim=c(2350000,2700000),lwd=2,border="grey",xlab="m",ylab="m")
plot(e2,add=T,col=lt.green,border="grey")
我有两个问题：
当我尝试使用proj4string(e) <- proj_coords分配CRS时，出现警告Warning message: In showSRID(uprojargs, format = "PROJ", multiline = "NO") : Discarded datum North_American_Datum_1983 in CRS definition, but +towgs84= values preserved这是什么意思，我应该做些不同的事情吗？
在我的真实情节中， e2是国家森林边界，其中一些边界超出了情节轴的范围。 当我尝试将它们添加到情节中时
lt.green <- adjustcolor("lightgreen", alpha.f = 0.5)
plot(e2,add=T,col=lt.green,border="grey")
所有国家森林边界图，但其中一些没有填充颜色。 当我在控制台中打印图形以及尝试将其导出为png时都是如此。
以下是实际png的示例。 具有空的，薄的灰色周长的多边形是无法填充颜色的国家森林。 任何帮助深表感谢。
声明：本站的技术帖子网页，遵循CC BY-SA 4.0协议，如果您需要转载，请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.