def test_get_output_xy(self):
"""Test calculation of output xy-coordinates."""
from pyresample.bilinear import _get_output_xy
from pyresample._spatial_mp import Proj
proj = Proj(self.target_def.proj_str)
out_x, out_y = _get_output_xy(self.target_def, proj)
self.assertTrue(out_x.all())
self.assertTrue(out_y.all())
def test_get_bounding_corners(self):
proj = Proj(self.target_def.proj4_string)
out_x, out_y = bil._get_output_xy(self.target_def, proj)
in_x, in_y = bil._get_input_xy(self.swath_def, proj, self.input_idxs,
self.idx_ref)
res = bil._get_bounding_corners(in_x, in_y, out_x, out_y,
self.neighbours, self.idx_ref)
for i in range(len(res) - 1):
pt_ = res[i]
for j in range(2):
# Only the sixth output location has four valid corners
self.assertTrue(np.isfinite(pt_[5, j]))
def test_get_bounding_corners(self):
proj = Proj(self.target_def.proj_str)
out_x, out_y = bil._get_output_xy(self.target_def, proj)
in_x, in_y = bil._get_input_xy(self.swath_def, proj,
self.input_idxs, self.idx_ref)
res = bil._get_bounding_corners(in_x, in_y, out_x, out_y,
self.neighbours, self.idx_ref)
for i in range(len(res) - 1):
pt_ = res[i]
for j in range(2):
# Only the sixth output location has four valid corners
self.assertTrue(np.isfinite(pt_[5, j]))
def test_get_bounding_corners(self):
"""Test calculation of bounding corners."""
from pyresample.bilinear import (_get_output_xy, _get_input_xy,
_get_bounding_corners)
from pyresample._spatial_mp import Proj
proj = Proj(self.target_def.proj_str)
out_x, out_y = _get_output_xy(self.target_def, proj)
in_x, in_y = _get_input_xy(self.swath_def, proj, self.input_idxs,
self.idx_ref)
res = _get_bounding_corners(in_x, in_y, out_x, out_y, self.neighbours,
self.idx_ref)
for i in range(len(res) - 1):
pt_ = res[i]
for j in range(2):
# Only the sixth output location has four valid corners
self.assertTrue(np.isfinite(pt_[5, j]))
def test_get_four_closest_corners(self):
"""Test calculation of bounding corners."""
from pyresample.bilinear import (_get_output_xy, _get_input_xy,
_get_four_closest_corners)
from pyresample._spatial_mp import Proj
proj = Proj(self.target_def.proj_str)
out_x, out_y = _get_output_xy(self.target_def)
in_x, in_y = _get_input_xy(self.source_def, proj, self.input_idxs,
self.idx_ref)
(pt_1, pt_2, pt_3,
pt_4), ia_ = _get_four_closest_corners(in_x, in_y, out_x, out_y,
self._neighbours, self.idx_ref)
self.assertTrue(pt_1.shape == pt_2.shape == pt_3.shape == pt_4.shape ==
(self.target_def.size, 2))
self.assertTrue(ia_.shape == (self.target_def.size, 4))
# Check which of the locations has four valid X/Y pairs by
# finding where there are non-NaN values
res = np.sum(pt_1 + pt_2 + pt_3 + pt_4, axis=1)
self.assertEqual(np.sum(~np.isnan(res)), 10)
def test_get_input_xy(self):
proj = Proj(self.target_def.proj4_string)
in_x, in_y = bil._get_output_xy(self.swath_def, proj)
self.assertTrue(in_x.all())
self.assertTrue(in_y.all())
def test_get_output_xy(self):
proj = Proj(self.target_def.proj4_string)
out_x, out_y = bil._get_output_xy(self.target_def, proj)
self.assertTrue(out_x.all())
self.assertTrue(out_y.all())
def test_get_input_xy(self):
proj = Proj(self.target_def.proj_str)
in_x, in_y = bil._get_output_xy(self.swath_def, proj)
self.assertTrue(in_x.all())
self.assertTrue(in_y.all())
def test_get_output_xy(self):
proj = Proj(self.target_def.proj_str)
out_x, out_y = bil._get_output_xy(self.target_def, proj)
self.assertTrue(out_x.all())
self.assertTrue(out_y.all())