Skip to content
Snippets Groups Projects
Commit 6859c838 authored by Huabin Zhou's avatar Huabin Zhou
Browse files

Add 'rlnTomoName' in the written starfile for the support of relion5

parent 99d9b416
No related merge requests found
...@@ -162,6 +162,7 @@ class ClashResolver: ...@@ -162,6 +162,7 @@ class ClashResolver:
def save_results(self, vol_array): def save_results(self, vol_array):
prefix = config.prefix prefix = config.prefix
tomoID = config.tomoID
# save xml file # save xml file
results_xml = save_xml_one_loop( results_xml = save_xml_one_loop(
self.ccc_save, self.rot_save, self.coord_save, self.model_save self.ccc_save, self.rot_save, self.coord_save, self.model_save
...@@ -181,6 +182,7 @@ class ClashResolver: ...@@ -181,6 +182,7 @@ class ClashResolver:
df.to_csv(self.output_path + prefix + ".csv", index=False) df.to_csv(self.output_path + prefix + ".csv", index=False)
# save star file # save star file
star = coord_to_relion(self.coord_save, self.rot_save, self.ccc_save) star = coord_to_relion(self.coord_save, self.rot_save, self.ccc_save)
star['rlnTomoName'] = tomoID
starfile.write(star, self.output_path + prefix + ".star", overwrite=True) starfile.write(star, self.output_path + prefix + ".star", overwrite=True)
# plot the annotation over the tomogram # plot the annotation over the tomogram
df2 = coord_to_relion(self.coord_save, self.rot_save) df2 = coord_to_relion(self.coord_save, self.rot_save)
......
...@@ -19,7 +19,7 @@ tomoID = split_path[-2] ...@@ -19,7 +19,7 @@ tomoID = split_path[-2]
prefix = tomoID + ".test1" prefix = tomoID + ".test1"
input_folder = "inputs/" input_folder = "inputs/"
tomogram = input_folder + tomoID + ".test1.lps25.mrc" tomogram = input_folder + tomoID + ".test1.lps25.mrc"
black_on_white = 1 # The constrast of tomogram, 0 if the particle has white contrast (positive values) black_on_white = 1 # The contrast of tomogram, 0 if the particle has white contrast (positive values)
templates = [input_folder + "mono-8Apx-lps30-box30-rot12-core.mrc"] templates = [input_folder + "mono-8Apx-lps30-box30-rot12-core.mrc"]
contour_level = [0.2] contour_level = [0.2]
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment