def test_CCI_SM_v33_025Img_img_reading_2D():
"""
2D test for the read function of the CCI_SM_v042_025Img class
"""
parameter = ['sm']
img_c = CCI_SM_025Img(os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v03.3", "combined", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20160101000000-fv03.3.nc"),
parameter=parameter)
image_c = img_c.read()
assert sorted(image_c.data.keys()) == sorted(parameter)
assert image_c.data['sm'].shape == (720, 1440)
assert image_c.lon[0, 0] == -179.875
assert image_c.lon[0, 1439] == 179.875
assert image_c.lat[0, 0] == 89.875
assert image_c.lat[719, 0] == -89.875
nptest.assert_almost_equal(image_c.data['sm'][203, 693], 0.23484, 5)
assert image_c.lon.shape == image_c.lat.shape == (720, 1440)
parameter = ['sm']
img_a = CCI_SM_025Img(os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v03.3", "active", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMS-ACTIVE-20160101000000-fv03.3.nc"),
parameter=parameter)
image_a = img_a.read()
assert sorted(image_a.data.keys()) == sorted(parameter)
assert image_a.data['sm'].shape == (720, 1440)
assert image_a.lon[0, 0] == -179.875
assert image_a.lon[0, 1439] == 179.875
assert image_a.lat[0, 0] == 89.875
assert image_a.lat[719, 0] == -89.875
nptest.assert_almost_equal(image_a.data['sm'][203, 693], 67.70157, 5)
assert image_a.lon.shape == image_a.lat.shape == (720, 1440)
parameter = ['sm']
img_p = CCI_SM_025Img(os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v03.3", "passive", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-PASSIVE-20160101000000-fv03.3.nc"),
parameter=parameter)
image_p = img_p.read()
assert sorted(image_p.data.keys()) == sorted(parameter)
assert image_p.data['sm'].shape == (720, 1440)
assert image_p.lon[0, 0] == -179.875
assert image_p.lon[0, 1439] == 179.875
assert image_p.lat[0, 0] == 89.875
assert image_p.lat[719, 0] == -89.875
nptest.assert_almost_equal(image_p.data['sm'][203, 693], 0.322685, 5)
assert image_p.lon.shape == image_p.lat.shape == (720, 1440)
def test_CCI_SM_v052_025Img_img_reading_2D():
"""
2D test for the read function of the CCI_SM_image class
"""
parameter = ['sm']
img_c = CCI_SM_025Img(os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v05.2", "combined", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20160607000000-fv05.2.nc"),
parameter=parameter)
image_c = img_c.read()
assert sorted(image_c.data.keys()) == sorted(parameter)
assert image_c.data['sm'].shape == (720, 1440)
assert image_c.lon[0, 0] == -179.875
assert image_c.lon[0, 1439] == 179.875
assert image_c.lat[0, 0] == 89.875
assert image_c.lat[719, 0] == -89.875
assert abs(image_c.data['sm'][273, 693] - 0.06405) <= 1e-5
assert image_c.lon.shape == image_c.lat.shape == (720, 1440)
parameter = ['sm']
img_a = CCI_SM_025Img(os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v05.2", "active", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMS-ACTIVE-20160606000000-fv05.2.nc"),
parameter=parameter)
image_a = img_a.read()
assert sorted(image_a.data.keys()) == sorted(parameter)
assert image_a.data['sm'].shape == (720, 1440)
assert image_a.lon[0, 0] == -179.875
assert image_a.lon[0, 1439] == 179.875
assert image_a.lat[0, 0] == 89.875
assert image_a.lat[719, 0] == -89.875
assert abs(image_a.data['sm'][273, 693] - 0.65197) <= 1e-5
assert image_a.lon.shape == image_a.lat.shape == (720, 1440)
parameter = ['sm']
img_p = CCI_SM_025Img(os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v05.2", "passive", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-PASSIVE-20160606000000-fv05.2.nc"),
parameter=parameter)
image_p = img_p.read()
assert sorted(image_p.data.keys()) == sorted(parameter)
assert image_p.data['sm'].shape == (720, 1440)
assert image_p.lon[0, 0] == -179.875
assert image_p.lon[0, 1439] == 179.875
assert image_p.lat[0, 0] == 89.875
assert image_p.lat[719, 0] == -89.875
assert abs(image_p.data['sm'][273, 693] - 0.0600) <= 1e-5
assert image_p.lon.shape == image_p.lat.shape == (720, 1440)
def test_CCI_SM_v022_025Img_img_reading_1D():
"""
1D test for the read function of the CCI_SM_image class
"""
parameter = ['sm']
img_c = CCI_SM_025Img(
os.path.join(os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v02.2", "combined", "2014",
"ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20140101000000-fv02.2.nc"),
parameter=parameter, array_1D=True)
image_c = img_c.read()
assert image_c.lon[901440] == -179.875
assert image_c.lat[901440] == 66.625
assert sorted(image_c.data.keys()) == sorted(parameter)
nptest.assert_almost_equal(image_c.data['sm'][1440 * (720 - 274) + 693],
0.142200, 5)
parameter = ['sm']
img_a = CCI_SM_025Img(
os.path.join(os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v02.2", "active", "2014",
"ESACCI-SOILMOISTURE-L3S-SSMS-ACTIVE-20140101000000-fv02.2.nc"),
parameter=parameter, array_1D=True)
image_a = img_a.read()
assert image_a.lon[901440] == -179.875
assert image_a.lat[901440] == 66.625
assert sorted(image_a.data.keys()) == sorted(parameter)
nptest.assert_almost_equal(image_a.data['sm'][1440 * (720 - 274) + 693],
13.030, 5)
parameter = ['sm']
img_p = CCI_SM_025Img(
os.path.join(os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v02.2", "passive", "2014",
"ESACCI-SOILMOISTURE-L3S-SSMV-PASSIVE-20140101000000-fv02.2.nc"),
parameter=parameter, array_1D=True)
image_p = img_p.read()
assert image_p.lon[901440] == -179.875
assert image_p.lat[901440] == 66.625
assert sorted(image_p.data.keys()) == sorted(parameter)
nptest.assert_almost_equal(image_p.data['sm'][1440 * (720 - 274) + 693],
0.0600, 5)
def test_CCI_SM_v042_025Img_img_reading_1D():
"""
1D test for the read function of the CCI_SM_v042_025Img class
_a for active
_p for passive
_c for combined
"""
parameter = ['sm']
img_a = CCI_SM_025Img(
os.path.join(os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v04.2", "active", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMS-ACTIVE-20160606000000-fv04.2.nc"),
parameter=parameter,
array_1D=True)
image_a = img_a.read()
assert image_a.lon[901440] == -179.875
assert image_a.lat[901440] == 66.625
assert sorted(image_a.data.keys()) == sorted(parameter)
assert abs(image_a.data['sm'][901440] - 100) <= 1e-3
img_p = CCI_SM_025Img(
os.path.join(os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v04.2", "passive", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-PASSIVE-20160606000000-fv04.2.nc"),
parameter=parameter,
array_1D=True)
image_p = img_p.read()
assert image_p.lon[911520] == -179.875
assert image_p.lat[911520] == 68.375
assert sorted(image_p.data.keys()) == sorted(parameter)
assert abs(image_p.data['sm'][911520] - 0.31) <= 1e-3
img_c = CCI_SM_025Img(
os.path.join(os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v04.2", "combined", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20160607000000-fv04.2.nc"),
parameter=parameter,
array_1D=True)
image_c = img_c.read()
assert image_c.lon[912960] == -179.875
assert image_c.lat[912960] == 68.625
assert sorted(image_c.data.keys()) == sorted(parameter)
assert abs(image_c.data['sm'][912960] - 0.1755) <= 1e-3
def test_CCI_SM_v052_025Img_img_reading_1D_combined():
"""
1D test for the read function of the CCI_SM_image class
"""
##### global grid
parameter = ['sm']
filename = os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v05.2", "combined", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20160607000000-fv05.2.nc")
img_c = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=None,
array_1D=True)
image_c = img_c.read()
assert img_c.grid.find_nearest_gpi(75.625, 14.625) == (602942, 0)
ref_lat = image_c.lat[434462]
ref_lon = image_c.lon[434462]
assert ref_lon == 75.625
assert ref_lat == 14.625
assert sorted(image_c.data.keys()) == sorted(parameter)
ref_sm = image_c.data['sm'][1440 * 301 + 1022]
nptest.assert_almost_equal(ref_sm, 0.37016, 5)
###### land grid
land_grid = CCILandGrid()
img_c = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=land_grid,
array_1D=True)
assert img_c.grid.find_nearest_gpi(75.625, 14.625) == (602942, 0)
image_c = img_c.read()
sm = image_c.data['sm'][177048]
lat = image_c.lat[177048]
lon = image_c.lon[177048]
assert ref_lat == lat
assert ref_lon == lon
nptest.assert_almost_equal(ref_sm, sm, 5)
def test_CCI_SM_v042_025Img_img_reading_1D_active():
"""
1D test for the read function of the CCI_SM_v042_025Img class
_a for active
_p for passive
_c for combined
"""
parameter = ['sm']
filename = os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v04.2", "active", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMS-ACTIVE-20160606000000-fv04.2.nc")
img_a = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=None,
array_1D=True)
image_a = img_a.read()
ref_lat = image_a.lat[1440 * 93] # = 133920
ref_lon = image_a.lon[1440 * 93]
assert ref_lon == -179.875
assert ref_lat == 66.625
assert sorted(image_a.data.keys()) == sorted(parameter)
ref_sm = image_a.data['sm'][133920]
assert abs(ref_sm - 100.) <= 1e-3
###land_grid
land_grid = CCILandGrid()
img_a = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=land_grid,
array_1D=True)
assert img_a.grid.find_nearest_gpi(-179.875, 66.625) == (901440, 0)
image_a = img_a.read()
sm = image_a.data['sm'][32830]
lat = image_a.lat[32830]
lon = image_a.lon[32830]
assert ref_lat == lat
assert ref_lon == lon
nptest.assert_almost_equal(ref_sm, sm, 5)
def test_CCI_SM_v042_025Img_img_reading_1D_combined():
"""
1D test for the read function of the CCI_SM_v042_025Img class
_a for active
_p for passive
_c for combined
"""
parameter = ['sm']
filename = os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v04.2", "combined", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20160607000000-fv04.2.nc")
img_c = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=None,
array_1D=True)
image_c = img_c.read()
ref_lat = image_c.lat[1440 * 85] # = 122400
ref_lon = image_c.lon[1440 * 85]
assert ref_lon == -179.875
assert ref_lat == 68.625
assert sorted(image_c.data.keys()) == sorted(parameter)
ref_sm = image_c.data['sm'][122400]
assert abs(ref_sm - 0.1755) <= 1e-3
###land_grid
land_grid = CCILandGrid()
img_c = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=land_grid,
array_1D=True)
assert img_c.grid.find_nearest_gpi(-179.875, 68.625) == (912960, 0)
image_c = img_c.read()
sm = image_c.data['sm'][24434]
lat = image_c.lat[24434]
lon = image_c.lon[24434]
assert ref_lat == lat
assert ref_lon == lon
nptest.assert_almost_equal(ref_sm, sm, 5)
def test_CCI_SM_v033_025Img_img_reading_1D_combined():
"""
1D test for the read function of the CCI_SM_v042_025Img class
"""
filename = os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v03.3", "combined", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20160101000000-fv03.3.nc")
parameter = ['sm']
img_c = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=None,
array_1D=True)
image_c = img_c.read()
ref_lat = image_c.lat[1440 * 273 + 693]
ref_lon = image_c.lon[1440 * 273 + 693]
assert ref_lon == -6.625
assert ref_lat == 21.625
assert sorted(image_c.data.keys()) == sorted(parameter)
ref_sm = image_c.data['sm'][1440 * 273 + 693]
nptest.assert_almost_equal(ref_sm, 0.142998, 5)
###### land grid
land_grid = CCILandGrid()
img_c = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=land_grid,
array_1D=True)
assert img_c.grid.find_nearest_gpi(-6.625, 21.625) == (642933, 0)
image_c = img_c.read()
sm = image_c.data['sm'][164759]
lat = image_c.lat[164759]
lon = image_c.lon[164759]
assert ref_lat == lat
assert ref_lon == lon
nptest.assert_almost_equal(ref_sm, sm, 5)
def test_CCI_SM_v033_025Img_img_reading_1D_active():
filename = os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v03.3", "active", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMS-ACTIVE-20160101000000-fv03.3.nc")
parameter = ['sm']
img_a = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=None,
array_1D=True)
image_a = img_a.read()
ref_lat = image_a.lat[1440 * 273 + 693]
ref_lon = image_a.lon[1440 * 273 + 693]
assert ref_lon == -6.625
assert ref_lat == 21.625
assert sorted(image_a.data.keys()) == sorted(parameter)
ref_sm = image_a.data['sm'][1440 * 273 + 693]
nptest.assert_almost_equal(ref_sm, 18.92771, 5)
image_a = img_a.read()
###### land grid
land_grid = CCILandGrid()
img_a = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=land_grid,
array_1D=True)
assert img_a.grid.find_nearest_gpi(-6.625, 21.625) == (642933, 0)
image_a = img_a.read()
sm = image_a.data['sm'][164759]
lat = image_a.lat[164759]
lon = image_a.lon[164759]
assert ref_lat == lat
assert ref_lon == lon
nptest.assert_almost_equal(ref_sm, sm, 5)
def test_CCI_SM_v052_025Img_img_reading_1D_passive():
parameter = ['sm']
filename = os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v05.2", "passive", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-PASSIVE-20160606000000-fv05.2.nc")
img_p = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=None,
array_1D=True)
image_p = img_p.read()
assert img_p.grid.find_nearest_gpi(-6.625, 21.625) == (642933, 0)
ref_lat = image_p.lat[393813]
ref_lon = image_p.lon[393813]
assert ref_lon == -6.625
assert ref_lat == 21.625
assert sorted(image_p.data.keys()) == sorted(parameter)
ref_sm = image_p.data['sm'][1440 * 273 + 693]
nptest.assert_almost_equal(ref_sm, 0.0600, 5)
###### land grid
land_grid = CCILandGrid()
img_p = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=land_grid,
array_1D=True)
assert img_p.grid.find_nearest_gpi(-6.625, 21.625) == (642933, 0)
image_p = img_p.read()
sm = image_p.data['sm'][164759]
lat = image_p.lat[164759]
lon = image_p.lon[164759]
assert ref_lat == lat
assert ref_lon == lon
nptest.assert_almost_equal(ref_sm, sm, 5)
def test_CCI_SM_v042_025Img_img_reading_1D_passive():
parameter = ['sm']
filename = os.path.join(
os.path.dirname(__file__), "esa_cci_sm-test-data",
"esa_cci_sm_dailyImages", "v04.2", "passive", "2016",
"ESACCI-SOILMOISTURE-L3S-SSMV-PASSIVE-20160606000000-fv04.2.nc")
img_p = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=None,
array_1D=True)
image_p = img_p.read()
ref_lat = image_p.lat[1440 * 86] # = 123840
ref_lon = image_p.lon[1440 * 86]
assert ref_lon == -179.875
assert ref_lat == 68.375
assert sorted(image_p.data.keys()) == sorted(parameter)
ref_sm = image_p.data['sm'][123840]
assert abs(ref_sm - 0.31) <= 1e-3
###land_grid
land_grid = CCILandGrid()
img_a = CCI_SM_025Img(filename=filename,
parameter=parameter,
subgrid=land_grid,
array_1D=True)
assert img_a.grid.find_nearest_gpi(-179.875, 68.375) == (911520, 0)
image_a = img_a.read()
sm = image_a.data['sm'][25410]
lat = image_a.lat[25410]
lon = image_a.lon[25410]
assert ref_lat == lat
assert ref_lon == lon
nptest.assert_almost_equal(ref_sm, sm, 5)