Skip to Content
Last modified by flenners on 2026-02-04 17:25
From version 16.2
edited by flenners
on 2026-02-04 17:15
Change comment: There is no comment for this version
To version 6.1
edited by flenners
on 2022-09-25 10:51
Change comment: There is no comment for this version

Summary

Details

Page properties
Content
... ... @@ -1,54 +1,6 @@
1 -== {{id name="03aReconstructionwithRecoGUI-Preperation:"/}}Preperation: ==
1 += {{id name="03aReconstructionwithRecoGUI-LoadingData"/}}Loading Data
2 +\\ =
2 2  
3 -Terminal:
4 -
5 -
6 -salloc ~-~-partition=psxcpu ~-~-nodes=1 ~-~-time=06:00:00
7 -\\(if you need gpu: ~-~-partition=psxgpu )
8 -
9 -//Answer~:// salloc: job [//number//] queued and waiting for resources
10 -
11 - salloc: job [number] has been allocated resources
12 -
13 - salloc: Granted job allocation [number]
14 -
15 - salloc: Waiting for resource configuration
16 -
17 - salloc: Nodes //max-wn035// are ready for job  (name of your node; code after “max-“ can vary, in this example wn035, if you do not get a specific node name, repeat salloc command)
18 -
19 -(% class="code" %)
20 -(((
21 -ssh max-wn035 (shell to your node)
22 -\\module load maxwell mamba
23 -)))
24 -
25 -(% class="code" %)
26 -(((
27 -. mamba-init (. is important! never do without, this can crash your FastX)
28 -\\mamba activate /asap3/petra3/gpfs/common/p05/nano/envs/mamba
29 -)))
30 -
31 -
32 -spyder &
33 -\\~-~-~-~-~-~-
34 -\\\\Spyder will open.
35 -
36 -On the right, choose "file explorer". Navigate to your beamtime folder and find processed/scripts/RecoGUI:
37 -
38 -[[image:attach:image2022-11-1_13-7-14.png||height="250"]]
39 -
40 -
41 -Important: you need to be in the right folder, else the script will not be able to find a file in the first cell of the script
42 -
43 -
44 -Open RecoGUI.py by double clicking the file and click "Run" (green arrow)
45 -
46 -
47 -
48 -= {{id name="03aReconstructionwithRecoGUI-LoadingData"/}}Loading Data =
49 -
50 -= =
51 -
52 52  **TXM DATA**
53 53  
54 54  Go to "Load' tab and select year, enter your beamtime ID (e.g. 11001234) and press "enter" on keyboard.
... ... @@ -63,18 +63,16 @@
63 63  
64 64  Continue with next step, "Preparation".
65 65  
18 +\\
66 66  
67 67  **HOLOTOMO DATA**
68 68  
69 -Hint: you have to perform the phase reconstruction of your projections first! This is usually done by beamline staff. Then you can continue:
22 +Hint: you have to perform the phase reconstruction of your projections first!
70 70  
71 -Go to "Load" tab, load processed data.
24 +Go to "Load" tab. Click on "Load normalized data" at the bottom of the page and select your phase reconstruced projections.
72 72  
73 -Check the subfolder path of the phase retievals. Standard path is holopipe/phase_retrieval
26 +The common path is: ///asap3/petra3/gpfs/p05/YEAR/data/BEAMTIMEID/processed/SCANNAME/reco_~#~#//
74 74  
75 -Click on "Load normalized data" at the bottom of the page and select your phase reconstruced projections.
76 -
77 -
78 78  You do not need to load any additional dark images.
79 79  
80 80  Skip "Normalize, Minus Log, Rotate" in the Preperation Tab. The optional steps are still optional. (see below)
... ... @@ -81,10 +81,11 @@
81 81  
82 82  In principle, you can directly proceed to the third tab "Reconstruction". 
83 83  
34 +\\
84 84  
85 -= {{id name="03aReconstructionwithRecoGUI-Preparationofdataforreconstruction"/}}Preparation of data for reconstruction =
36 += {{id name="03aReconstructionwithRecoGUI-Preparationofdataforreconstruction"/}}Preparation of data for reconstruction =
86 86  
87 -(**Can be skipped for Holotomography**)
38 +(Can be skipped for Holotomography)
88 88  
89 89  Go to "Prep" tab.
90 90  
... ... @@ -96,10 +96,11 @@
96 96  
97 97  Click "Rotate" button.
98 98  
50 +\\
99 99  
100 100  The following steps are **optional:**
101 101  
102 -**Cropping:** If your samples are significantly (% style="color:#003366" %)__smaller__(%%) than the field of view, you can crop the data to reduce computation time and storage. Draw a rectangle around your sample by keeping the left mouse button pressed. Check that your sample stays inside that rectangle at all angles. If you are happy, press the "crop" button to crop your data.
54 +**Cropping:** If your samples are significantly (% style="color: rgb(0,51,102);" %)__smaller__(%%) than the field of view, you can crop the data to reduce computation time and storage. Draw a rectangle around your sample by keeping the left mouse button pressed. Check that your sample stays inside that rectangle at all angles. If you are happy, press the "crop" button to crop your data.
103 103  
104 104  **Binning:** You can bin your data before reconstruction. Since the detector has a point spread function of 2-3 pixels, it is save to bin by a factor of 2 without loosing spatial resolution. This reduces computation time and storage space needed.
105 105  
... ... @@ -107,103 +107,39 @@
107 107  
108 108  **Filter Projections before reconstruction:** You can filter the data before reconstruction. This reduces the noise, but also can induce blurring in the data. Only recommended for very noisy data. Filtering after reconstruction is usually better.
109 109  
62 +\\
110 110  
111 111  **Linogram alignment:** Linogram alignment if sample moved. Ask your local contact if this is necessary.
112 112  
66 +\\
113 113  
114 114  **Save current stack.** Writes the normalized projections to your beamtime folder.
115 115  
70 +\\
116 116  
117 -= {{id name="03aReconstructionwithRecoGUI-Reconstruction"/}}Reconstruction =
72 += {{id name="03aReconstructionwithRecoGUI-Reconstruction"/}}Reconstruction =
118 118  
74 +\\
119 119  
120 120  **Finding the correct center of rotation. **
121 121  
122 -First, you have to find the correct center of rotation. For this, a single slice is reconstructed with different centers of rotation
78 +First, you have to find the correct center of rotation. For this, a single slice is reconstructed with different centers of rotation.
123 123  
124 -**~1. Rot center: **Defines the rotation center in pixel coordinates. For the start, enter half of your image size (1024 for unbinned data, 512 for data with binning 2.)
80 +**Rot center:** For the start, enter half of your image size (1024 for unbinned data, 512 for data with binning 2.)
125 125  
126 -**2. Delta: **Range of pixels around the in step 1 given rot center for which the test slice is reconstructed. Good starting value is 50. For fine rot center, choose 10. 
82 +**Delta:** Range which is reconstructed. Good starting value is 50. For fine rot center, choose 10. 
127 127  
128 -**3. Stepsize:** Stepsize beween the different rotation centers. Good start value is 5 and 1 for fine rot center. (for example [ ... , 507, 512, 517, ... ] for a stepsize of 5)
84 +**Stepsize:** Stepsize beween the different rotation centers. Good start value is 5 and 1 for fine rot center.
129 129  
130 -**4. Slice:** Defines the pixel row which is used to reconstruct the test slice. Tip: Choose a region where you expect to see distinctive structures. (You can also check different slices for checking the rotation center.)
86 +Slice: Slice which will be reconstructed. You can also check different slices for checking the rotation center.
131 131  
132 -**5. Check rotation center: **Click on the "Check rotation center" button.** **In the preview window move the slider around until you find the positions with the least artifacts. Remember the position number and check the Spyder Console to get the new rot center coordinates in pixel.
88 +\\
133 133  
134 -[[image:attach:image2022-11-14_14-4-49.png||thumbnail="true" height="250"]]
135 -
136 -in blue the coordinate of the center of rotation; in yellow the position number of the preview slider
137 -
138 -**6. Update and Repeat:** Replace the Rot Center (step 1) with the new found coordinate, lower the Delta (step 2) and Stepsize (step 3) and repeat the process until you are satisfied with the result.
139 -
140 -
141 -How does a good rotation center look like?
142 -
143 -(% class="relative-table wrapped" style="width:37.5112%" %)
144 -|=(((
145 -good
146 -)))|=(((
147 -bad
148 -)))|=(((
149 -comment
150 -)))
151 -|(((
152 -(% class="content-wrapper" %)
153 -(((
154 -[[image:attach:image2022-11-14_14-17-45.png||height="250"]]
155 -)))
156 -)))|(((
157 -(% class="content-wrapper" %)
158 -(((
159 -[[image:attach:image2022-11-14_14-19-33.png||height="250"]]
160 -)))
161 -)))|(((
162 -
163 -)))
164 -|(((
165 -(% class="content-wrapper" %)
166 -(((
167 -[[image:attach:image2022-11-14_15-16-34.png||height="250"]]
168 -)))
169 -)))|(((
170 -(% class="content-wrapper" %)
171 -(((
172 -[[image:attach:image2022-11-14_15-7-13.png||height="250"]]
173 -)))
174 -)))|(((
175 -The ghosting effect on the edges is due to sample movement during the scan.
176 -)))
177 -|(((
178 -
179 -)))|(((
180 -
181 -)))|(((
182 -
183 -)))
184 -|(((
185 -(% class="content-wrapper" %)
186 -(((
187 -[[image:attach:image2022-11-14_15-29-47.png||height="250"]]
188 -)))
189 -)))|(((
190 -(% class="content-wrapper" %)
191 -(((
192 -**[[image:attach:image2022-11-14_15-28-50.png||height="250"]]**
193 -)))
194 -)))|(((
195 -
196 -)))
197 -
198 -[[image:attach:image2023-1-23_11-39-57.png||height="400"]]
199 -
200 -
201 -
202 202  **Reconstruction Parameters. **
203 203  
204 204  When you found the correct rotation center, you can continue with the reconstruction.
205 205  
206 -Enter the final roation center from the above step.
94 +Enter the final roation center from the above step.
207 207  
208 208  You can select different reconstruction algorithms and filters. Standard is gridrec and shepp.
209 209  
... ... @@ -213,6 +213,7 @@
213 213  
214 214  before the final reconstruction, you can reconstruct a test slice to check your reconstruction parameters. Enter the slice number you want to reconstruct and press " Reconstruct slice".
215 215  
104 +\\
216 216  
217 217  When you are happy with your result, press "Reconstruct full stack".
218 218  
... ... @@ -220,7 +220,10 @@
220 220  
221 221  When you are finished with one scan, please press the "clear all" button! This saves the log data you see on the right, so you can later check what you did. The data are also removed from the memory.
222 222  
112 +\\
223 223  
114 +\\
224 224  
116 +\\
225 225  
226 -
118 +\\
image2022-11-14_14-1-39.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -170.0 KB
Content
image2022-11-14_14-17-45.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -473.7 KB
Content
image2022-11-14_14-19-33.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -473.7 KB
Content
image2022-11-14_14-4-49.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -172.4 KB
Content
image2022-11-14_15-16-34.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -463.7 KB
Content
image2022-11-14_15-28-50.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -455.4 KB
Content
image2022-11-14_15-29-47.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -428.5 KB
Content
image2022-11-14_15-7-13.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.wirtenss
Size
... ... @@ -1,1 +1,0 @@
1 -431.4 KB
Content
image2022-11-1_13-7-14.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.flenners
Size
... ... @@ -1,1 +1,0 @@
1 -165.5 KB
Content
image2023-1-23_11-39-57.png
Author
... ... @@ -1,1 +1,0 @@
1 -XWiki.flenners
Size
... ... @@ -1,1 +1,0 @@
1 -1.4 MB
Content
Confluence.Code.ConfluencePageClass[0]
Id
... ... @@ -1,1 +1,1 @@
1 -236814957
1 +284691808
URL
... ... @@ -1,1 +1,1 @@
1 -https://confluence.desy.de/spaces/P5I/pages/236814957/03a Reconstruction with Reco GUI
1 +https://confluence.desy.de/spaces/P5I/pages/284691808/03a Reconstruction with Reco GUI