def METHOD_NAME(coord_system):
    r"""Run the xyz extraction test sequence for a given
    xyz extraction function"""
    expected_xyz_path = download.get_file_within_data_zip(
        "dicom_dose_test_data.zip", "expected_{}_xyz.json".format(coord_system.lower())
    )
    with open(expected_xyz_path) as fp:
        expected_xyz = json.load(fp)
    assert set(expected_xyz.keys()) == set(ORIENTATIONS_SUPPORTED)
    test_ds_dict = {
        key: pydicom.dcmread(get_data_file(key)) for key in ORIENTATIONS_SUPPORTED
    }
    for orient, dicom in test_ds_dict.items():
        test_xyz = coords.xyz_axes_from_dataset(dicom, coord_system)
        expected_xyz[orient] = np.array(expected_xyz[orient], dtype=object)
        # These tests were being skipped in the previous code
        assert np.array_equal(test_xyz[0], expected_xyz[orient][0])
        assert np.array_equal(test_xyz[1], expected_xyz[orient][1])
        assert np.array_equal(test_xyz[2], expected_xyz[orient][2])