[GRASS-user] Clipping rasters

[Hanlie Pretorius]

Hi,

I thought I created a clipped DEM by rasterising an irregular vector
boundary and then using a mask to create the clipped DEM from a
rectangular DEM.

For example:
-----
# convert vector border to raster
v.to.rast input="border" output="border" use="cat" type="area" layer=1
value=1 rows=4096                                |

# set the mask to the raster border
r.mask border

# create the clipped DEM
r.mapcalc
>dem_clipped=dem_full

# remove the mask
r.mask -r
-----

But it turns out I was wrong, the 'clipped' raster just contains NULL
values outside of the border raster. What did I do wrong? Is there a
better way to clip a raster to a vector boundary?

Thanks
Hanlie