surf = plot.surf_points
inter_z = interpolate.LinearNDInterpolator(nodes, nodeDisplacementZ, fill_value=0.0, rescale=True)
extrapo_z = inter_z.__call__(surf)
plot.plot3dpoint(extrapo_z)
#####
plot3 = Visualization.DataVisulization(ndimage.gaussian_filter(newim3_map,5), 110)
plot3.contour3d()
plot4 = Visualization.DataVisulization(ndimage.gaussian_filter(newim3,5), 110)
plot4.contour3d()
#src3 = mlab.pipeline.scalar_field(ndimage.gaussian_filter(newim3,3))
#mlab.pipeline.iso_surface(src3, contours=[70], colormap='Oranges')
save = SaveLoad.saveDataBase("/Users/junchaowei/Desktop/SpaceRegistration_000_125/")
datadict={}
datadict['volume000'] = orig_volume1
datadict['volume125'] = orig_volume2
datadict['volume125_regi'] = newimage2
datadict['target_point'] = target
datadict['source_point'] = source
save.save(datadict)
newimage2[newimage2>255.0]=255
newimage2[newimage2==0]=20
###### multi-channel color images ####
sy,sx = im1[:,:,0].shape
registration_check = np.zeros((sy,sx,3))
registration_check[:,:,0] = 255-im1[:,:,70]
nodeDisplacementZ[inode] = new_cord[2]-nodes[inode][2]
surf = plot.surf_points
inter_z = interpolate.LinearNDInterpolator(nodes, nodeDisplacementY, fill_value=0.0, rescale=True)
extrapo_z = inter_z.__call__(surf)
plot.plot3dpoint(extrapo_z)
mlab.colorbar(object=None, title="test", orientation="vertical",nb_labels=20, nv_colors=None, label_fmt=None)
#####
plot3 = Visualization.DataVisulization(ndimage.gaussian_filter(newim3_map,5), 110)
plot3.contour3d()
plot4 = Visualization.DataVisulization(ndimage.gaussian_filter(newim3,5), 110)
plot4.contour3d()
save = SaveLoad.saveDataBase("C:\Users\HT\Desktop/test run/")
datadict={}
datadict['volume000'] = nodeDisplacementX
datadict['volume125'] = nodeDisplacementX
datadict['volume125_regi'] = nodeDisplacementX
datadict['target_point'] = im1
datadict['source_point'] = im2_matched
datadict['source_point'] = newim3
datadict['target_point'] = target
datadict['source_point'] = source
save.save(datadict)
###### multi-channel color images ####
sy,sx = im1[:,:,0].shape
registration_check = np.zeros((sy,sx,3))
