Creating an array to use in: volumeNode.GetIJKToRASDirections (array)

[ihnorton (Isaiah Norton)]

Here is a post with some example code that should help:

ACPC transform question

Hi Alexis, You can try this code on the Python console: import numpy as np import SampleData def getSampleVolume(): # Download MRHead sampleDataLogic = SampleData.SampleDataLogic() volumeNode = sampleDataLogic.downloadMRHead() return volumeNode def getMatrixToACPC(ac, pc, ih): # Anteroposterior axis pcAc = ac - pc yAxis = pcAc / np.linalg.norm(pcAc) # Lateral axis acIhDir = ih - ac xAxis = np.cross(yAxis, acIhDir) xAxis /= np.linalg.norm(xAxis) …