Edzer Pebesma

> but it doesn't seem to make sense to me why this function should use s2 as a default when s2 can't actually produce the same results. because GEOS can...

Units can still be in degrees: ```r library(sf) library(tidyverse) tibble(id = 1, lon = 0, lat = 50) %>% st_as_sf(coords = c("lon", "lat")) %>% st_set_crs(4326) %>% #st_buffer(111000, nQuadSegs = 1,...

> so i think my PR is all that can be done. I agree, but you didn't finish it.

Could you specify "would not work"? I downloaded the map you point to, which obviously is not valid on S2, but `st_make_valid()` makes it valid without issue.

It would be great if you could provide the data set that generated this, if needed off-line (email), to help us understand what causes this. The solution to this is...

Thanks! As an alternative you might try ```r p

> However, when I load sf, it says: "Linking to GEOS 3.9.1, GDAL 3.2.3, PROJ 7.2.1" This may come from a _binary_ package install, which includes these three libraries (ignoring...

The `gdalinfo` command line utility doesn't return anything useful when given this data source, so I guess it can't be ready by GDAL.

Right, I see ``` ogrinfo -so "/vsizip/vsicurl/https://prd-tnm.s3.amazonaws.com/StagedProducts/Hydrography/NHD/National/HighResolution/GDB/NHD_H_National_GDB.zip" NHDFlowline INFO: Open of `/vsizip/vsicurl/https://prd-tnm.s3.amazonaws.com/StagedProducts/Hydrography/NHD/National/HighResolution/GDB/NHD_H_National_GDB.zip' using driver `OpenFileGDB' successful. Layer name: NHDFlowline Geometry: 3D Measured Multi Line String Feature Count: 29321407 Extent: (-179.150067,...

Here is an example where `wkt_filter` does seem to work: ```r library(sf) # Linking to GEOS 3.9.1, GDAL 3.3.2, PROJ 7.2.1; sf_use_s2() is TRUE demo(nc, echo = FALSE) # Reading...