def test_fill_col_borders(self):
lons = np.arange(10).reshape((2, 5), order="F")
lats = np.arange(10).reshape((2, 5), order="C")
lines = np.array([2, 7])
cols = np.array([2, 7, 12, 17, 22])
hlines = np.arange(10)
hcols = np.arange(24)
satint = GeoInterpolator((lons, lats), (lines, cols), (hlines, hcols))
satint._fill_col_borders()
np.testing.assert_allclose(satint.tie_data[0], TIES_EXP7)
np.testing.assert_allclose(satint.col_indices,
np.array([0, 2, 7, 12, 17, 22, 23]))
def test_fill_col_borders(self):
lons = np.arange(10).reshape((2, 5), order="F")
lats = np.arange(10).reshape((2, 5), order="C")
lines = np.array([2, 7])
cols = np.array([2, 7, 12, 17, 22])
hlines = np.arange(10)
hcols = np.arange(24)
satint = GeoInterpolator((lons, lats), (lines, cols), (hlines, hcols))
satint._fill_col_borders()
np.testing.assert_allclose(satint.tie_data[0], TIES_EXP7)
np.testing.assert_allclose(satint.col_indices,
np.array([0, 2, 7, 12, 17, 22, 23]))
def test_fill_col_borders(self):
lons = np.arange(10).reshape((2, 5), order="F")
lats = np.arange(10).reshape((2, 5), order="C")
lines = np.array([2, 7])
cols = np.array([2, 7, 12, 17, 22])
hlines = np.arange(10)
hcols = np.arange(24)
satint = GeoInterpolator((lons, lats), (lines, cols), (hlines, hcols))
satint._fill_col_borders()
self.assertTrue(np.allclose(satint.tie_data[0],
np.array([[ 6372937.31273379, 6370997. , 6366146.21816553,
6351605.98629588, 6327412.61244969, 6293626.50067273,
6286869.27831734],
[ 6353136.46335726, 6345786.79166939, 6327412.61244969,
6299445.69529922, 6261968.60390423, 6215087.60607344,
6205711.40650728]])))
self.assertTrue(np.allclose(satint.col_indices,
np.array([ 0, 2, 7, 12, 17, 22, 23])))
