Fix geos proj parameters for insat 3d satellites

mraspaud · Nov 6, 2024 · b7f0738 · b7f0738
1 parent 2424541
commit b7f0738
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/satpy/readers/insat3d_img_l1b_h5.py b/satpy/readers/insat3d_img_l1b_h5.py
@@ -179,7 +179,10 @@ def get_area_def(self, ds_id):
         #fov = self.datatree.attrs["Field_of_View(degrees)"]
         fov = 18
         cfac = 2 ** 16 / (fov / cols)
-        lfac = 2 ** 16 / (fov / lines)
+
+        # From reverse engineering metadata from a netcdf file, we discovered
+        # the lfac is actually the same as cfac, ie dependend on cols, not lines!
+        lfac = 2 ** 16 / (fov / cols)
 
         h = self.datatree.attrs["Observed_Altitude(km)"] * 1000
         # WGS 84

diff --git a/satpy/tests/reader_tests/test_insat3d_img_l1b_h5.py b/satpy/tests/reader_tests/test_insat3d_img_l1b_h5.py
@@ -280,6 +280,8 @@ def test_filehandler_returns_area(insat_filehandler):
     area_def = fh.get_area_def(ds_id)
     _ = area_def.get_lonlats(chunks=1000)
     assert subsatellite_longitude == area_def.crs.to_cf()["longitude_of_projection_origin"]
+    np.testing.assert_allclose(area_def.area_extent, [-5618068.510660236, -5640108.009097205,
+                                                      5622075.692194229, 5644115.1906312])
 
 
 def test_filehandler_has_start_and_end_time(insat_filehandler):