Add some additional options to make_normalised_raster()

[scabecks]

• add default resampling method defaults based on comparing resolutions of both PU and input raster.
• add manual override using terra default methods
• add conditional_expression to allow user to do simple single raster algebrda using defined function, (e.g., function (r) ifel(r > 4, 1, 0) )
• add fill_na option (default = TRUE) to fill in NA values to 0 before masking to PU layer
• add default output datatype (INT1S or FLT4S) and NAflag when writing to file based on raster number type.

[sandra-neubert]

Looks good @scabecks, just some small things: I think it would be good to generalise the structure of functions and use %>% instead of |> (just because I used this in all other functions so far). I only saw this once now in line 96 in make_normalised_raster(). And when writing the roxygen skeleton I think it would also be good to have a generic way of doing it - I don't mind which way, but I saw that in make_normalised_raster() you give the inputs with [] (e.g. raster_in [SpatRaster] ), but not in the other new function you added. I really like the idea of including the conditional expression!

[scabecks]

@sandra, the roxygen formatting is just me being lazy and having chatgpt fill in the commenting - and it's never constant. Good idea though - we'll stick with the generic way that exists already.