def test_cube_resampling(loadsfile1):
"""Import a cube, then make a smaller and resample, then export the new"""
logger.info("Import SEGY format via SEGYIO")
incube = loadsfile1
newcube = Cube(
xori=460500,
yori=5926100,
zori=1540,
xinc=40,
yinc=40,
zinc=5,
ncol=200,
nrow=100,
nlay=100,
rotation=incube.rotation,
yflip=incube.yflip,
)
newcube.resample(incube, sampling="trilinear", outside_value=10.0)
assert newcube.values.mean() == pytest.approx(5.3107, 0.0001)
assert newcube.values[20, 20, 20] == pytest.approx(10.0, 0.0001)
def test_cube_resampling(loadsfile1):
"""Import a cube, then make a smaller and resample, then export the new"""
logger.info("Import SEGY format via SEGYIO")
incube = loadsfile1
newcube = Cube(
xori=460500,
yori=5926100,
zori=1540,
xinc=40,
yinc=40,
zinc=5,
ncol=200,
nrow=100,
nlay=100,
rotation=incube.rotation,
yflip=incube.yflip,
)
newcube.resample(incube, sampling="trilinear", outside_value=10.0)
tsetup.assert_almostequal(newcube.values.mean(), 5.3107, 0.0001)
tsetup.assert_almostequal(newcube.values[20, 20, 20], 10.0, 0.0001)
newcube.to_file(join(TMD, "cube_resmaple1.segy"))
