Density sampling and EM projections not working as expected #915
Comments
|
Even though this issue is mostly solved with recent changes in IMP::em, I'm still finding errors in simulated maps. They still have a slight offset from the structure they're simulated from. |
|
What do you mean by offset? Resolution-wise or coordinate-wise? |
|
If you mean coordinate-wise, it would be simple to write a (failing) test that simulated a map from a single atom, then compared the map centroid to the atom coordinates, no? |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
The current version of the EM and EM2D modules uses the KernelParameters.h file to sample the density at a given resolution and angstrom per pixel (apix).
However both uses the RadiusDependentKernelParameters (RDKP) to define the volume and sigma at which the density will be sampled.
At low resolution, the models looked okay given the resolution. However, upon using high resolution data, the simulated maps (respectively the projection) are not correct. This is easily shown by using standard EM software such as EMAN2.
Disabling the RadiusDependentKernelParameters seems to solve one of the issue at high resolution (proper density location but wrong voxel value even when using NumberOfElectron rather than mass when calculating the KP) but low resolution is still not working as expected.
Ubiquitin looks like a ball at 15A losing all the features (EMAN2 actually give proper features even at 15A resolution).
Obviously, disabling the RDKP is not ideal. I am not sure why it was implemented but I am sure there is a good reason for it to be there.
TL/TR: Incorrect density sampling (resolution and voxel intensity) from PDB when using RDKP.
The text was updated successfully, but these errors were encountered: