def setUpClass(cls):
"""
Creates GeoTIFF test, multiple single-time and single-variable NetCDF and single multi-time and
multi-variable NetCDF data.
"""
setup_gt_test_data()
setup_nc_multi_test_data()
setup_nc_single_test_data()
def setUp(self):
"""
Retrieves test data filepaths and auxiliary data and creates temporary reference data as NumPy arrays,
xarray arrays and Pandas data frames.
"""
self.gt_filepaths, self.timestamps = setup_gt_test_data()
self.nc_filepaths, _ = setup_nc_multi_test_data()
self.nc_filepath, _ = setup_nc_single_test_data()
self.lon = 5.
self.lat = 44.
sref = osr.SpatialReference()
sref.ImportFromEPSG(4326)
self.sref = sref
row = 970
col = 246
self.x = 4323250.
self.y = 1314750.
self.ref_np_ar = (np.array([[[row + col] * 4]]).T +
np.arange(0, 4)[:, None, None]).astype(float)
xr_ar = xr.DataArray(data=da.array(
self.ref_np_ar.astype(float)).rechunk((1, 1, 1)),
coords={
'time': self.timestamps,
'y': [self.y],
'x': [self.x]
},
dims=['time', 'y', 'x'])
self.ref_xr_ds = xr.Dataset(data_vars={'1': xr_ar})
self.ref_pd_df = self.ref_xr_ds.to_dataframe()
def setUp(self):
"""
Retrieves test data filepaths and auxiliary data and creates temporary reference data as NumPy arrays,
xarray arrays and Pandas data frames.
"""
self.gt_filepaths, self.timestamps = setup_gt_test_data()
self.nc_filepaths, _ = setup_nc_multi_test_data()
self.nc_filepath, _ = setup_nc_single_test_data()
row = 246
col = 970
x = 4323250.
y = 1314750.
row_size = 10
col_size = 16
sres = 500.
# define bounding boxes
# general bounding box for reading data
x_min = x
x_max = x + col_size * sres
y_min = y - row_size * sres
y_max = y
self.bbox = [(x_min, y_min), (x_max, y_max)]
# defines bounding box being partially outside the tile (lower right corner)
x_lr = 4800000.0
y_lr = 1200000.0
x_min = x_lr - col_size * sres
x_max = x_lr + 2 * col_size * sres
y_min = y_lr - 2 * row_size * sres
y_max = y_lr + row_size * sres
self.partial_outside_bbox = [(x_min, y_min), (x_max, y_max)]
self.outside_bbox = [(0, 0), (1, 1)]
rows, cols = np.meshgrid(np.arange(row, row + row_size),
np.arange(col, col + col_size),
indexing='ij')
xs = np.arange(x, x + col_size * sres, sres)
ys = np.arange(y, y - row_size * sres, -sres)
base_np_ar_2D = rows + cols
base_np_ar = np.stack([base_np_ar_2D] * 4, axis=0)
self.ref_np_ar_area = (base_np_ar +
np.arange(0, 4)[:, None, None]).astype(float)
xr_ar = xr.DataArray(data=da.array(
self.ref_np_ar_area.astype(float)).rechunk(
(1, row_size, col_size)),
coords={
'time': self.timestamps,
'y': ys,
'x': xs
},
dims=['time', 'y', 'x'])
self.ref_xr_ds_area = xr.Dataset(data_vars={'1': xr_ar})
self.ref_pd_df_area = self.ref_xr_ds_area.to_dataframe()
def setUp(self):
"""
Retrieves test data filepaths and auxiliary data and creates temporary reference data as NumPy arrays,
xarray arrays and Pandas data frames.
"""
self.gt_filepaths, self.timestamps = setup_gt_test_data()
self.nc_filepaths, _ = setup_nc_multi_test_data()
self.nc_filepath, _ = setup_nc_single_test_data()
self.row = 970
self.col = 246
self.row_size = 10
self.col_size = 16
x = 4323250.
y = 1314750.
sres = 500.
self.ref_np_ar = (np.array([[[self.row + self.col] * 4]]).T +
np.arange(0, 4)[:, None, None]).astype(float)
xr_ar = xr.DataArray(data=self.ref_np_ar.astype(float),
coords={
'time': self.timestamps,
'y': [y],
'x': [x]
},
dims=['time', 'y', 'x'])
self.ref_xr_ds = xr.Dataset(data_vars={'1': xr_ar})
self.ref_pd_df = self.ref_xr_ds.to_dataframe()
rows, cols = np.meshgrid(np.arange(self.row, self.row + self.row_size),
np.arange(self.col, self.col + self.col_size),
indexing='ij')
xs = np.arange(x, x + self.col_size * sres, sres)
ys = np.arange(y, y - self.row_size * sres, -sres)
base_np_ar_2D = rows + cols
base_np_ar = np.stack([base_np_ar_2D] * 4, axis=0)
self.ref_np_ar_area = (base_np_ar +
np.arange(0, 4)[:, None, None]).astype(float)
xr_ar = xr.DataArray(data=da.array(
self.ref_np_ar_area.astype(float)).rechunk(
(1, self.row_size, self.col_size)),
coords={
'time': self.timestamps,
'y': ys,
'x': xs
},
dims=['time', 'y', 'x'])
self.ref_xr_ds_area = xr.Dataset(data_vars={'1': xr_ar})
self.ref_pd_df_area = self.ref_xr_ds_area.to_dataframe()
def setUp(self):
""" Retrieves test data filepaths and auxiliary data. """
self.filepaths, self.timestamps = setup_gt_test_data()
self.data_dirpath = os.path.join(dirpath_test(), 'data',
'Sentinel-1_CSAR')
def setUpClass(cls):
""" Creates GeoTIFF test data. """
setup_gt_test_data()
def setUp(self):
""" Retrieves test data filepaths and auxiliary data. """
self.gt_filepaths, self.timestamps = setup_gt_test_data()