def METHOD_NAME(lt_ctx):
    data = _mk_random(size=(16, 16, 16, 16))
    dataset = MemoryDataSet(
        data=data,
        tileshape=(1, 16, 16),
        num_partitions=2,
        sig_dims=2
    )
    analysis = lt_ctx.create_pick_analysis(dataset=dataset, x=8, y=7)
    result = lt_ctx.run(analysis)
    assert result.intensity.raw_data.shape == (16, 16)
    assert np.allclose(result.intensity.raw_data, data[7, 8])
    assert np.allclose(result.intensity_lin.raw_data, data[7, 8])