Add reader for FIDUCEO MVIRI FCDR data

doc/source/api/satpy.readers.rst

@@ -364,6 +364,14 @@ satpy.readers.msi\_safe module
    :undoc-members:
    :show-inheritance:
 
+satpy.readers.mviri\_l1b\_fiduceo\_nc module
+--------------------------------------------
+
+.. automodule:: satpy.readers.mviri_l1b_fiduceo_nc
+   :members:
+   :undoc-members:
+   :show-inheritance:
+
 satpy.readers.netcdf\_utils module
 ----------------------------------
 

doc/source/index.rst

@@ -93,6 +93,9 @@ the base Satpy installation.
     * - MSG (Meteosat 8 to 11) L2 products in GRIB2 format
       - `seviri_l2_grib`
       - In development, CLM, OCA and FIR products supported
+    * - MFG (Meteosat 2 to 7) MVIRI data in netCDF format (FIDUCEO FCDR)
+      - `mviri_l1b_fiduceo_nc`
+      - Beta
     * - Himawari 8 and 9 AHI data in HSD format
       - `ahi_hsd`
       - Nominal

satpy/etc/readers/mviri_l1b_fiduceo_nc.yaml

@@ -0,0 +1,135 @@
+# References:
+#    - MFG User Handbook
+#    - FIDUCEO MVIRI FCDR Product User Guide
+
+reader:
+  name: mviri_l1b_fiduceo_nc
+  short_name: FIDUCEO MVIRI FCDR
+  long_name: >
+    Fundamental Climate Data Record of re-calibrated Level 1.5 Infrared, Water Vapour, and Visible radiances
+    from the Meteosat Visible Infra-Red Imager (MVIRI) instrument onboard the Meteosat First Generation satellites
+  description: >
+    Reader for FIDUCEO MVIRI FCDR data in netCDF format. For documentation see
+    http://doi.org/10.15770/EUM_SEC_CLM_0009 .
+  sensors: [mviri]
+  default_channels: [VIS, WV, IR]
+  reader: !!python/name:satpy.readers.yaml_reader.GEOFlippableFileYAMLReader
+
+file_types:
+    nc_easy:
+        file_reader: !!python/name:satpy.readers.mviri_l1b_fiduceo_nc.FiduceoMviriEasyFcdrFileHandler
+        file_patterns: [
+          'FIDUCEO_FCDR_{level}_{sensor}_{platform}-{projection_longitude:f}_{start_time:%Y%m%d%H%M}_{end_time:%Y%m%d%H%M}_EASY_{processor_version}_{format_version}.nc'
+          # Example: FIDUCEO_FCDR_L15_MVIRI_MET7-57.0_201701201000_201701201030_EASY_v2.6_fv3.1.nc
+        ]
+    nc_full:
+        file_reader: !!python/name:satpy.readers.mviri_l1b_fiduceo_nc.FiduceoMviriFullFcdrFileHandler
+        file_patterns: [
+          'FIDUCEO_FCDR_{level}_{sensor}_{platform}-{projection_longitude:f}_{start_time:%Y%m%d%H%M}_{end_time:%Y%m%d%H%M}_FULL_{processor_version}_{format_version}.nc'
+          # Example: FIDUCEO_FCDR_L15_MVIRI_MET7-57.0_201701201000_201701201030_FULL_v2.6_fv3.1.nc
+        ]
+
+datasets:
+  VIS:
+    name: VIS
+    resolution: 2250
+    wavelength: [0.5, 0.7, 0.9]
+    calibration:
+      reflectance:
+        standard_name: toa_bidirectional_reflectance
+        units: "%"
+      radiance:
+        # Confirmed by EUM: No (1/wavenumber) here. Hence no standard name.
+        units: W m-2 sr-1
+      counts:
+        standard_name: counts
+        units: count
+    file_type: [nc_easy, nc_full]
+
+  WV:
+    name: WV
+    resolution: 4500
+    wavelength: [5.7, 6.4, 7.1]
+    calibration:
+      brightness_temperature:
+        standard_name: toa_brightness_temperature
+        units: K
+      radiance:
+        standard_name: toa_outgoing_radiance_per_unit_wavenumber
+        units: mW m-2 sr-1 (cm-1)-1
+      counts:
+        standard_name: counts
+        units: count
+    file_type: [nc_easy, nc_full]
+
+  IR:
+    name: IR
+    resolution: 4500
+    wavelength: [10.5, 11.5, 12.5]
+    calibration:
+      brightness_temperature:
+        standard_name: toa_brightness_temperature
+        units: K
+      radiance:
+        standard_name: toa_outgoing_radiance_per_unit_wavenumber
+        units: mW m-2 sr-1 (cm-1)-1
+      counts:
+        standard_name: counts
+        units: count
+    file_type: [nc_easy, nc_full]
+
+  quality_pixel_bitmask:
+    name: quality_pixel_bitmask
+    resolution: 2250
+    file_type: [nc_easy, nc_full]
+
+  data_quality_bitmask:
+    name: data_quality_bitmask
+    resolution: 2250
+    file_type: [nc_easy, nc_full]
+
+  u_independent_toa_bidirectional_reflectance:
+    name: u_independent_toa_bidirectional_reflectance
+    long_name: "independent uncertainty per pixel"
+    units: "%"
+    resolution: 2250
+    file_type: [nc_easy]
+
+  u_structured_toa_bidirectional_reflectance:
+    name: u_structured_toa_bidirectional_reflectance
+    long_name: "structured uncertainty per pixel"
+    units: "%"
+    resolution: 2250
+    file_type: [nc_easy]
+
+  solar_zenith_angle:
+    name: solar_zenith_angle
+    standard_name: solar_zenith_angle
+    long_name: "Solar zenith angle"
+    units: degree
+    resolution: [2250, 4500]
+    file_type: [nc_easy, nc_full]
+
+  solar_azimuth_angle:
+    name: solar_azimuth_angle
+    standard_name: solar_azimuth_angle
+    long_name: "Solar azimuth angle"
+    units: degree
+    resolution: [2250, 4500]
+    file_type: [nc_easy, nc_full]
+
+  satellite_zenith_angle:
+    name: satellite_zenith_angle
+    standard_name: sensor_zenith_angle
+    long_name: "Satellite zenith angle"
+    units: degree
+    resolution: [2250, 4500]
+    file_type: [nc_easy, nc_full]
+
+  satellite_azimuth_angle:
+    name: satellite_azimuth_angle
+    standard_name: sensor_azimuth_angle
+    long_name: "Satellite azimuth angle"
+    units: degree
+    resolution: [2250, 4500]
+    file_type: [nc_easy, nc_full]

satpy/readers/_geos_area.py

@@ -145,3 +145,24 @@ def get_area_definition(pdict, a_ext):
         a_ext)
 
     return a_def
+
+
+def sampling_to_lfac_cfac(sampling):
+    """Convert angular sampling to line/column scaling factor (aka LFAC/CFAC).
+
+    Reference: `MSG Ground Segment LRIT HRIT Mission Specific Implementation`_,
+    Appendix E.2.
+
+
+    .. _MSG Ground Segment LRIT HRIT Mission Specific Implementation: http://www.eumetsat.int/\
+website/wcm/idc/idcplg?IdcService=GET_FILE&dDocName=PDF_TEN_05057_SPE_MSG_LRIT_HRI\
+&RevisionSelectionMethod=LatestReleased&Rendition=Web
+
+    Args:
+        sampling: float
+            Angular sampling (rad)
+
+    Returns:
+        Line/column scaling factor (deg-1)
+    """
+    return 2.0 ** 16 / np.rad2deg(sampling)