def test_decimate(self): side1_lons = np.arange(8) side1_lats = np.arange(8) + 30 side2_lons = np.arange(8) + 7 side2_lats = np.arange(8) + 30 + 7 side3_lons = np.arange(8) + 14 side3_lats = np.arange(8) + 30 + 14 side4_lons = np.arange(8) + 21 side4_lats = np.arange(8) + 30 + 21 bond = AreaBoundary((side1_lons, side1_lats), (side2_lons, side2_lats), (side3_lons, side3_lats), (side4_lons, side4_lats)) bond.decimate(5) lons, lats = bond.contour() self.assertTrue(np.allclose(lons, np.array([0, 1, 6, 7, 8, 13, 14, 15, 20, 21, 22, 27]))) self.assertTrue(np.allclose(lats, np.array([30, 31, 36, 37, 38, 43, 44, 45, 50, 51, 52, 57])))
def test_contour(self): side1_lons = np.arange(4) side1_lats = np.arange(4) + 20 side2_lons = np.arange(4) + 3 side2_lats = np.arange(4) + 20 + 3 side3_lons = np.arange(4) + 6 side3_lats = np.arange(4) + 20 + 6 side4_lons = np.arange(4) + 9 side4_lats = np.arange(4) + 20 + 9 bond = AreaBoundary((side1_lons, side1_lats), (side2_lons, side2_lats), (side3_lons, side3_lats), (side4_lons, side4_lats)) lons, lats = bond.contour() self.assertTrue(np.allclose(lons, np.arange(12))) self.assertTrue(np.allclose(lats, np.arange(12) + 20))