一些多边形无法填充颜色

编程问答 2022-05-28

问题描述

我在 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 (将#修改为@）

geospatialplotpolygonrrraster