

In [1]:

    
import NiftyCore
import numpy
import occiput; occiput.initialise()
import ctypes
from numpy import pi



In [2]:

    
imA = occiput.DataSources.FileSource('/Users/spedemon/Desktop/notebooks/Ralpfh2.nii')
imA.data = imA.data[range(0,imA.data.shape[0],2),0:200,:]

Rigid Transformation:



In [18]:

    
r = NiftyCore.NiftyReg.resample_image_rigid(imA.data, [0,0,0], [0.3,0.0,0.0], 0.5*numpy.asarray(imA.shape), None )



In [19]:

    
occiput.Visualization.MultipleVolumes((imA.data,r))

Spatial gradient:



In [20]:

    
g = NiftyCore.NiftyReg.deriv_intensity_wrt_space_rigid(imA.data, [0,0,0], [0.3,0,0], 0.5*numpy.asarray(imA.shape), None )



In [17]:

    
occiput.Visualization.MultipleVolumes([rx,g[:,:,:,0],g[:,:,:,1],g[:,:,:,2]],scales=(1,3,3,3))

Derivatives of rigid tranformation parameters



In [21]:

    
gt = NiftyCore.NiftyReg.deriv_intensity_wrt_transformation_rigid(imA.data, [0,0,0], [0.3,0,0], 0.5*numpy.asarray(imA.shape), None )









    



---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-21-c0e274f2154b> in <module>()
----> 1 gt = NiftyCore.NiftyReg.deriv_intensity_wrt_transformation_rigid(imA.data, [0,0,0], [0.3,0,0], 0.5*numpy.asarray(imA.shape), None )

TypeError: deriv_intensity_wrt_transformation_rigid() takes at most 2 arguments (5 given)



In [22]:

    
occiput.Visualization.MultipleVolumes((r,(gt[:,:,:,3])),scales=(1,10))









    



---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-22-1feaedc95fa6> in <module>()
----> 1 occiput.Visualization.MultipleVolumes((r,(gt[:,:,:,3])),scales=(1,10))

NameError: name 'gt' is not defined



In [9]:



In [9]:



In [9]:



In [9]:



In [9]:



In [ ]: