def test_reduce(self):
data = numpy.fromfunction(lambda y, x: (y + x), (1000, 1000))
lons = numpy.fromfunction(lambda y, x: -180 + (360.0 / 1000) * x, (1000, 1000))
lats = numpy.fromfunction(lambda y, x: -90 + (180.0 / 1000) * y, (1000, 1000))
grid_lons, grid_lats = self.area_def.get_lonlats()
lons, lats, data = data_reduce.swath_from_lonlat_grid(grid_lons, grid_lats, lons, lats, data, 7000)
cross_sum = data.sum()
expected = 20514375.0
self.assertAlmostEqual(cross_sum, expected, msg="Reduce data failed")
def test_reduce(self):
data = np.fromfunction(lambda y, x: (y + x), (1000, 1000))
lons = np.fromfunction(lambda y, x: -180 + (360.0 / 1000) * x,
(1000, 1000))
lats = np.fromfunction(lambda y, x: -90 + (180.0 / 1000) * y,
(1000, 1000))
grid_lons, grid_lats = self.area_def.get_lonlats()
lons, lats, data = swath_from_lonlat_grid(grid_lons, grid_lats, lons,
lats, data, 7000)
cross_sum = data.sum()
expected = 20685125.0
self.assertAlmostEqual(cross_sum, expected)
def test_reduce_north_pole(self):
"""Test reducing around the poles."""
from pyresample import utils
area_id = 'ease_sh'
description = 'Antarctic EASE grid'
proj_id = 'ease_sh'
projection = '+proj=laea +lat_0=-90 +lon_0=0 +a=6371228.0 +units=m'
x_size = 425
y_size = 425
area_extent = (-5326849.0625, -5326849.0625, 5326849.0625,
5326849.0625)
area_def = utils.get_area_def(area_id, description, proj_id,
projection, x_size, y_size, area_extent)
grid_lons, grid_lats = area_def.get_lonlats()
area_id = 'ease_sh'
description = 'Antarctic EASE grid'
proj_id = 'ease_sh'
projection = '+proj=laea +lat_0=-90 +lon_0=0 +a=6371228.0 +units=m'
x_size = 1000
y_size = 1000
area_extent = (-532684.0625, -532684.0625, 532684.0625, 532684.0625)
smaller_area_def = utils.get_area_def(area_id, description, proj_id,
projection, x_size, y_size,
area_extent)
data = np.fromfunction(lambda y, x: (y + x), (1000, 1000))
lons, lats = smaller_area_def.get_lonlats()
lons, lats, data = swath_from_lonlat_grid(grid_lons, grid_lats, lons,
lats, data, 7000)
cross_sum = data.sum()
expected = 999000000.0
self.assertAlmostEqual(cross_sum, expected)
def test_reduce_north_pole(self):
"""Test reducing around the poles."""
from pyresample import utils
area_id = 'ease_sh'
description = 'Antarctic EASE grid'
proj_id = 'ease_sh'
projection = '+proj=laea +lat_0=-90 +lon_0=0 +a=6371228.0 +units=m'
x_size = 425
y_size = 425
area_extent = (-5326849.0625, -5326849.0625,
5326849.0625, 5326849.0625)
area_def = utils.get_area_def(area_id, description, proj_id,
projection, x_size, y_size, area_extent)
grid_lons, grid_lats = area_def.get_lonlats()
area_id = 'ease_sh'
description = 'Antarctic EASE grid'
proj_id = 'ease_sh'
projection = '+proj=laea +lat_0=-90 +lon_0=0 +a=6371228.0 +units=m'
x_size = 1000
y_size = 1000
area_extent = (-532684.0625, -532684.0625, 532684.0625, 532684.0625)
smaller_area_def = utils.get_area_def(area_id, description, proj_id,
projection, x_size, y_size,
area_extent)
data = np.fromfunction(lambda y, x: (y + x), (1000, 1000))
lons, lats = smaller_area_def.get_lonlats()
lons, lats, data = swath_from_lonlat_grid(grid_lons, grid_lats,
lons, lats, data, 7000)
cross_sum = data.sum()
expected = 999000000.0
self.assertAlmostEqual(cross_sum, expected)
