def test_cube_slice_auto4d_data():
"""Get cube slice aka Auto4D input, with synthetic/scrambled data"""
xs1 = RegularSurface(XTOP1, fformat="gri")
xs1.describe()
xs1out = ojn(td, "XTOP1.ijxyz")
xs1.to_file(xs1out, fformat="ijxyz")
xs2 = RegularSurface(xs1out, fformat="ijxyz")
assert xs1.values.mean() == pytest.approx(xs2.values.mean(), abs=0.0001)
kube1 = Cube(XCUB1)
kube1.describe()
assert xs2.nactive == 10830
xs2.slice_cube_window(kube1,
sampling="trilinear",
mask=True,
attribute="max")
xs2out1 = ojn(td, "XTOP2_sampled_from_cube.ijxyz")
xs2out2 = ojn(td, "XTOP2_sampled_from_cube.gri")
xs2out3 = ojn(td, "XTOP2_sampled_from_cube.png")
xs2.to_file(xs2out1, fformat="ijxyz")
xs2.to_file(xs2out2)
assert xs2.nactive == 3275 # 3320 # shall be fewer cells
xs2.quickplot(
filename=xs2out3,
colortable="seismic",
title="Auto4D Test",
minmax=(0, 12000),
infotext="Method: max",
)
def test_cube_slice_auto4d_data():
"""Get cube slice aka Auto4D input, with synthetic/scrambled data"""
xs1 = RegularSurface(xtop1, fformat='gri')
xs1.describe()
xs1out = ojn(td, 'xtop1.ijxyz')
xs1.to_file(xs1out, fformat='ijxyz')
xs2 = RegularSurface(xs1out, fformat='ijxyz')
assert xs1.values.mean() == pytest.approx(xs2.values.mean(), abs=0.0001)
kube1 = Cube(xcub1)
kube1.describe()
assert xs2.nactive == 10830
xs2.slice_cube_window(kube1,
sampling='trilinear',
mask=True,
attribute='max')
xs2out1 = ojn(td, 'xtop2_sampled_from_cube.ijxyz')
xs2out2 = ojn(td, 'xtop2_sampled_from_cube.gri')
xs2out3 = ojn(td, 'xtop2_sampled_from_cube.png')
xs2.to_file(xs2out1, fformat='ijxyz')
xs2.to_file(xs2out2)
assert xs2.nactive == 3320 # shall be fewer cells
xs2.quickplot(filename=xs2out3,
colortable='seismic',
title='Auto4D Test',
minmax=(0, 12000),
infotext='Method: max')
