From 55a040e233a1edb0f8aeb23c84236c424444d411 Mon Sep 17 00:00:00 2001 From: Mihai Cara Date: Wed, 18 Dec 2024 18:56:38 -0500 Subject: [PATCH] Fix inequality sign and var name in outside bbox test --- gwcs/wcs.py | 18 +++++++++--------- 1 file changed, 9 insertions(+), 9 deletions(-) diff --git a/gwcs/wcs.py b/gwcs/wcs.py index e0cbf191..2953fc85 100644 --- a/gwcs/wcs.py +++ b/gwcs/wcs.py @@ -508,19 +508,19 @@ def outside_footprint(self, world_arrays): axis_range = footprint[:, idim] else: axis_range = footprint - range = [axis_range.min(), axis_range.max()] + min_ax = axis_range.min() + max_ax = axis_range.max() if (axtyp == 'SPATIAL' and str(phys).endswith((".ra", ".lon")) - and range[1] - range[0] > 180): + and max_ax - min_ax > 180): # most likely this coordinate is wrapped at 360 - d = np.mean(range) - range = [ - axis_range[axis_range < d].max(), - axis_range[axis_range > d].min() - ] - outside = (coord >= range[0]) & (coord < range[1]) + d = 0.5 * (min_ax + max_ax) + m = axis_range <= d + min_ax = axis_range[m].max() + max_ax = axis_range[~m].min() + outside = (coord > min_ax) & (coord < max_ax) else: - outside = (coord < range[0]) | (coord > range[1]) + outside = (coord < min_ax) | (coord > max_ax) if np.any(outside): if np.isscalar(coord): coord = np.nan