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