def test_hncc_dnb(self, area, c01, c02, c03):
"""Test the 'hncc_dnb' compositor."""
from satpy.composites.viirs import NCCZinke
comp = NCCZinke('hncc_dnb',
prerequisites=('dnb', ),
standard_name='toa_outgoing_radiance_per_'
'unit_wavelength')
mif = xr.DataArray(da.zeros((5, ), chunks=5) + 0.1,
dims=('y', ),
attrs={
'name': 'moon_illumination_fraction',
'area': area
})
res = comp((c01, c02, c03, mif))
assert isinstance(res, xr.DataArray)
assert isinstance(res.data, da.Array)
assert res.attrs['name'] == 'hncc_dnb'
assert res.attrs['standard_name'] == 'ncc_radiance'
data = res.compute()
unique = np.unique(data)
np.testing.assert_allclose(unique, [
3.48479712e-04, 6.96955799e-04, 1.04543189e-03, 4.75394738e-03,
9.50784532e-03, 1.42617433e-02, 1.50001560e+03, 3.00001560e+03,
4.50001560e+03
])
def test_hncc_dnb(self):
"""Test the 'hncc_dnb' compositor."""
import dask.array as da
import numpy as np
import xarray as xr
from pyresample.geometry import AreaDefinition
from satpy.composites.viirs import NCCZinke
rows = 5
cols = 10
area = AreaDefinition(
'test', 'test', 'test',
{'proj': 'eqc', 'lon_0': 0.0,
'lat_0': 0.0},
cols, rows,
(-20037508.34, -10018754.17, 20037508.34, 10018754.17))
comp = NCCZinke('hncc_dnb', prerequisites=('dnb',),
standard_name='toa_outgoing_radiance_per_'
'unit_wavelength')
dnb = np.zeros((rows, cols)) + 0.25
dnb[3, :] += 0.25
dnb[4:, :] += 0.5
dnb = da.from_array(dnb, chunks=25)
c01 = xr.DataArray(dnb,
dims=('y', 'x'),
attrs={'name': 'DNB', 'area': area})
sza = np.zeros((rows, cols)) + 70.0
sza[:, 3] += 20.0
sza[:, 4:] += 45.0
sza = da.from_array(sza, chunks=25)
c02 = xr.DataArray(sza,
dims=('y', 'x'),
attrs={'name': 'solar_zenith_angle', 'area': area})
lza = np.zeros((rows, cols)) + 70.0
lza[:, 3] += 20.0
lza[:, 4:] += 45.0
lza = da.from_array(lza, chunks=25)
c03 = xr.DataArray(lza,
dims=('y', 'x'),
attrs={'name': 'lunar_zenith_angle', 'area': area})
mif = xr.DataArray(da.zeros((5,), chunks=5) + 0.1,
dims=('y',),
attrs={'name': 'moon_illumination_fraction', 'area': area})
res = comp((c01, c02, c03, mif))
self.assertIsInstance(res, xr.DataArray)
self.assertIsInstance(res.data, da.Array)
self.assertEqual(res.attrs['name'], 'hncc_dnb')
self.assertEqual(res.attrs['standard_name'],
'ncc_radiance')
data = res.compute()
unique = np.unique(data)
np.testing.assert_allclose(
unique, [3.48479712e-04, 6.96955799e-04, 1.04543189e-03, 4.75394738e-03,
9.50784532e-03, 1.42617433e-02, 1.50001560e+03, 3.00001560e+03,
4.50001560e+03])