volume morph reverts orientation

[jaremko@…]

The following bug report has been submitted:
Platform:        macOS-10.16-x86_64-i386-64bit
ChimeraX Version: 1.2.dev202103292307 (2021-03-29 23:07:54 UTC)
Description
volume morph reverts mrc realigned and moved maps back to original orientation

Log:
> ui mousemode right zoom

> set bgColor white

> graphics silhouettes true

> graphics silhouettes width 3

> hide all atoms

> style stick

Changed 0 atom styles  

> style ions sphere

Changed 0 atom styles  

> style ions sphere

Changed 0 atom styles  

> show all cartoons

> show ligand target ab

> show ions atoms

> show sidechain & (ligand | ions) :< 3.5 target ab

> color /a #55c087

> color /b #ffe255

> color /c #f6986c

> color /d #6b80bc

> color /e #aa7fba

> color /g #FFC0CB

> color /h #b2b2b2

> color /i red

> color byhetero

> lighting soft

> lighting depthCue false

> cartoon suppressBackboneDisplay false

> volume all color #b2b2b2

Error running startup command 'vol all color #b2b2b2': No volumes specified  
UCSF ChimeraX version: 1.2.dev202103292307 (2021-03-29)  
© 2016-2020 Regents of the University of California. All rights reserved.  
How to cite UCSF ChimeraX  

> open
> /Users/jaremko/Documents/cryoem/2021-03-25_x139-Hs128-DNA/Refine3D/job022/run_class001.mrc

Opened run_class001.mrc as #1, grid size 320,320,320, pixel 1.1, shown at
level 0.00546, step 2, values float32  

> volume #1 step 1

> volume #1 level 0.007836

> volume #1 level 0.00654

> volume #1 level 0.007296

> open
> /Users/jaremko/Documents/cryoem/2021-03-26_x138-Hs145-128-60bpDNA/Refine3D/job055/run_class001.mrc
> format mrc

Opened run_class001.mrc as #2, grid size 320,320,320, pixel 1.1, shown at
level 0.00375, step 2, values float32  

> select #2

2 models selected  

> ui mousemode right "rotate selected models"

> view matrix models
> #2,-0.53137,-0.83764,-0.12652,437.86,-0.82872,0.54496,-0.12741,247.13,0.17567,0.037148,-0.98375,294.55

> view matrix models
> #2,-0.99188,-0.071703,0.10499,346.83,-0.069815,0.99733,0.021556,9.1917,-0.10626,0.014051,-0.99424,350.14

> view matrix models
> #2,-0.97404,-0.050809,0.22059,320.68,-0.093099,0.97817,-0.18579,51.356,-0.20633,-0.2015,-0.95751,399.28

> view matrix models
> #2,-0.96225,0.067394,0.2637,290.75,0.022339,0.98515,-0.17026,27.145,-0.27126,-0.15794,-0.94946,401.8

> ui mousemode right select

> ui mousemode right zoom

> ui tool show "Fit in Map"

Fit map run_class001.mrc in map run_class001.mrc using 40955 points  
correlation = 0.6429, correlation about mean = 0.09574, overlap = 4.262  
steps = 120, shift = 19, angle = 11.1 degrees  
  
