def __get_inside_cells(self, geometry: Geometry): """ :param geometry: Geometry to check :return: matrix with 1 for cells inside the geometry, 0 otherwise """ inside = np.zeros_like(self.gridX) for i in range(self.dimX): for j in range(self.dimY): x, y = self.get_coordinates(i, j) if geometry.is_in_geometry(x, y): inside[i][j] = 1 inside = inside + self.__get_entrance_cells(geometry) return inside