Получить CRS из prj4string - PullRequest
Новая
       28

Получить CRS из prj4string

0 голосов
/ 10 мая 2019

У меня есть файл формы, и я хочу объединить его с другим фреймом данных. Но для этого мне нужно знать Crs. Но у файла формы нет crs, только это описание: 

 Simple feature collection with 96 features and 9 fields
    geometry type:  POLYGON
    dimension:      XY
    bbox:           xmin: 313434.8 ymin: 7343789 xmax: 360663.2 ymax: 7416202
    epsg (SRID):    NA
    proj4string:    +proj=utm +zone=23 +south +ellps=aust_SA +units=m +no_defs

Я пытался найти crs из строки proj4, но не нашел. Кто-нибудь знает, как я могу получить это?

1 Ответ

0 голосов
/ 10 мая 2019

Обычно, если вы устанавливаете пакет rgdal, вы можете найти имя системы координат с помощью proj4string в предоставленных вами метаданных шейп-файла: 

library(rgdal)
showEPSG('+proj=utm +zone=23 +south +ellps=aust_SA +units=m +no_defs')

Но в этом случае ничего полезного не возвращается: "OGRERR_UNSUPPORTED_SRS" Из строки proj4 у вас есть UTM-проекция в южной зоне 23. Поиск в строке proj4 привел меня сюда: https://wiki.osgeo.org/wiki/Brazilian_Coordinate_Reference_Systems

Кроме того, здесь приведена дополнительная справка, объясняющая использование систем координат в R: https://www.nceas.ucsb.edu/~frazier/RSpatialGuides/OverviewCoordinateReferenceSystems.pdf

...