Finding indices and number pixels for rectangle shapes

This module will find the indices of rectangle or square shape and count the
  number of pixels in that shape.

and you can assign a scalar
```python
mesh[slc1, slc2] = (j + 1)
```