Hello, I am having some problems using/adapting a script I found by searching the mailing list. The script is pasted below but to summarize it should print out the transformation matrix used to fit one map to another. In principle I would like to do the same but in addition to writing out the matrix I would also like to apply it and save the resulting map. I am currently using chimera 1.4 build 27496 (Aqua on OSX 10.5.7) the first problem I have is just invoking the python interpreter..... if I use: Chimera.app/Contents/MacOS/chimera --nogui fittomap.py Executing fittomap.py... Traceback (most recent call last): File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ __main__.py", line 65, in <module> value = chimeraInit.init(sys.argv) File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ chimeraInit.py", line 407, in init chimera.openModels.open(a, prefixableType=1) File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ chimera/__init__.py", line 1357, in open models = func(filename, *args, **kw) File "/Users/src/Downloads/Chimera.app/Contents/Resources/share/ chimera/__init__.py", line 794, in _openPython loadFunc(sandboxName, fileName, f) File "fittomap.py", line 70, in <module> fit_map_in_map(p1, p2) File "fittomap.py", line 24, in fit_map_in_map from VolumeViewer import Data_Region, full_region, Rendering_Options ImportError: cannot import name Data_Region Error while processing fittomap.py: ImportError: cannot import name Data_Region (see reply log for Python traceback info) # ----------------------------------------------------------------------------- # Example script for fitting one map in another without the graphical user # interface. # # chimera --nogui fitnogui.py # # It can also be run using the graphical Chimera interface using File/ Open. # def fit_map_in_map(map1_path, map2_path, initial_map1_transform = None, map1_threshold = None, ijk_step_size_min = 0.01, # Grid index units ijk_step_size_max = 0.5, # Grid index units max_steps = 100, optimize_translation = True, optimize_rotation = True): # Files have to have file suffix indicating volume format. from VolumeData import open_file data1 = open_file(map1_path) data2 = open_file(map2_path) from VolumeViewer import Data_Region, full_region, Rendering_Options ro = Rendering_Options() map1 = Data_Region(data1, full_region(data1.size), '', ro) map2 = Data_Region(data2, full_region(data2.size), '', ro) if initial_map1_transform: from PDBmatrices.matrices import chimera_xform xf = chimera_xform(initial_map1_transform) map1.surface_model().openState.globalXform(xf) use_threshold = (map1_threshold != None) from FitMap import map_points_and_weights, motion_to_maximum points, point_weights = map_points_and_weights(map1, use_threshold) if len(points) == 0: if use_threshold: print 'No grid points above map threshold.' else: print 'Map has no non-zero values.' return move_tf, stats = motion_to_maximum(points, point_weights, map2, max_steps, ijk_step_size_min, ijk_step_size_max, optimize_translation, optimize_rotation) if initial_map1_transform: from PDBmatrices.matrices import multiply_matrices move_tf = multiply_matrices(move_tf, initial_map1_transform) header = ('\nFit map %s in map %s using %d points\n' % (map1.full_name(), map2.full_name(), stats['points']) + ' correlation = %.4g, overlap = %.4g\n' % (stats['correlation'], stats['overlap']) + ' steps = %d, shift = %.3g, angle = %.3g degrees\n' % (stats['steps'], stats['shift'], stats['angle'])) print header from FitMap.gui import transformation_description tfs = transformation_description(move_tf) print tfs # ----------------------------------------------------------------------------- # Example call fitting map into itself. # p1 = p2 = '/usr/local/src/chimera-demos/volume/examples/groel.mrc' fit_map_in_map(p1, p2)