@property
    def band_no(self):
        level_name_dict = self._get_level_name_dict()
        return {band_tag: level_name_dict[band_tag]['tile_layer'] for band_tag in level_name_dict}
Only pixels not masked out should be converted to features
dataset/ETHPy150Open GeoscienceAustralia/agdc/src/band_lookup.py/BandLookup.band_no--crossover with--dataset/ETHPy150Open mapbox/rasterio/tests/test_features.py/test_shapes_mask