r [sf] Как включить поддержку GEOS?

Question

Я пытаюсь преобразовать CRS полигонов, пересекающих линию даты, используя пакет sf, следующий примеру здесь , но получаю предупреждение

Предупреждающее сообщение:

In CPL_wrap_dateline(x, options, quiet) :
  GDAL Error 6: GEOS support not enabled.

, что приводит к сбою последней части кода:

library(sf)
#Linking to GEOS 3.6.1, GDAL 2.1.3, proj.4 4.9.3

point1 <- st_point(c(-179.5, -50))
point2 <- st_point(c(0, -30))
points <- st_sfc(point1, point2, crs = 4236)

points_trans <- st_transform(points, 32621)
circles <- st_buffer(points_trans, 100000)
circles_latlon <- st_transform(circles, 4326)
circles_latlon2 <- st_wrap_dateline(circles_latlon, options = c("WRAPDATELINE=YES"))

Я посмотрел в Интернете, как включить поддержку GEOS (или, возможно, просто обновить GDAL), но не могу найти простое объяснение. Есть предложения?

Я использую MACOS 10.13.6; R 3.5.0 (проверено также с 3.5.1, но та же проблема).

Answer 1

Я наконец решил это с помощью

1. загрузка и установка обновленных библиотек GDAL и GEOS из здесь
2. загрузка исходного файла sf из cran
3. (пере) установка sf с использованием R CMD INSTALL sf_0.6-3.tar.gz --configure-args='--with-gdal-config=/Library/Frameworks/GDAL.framework/Versions/2.2/unix/bin/gdal-config -with-geos-config=/Library/Frameworks/GEOS.framework/Versions/3B/unix/bin/geos-config'. Обратите внимание, что мне нужно было (повторно) добавить расширение .gz, так как Apple его удобно скрывает.