Position of run_class001.mrc (#2) relative to run_class001.mrc (#1)
coordinates:  
Matrix rotation and translation  
-0.98677979 -0.03267663 0.15873844 323.65384861  
-0.03788562 0.99883467 -0.02989961 7.74393040  
-0.15757644 -0.03551823 -0.98686783 383.60570325  
Axis -0.01775755 0.99970678 -0.01646292  
Axis point 177.16545948 0.00000000 178.98395679  
Rotation angle (degrees) 170.89734474  
Shift along axis -4.32090968  
  

> ui mousemode right "rotate selected models"

> view matrix models
> #2,-0.97789,-0.024501,0.20767,312.47,-0.049843,0.9918,-0.11769,25.78,-0.20308,-0.12544,-0.97109,404.69

> view matrix models
> #2,0.72007,0.69164,0.055882,-87.044,0.69032,-0.7222,0.043408,167.11,0.070381,0.007319,-0.99749,337.64

> view matrix models
> #2,0.98955,-0.08897,0.11345,-8.0841,-0.086052,-0.99583,-0.030374,364.36,0.11567,0.020295,-0.99308,326.63

Fit map run_class001.mrc in map run_class001.mrc using 40955 points  
correlation = 0.9713, correlation about mean = 0.9018, overlap = 8.325  
steps = 120, shift = 9.88, angle = 12.5 degrees  
  
Position of run_class001.mrc (#2) relative to run_class001.mrc (#1)
coordinates:  
Matrix rotation and translation  
0.99955879 0.02796161 -0.01001856 -5.22990750  
0.02890482 -0.99335211 0.11142757 327.84608404  
-0.00683627 -0.11166800 -0.99372206 372.43796645  
Axis -0.99988935 -0.01426269 0.00422737  
Axis point 0.00000000 174.38107268 177.03379747  
Rotation angle (degrees) 173.59474201  
Shift along axis 2.12779371  
  

> ui mousemode right zoom

> select clear

> volume #2 step 1

> volume #2 level 0.008549

> volume #2 level 0.01041

> volume #1 level 0.008408

> volume #2 level 0.01319

> volume morph

Missing or invalid "volumes" argument: empty atom specifier  

> volume #1 level 0.007407

> volume #1 level 0.006295

> volume #1 level 0.006295

> volume #1 level 0.007518

> volume morph

Missing or invalid "volumes" argument: empty atom specifier  

> volume morph #1-2

Opened morph as #3, grid size 320,320,320, pixel 1.1, shown at step 1, values
float32  

> close #3

> show #!2 models

> show #!1 models

> hide #!2 models

> hide #!1 models

> volume morph #1-2

Opened morph as #3, grid size 320,320,320, pixel 1.1, shown at step 1, values
float32  

> close #3

> show #!2 models

> show #!1 models

> save

Missing or invalid "fileName" argument: Expected a file name  




OpenGL version: 4.1 ATI-4.2.15
OpenGL renderer: AMD Radeon Pro 5500M OpenGL Engine
OpenGL vendor: ATI Technologies Inc.Hardware:

    Hardware Overview:

      Model Name: MacBook Pro
      Model Identifier: MacBookPro16,1
      Processor Name: 8-Core Intel Core i9
      Processor Speed: 2.4 GHz
      Number of Processors: 1
      Total Number of Cores: 8
      L2 Cache (per Core): 256 KB
      L3 Cache: 16 MB
      Hyper-Threading Technology: Enabled
      Memory: 32 GB
      System Firmware Version: 1554.80.3.0.0 (iBridge: 18.16.14347.0.0,0)

Software:

    System Software Overview:

      System Version: macOS 11.2.3 (20D91)
      Kernel Version: Darwin 20.3.0
      Time since boot: 20:45

Graphics/Displays:

    Intel UHD Graphics 630:

      Chipset Model: Intel UHD Graphics 630
      Type: GPU
      Bus: Built-In
      VRAM (Dynamic, Max): 1536 MB
      Vendor: Intel
      Device ID: 0x3e9b
      Revision ID: 0x0002
      Automatic Graphics Switching: Supported
      gMux Version: 5.0.0
      Metal Family: Supported, Metal GPUFamily macOS 2

    AMD Radeon Pro 5500M:

      Chipset Model: AMD Radeon Pro 5500M
      Type: GPU
      Bus: PCIe
      PCIe Lane Width: x8
      VRAM (Total): 8 GB
      Vendor: AMD (0x1002)
      Device ID: 0x7340
      Revision ID: 0x0040
      ROM Revision: 113-D3220E-190
      VBIOS Version: 113-D32206U1-020
      Option ROM Version: 113-D32206U1-020
      EFI Driver Version: 01.A1.190
      Automatic Graphics Switching: Supported
      gMux Version: 5.0.0
      Metal Family: Supported, Metal GPUFamily macOS 2
      Displays:
        Color LCD:
          Display Type: Built-In Retina LCD
          Resolution: 3072 x 1920 Retina
          Framebuffer Depth: 24-Bit Color (ARGB8888)
          Main Display: Yes
          Mirror: Off
          Online: Yes
          Automatically Adjust Brightness: Yes
          Connection Type: Internal

Locale: (None, 'UTF-8')
PyQt5 5.15.2, Qt 5.15.2
Installed Packages:
    alabaster: 0.7.12
    appdirs: 1.4.4
    appnope: 0.1.2
    Babel: 2.9.0
    backcall: 0.2.0
    blockdiag: 2.0.1
    certifi: 2020.12.5
    cftime: 1.4.1
    chardet: 3.0.4
    ChimeraX-AddCharge: 1.0
    ChimeraX-AddH: 2.1.4
    ChimeraX-AlignmentAlgorithms: 2.0
    ChimeraX-AlignmentHdrs: 3.2
    ChimeraX-AlignmentMatrices: 2.0
    ChimeraX-Alignments: 2.1
    ChimeraX-AmberInfo: 1.0
    ChimeraX-Arrays: 1.0
    ChimeraX-Atomic: 1.12
    ChimeraX-AtomicLibrary: 3.0
    ChimeraX-AtomSearch: 2.0
    ChimeraX-AtomSearchLibrary: 1.0
    ChimeraX-AxesPlanes: 2.0
    ChimeraX-BasicActions: 1.1
    ChimeraX-BILD: 1.0
    ChimeraX-BlastProtein: 1.1
    ChimeraX-BondRot: 2.0
    ChimeraX-BugReporter: 1.0
    ChimeraX-BuildStructure: 2.5.1
    ChimeraX-Bumps: 1.0
    ChimeraX-BundleBuilder: 1.1
    ChimeraX-ButtonPanel: 1.0
    ChimeraX-CageBuilder: 1.0
    ChimeraX-CellPack: 1.0
    ChimeraX-Centroids: 1.1
    ChimeraX-ChemGroup: 2.0
    ChimeraX-Clashes: 2.1
    ChimeraX-ColorActions: 1.0
    ChimeraX-ColorGlobe: 1.0
    ChimeraX-ColorKey: 1.2.1
    ChimeraX-CommandLine: 1.1.3
    ChimeraX-ConnectStructure: 2.0
    ChimeraX-Contacts: 1.0
    ChimeraX-Core: 1.2.dev202103292307
    ChimeraX-CoreFormats: 1.0
    ChimeraX-coulombic: 1.1
    ChimeraX-Crosslinks: 1.0
    ChimeraX-Crystal: 1.0
    ChimeraX-CrystalContacts: 1.0
    ChimeraX-DataFormats: 1.1
    ChimeraX-Dicom: 1.0
    ChimeraX-DistMonitor: 1.1
    ChimeraX-DistUI: 1.0
    ChimeraX-Dssp: 2.0
    ChimeraX-EMDB-SFF: 1.0
    ChimeraX-ExperimentalCommands: 1.0
    ChimeraX-FileHistory: 1.0
    ChimeraX-FunctionKey: 1.0
    ChimeraX-Geometry: 1.1
    ChimeraX-gltf: 1.0
    ChimeraX-Graphics: 1.0
    ChimeraX-Hbonds: 2.1
    ChimeraX-Help: 1.1
    ChimeraX-HKCage: 1.3
    ChimeraX-IHM: 1.0
    ChimeraX-ImageFormats: 1.1
    ChimeraX-IMOD: 1.0
    ChimeraX-IO: 1.0.1
    ChimeraX-Label: 1.0
    ChimeraX-ListInfo: 1.1
    ChimeraX-Log: 1.1.2
    ChimeraX-LookingGlass: 1.1
    ChimeraX-Maestro: 1.8.1
    ChimeraX-Map: 1.0.2
    ChimeraX-MapData: 2.0
    ChimeraX-MapEraser: 1.0
    ChimeraX-MapFilter: 2.0
    ChimeraX-MapFit: 2.0
    ChimeraX-MapSeries: 2.0
    ChimeraX-Markers: 1.0
    ChimeraX-Mask: 1.0
    ChimeraX-MatchMaker: 1.1
    ChimeraX-MDcrds: 2.2
    ChimeraX-MedicalToolbar: 1.0.1
    ChimeraX-Meeting: 1.0
    ChimeraX-MLP: 1.1
    ChimeraX-mmCIF: 2.3
    ChimeraX-MMTF: 2.1
    ChimeraX-Modeller: 1.0
    ChimeraX-ModelPanel: 1.0.1
    ChimeraX-ModelSeries: 1.0
    ChimeraX-Mol2: 2.0
    ChimeraX-Morph: 1.0
    ChimeraX-MouseModes: 1.1
    ChimeraX-Movie: 1.0
    ChimeraX-Neuron: 1.0
    ChimeraX-Nucleotides: 2.0
    ChimeraX-OpenCommand: 1.5
    ChimeraX-PDB: 2.4
    ChimeraX-PDBBio: 1.0
    ChimeraX-PDBLibrary: 1.0
    ChimeraX-PDBMatrices: 1.0
    ChimeraX-PickBlobs: 1.0
    ChimeraX-Positions: 1.0
    ChimeraX-PresetMgr: 1.0.1
    ChimeraX-PubChem: 2.0.1
    ChimeraX-ReadPbonds: 1.0
    ChimeraX-Registration: 1.1
    ChimeraX-RemoteControl: 1.0
    ChimeraX-ResidueFit: 1.0
    ChimeraX-RestServer: 1.1
    ChimeraX-RNALayout: 1.0
    ChimeraX-RotamerLibMgr: 2.0
    ChimeraX-RotamerLibsDunbrack: 2.0
    ChimeraX-RotamerLibsDynameomics: 2.0
    ChimeraX-RotamerLibsRichardson: 2.0
    ChimeraX-SaveCommand: 1.4
    ChimeraX-SchemeMgr: 1.0
    ChimeraX-SDF: 2.0
    ChimeraX-Segger: 1.0
    ChimeraX-Segment: 1.0
    ChimeraX-SeqView: 2.3
    ChimeraX-Shape: 1.0.1
    ChimeraX-Shell: 1.0
    ChimeraX-Shortcuts: 1.0
    ChimeraX-ShowAttr: 1.0
    ChimeraX-ShowSequences: 1.0
    ChimeraX-SideView: 1.0
    ChimeraX-Smiles: 2.0.1
    ChimeraX-SmoothLines: 1.0
    ChimeraX-SpaceNavigator: 1.0
    ChimeraX-StdCommands: 1.3
    ChimeraX-STL: 1.0
    ChimeraX-Storm: 1.0
    ChimeraX-Struts: 1.0
    ChimeraX-Surface: 1.0
    ChimeraX-SwapAA: 2.0
    ChimeraX-SwapRes: 2.1
    ChimeraX-TapeMeasure: 1.0
    ChimeraX-Test: 1.0
    ChimeraX-Toolbar: 1.0.1
    ChimeraX-ToolshedUtils: 1.1
    ChimeraX-Tug: 1.0
    ChimeraX-UI: 1.7.1
    ChimeraX-uniprot: 2.1
    ChimeraX-UnitCell: 1.0
    ChimeraX-ViewDockX: 1.0
    ChimeraX-Vive: 1.1
    ChimeraX-VolumeMenu: 1.0
    ChimeraX-VTK: 1.0
    ChimeraX-WavefrontOBJ: 1.0
    ChimeraX-WebCam: 1.0
    ChimeraX-WebServices: 1.0
    ChimeraX-Zone: 1.0
    colorama: 0.4.3
    comtypes: 1.1.7
    cxservices: 1.0
    cycler: 0.10.0
    Cython: 0.29.21
    decorator: 4.4.2
    distlib: 0.3.1
    docutils: 0.16
    filelock: 3.0.12
    funcparserlib: 0.3.6
    grako: 3.16.5
    h5py: 2.10.0
    html2text: 2020.1.16
    idna: 2.10
    ihm: 0.17
    imagecodecs: 2020.5.30
    imagesize: 1.2.0
    ipykernel: 5.3.4
    ipython: 7.18.1
    ipython-genutils: 0.2.0
    jedi: 0.17.2
    Jinja2: 2.11.2
    jupyter-client: 6.1.7
    jupyter-core: 4.7.1
    kiwisolver: 1.3.1
    line-profiler: 2.1.2
    lxml: 4.6.2
    lz4: 3.1.0
    MarkupSafe: 1.1.1
    matplotlib: 3.3.2
    msgpack: 1.0.0
    netCDF4: 1.5.4
    networkx: 2.5
    numexpr: 2.7.3
    numpy: 1.19.2
    numpydoc: 1.1.0
    openvr: 1.14.1501
    packaging: 20.9
    ParmEd: 3.2.0
    parso: 0.7.1
    pexpect: 4.8.0
    pickleshare: 0.7.5
    Pillow: 7.2.0
    pip: 21.0.1
    pkginfo: 1.5.0.1
    prompt-toolkit: 3.0.18
    psutil: 5.7.2
    ptyprocess: 0.7.0
    pycollada: 0.7.1
    pydicom: 2.0.0
    Pygments: 2.7.1
    PyOpenGL: 3.1.5
    PyOpenGL-accelerate: 3.1.5
    pyparsing: 2.4.7
    PyQt5-commercial: 5.15.2
    PyQt5-sip: 12.8.1
    PyQtWebEngine-commercial: 5.15.2
    python-dateutil: 2.8.1
    pytz: 2021.1
    pyzmq: 22.0.3
    qtconsole: 4.7.7
    QtPy: 1.9.0
    RandomWords: 0.3.0
    requests: 2.24.0
    scipy: 1.5.2
    setuptools: 50.3.2
    sfftk-rw: 0.6.7.dev1
    six: 1.15.0
    snowballstemmer: 2.1.0
    sortedcontainers: 2.2.2
    Sphinx: 3.2.1
    sphinxcontrib-applehelp: 1.0.2
    sphinxcontrib-blockdiag: 2.0.0
    sphinxcontrib-devhelp: 1.0.2
    sphinxcontrib-htmlhelp: 1.0.3
    sphinxcontrib-jsmath: 1.0.1
    sphinxcontrib-qthelp: 1.0.3
    sphinxcontrib-serializinghtml: 1.1.4
    suds-jurko: 0.6
    tables: 3.6.1
    tifffile: 2020.9.3
    tinyarray: 1.2.3
    tornado: 6.1
    traitlets: 5.0.5
    urllib3: 1.25.11
    wcwidth: 0.2.5
    webcolors: 1.11.1
    wheel: 0.36.0
    wheel-filename: 1.3.0

[Eric Pettersen]

Reported by Matt Jaremko

[Tom Goddard]

I think this is a misunderstanding, not a bug. The volume morph command just linearly interpolates the two maps at corresponding grid points. If you move one map relative to the other for example to align them that makes no difference, it ignores the alignment. It would certainly be nicer if the command resampled the second map on the grid of the first map if the grids are not aligned, but it currently does not do that. Unfortunately it does not warn you that the grids are not aligned. So you have to use "volume resample" beforehand to interpolate the second map on the grid of the first map. This is documented in the volume morph documentation

​https://www.cgl.ucsf.edu/chimerax/docs/user/commands/volume.html#morph

I have made the volume morph command check that the grids are aligned and issue an error if they are not aligned to avoid this confusion.

I will make a separate ticket to improve the volume morph command so it can do the interpolation of the second map when the grids are not aligned.

[jaremko@…]

Hi Tom,

I did end up using the volume resample. Thanks for the reply.

Best,
Matt

[goddard@…]

Glad you got it to work.  I appreciate the bug report, the better error reporting I added will help others, and hopefully I will find time to make volume morph do the interpolation itself in the future.  I think your other  request for a volume morph slider is more important and long overdue.