code: def cutout_shape(self, shape_obj):
          
  
          view, mask = self.get_shape_view(shape_obj)
  
          # cutout our enclosing (possibly shortened) bbox
          data = self._slice(view)
  
          # mask non-containing members
          mdata = np.ma.array(data, mask=np.logical_not(mask))
          return mdata
docs: Cut out and return a portion of the data corresponding to `shape_obj`.
          A masked numpy array is returned, where the pixels not enclosed in
          the shape are masked out.