make test_from_shape_zurich_pass pass

Author: emanuel-schmid

climada/entity/exposures/litpop/litpop.py

@@ -1285,8 +1285,8 @@ def reproject_input_data(
         ...     'height': 1939,
         ...     'count': 1,
         ...     'crs': CRS.from_epsg(4326),
-        ...     'transform': Affine(0.00833333333333333, 0.0, -18.175000000000068,
-        ...                         0.0, -0.00833333333333333, 43.79999999999993),
+        ...     'transform': Affine(30/3600, 0.0, -18.175,
+        ...                         0.0, -30/3600, 43.8),
         ... }
 
         The meta data with the reference grid used to define the global destination

climada/test/test_litpop_integr.py

@@ -169,24 +169,30 @@ def test_from_shape_zurich_pass(self):
         )
         self.assertAlmostEqual(ent.value.sum(), 1000.0)
         self.assertEqual(ent.value.min(), 0.0)
+        self.assertAlmostEqual(ent.value.max(), 5.058, places=4)
+        # Note: up to climada 6.0 this value used to be 5.05792616746308, then 5.058035160766561.
+        # The discrepancy assumedly caused by rounding differences of the slightly changed
+        # calculation procedure.
         self.assertEqual(ent.region_id.min(), 756)
         self.assertEqual(ent.region_id.max(), 756)
-        self.assertAlmostEqual(ent.latitude.min(), 47.20416666666661)
+        self.assertAlmostEqual(ent.latitude.min(), 47.2 + 15 / 3600)
         # index and coord. of largest value:
         self.assertEqual(
-            ent.gdf.loc[ent.gdf["value"] == ent.gdf["value"].max()].index[0], 482
+            ent.gdf.loc[ent.gdf["value"] == ent.gdf["value"].max()].index[0], 434
         )
+        # Note: up to climada 6.0 this index used to be 482, because of its then irregular order.
+        # Specifically: the dataframe started at 36 and then skipped 71, 107, 143, etc.
         self.assertAlmostEqual(
             ent.gdf.loc[ent.gdf["value"] == ent.gdf["value"].max()].geometry.y.values[
                 0
             ],
-            47.34583333333325,
+            5681.5 * 30 / 3600,
         )
         self.assertAlmostEqual(
             ent.gdf.loc[ent.gdf["value"] == ent.gdf["value"].max()].geometry.x.values[
                 0
             ],
-            8.529166666666658,
+            1023.5 * 30 / 3600,
         )
 
     def test_from_shape_and_countries_zurich_pass(self):