This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Last revision Both sides next revision | ||
pipeline:utilities:sxpipe_moon_eliminator [2018/05/14 16:19] 127.0.0.1 external edit |
pipeline:utilities:sxpipe_moon_eliminator [2019/04/02 09:30] lusnig |
||
---|---|---|---|
Line 1: | Line 1: | ||
~~NOTOC~~ | ~~NOTOC~~ | ||
- | ===== sxpipe | + | ===== sp_pipe |
- | Moon eliminator: Eliminate moons or remove dusts from the background of a 3D density map based on the expected molecular mass. Optionally, create 3D mask also. | + | Volume Adjustment: Eliminate moons or remove dusts from the background of a 3D density map based on the expected molecular mass. |
\\ | \\ | ||
Line 9: | Line 9: | ||
Usage in command line | Usage in command line | ||
- | | + | |
+ | |||
\\ | \\ | ||
===== Typical usage ===== | ===== Typical usage ===== | ||
- | sxpipe | + | sp_pipe |
- | Create reference 3D volume and 3D mask from R-VIPER 3D model with the resample ratio used in ISAC2 using an expected molecular mass [kDa]. | + | There are two modes to run the program: |
- | sxpipe.py moon_eliminator ' | + | \\ |
+ | __1. Single Volume Mode:__ | ||
- | Create reference 3D volume | + | Create reference 3D structure |
- | | + | |
- | Create reference 3D volume | + | Create reference 3D structure |
- | | + | |
+ | |||
+ | Create reference 3D structure and 3D mask from post-refined MERIDIEN 3D model using the expected molecular mass [kDa]. | ||
+ | |||
+ | sp_pipe.py moon_eliminator ' | ||
+ | |||
+ | \\ | ||
+ | __2. Halfset Volumes Mode:__ | ||
+ | |||
+ | Create reference 3D structure and 3D mask from halfset unfiltered maps produced by MERIDIEN, using the expected molecular mass [kDa]. | ||
+ | |||
+ | sp_pipe.py moon_eliminator ' | ||
Line 33: | Line 45: | ||
=== Main Parameters === | === Main Parameters === | ||
; input_volume_path : Input volume path: Path to input volume file containing the 3D density map. (default required string) | ; input_volume_path : Input volume path: Path to input volume file containing the 3D density map. (default required string) | ||
- | ; output_directory : Output directory: The results will be written here. This directory will be created automatically and it must not exist previously. (default required string) | + | ; output_directory : Output directory: The results will be written here. It cannot |
- | | + | ; %%--%%pixel_size : Output pixel size [A]: The original pixel size of dataset. This must be the pixel size after resampling when resample_ratio != 1.0. That is, it will be the pixel size of the output |
- | | + | ; %%--%%use_mol_mass |
- | ; %%--%%moon_distance | + | ; %%--%%mol_mass |
- | ; %%--%%falloff_speed | + | ; %%--%%use_density_threshold : Use ad-hoc density threshold: Use user-provided ad-hoc density threshold, instead of computing the value from the molecular mass. Below this density value, the data is assumed not to belong to the main body of the particle density. (default none) : %%--%%use_mol_mass==False |
- | ; %%--%%use_density_threshold : Use ad-hoc density threshold: Use user-provided ad-hoc density threshold, instead of computing the value from the molecular mass. Below this density value, the data is assumed not to belong to the main body of the particle density. (default none) | + | ; %%--%%moon_distance : Distance to the nearest moon [Pixels]: The moons further than this distance from the density surface will be elminated. The value smaller than the default is not recommended because it is difficult to avoid the stair-like gray level change at the edge of the density surface. (default 3.0) |
- | ; %%--%%resample_ratio : Resample ratio: Specify a value larger than 0.0. By default, the program does not resmaple | + | ; %%--%%resample_ratio : Resample ratio: Specify a value larger than 0.0. By default, the program does not resample |
; %%--%%box_size : Output box size [Pixels]: The x, y, and z dimensions of cubic area to be windowed from input 3D volume for output 3D volumes. This must be the box size after resampling when resample_ratio != 1.0. (default none) | ; %%--%%box_size : Output box size [Pixels]: The x, y, and z dimensions of cubic area to be windowed from input 3D volume for output 3D volumes. This must be the box size after resampling when resample_ratio != 1.0. (default none) | ||
- | | + | ; %%--%%invert_handedness : Invert handedness: Invert the handedness of the 3D map. (default False) |
- | ; %%--%%shift3d_x : 3D x-shift [Pixels]: Provide 3D x-shift corresponding to shifting the 3D volume along x-axis. (default 0) | + | |
- | ; %%--%%shift3d_y : 3D y-shift [Pixels]: Provide 3D y-shift corresponding to shifting the 3D volume along y-axis. (default 0) | + | |
- | ; %%--%%shift3d_z : 3D z-shift [Pixels]: Provide 3D z-shift corresponding to shifting the 3D volume along z-axis. (default 0) | + | |
- | | + | |
; %%--%%fl : Low-pass filter resolution [A]: >0.0: low-pass filter to the value in Angstrom; =-1.0: no low-pass filter. The program applies this low-pass filter before the moon elimination. (default -1.0) | ; %%--%%fl : Low-pass filter resolution [A]: >0.0: low-pass filter to the value in Angstrom; =-1.0: no low-pass filter. The program applies this low-pass filter before the moon elimination. (default -1.0) | ||
- | ; %%--%%aa : Low-pass filter fall-off [1/Pixels]: Low-pass filter fall-off in absolute frequency. The program applies this low-pass filter before the moon elimination. Effective only when %%--%%fl > 0.0. (default 0.1) | ||
- | ; %%--%%generate_mask : Generate soft-edged mask: Generate a soft-edged mask with specified method from the 3D binary volume corresponding to the specified molecular mass or density threshold, using the values provided through %%--%%gm_edge_width and/or %%--%%gm_falloff_speed. Available methods are (1) ' | ||
\\ | \\ | ||
=== Advanced Parameters === | === Advanced Parameters === | ||
- | ; %%--%%gm_edge_width | + | |
- | ; %%--%%gm_falloff_speed | + | ; %%--%%ndilation : Dilation width [Pixels]: The pixel width to dilate the 3D binary volume corresponding to the specified molecular mass or density threshold prior to softening the edge. By default, it is set to half of --moon_distance so that the voxels with 1.0 values in the mask are same as the hard-edged molecular-mass binary volume. (default -1.0) |
+ | | ||
+ | ; %%--%%edge_type | ||
; %%--%%outputs_root : Root name of outputs: Specify the root name of all outputs. It cannot be empty string or only white spaces. (default vol3d) | ; %%--%%outputs_root : Root name of outputs: Specify the root name of all outputs. It cannot be empty string or only white spaces. (default vol3d) | ||
+ | ; %%--%%resampled_shift3d : Providing resampled 3D shifts: Use this option when you are providing the resampled 3D shifts (using pixel size of outputs) when %%--%%resample_ratio!=1.0. By default, the program assums the provided shifts are not resampled. (default False) | ||
+ | ; %%--%%shift3d_x : 3D x-shift [Pixels]: 3D x-shift value. (default 0) | ||
+ | ; %%--%%shift3d_y : 3D y-shift [Pixels]: 3D y-shift value. (default 0) | ||
+ | ; %%--%%shift3d_z : 3D z-shift [Pixels]: 3D z-shift value. (default 0) | ||
+ | ; %%--%%aa : Low-pass filter fall-off [1/Pixels]: Low-pass filter fall-off in absolute frequency. The program applies this low-pass filter before the moon elimination. Effective only when %%--%%fl > 0.0. (default 0.1) : %%--%%fl!=-1.0 | ||
; %%--%%debug : Run with debug mode: Mainly for developer. (default False) | ; %%--%%debug : Run with debug mode: Mainly for developer. (default False) | ||
Line 62: | Line 75: | ||
=== List of output Files === | === List of output Files === | ||
|| **File Name** || **Discription** || | || **File Name** || **Discription** || | ||
+ | || *_ref_before_moon_elimination.hdf || File containing the 3D reference map before moon elimination (i.e., the 3D map just before applying moon elimination.). || | ||
|| *_ref_moon_eliminated.hdf || File containing the moon eliminated 3D reference. || | || *_ref_moon_eliminated.hdf || File containing the moon eliminated 3D reference. || | ||
|| *_mask_moon_elminator.hdf || File containing the moon elminator 3D mask. || | || *_mask_moon_elminator.hdf || File containing the moon elminator 3D mask. || | ||
Line 71: | Line 85: | ||
This command executes the following processes: | This command executes the following processes: | ||
- Extract resample ratio from ISAC run directory if necessary (mainly for R-VIPER models). | - Extract resample ratio from ISAC run directory if necessary (mainly for R-VIPER models). | ||
- | - Resample and window the volume | + | - Resample and window the map if necessary (mainly for R-VIPER models) |
- | - Shift 3D volume | + | - Shift 3D map if necessary. |
- Invert the handedness if necessary. | - Invert the handedness if necessary. | ||
- | - Apply low-pass filter to the input volume | + | - Apply low-pass filter to the input map before the moon elimination if necessary. |
- | - Create reference 3D volume | + | - Save reference 3D map before eliminating the moons. |
- | - Create 3D mask from the 3D bainary | + | - Create reference 3D map by eliminating the moons from the input map and save the results. |
+ | - Create 3D mask from the 3D binary | ||
\\ | \\ | ||
Line 86: | Line 101: | ||
\\ | \\ | ||
==== Developer Notes ==== | ==== Developer Notes ==== | ||
- | === 2018/04/12 Toshio Moriya === | + | === 2018/06/18 Toshio Moriya === |
- | Tips about balancing settings of moon_distance and falloff_speed | + | |
+ | Wish | ||
+ | * Add options for 3D rotation of the map. | ||
+ | |||
+ | === 2018/06/18 Toshio Moriya === | ||
+ | |||
+ | Tips about balancing settings of moon_distance, dilation, | ||
* moon_distance | * moon_distance | ||
- | * A shorter | + | * In principle, |
- | * If it is too long, the moons will be connected and creates strange low density shape at the edge of the monn distance. | + | * If moon_distance |
* On the other hand, if it is too short, soft-edge will have the stair-like gray level change because of quantization or digitization. | * On the other hand, if it is too short, soft-edge will have the stair-like gray level change because of quantization or digitization. | ||
- | * falloff_speed | + | * ndilation |
- | * The falloff speed seems to have be smaller (or slower) than moon_distance, so that Gaussian sigma will be at least larger than 1[pixel]. | + | * Setting dilation |
- | * If it is larger (or faster), | + | * Setting dilation to smaller than half of moon_distance |
- | * In addition, moon-eliminated reference 3D volume will have a strange dent near zero. | + | * edge_sigma |
- | * Tested combinations with a density | + | * In principle, smaller edge_sigma is better. |
- | * moon_distance=3 & falloff_speed=3.0 | + | * However, |
- | * Very nice balance! | + | * If not, the density distribution of moon eliminator |
- | * moon_distance=3 & falloff_speed=5.0 (gauss_sigma = 0.6[Pixels]) | + | * In addition, moon-eliminated reference 3D map will have a strange dent near zero. |
- | * NG. The density distribution of moon eliminator 3D mask is not smooth. | + | |
- | * The moon-eliminated reference 3D volume has a strange dent near zero. | + | |
- | * That is, falloff_speed must be too fast for this moon_distance! | + | |
- | * moon_distance=6 & falloff_speed=3.0 (gauss_sigma = 2.0[Pixels]) | + | |
- | * OK but created strange low density shape at the edge of the monn distance. | + | |
- | * moon_distance=6 & falloff_speed=5.0 (gauss_sigma = 1.2[Pixels]) | + | |
- | * OK but created strange low density shape at the edge of the monn distance. | + | |
- | + | ||
- | Tips about balancing settings of gm_edge_width and gm_falloff_speed for Gaussian soft-edge. | + | |
- | * gm_edge_width | + | |
- | * A wider edge seems to be better. If it is too narrow, soft-edge will have the stair-like gray level change because of quantization or digitization. | + | |
- | * However, | + | |
- | * gm_falloff_speed | + | |
- | * The falloff speed seems to have be smaller than the edge width, so that Gaussian sigma will be at least larger than 1[pixel]. | + | |
- | * If it is larger, the density distribution of moon-eliminated | + | |
\\ | \\ | ||
Line 126: | Line 132: | ||
\\ | \\ | ||
==== Files ==== | ==== Files ==== | ||
- | sparx/bin/sxpipe.py | + | sparx/bin/sp_pipe.py |
\\ | \\ | ||
==== See also ==== | ==== See also ==== | ||
- | [[pipeline: | + | [[pipeline: |
\\ | \\ |