R: "индекс вне границ" при обрезке сетки - PullRequest
Новая
       63

R: "индекс вне границ" при обрезке сетки

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

Я собираюсь сделать пространственную интерполяцию с помощью метода кригинга на моих данных.Для продолжения мне нужна модель grid на карте Франции.

Для этого я следую главе 4 этого курса: https://www.datacamp.com/courses/spatial-statistics-in-r

К сожалению, когда я хочу обрезать сетку в области исследования, у меня появляется эта ошибка: "Error in x@coords[i, , drop = FALSE] : subscript out of bounds".

Даже если я сделаю сетку большего или меньшего размера, ошибка все равно появится.

Вот мой код R: 

#Read the France map (spatial polygon dataframe)
communes <- readOGR(dsn=paste("donnees/GEOFLA_2-2_COMMUNE_SHP_LAMB93_FXX_2016-06-28/GEOFLA/1_DONNEES_LIVRAISON_2016-06-00236/GEOFLA_2-2_SHP_LAMB93_FR-ED161/COMMUNE", sep = ""), layer="COMMUNE")

# Find the corners of the boundary
bbox(communes)
        min     max
x   99217.1 1242417
y 6049646.3 7110480

# Define a 10km square grid over the polygon extent. The first parameter is the bottom left corner.
grid <- GridTopology(c(99216,6049645), c(10000, 10000), c(120, 120))

# Create points with the same coordinate system as the boundary
gridpoints <- SpatialPoints(grid, proj4string = CRS(projection(communes)))

# Crop out the points outside the boundary
cropped_gridpoints <- crop(gridpoints, communes)

Error in x@coords[i, , drop = FALSE] : subscript out of bounds

Споследняя строка, сетка должна соответствовать краю карты, но появляется ошибка.Вот ссылка на dl шейп-файл: https://www.dropbox.com/sh/g2gm30zxpdgr0ok/AABrGMDrA3KX8LSWLzBmFDb7a?dl=0

...