Экстент файла формы с R для перекрытия

Question

У меня есть 2 файла формы: данные (dataframe пространственной точки) и polys (dataframe полигонов). Я хочу сделать наложение, но кажется, что это не работает.

Вот данные и данные:

> data
class       : SpatialPointsDataFrame 
features    : 12527 
extent      : 10.20075, 20.6108, 54.08669, 57.75905  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs 
variables   : 3
names       :    timestamp_pretty,     timestamp,     imo 
min values  : 01/04/2006 00:00:55, 1143849655232, 9048392 
max values  : 30/04/2006 23:59:36, 1146441576823, 9191541 
> polys
class       : SpatialPolygonsDataFrame 
features    : 436375 
extent      : 4210000, 5441000, 3395000, 4813000  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs 
variables   : 2
names       :    Id, Count 
min values  :     0,     0 
max values  : 99999,     9

чтобы управлять перекрытием, я использую

proj4string(data) <- proj4string(polys)                            # to confirm the same reference system
inside <- !is.na(over(data, as(polys, "SpatialPolygons")))    # overlapping shape file and data

а потом mean(inside) проверить среднее количество баллов в полисе.

Но ничего не происходит, значит всегда 0. Я использовал это много раз раньше, и это всегда работает, я думаю, это не работает не потому, что размер двух sph-файлов различен. Есть ли способ отредактировать это?

Спасибо!

Answer 1

По моему опыту, проблема, вероятно, заключается в том, как вы назначаете проекцию.

require(geosphere)
require(rgeos)
require(rgdal)
proj4string(data) <- spTransform(polys, "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs")   #this is the correct way to really project the spatial object on the same projection of polys
over(data, polys)   #this is just to check whether the two sp objects interact

попробуйте запустить этот код и скажите мне, если что-то случится.