我刚刚从以下位置下载了空间信息国家湿地清查 https://www.fws.gov/wetlands/data/Mapper.html(西北地区)。我正在尝试操纵 R 中的几何图形(st_transform、st_cast、st_crop 等),但我从未见过如此复杂的几何图形。
library(tidyverse)
library(sf)
wetlands <- st_read('./CA_geodatabase_wetlands.gdb', layer = 'CA_Wetlands')
wetlands[1,] %>% st_geometry %>% st_as_text
[1]“多曲面(曲线多边形(复合曲线)(线串(-2193424 1485786、-2193424 1485786、-2193424 1485787、-2193424 1485787、-2193426 1485793、-219342) 5 1485795, -2193419 1485806, -2193418 1485806, -2193418 1485807, -2193405 1485820, -2193404 1485820, -2193404 1485821, -2193393 1485828, -2193393 1485828, -2193391 1485829, -2193388 1485831, -2193368 14 85843)、环形线 (-2193368 1485843、-2193367 1485843、-2193366 1485843)、线线 (-2193366 1485843, -2193363 1485843, -2193359 1485845, -2193359 1485845, -2193359 1485846, -2193345 1485848, -2193344 1485848, -2193341 14 85849)、圆形线 (-2193341 1485849、-2193340 1485850、-2193339 1485850)、线线 (-2193339 1485850, -2193295 1485854, -2193294 1485855, -2193282 1485857, -2193275 1485860, -2193275 1485860, -2193273 1485861, -2193266 14 85866、-2193245 1485884、-2193239 1485889、-2193237 1485892、-2193237 1485892、-2193234 1485896、 -2193234 1485896、-2193221 1485910)、圆线(-2193221 1485910、-2193220 1485910、-2193220 1485910)、....."
由于长度的原因,我剪掉了整个打印输出,但你明白了:我有一系列嵌套的几何图形,包括 MULTISURFACE、CURVEPOLYGON、COMPOUNDCURVE 和 CIRCULARSTRING,它们都不符合标准sf
功能。我什至无法绘制它来看看它是什么样子。如何将其简化为更可行的几何形状(即多多边形或多线)?