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"))