Отсутствует файл проекции в пакете Rgdal для геопространственного анализа.

У меня проблема, которую я не могу решить с помощью rgdal, из-за которой отсутствует файл проекции. Это ошибка, которую я получаю:

proj4string(d) <- CRS("+init=epsg:28992")

Error in CRS("+init=epsg:28992") : no system list, errno: 2

Это сбивает с толку, потому что, когда я загружаю пакет, он, кажется, успешно автоматически определяет файлы PROJ.4:

> library(rgdal)
rgdal: version: 0.9-2, (SVN revision 526)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.11.2, released 2015/02/10
Path to GDAL shared files:
Loaded PROJ.4 runtime: Rel. 4.9.1, 04 March 2015, [PJ_VERSION: 491]
Path to PROJ.4 shared files: (autodetected)

Но когда я это делаю, он возвращает ноль:

.Call("PROJcopyEPSG", tempfile(), PACKAGE = "rgdal")

Единственный ответ, который я нашел при поиске в Google и поиске SO, заключается в том, что epsg должен быть строчным, как в моем коде.

Я работаю на Mac с Yosemite, R 3.2.0 и работаю в RStudio (версия 0.99.441). Другие пакеты и версии из моей информации о сеансе:

rgeos_0.3-8          ggmap_2.4            ggplot2_1.0.1        rgdal_0.9-2
leaflet_0.0.15       tidyr_0.2.0          dplyr_0.4.1
UScensus2010blk_1.00 UScensus2010_0.11    foreign_0.8-63
maptools_0.8-36      sp_1.1-0

rstudio r rgdal gdal
person eamcvey    schedule 08.06.2015    source источник
comment
И CRS("+init=epsg:28992"), и .Call("PROJcopyEPSG", tempfile(), PACKAGE = "rgdal"), но у меня более старые среды выполнения GDAL и PROJ.4 в OS X (через доморощенный)   -  person hrbrmstr    schedule 08.06.2015
comment
Когда вы загружаете rgdal, он должен сообщить вам, где автоопределение находит общие файлы PROJ.4. Я думаю, что (autodetected) означает, что он на самом деле не находит ваши установленные файлы, что мне кажется, что ваша установка или разрешения могут быть неправильными. Тем не менее, я использую rgdal 0.9-3 в Linux и Windows, так что это может быть просто ошибка?   -  person Forrest R. Stevens    schedule 08.06.2015

Ответы (1)


arrow_upward
1
arrow_downward

Я быстро погуглил и наткнулся на этот пост что предполагает, что это была ошибка OS X, когда пакеты были созданы неправильно.

В моей системе Mac мне не удалось воспроизвести вашу ошибку, поэтому я установил свои версии rgdal и rgeos из http://www.kyngchaos.com/software:frameworks Может быть, вы можете переустановить свои версии этих пакетов и посмотреть, решит ли это проблему? Я бы также обновил sp одновременно.

В качестве обходного пути вы можете указать проекцию вручную. Строку proj для вашего прогноза можно получить по адресу: http://spatialreference.org/ref/epsg/amersfoort-rd-new/

Так, например, вместо:

proj4string(d) <- CRS("+init=epsg:28992")

вы можете указать:

proj4string(d) <- CRS("+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs ")

Не ответ, может быть, но стоит попробовать? Удачи.

person Phil    schedule 08.06.2015
comment
Это сработало - необходимо было обновить sp, просто переустановить было недостаточно. Спасибо! - person eamcvey; 22.06.2015