微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

一些多边形无法填充颜色

如何解决一些多边形无法填充颜色

我在 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"

# 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")

我有两个问题:

  1. 当我尝试使用 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 这是什么意思,我应该做些不同的事情吗?

  2. 在我的真实情节中,e2 是国家森林边界,其中一些超出了情节轴的范围。当我尝试使用

    将它们添加到绘图中时
lt.green <- adjustcolor("lightgreen",alpha.f = 0.5) 
plot(e2,border="grey")

所有国家森林边界图,但其中一些没有填充颜色。无论是在控制台中打印图形还是尝试将其导出为 png 时,都是如此。

以下是实际 png 的示例。具有空的、细长的灰色周长的多边形是无法填充颜色的国家森林。任何帮助深表感谢。

States,NFs,NPs,Fires

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。