|
17 | 17 | from ..tractogram_file import HeaderError, HeaderWarning
|
18 | 18 |
|
19 | 19 | from .. import trk as trk_module
|
20 |
| -from ..trk import TrkFile, encode_value_in_name, decode_value_from_name |
| 20 | +from ..trk import TrkFile, encode_value_in_name, decode_value_from_name, get_affine_trackvis_to_rasmm |
21 | 21 | from ..header import Field
|
22 | 22 |
|
23 | 23 | DATA = {}
|
@@ -110,6 +110,17 @@ def trk_with_bytes(self, trk_key='simple_trk_fname', endian='<'):
|
110 | 110 | return trk_struct, trk_bytes
|
111 | 111 |
|
112 | 112 | def test_load_file_with_wrong_information(self):
|
| 113 | + # Simulate a TRK file where `voxel_order` is lowercase. |
| 114 | + trk_struct1, trk_bytes1 = self.trk_with_bytes() |
| 115 | + trk_struct1[Field.VOXEL_ORDER] = b'LAS' |
| 116 | + trk1 = TrkFile.load(BytesIO(trk_bytes1)) |
| 117 | + trk_struct2, trk_bytes2 = self.trk_with_bytes() |
| 118 | + trk_struct2[Field.VOXEL_ORDER] = b'las' |
| 119 | + trk2 = TrkFile.load(BytesIO(trk_bytes2)) |
| 120 | + trk1_aff2rasmm = get_affine_trackvis_to_rasmm(trk1.header) |
| 121 | + trk2_aff2rasmm = get_affine_trackvis_to_rasmm(trk2.header) |
| 122 | + assert_array_equal(trk1_aff2rasmm,trk2_aff2rasmm) |
| 123 | + |
113 | 124 | # Simulate a TRK file where `count` was not provided.
|
114 | 125 | trk_struct, trk_bytes = self.trk_with_bytes()
|
115 | 126 | trk_struct[Field.NB_STREAMLINES] = 0
|
|
0 commit comments