Scripts: mapdiff.py

File mapdiff.py, 408 bytes (added by goddard, 8 years ago)
Line 
1# Report the minimum and maximum difference between two maps.
2path1 = '/Users/goddard/Downloads/Chimera/EMDB/emd_1733.map'
3path2 = '/Users/goddard/Downloads/Chimera/EMDB/emd_1740.map'
4
5from chimera import openModels as om
6m1 = om.open(path1)[0]
7m2 = om.open(path2)[0]
8
9mdiff = m1.full_matrix() - m2.full_matrix()
10print('Difference %s - %s, max %.3g, min %.3g' % (m1.name, m2.name, mdiff.max(), mdiff.min()))