Crop to shape not working properly

[Ludovic Räss]

Trying to crop bedrock elevation to shape (glacier outline) does not work as expected in here

GlacioTools.jl/src/Alpine_glaciers.jl

Line 92 in 4650bc8

 push!(BedrElev_stack, GlacioTools.crop_padded(BedrElev, shape, padding)) 

Shape is following

But cropped output is following

Exclude the plotting distortion.

[Mauro]

Can you post a minimal example?

[Ludovic Räss]

using

mask_trim(raster, poly; pad=10) = trim(mask(raster; with=poly); pad=pad)

instead of GlacioTools.crop_padded such as

push!(IceThick_stack, mask_trim(IceThick, shape; pad=padding))

in

GlacioTools.jl/src/Alpine_glaciers.jl

Line 107 in 3f8287d

 push!(IceThick_stack, GlacioTools.crop_padded(IceThick, shape, padding)) 

seems to fix the thing now.

[Ludovic Räss]

Should now be implemented in #14 and close this upon merge.

[Mauro]

Closed with #14