feat: Enable to read for the Q4 coverage and the IQTI files for the fci l1c data

satpy/readers/fci_l1c_nc.py

@@ -216,18 +216,27 @@
         logger.debug("Start: {}".format(self.start_time))
         logger.debug("End: {}".format(self.end_time))
 
+        if self.filename_info["coverage"] == "Q4":
+            # change number of chunk so that padding gets activated correctly on missing chunks
+            self.filename_info["count_in_repeat_cycle"] += 28
+
+        if self.filename_info["facility_or_tool"] == "IQTI":
+            self.is_iqt = True
+        else :
+            self.is_iqt = False
+
         self._cache = {}
 
     @property
     def rc_period_min(self):
-        """Get nominal repeat cycle duration.
-
-        As RSS is not yet implemented an error will be raised if RSS are tried to be read.
-        """
-        if self.filename_info["coverage"] not in ["FD","AF"]:
-            raise NotImplementedError(f"coverage for {self.filename_info['coverage']} not supported by this reader")
+        """Get nominal repeat cycle duration."""
+        if "Q4" in self.filename_info["coverage"]:
             return 2.5
-        return 10
+        elif self.filename_info["coverage"] in ["FD","AF"]:
+            return 10
+        else:
+            raise NotImplementedError(f"coverage for {self.filename_info['coverage']}"
+                                      " not supported by this reader")
 
     @property
     def nominal_start_time(self):
@@ -378,11 +387,12 @@
             self["attr/platform"], self["attr/platform"])
 
         # remove unpacking parameters for calibrated data
-        if key["calibration"] in ["brightness_temperature", "reflectance"]:
+        if key["calibration"] in ["brightness_temperature", "reflectance","radiance"]:
             res.attrs.pop("add_offset")
             res.attrs.pop("warm_add_offset")
             res.attrs.pop("scale_factor")
             res.attrs.pop("warm_scale_factor")
+            res.attrs.pop("valid_range")
 
         # remove attributes from original file which don't apply anymore
         res.attrs.pop("long_name")
@@ -400,9 +410,21 @@
     @cached_property
     def orbital_param(self):
         """Compute the orbital parameters for the current segment."""
-        actual_subsat_lon = float(np.nanmean(self._get_aux_data_lut_vector("subsatellite_longitude")))
-        actual_subsat_lat = float(np.nanmean(self._get_aux_data_lut_vector("subsatellite_latitude")))
-        actual_sat_alt = float(np.nanmean(self._get_aux_data_lut_vector("platform_altitude")))
+        if self.is_iqt:
+            actual_subsat_lon = float(self.get_and_cache_npxr("data/mtg_geos_projection/attr/"
+                                                              "longitude_of_projection_origin"))
+            #actual_subsat_lon = -3.4
+            actual_subsat_lat = 0.0
+            actual_sat_alt = float(self.get_and_cache_npxr("data/mtg_geos_projection/attr/perspective_point_height"))
+            #actual_sat_alt = 0.0
+            logger.info("IQT data the following parameters are hardcoded "
+                        f"satellite_actual_longitude = {actual_subsat_lon} ,"
+                        f" satellite_actual_latitude = {actual_subsat_lat} ,"
+                        f"satellite_sat_alt = {actual_sat_alt}")
+        else:
+             actual_subsat_lon = float(np.nanmean(self._get_aux_data_lut_vector("subsatellite_longitude")))
+             actual_subsat_lat = float(np.nanmean(self._get_aux_data_lut_vector("subsatellite_latitude")))
+             actual_sat_alt = float(np.nanmean(self._get_aux_data_lut_vector("platform_altitude")))
         # The "try" is a temporary part of the code as long as the AF data are not modified
         try :
             nominal_and_proj_subsat_lon = float(
@@ -497,16 +519,6 @@
         for coord in "xy":
             coord_radian = self.get_and_cache_npxr(measured + "/{:s}".format(coord))
 
-            # TODO remove this check when old versions of IDPF test data (<v4) are deprecated.
-            if coord == "x" and coord_radian.attrs["scale_factor"] > 0:
-                coord_radian.attrs["scale_factor"] *= -1
-
-            # TODO remove this check when old versions of IDPF test data (<v5) are deprecated.
-            if type(coord_radian.attrs["scale_factor"]) is np.float32:
-                coord_radian.attrs["scale_factor"] = coord_radian.attrs["scale_factor"].astype("float64")
-            if type(coord_radian.attrs["add_offset"]) is np.float32:
-                coord_radian.attrs["add_offset"] = coord_radian.attrs["add_offset"].astype("float64")
-
             coord_radian_num = coord_radian[:] * coord_radian.attrs["scale_factor"] + coord_radian.attrs["add_offset"]
 
             # FCI defines pixels by centroids (see PUG), while pyresample
@@ -687,15 +699,8 @@
                 "cannot produce reflectance for {:s}.".format(measured))
             return radiance * np.float32(np.nan)
 
-        sun_earth_distance = np.mean(
-            self.get_and_cache_npxr("state/celestial/earth_sun_distance")) / 149597870.7  # [AU]
-
-        # TODO remove this check when old versions of IDPF test data (<v5) are deprecated.
-        if sun_earth_distance < 0.9 or sun_earth_distance > 1.1:
-            logger.info("The variable state/celestial/earth_sun_distance contains unexpected values"
-                        "(mean value is {} AU). Defaulting to 1 AU for reflectance calculation."
-                        "".format(sun_earth_distance))
-            sun_earth_distance = 1
+        sun_earth_distance = np.nanmean(
+            self._get_aux_data_lut_vector("earth_sun_distance")) / 149597870.7  # [AU]
 
         res = 100 * radiance * np.float32(np.pi) * np.float32(sun_earth_distance) ** np.float32(2) / cesi
         return res