##all layers are resampled to 30 meters with 3 x 3 neighborhood recalculations of pixels (ask Åshild about this)

###quick overview of neighborhood layers and what they are

## biomass
biomass<-raster(paste(raster.90m.dir, "biomass30_3_na.tif", sep = "/"))
proj4string(biomass) <- proj.rein  ## check projections

## cost distance to nearest bird colony
bird_cost_distance<-raster(paste(raster.90m.dir, "cost_air_3_na.tif", sep = "/"))
proj4string(bird_cost_distance) <- proj.rein  ## check projections

## DEM
DEM<-raster(paste(raster.90m.dir, "DEM30_3_na.tif", sep = "/"))
proj4string(DEM) <- proj.rein  ## check projections

## NDVI
NDVI<-raster(paste(raster.90m.dir, "NDVI_30_3_na.tif", sep = "/"))
proj4string(NDVI) <- proj.rein  ## check projections

## slope in degrees
slope<-raster(paste(raster.90m.dir, "slope_deg30_3_na.tif", sep = "/"))
proj4string(slope) <- proj.rein  ## check projections

## terrain ruggedness
ruggedness<-raster(paste(raster.90m.dir, "VRM30_3_na.tif", sep = "/"))
proj4string(ruggedness) <- proj.rein  ## check projections

## barren ground cover
barren<-raster(paste(raster.90m.dir, "VEG30_barren_na_3.tif", sep = "/"))
proj4string(barren) <- proj.rein  ## check projections

## ground cover "heath"
heath<-raster(paste(raster.90m.dir, "VEG30_heath_na_3.tif", sep = "/"))
proj4string(heath) <- proj.rein  ## check projections

## ground cover "moss"
moss<-raster(paste(raster.90m.dir, "VEG30_moss_na_3.tif", sep = "/"))
proj4string(moss) <- proj.rein  ## check projections