def compute_gradient_on_grid(self, grid, affine_grid_to_world=None, verify_mapping=True): resampled_data = resample_image_on_grid(self, grid, affine_grid_to_world, verify_mapping, self.background, self.use_gpu) # create new Image3D object #print resampled_data.max() gradient_data = numpy.gradient(resampled_data) #FIXME: use NiftyPy return gradient_data
def compute_resample_on_grid(self,grid,affine_grid_to_world=None, verify_mapping=True,interpolation_mode=INTERPOLATION_LINEAR): resampled_data = resample_image_on_grid(self, grid, affine_grid_to_world, verify_mapping, self.background, self.use_gpu, interpolation_mode) # create new Image3D object return Image3D(data=resampled_data) #FIXME: perhaps just return the raw resampled data