Merge branch 'sabine_issue64-services-provide-unified-missing-values' into 'master'

services provide harmonized missing values (except for those, that have their own missing value labeled)

See merge request !87

services/stable_nightlights/views.py

@@ -30,7 +30,7 @@ class StableNightlightsView(YearMixin, CommonView):
 
     def __init__(self):
 
-        opts = dict(default_value=-99.,
+        opts = dict(default_value=-999.,
                     out_of_bounds_value=-888.,
                     min_valid=0.,
                     max_valid=2.e6,

services/topography_srtm/views.py

@@ -27,7 +27,7 @@ class TopographySRTMView(RelativeMixin, CommonView):
 
     def __init__(self):
 
-        opts = dict(default_value=-9999.,
+        opts = dict(default_value=-999.,
                     out_of_bounds_value=0.,
                     min_valid=-500.,
                     max_valid=9000.,

utils/extraction_tools.py

@@ -96,8 +96,8 @@ def extract_value(lonvec, latvec, data, lon=-999., lat=-999.,
             res = default_value
         else:
             if relative:
-                tma = extract_value_stats(lonvec, latvec, data, lon, lat, default_value=-9999.,
-                                          out_of_bounds_value=-9999., min_valid=0., max_valid=8050.,
+                tma = extract_value_stats(lonvec, latvec, data, lon, lat, default_value=-999.,
+                                          out_of_bounds_value=-888., min_valid=0., max_valid=8050.,
                                           radius=5., agg=np.min)
                 tma = 0. if tma < -9000. else tma
                 res = 0. if res < 0. else res-tma
@@ -238,8 +238,8 @@ def extract_value_stats(lonvec, latvec, data, lon=-999., lat=-999.,
         if subset is not None:
             res = agg(subset)
             if relative:
-                tma = extract_value_stats(lonvec, latvec, data, lon, lat, default_value=-9999.,
-                                          out_of_bounds_value=-9999., min_valid=0., max_valid=8050.,
+                tma = extract_value_stats(lonvec, latvec, data, lon, lat, default_value=-999.,
+                                          out_of_bounds_value=-888., min_valid=0., max_valid=8050.,
                                           radius=5., agg=np.min)
                 tma = 0. if tma < -9000. else tma
                 res = 0. if res < 0. else res-tma