creating the convolution matrix for a region

[keewis]

Subtask of #18, towards #19: So far, we're using the cell ids as coordinates for the sparse matrix. However, this only makes sense for global coverage, as in that case the cell ids are also the cell number.

For regional matrices, we need to find a way to translate cell ids to the regional cell number, while still getting the same sparse matrix.