Make use of the fact that longitudes are an xarray data array, and us…

…e 'y' dimension for rows/scanlines This makes the code resillient towards any other (awkward) data layout than the standard (where first dimension is usually the rows. Signed-off-by: Adam.Dybbroe <[email protected]>
adybbroe · Nov 17, 2023 · 05d7e37 · 05d7e37
1 parent 6a8afc0
commit 05d7e37
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/pyresample/geometry.py b/pyresample/geometry.py
@@ -1,7 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
 #
-# Copyright (C) 2010-2020 Pyresample developers
+# Copyright (C) 2010-2023 Pyresample developers
 #
 # This program is free software: you can redistribute it and/or modify it under
 # the terms of the GNU Lesser General Public License as published by the Free
@@ -693,16 +693,17 @@ def geocentric_resolution(self, ellps='WGS84', radius=None, nadir_factor=2):
         if self.ndim == 1:
             raise RuntimeError("Can't confidently determine geocentric "
                                "resolution for 1D swath.")
-        rows = self.shape[0]
+
+        rows = self.lons['y'].shape[0]
         start_row = rows // 2  # middle row
         src = CRS('+proj=latlong +datum=WGS84')
         if radius:
             dst = CRS("+proj=cart +a={} +b={}".format(radius, radius))
         else:
             dst = CRS("+proj=cart +ellps={}".format(ellps))
         # simply take the first two columns of the middle of the swath
-        lons = self.lons[start_row: start_row + 1, :2]
-        lats = self.lats[start_row: start_row + 1, :2]
+        lons = self.lons.sel(y=start_row)[:2]
+        lats = self.lats.sel(y=start_row)[:2]
         if hasattr(lons.data, 'compute'):
             # dask arrays, compute them together
             import dask.array as da