Basics: elevation data
In this example, we will look at elevation data from the worldclim 2 data, crop it for Western Europe, and then change the resolution to aggregate the data. The first step is to get the worldclim layer for elevation:
using SimpleSDMLayers
using StatsPlots
import Statistics
elevation = convert(Float32, SimpleSDMPredictor(WorldClim, Elevation))
SDM response → 1080×2160 grid with 808053 Float32-valued cells
Latitudes -90.0 ⇢ 90.0
Longitudes -180.0 ⇢ 180.0
Thanks to the integration with Plots and StatsPlots, we can very rapidly visualize these data:
heatmap(elevation, c=:cividis, frame=:box)
xaxis!("Longitude")
yaxis!("Latitude")
Let's also have a look at the density while we're at it:
density(elevation, frame=:zerolines, c=:grey, fill=(0, :grey, 0.2), leg=false)
xaxis!("Elevation")
The next step is to clip the data to the region of interest. This requires a the coordinates of the bounding box as two tuples (for longitude and latitude) – we can also make a quick heatmap to see what the region looks like:
elevation_europe = clip(elevation; left=-11.0, right=31.5, bottom=29.0, top=71.5)
SDM response → 255×255 grid with 38507 Float32-valued cells
Latitudes 29.0 ⇢ 71.5
Longitudes -11.0 ⇢ 31.5
heatmap(elevation_europe, c=:cividis, aspectratio=1, frame=:box)