added get_valid_midpoints to gridded forecast

csep/core/forecasts.py

@@ -104,6 +104,20 @@ def get_longitudes(self):
         """ Returns the lognitude of the lower left node of the spatial grid """
         return self.region.origins()[:,0]
 
+    def get_valid_midpoints(self):
+        """ Returns the midpoints of the valid testing region
+
+            Returns:
+                lons (numpy.array), lats (numpy.array): two numpy arrays containing the valid midpoints from the forecast
+        """
+        latitudes = []
+        longitudes = []
+        for idx in range(self.region.num_nodes):
+            if self.region.bbox_max[idx] == 0:
+                latitudes.append(self.region.midpoints()[idx,1])
+                longitudes.append(self.region.midpoints()[idx,0])
+        return numpy.array(longitudes), numpy.array(latitudes)
+
     @property
     def polygons(self):
         return self.region.polygons

csep/core/regions.py

@@ -548,7 +548,7 @@ def __init__(self, polygons, dh, name='cartesian2d', mask=None):
         self.name = name
         a, xs, ys = self._build_bitmask_vec()
         # in mask, True = bad value and False = good value
-        self.bbox_mask = a[:, :, 0]
+        self.bbox_mask = a[:,:,0]
         # contains the mapping from polygon_index to the mask
         self.idx_map = a[:,:,1]
         # index values of polygons array into the 2d cartesian grid, based on the midpoint.