library(sf)

## Linking to GEOS 3.6.2, GDAL 2.2.3, proj.4 4.9.3

dp <- getOption("default.datadir")
f <- file.path(dp, "data_staging/maps.nyc.gov/download/3dmodel/DA_WISE_Multipatch/DA_WISE_Multipatch/DA10_3D_Buildings_Multipatch.gdb")
d <- read_sf(f)
## I did some hacky to find some nearby polygons from this huge
## data set
midx <- c(4021792L, 4021783L, 4021790L, 4021791L, 4021795L, 4021797L, 
          4021788L, 4021789L, 4021796L, 4021794L, 4021782L, 4437054L, 4469157L)
idx <- match(midx, d$BIN)

library(anglr)
library(silicate)
rgl::rgl.useNULL()

## [1] TRUE

## plot as wiremesh
plot3d(SC(d[idx, ]))
rgl::rglwidget()

rgl::rgl.clear()
## plot as planar triangulate polygon layers
plot3d(TRI(d[idx, ]))
rgl::rglwidget()