From 36af62233c9e03e9a064d2f0f68e957f1a979de9 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Fri, 1 May 2020 21:14:04 +0200 Subject: check_gtiff_grid.py: fix extent comparison when geotransform has gt[5] > 0 --- grid_tools/check_gtiff_grid.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'grid_tools') diff --git a/grid_tools/check_gtiff_grid.py b/grid_tools/check_gtiff_grid.py index 3983d02..55b4265 100755 --- a/grid_tools/check_gtiff_grid.py +++ b/grid_tools/check_gtiff_grid.py @@ -522,10 +522,10 @@ class GlobalInfo(object): def get_extent(ds): gt = ds.GetGeoTransform() xmin = gt[0] + gt[1] / 2 - ymax = gt[3] + gt[5] / 2 + ytop = gt[3] + gt[5] / 2 xmax = gt[0] + gt[1] * ds.RasterXSize - gt[1] / 2 - ymin = gt[3] + gt[5] * ds.RasterYSize - gt[5] / 2 - return xmin, ymin, xmax, ymax + ybottom = gt[3] + gt[5] * ds.RasterYSize - gt[5] / 2 + return xmin, min(ytop, ybottom), xmax, max(ytop, ybottom) def validate_ifd(global_info, ds, is_first_subds, first_subds): -- cgit v1.2.3