一些多边形无法填充颜色

时间:2021-01-25 20:34:41

标签: r plot geospatial polygon raster

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

我有两个问题:

  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,add=T,col=lt.green,border="grey")

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

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

0 个答案:

没有答案