Support API for CHM calculation from DEM

API like:

`roi.get_chm(DTM, DEM)`

Currently, `roi.crop` only can return one layer with nodata inside

it will cause problem on the plot boundary.

The new function will first crop roi with bbox, and calculate the DSM-DTM, then crop by ROI polygon, will solve this problem