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))