Намиране на граничните точки във файл с растерни данни

Ако вашите данни са от клас RasterLayer, тогава extent ще ви даде степента на растера и xmin, 'min, ymax и xmax за достъп до различните слотове.

eg

# create a dummy raster
r1 <- raster(nrows=108, ncols=21, xmn=0, xmx=10)
r1[] <-1

extent(r1)
## class       : Extent 
## xmin        : 0 
## xmax        : 10 
## ymin        : -90 
## ymax        : 90 

Можете да получите достъп до различните слотове, като използвате

xmin(r1)
## [1] 0
xmax(r1)
##[1] 10
ymin(r1)
## [1] -90
ymax(r1)
## [1] 90

Ако вашите данни са SpatialGridDataFrame, тогава bbox ще върне ограничителната кутия

.grid <- as(r1,'SpatialGridDataFrame')

bbox(.grid)
##    min max
## s1   0  10
##  2 -90  90

Моята xy координата лежи ли в границата на растера

можете да използвате cellFromXY, за да намерите клетката id и тя ще върне NA, ако е извън текста

eg

 # some data
 .points <- rbind(c(1,1),c(-4,1))
 # the first point lies within the raster, the second not

 # cell from XY will tell you this.
 cellFromXY(r1,.points)
 ## [1] 1116   NA

РЕДАКТИРАНЕ за ggmap

ако имате карта, придобита от get_map, тя е ggmap обект и няма да работи с пакета raster без малко помощ от вас .

можете да получите ограничителната кутия като атрибут bb.

  hdf <- get_map()
  attr(hdf,'bb')
##    ll.lat    ll.lon   ur.lat    ur.lon
## 1 29.38048 -95.80204 30.14344 -94.92313

Помощна функция, която ще създаде RasterStack от обект ggmap

ggmap_rasterlayer <- function(map){
  map_bbox <- attr(map, 'bb') 
  .extent <- extent(as.numeric(map_bbox[c(2,4,1,3)]))
  my_map <- raster(.extent, nrow= nrow(map), ncol = ncol(map))
  rgb_cols <- setNames(as.data.frame(t(col2rgb(map))), c('red','green','blue'))
  red <- my_map
  values(red) <- rgb_cols[['red']]
  green <- my_map
  values(green) <- rgb_cols[['green']]
  blue <- my_map
  values(blue) <- rgb_cols[['blue']]
  stack(red,green,blue)

}

my_map <- ggmap_rasterlayer(hdf)