From 4d1664711419d2fb38c4b932b88588e18514f3ae Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Tue, 3 May 2022 12:45:43 -0500
Subject: [PATCH] snapshot...
---
modules/util/util.py | 24 ++++++++----------------
1 file changed, 8 insertions(+), 16 deletions(-)
diff --git a/modules/util/util.py b/modules/util/util.py
index 6136268b..53c8c55d 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -1067,7 +1067,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
dim_0_name = 'x'
dim_1_name = 'y'
time_dim_name = 'time'
- geo_coords = 'time y x'
+ geo_coords = 'longitude latitude'
dim_1_len, dim_0_len = lons.shape
@@ -1092,7 +1092,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_pred_ds = rootgrp.createVariable('icing_prediction_level_'+flt_level_ranges_str[flvl], 'i2', var_dim_list)
icing_pred_ds.setncattr('coordinates', geo_coords)
- #icing_pred_ds.setncattr('grid_mapping', 'Projection')
+ icing_pred_ds.setncattr('grid_mapping', 'Projection')
icing_pred_ds.setncattr('missing', -1)
if has_time:
preds = preds.reshape((1, dim_1_len, dim_0_len))
@@ -1104,7 +1104,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_prob_ds = rootgrp.createVariable('icing_probability_level_'+flt_level_ranges_str[flvl], 'f4', var_dim_list)
icing_prob_ds.setncattr('coordinates', geo_coords)
- #icing_prob_ds.setncattr('grid_mapping', 'Projection')
+ icing_prob_ds.setncattr('grid_mapping', 'Projection')
if not use_nan:
icing_prob_ds.setncattr('missing', -1.0)
else:
@@ -1124,7 +1124,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_prob_ds = rootgrp.createVariable('max_icing_probability_column', 'f4', var_dim_list)
icing_prob_ds.setncattr('coordinates', geo_coords)
- #icing_prob_ds.setncattr('grid_mapping', 'Projection')
+ icing_prob_ds.setncattr('grid_mapping', 'Projection')
if not use_nan:
icing_prob_ds.setncattr('missing', -1.0)
else:
@@ -1140,14 +1140,14 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_pred_ds = rootgrp.createVariable('max_icing_probability_level', 'i2', var_dim_list)
icing_pred_ds.setncattr('coordinates', geo_coords)
- #icing_pred_ds.setncattr('grid_mapping', 'Projection')
+ icing_pred_ds.setncattr('grid_mapping', 'Projection')
icing_pred_ds.setncattr('missing', -1)
icing_pred_ds[:,] = max_lvl
if bt_10_4 is not None:
bt_ds = rootgrp.createVariable('bt_10_4', 'f4', var_dim_list)
bt_ds.setncattr('coordinates', geo_coords)
- #bt_ds.setncattr('grid_mapping', 'Projection')
+ bt_ds.setncattr('grid_mapping', 'Projection')
bt_ds[:,] = bt_10_4
lon_ds = rootgrp.createVariable('longitude', 'f4', [dim_1_name, dim_0_name])
@@ -1158,15 +1158,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
lat_ds.units = 'degrees_north'
lat_ds[:,] = lats
- # proj_ds = rootgrp.createVariable('Projection', 'b')
- # proj_ds.setncattr('long_name', long_name)
- # proj_ds.setncattr('grid_mapping_name', 'geostationary')
- # proj_ds.setncattr('sweep_angle_axis', cf_nav_dct['sweep_angle_axis'])
- # proj_ds.setncattr('semi_major_axis', cf_nav_dct['semi_major_axis'])
- # proj_ds.setncattr('semi_minor_axis', cf_nav_dct['semi_minor_axis'])
- # proj_ds.setncattr('inverse_flattening', cf_nav_dct['inverse_flattening'])
- # proj_ds.setncattr('perspective_point_height', cf_nav_dct['perspective_point_height'])
- # proj_ds.setncattr('latitude_of_projection_origin', cf_nav_dct['latitude_of_projection_origin'])
- # proj_ds.setncattr('longitude_of_projection_origin', cf_nav_dct['longitude_of_projection_origin'])
+ proj_ds = rootgrp.createVariable('Projection', 'b')
+ proj_ds.setncattr('grid_mapping_name', 'latitude_longitude')
rootgrp.close()
--
GitLab