Установите правильные crs на объекте sf для построения координатных точек - PullRequest
Новая
       104

Установите правильные crs на объекте sf для построения координатных точек

0 голосов
/ 08 мая 2018

Я пытаюсь определить правильный CRS для моего sf объекта. Я хочу нанести точки на следующий слой (страна: Нидерланды): 

Simple feature collection with 380 features and 3 fields
geometry type:  MULTIPOLYGON
dimension:      XY
bbox:           xmin: 13565.4 ymin: 306846.2 xmax: 278026.1 ymax: 619232.6
epsg (SRID):    NA
proj4string:    +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

Выход:

enter image description here

Этот слой имеет правильную проекцию.

Но слой POINT не имеет правильного проекта CRS, предположительно, потому что у него нет proj4string? 

Simple feature collection with 566 features and 5 fields
geometry type:  POINT
dimension:      XY
bbox:           xmin: 3.5837 ymin: 50.86487 xmax: 7.120998 ymax: 53.44835
epsg (SRID):    NA
proj4string:    NA

enter image description here

Как установить ту же проекцию, что и на предыдущей карте, чтобы можно было наносить на нее точки с координатами?

Ответы [ 2 ]

0 голосов
/ 08 мая 2018

Кроме того, есть функция st_set_crs(), которую можно использовать в трубе. Например,

points %>% st_set_crs(st_crs(polygons))

0 голосов
/ 08 мая 2018

Вы можете использовать st_crs.

Как сказано в примерах: 

st_crs(x) <- value

В вашем случае это, вероятно, будет 

st_crs(points) <- st_crs(polygons)

Примечание: точки должны быть зарегистрированы в одной и той же системе координат, конечно. Если нет, вам нужно будет найти правильную систему координат для них и затем использовать st_transform, чтобы поместить их в ту же систему координат, что и полигоны.

...