calculate the domain cells for a given band

[MAfarrag]

Currently the method assumes that the dataset has only one band, but if the dataset is a multi-band the method will miscalculate or give an error.

add a band parameter to the method and use the read_array method to retrieve the array values for each band.

[pyramids] pyramids/dataset.py (Lines 1054-1067)

def count_domain_cells(self):
        """Count cells inside the domain

        Returns
        -------
        int:
            Number of cells
        """
        # count cells inside the domain
        arr = self.raster.ReadAsArray()
        domain_count = np.size(arr[:, :]) - np.count_nonzero(
            (arr[np.isclose(arr, self.no_data_value[0], rtol=0.001)])
        )
        return domain_count

Open in IDE · Open on GitHub

^{Created from JetBrains using CodeStream}