pytroll / satpy

Python package for earth-observing satellite data processing
http://satpy.readthedocs.org/en/latest/
GNU General Public License v3.0
1.05k stars 288 forks source link

Highlights blacked-out in MODIS false-color image #2455

Closed yukaribbba closed 1 year ago

yukaribbba commented 1 year ago

Describe the bug I use exactly the same composite and enhancement as Polar2Grid.

Satpy image

P2G image

There's no problem with true_color

To Reproduce

files = find_files_and_readers(base_dir="/aqua_2008034_0515", reader="modis_l1b")
scn = Scene(filenames=files)
scn.load(["false_color"])
area_id = "lcc_cn_250"
projection = "+proj=lcc +lon_0=105 +lat_1=25 +lat_2=47 +datum=WGS84 +ellps=WGS84"
area_def = create_area_def(area_id=area_id, projection=projection, resolution=250)
scn = scn.resample(area_def, resampler="ewa", weight_delta_max=10, weight_distance_max=1)
scn.save_dataset("false_color", writer="geotiff", fill_value=0, filename="aqua_modis_false_color_20080203_051500_lcc_cn_250_temp.tif", compress=None)

Expected behavior A clear and concise description of what you expected to happen.

Actual results C:\Users\45107\miniconda3\lib\site-packages\satpy_config.py:129: DeprecationWarning: SelectableGroups dict interface is deprecated. Use select. for entry_point in cached_entry_points().get(name, []): [DEBUG: 2023-04-25 22:57:02 : satpy.readers.yaml_reader] Reading ('C:\Users\45107\miniconda3\Lib\site-packages\satpy\etc\readers\modis_l1b.yaml',) [DEBUG: 2023-04-25 22:57:02 : satpy.readers.yaml_reader] Reading ('C:\Users\45107\miniconda3\Lib\site-packages\satpy\etc\readers\modis_l1b.yaml',) [DEBUG: 2023-04-25 22:57:02 : satpy.readers.yaml_reader] Assigning to modis_l1b: ['C:/Users/45107/Downloads/Sat/Polar\aqua_2008034_0515\MYD02QKM.A2008034.0515.061.2018031165158.hdf', 'C:/Users/45107/Downloads/Sat/Polar\aqua_2008034_0515\MYD02HKM.A2008034.0515.061.2018031165158.hdf', 'C:/Users/45107/Downloads/Sat/Polar\aqua_2008034_0515\MYD021KM.A2008034.0515.061.2018031165158.hdf', 'C:/Users/45107/Downloads/Sat/Polar\aqua_2008034_0515\MYD03.A2008034.0515.061.2018031042907.hdf'] [DEBUG: 2023-04-25 22:57:02 : satpy.composites.config_loader] Looking for composites config file modis.yaml [DEBUG: 2023-04-25 22:57:03 : satpy.composites.config_loader] Looking for composites config file visir.yaml [DEBUG: 2023-04-25 22:57:03 : satpy.readers.hdfeos_base] Loading latitude [DEBUG: 2023-04-25 22:57:04 : satpy.readers.hdfeos_base] Interpolating from 1000 to 500 [DEBUG: 2023-04-25 22:57:04 : satpy.readers.yaml_reader] No coordinates found for DataID(name='latitude', resolution=500, modifiers=()) [DEBUG: 2023-04-25 22:57:04 : satpy.readers.hdfeos_base] Loading latitude [DEBUG: 2023-04-25 22:57:05 : satpy.readers.hdfeos_base] Interpolating from 1000 to 250 [DEBUG: 2023-04-25 22:57:05 : satpy.readers.yaml_reader] No coordinates found for DataID(name='latitude', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:05 : satpy.readers.hdfeos_base] Loading longitude [DEBUG: 2023-04-25 22:57:05 : satpy.readers.yaml_reader] No coordinates found for DataID(name='longitude', resolution=500, modifiers=()) [DEBUG: 2023-04-25 22:57:05 : satpy.readers.hdfeos_base] Loading longitude [DEBUG: 2023-04-25 22:57:05 : satpy.readers.yaml_reader] No coordinates found for DataID(name='longitude', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:05 : satpy.readers.hdfeos_base] Loading solar_zenith_angle [DEBUG: 2023-04-25 22:57:06 : satpy.readers.hdfeos_base] Interpolating from 1000 to 250 [DEBUG: 2023-04-25 22:57:06 : satpy.readers.hdfeos_base] Loading solar_azimuth_angle [DEBUG: 2023-04-25 22:57:06 : satpy.readers.hdfeos_base] Loading satellite_zenith_angle [DEBUG: 2023-04-25 22:57:06 : satpy.readers.hdfeos_base] Interpolating from 1000 to 250 [DEBUG: 2023-04-25 22:57:07 : satpy.readers.hdfeos_base] Loading satellite_azimuth_angle [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) because of incompatible areas [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of dependency's delayed generation: DataID(name='solar_zenith_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delayed prerequisite for 'DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected'))': 'DataID(name='solar_zenith_angle', resolution=250, modifiers=())' [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) because of incompatible areas [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of dependency's delayed generation: DataID(name='solar_zenith_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delayed prerequisite for 'DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected'))': 'DataID(name='solar_zenith_angle', resolution=250, modifiers=())' [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) because of incompatible areas [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of dependency's delayed generation: DataID(name='solar_zenith_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delayed prerequisite for 'DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected'))': 'DataID(name='solar_zenith_angle', resolution=250, modifiers=())' [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='false_color') because of dependency's delayed generation: DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delayed prerequisite for 'DataID(name='false_color')': 'DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected'))' [DEBUG: 2023-04-25 22:57:07 : satpy.modifiers.geometry] Applying sun zen correction [DEBUG: 2023-04-25 22:57:07 : satpy.modifiers.geometry] Apply the standard sun-zenith correction [1/cos(sunz)] [DEBUG: 2023-04-25 22:57:07 : satpy.modifiers._crefl] Loading CREFL averaged elevation information from: modifiers/tbase.hdf [DEBUG: 2023-04-25 22:57:07 : satpy.modifiers._crefl_utils] Using average elevation information provided to CREFL [DEBUG: 2023-04-25 22:57:07 : satpy.modifiers._crefl_utils] Using MODIS CREFL algorithm [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of dependency's delayed generation: DataID(name='solar_zenith_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delayed prerequisite for 'DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected'))': 'DataID(name='solar_zenith_angle', resolution=250, modifiers=())' [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) because of incompatible areas [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) because of incompatible areas [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) because of incompatible areas [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delaying generation of DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) because of dependency's delayed generation: DataID(name='solar_zenith_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Delayed prerequisite for 'DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected'))': 'DataID(name='solar_zenith_angle', resolution=250, modifiers=())' [WARNING: 2023-04-25 22:57:07 : satpy.scene] The following datasets were not created and may require resampling to be generated: DataID(name='false_color') [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Unloading dataset: DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=250, calibration=, modifiers=()) [DEBUG: 2023-04-25 22:57:07 : satpy.scene] Unloading dataset: DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=250, calibration=, modifiers=('sunz_corrected',)) C:\Users\45107\miniconda3\lib\site-packages\pyproj\crs\crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems proj = self._crs.to_proj4(version=version) [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=()) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=250, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=()) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=()) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='satellite_azimuth_angle', resolution=250, modifiers=()) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='satellite_zenith_angle', resolution=250, modifiers=()) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='solar_azimuth_angle', resolution=250, modifiers=()) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.scene] Resampling DataID(name='solar_zenith_angle', resolution=250, modifiers=()) [INFO: 2023-04-25 22:57:20 : satpy.scene] Not reducing data before resampling. [DEBUG: 2023-04-25 22:57:20 : satpy.modifiers.geometry] Applying sun zen correction [DEBUG: 2023-04-25 22:57:20 : satpy.modifiers.geometry] Apply the standard sun-zenith correction [1/cos(sunz)] [DEBUG: 2023-04-25 22:57:20 : satpy.modifiers._crefl] Loading CREFL averaged elevation information from: modifiers/tbase.hdf [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl_utils] Using average elevation information provided to CREFL [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl_utils] Using MODIS CREFL algorithm [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers.geometry] Applying sun zen correction [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers.geometry] Apply the standard sun-zenith correction [1/cos(sunz)] [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl] Loading CREFL averaged elevation information from: modifiers/tbase.hdf [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl_utils] Using average elevation information provided to CREFL [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl_utils] Using MODIS CREFL algorithm [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers.geometry] Applying sun zen correction [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers.geometry] Apply the standard sun-zenith correction [1/cos(sunz)] [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl] Loading CREFL averaged elevation information from: modifiers/tbase.hdf [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl_utils] Using average elevation information provided to CREFL [DEBUG: 2023-04-25 22:57:21 : satpy.modifiers._crefl_utils] Using MODIS CREFL algorithm [DEBUG: 2023-04-25 22:57:21 : satpy.composites] Sharpening image with high resolution green band [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=()) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=250, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=()) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=()) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='satellite_azimuth_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='satellite_zenith_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='solar_azimuth_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='solar_zenith_angle', resolution=250, modifiers=()) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='7', wavelength=WavelengthRange(min=2.105, central=2.13, max=2.155, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='2', wavelength=WavelengthRange(min=0.841, central=0.8585, max=0.876, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected',)) [DEBUG: 2023-04-25 22:57:21 : satpy.scene] Unloading dataset: DataID(name='1', wavelength=WavelengthRange(min=0.62, central=0.645, max=0.67, unit='µm'), resolution=500, calibration=, modifiers=('sunz_corrected', 'rayleigh_corrected')) [DEBUG: 2023-04-25 22:57:21 : satpy.writers] Reading ['C:\Users\45107\miniconda3\Lib\site-packages\satpy\etc\writers\geotiff.yaml'] [DEBUG: 2023-04-25 22:57:21 : rasterio.session] Could not import boto3, continuing with reduced functionality. [DEBUG: 2023-04-25 22:57:21 : rasterio.env] PROJ data files are available at built-in paths. C:\Users\45107\miniconda3\lib\site-packages\satpy_config.py:129: DeprecationWarning: SelectableGroups dict interface is deprecated. Use select. for entry_point in cached_entry_points().get(name, []): [DEBUG: 2023-04-25 22:57:22 : satpy.writers] Adding enhancement configuration from file: C:\Users\45107\miniconda3\Lib\site-packages\satpy\etc\enhancements\generic.yaml [DEBUG: 2023-04-25 22:57:22 : satpy.writers] Adding enhancement configuration from file: D:\satpy_config\enhancements\generic.yaml [DEBUG: 2023-04-25 22:57:22 : satpy.writers] Data for DataID(name='false_color', resolution=250) will be enhanced with options: [{'name': 'reflectance_range', 'method': <function stretch at 0x0000024518DE5750>, 'kwargs': {'stretch': 'crude', 'min_stretch': -1.0, 'max_stretch': 110.0}}, {'name': 'Linear interpolation', 'method': <function piecewise_linear_stretch at 0x000002451B049C60>, 'kwargs': {'xp': [0.0, 25.0, 55.0, 100.0, 255.0], 'fp': [0.0, 90.0, 140.0, 175.0, 255.0], 'reference_scale_factor': 255}}] [DEBUG: 2023-04-25 22:57:22 : trollimage.xrimage] Applying stretch crude with parameters {'min_stretch': -1.0, 'max_stretch': 110.0} [DEBUG: 2023-04-25 22:57:22 : satpy.enhancements] Applying the piecewise_linear_stretch [DEBUG: 2023-04-25 22:57:22 : rasterio.env] Entering env context: <rasterio.env.Env object at 0x0000024519487250> [DEBUG: 2023-04-25 22:57:22 : rasterio.env] Starting outermost env [DEBUG: 2023-04-25 22:57:22 : rasterio.env] No GDAL environment exists [DEBUG: 2023-04-25 22:57:22 : rasterio.env] New GDAL environment <rasterio._env.GDALEnv object at 0x000002451968F6D0> created [DEBUG: 2023-04-25 22:57:22 : rasterio._filepath] Installing FilePath filesystem handler plugin... [DEBUG: 2023-04-25 22:57:22 : rasterio._env] GDAL_DATA found in environment. [DEBUG: 2023-04-25 22:57:22 : rasterio._env] PROJ data files are available at built-in paths. [DEBUG: 2023-04-25 22:57:22 : rasterio._env] Started GDALEnv: self=<rasterio._env.GDALEnv object at 0x000002451968F6D0>. [DEBUG: 2023-04-25 22:57:22 : rasterio.env] Entered env context: <rasterio.env.Env object at 0x0000024519487250> [DEBUG: 2023-04-25 22:57:22 : rasterio._io] Path: _ParsedPath(path='C:/Users/45107/Downloads/Sat/Polar/aqua_modis_false_color_20080203_051500_lcc_cn_250_temp.tif', archive=None, scheme=None), mode: w, driver: GTiff [DEBUG: 2023-04-25 22:57:22 : rasterio._io] Skipped delete for overwrite, dataset does not exist: 'C:/Users/45107/Downloads/Sat/Polar/aqua_modis_false_color_20080203_051500_lcc_cn_250_temp.tif' [DEBUG: 2023-04-25 22:57:22 : rasterio._base] Nodata success: 1, Nodata value: 0.000000 [DEBUG: 2023-04-25 22:57:22 : rasterio._base] Nodata success: 1, Nodata value: 0.000000 [DEBUG: 2023-04-25 22:57:22 : rasterio._base] Nodata success: 1, Nodata value: 0.000000 [DEBUG: 2023-04-25 22:57:22 : rasterio.env] Exiting env context: <rasterio.env.Env object at 0x0000024519487250> [DEBUG: 2023-04-25 22:57:22 : rasterio.env] Cleared existing <rasterio._env.GDALEnv object at 0x000002451968F6D0> options [DEBUG: 2023-04-25 22:57:22 : rasterio._env] Stopped GDALEnv <rasterio._env.GDALEnv object at 0x000002451968F6D0>. [DEBUG: 2023-04-25 22:57:22 : rasterio.env] Exiting outermost env [DEBUG: 2023-04-25 22:57:22 : rasterio.env] Exited env context: <rasterio.env.Env object at 0x0000024519487250>

Environment Info:

djhoese commented 1 year ago

This is the mask_saturated handling which P2G defaults to the opposite of Satpy:

https://github.com/pytroll/satpy/blob/fb9bb8dcbb847b4e3e7431cef9c4f00659cb82a9/satpy/readers/modis_l1b.py#L36-L60

yukaribbba commented 1 year ago

@djhoese I realy shoud've read that😳. Closing now