From 7c5af871a5eae24d741d1253288392bfbe70fe4a Mon Sep 17 00:00:00 2001 From: "Morten W. Hansen" Date: Sun, 9 Jun 2024 11:49:41 +0000 Subject: [PATCH] Use lat/lon grids for initialization --- nansat/mappers/mapper_netcdf_cf.py | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/nansat/mappers/mapper_netcdf_cf.py b/nansat/mappers/mapper_netcdf_cf.py index 2551ab7b..637fb938 100644 --- a/nansat/mappers/mapper_netcdf_cf.py +++ b/nansat/mappers/mapper_netcdf_cf.py @@ -497,6 +497,7 @@ def _create_empty_from_subdatasets(self, gdal_dataset, metadata): no_projection = False break if no_projection: + #raise NansatMissingProjectionError sub_datasets = gdal_dataset.GetSubDatasets() filenames = [f[0] for f in sub_datasets] for _, filename in enumerate(filenames): @@ -504,6 +505,10 @@ def _create_empty_from_subdatasets(self, gdal_dataset, metadata): xDatasetSource = filename if 'latitude' in filename: yDatasetSource = filename + if not "xDatasetSource" in locals(): + raise NansatMissingProjectionError + if not "yDatasetSource" in locals(): + raise NansatMissingProjectionError lon_dataset = gdal.Open(xDatasetSource) lon = lon_dataset.GetRasterBand(1).ReadAsArray() lat_dataset = gdal.Open(yDatasetSource)