Skip to content

Validation

This chapter compares our output with that of gdaldem.

We test with the same elevation model as that in the Usage page, and also include an upside down version of it, to ensure that the algorithms are robust to different axes conventions.

julia
@info Geomorphometry.cellsize(dem)  # default negative y spacing
@info Geomorphometry.cellsize(ndem)  # reversed y spacing
┌ Error: Unknown CRS type, please report this issue for the given crs/file
└ @ GeoArrays ~/.julia/packages/GeoArrays/X63Mm/src/geointerface.jl:11
[ Info: (5.0, -5.0)
┌ Error: Unknown CRS type, please report this issue for the given crs/file
└ @ GeoArrays ~/.julia/packages/GeoArrays/X63Mm/src/geointerface.jl:11
[ Info: (5.0, 5.0)

Hillshade

julia
heatmap(hillshade(dem))

Multihillshade

julia
heatmap(multihillshade(dem))

Slope

Horn

julia
heatmap(slope(dem, method=Horn()))

ZevenbergenThorne

julia
heatmap(slope(dem, method=ZevenbergenThorne()))

Aspect

Horn

julia
heatmap(aspect(dem, method=Horn()))

ZevenbergenThorne

julia
heatmap(aspect(dem, method=ZevenbergenThorne()))