mx/Jungfraujoch
0
0
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases 37 Wiki Activity

jfjoch_process: Clean previous image information (to be optimized) + option to run rotation indexing

Browse Source
This commit is contained in:
leonarski_f
2026-01-19 12:09:40 +01:00
parent af2ec340db
commit 05da6bc2dd
1 changed files with 15 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
tools/jfjoch_process.cpp
View File

@@ -44,6 +44,7 @@ int main(int argc, char **argv) {
int start_image = 0;
int end_image = -1; // -1 indicates process until end
bool verbose = false;
bool rotation_index = false;
if (argc == 1) {
print_usage(logger);
@@ -51,7 +52,7 @@ int main(int argc, char **argv) {
}
int opt;
while ((opt = getopt(argc, argv, "o:N:s:e:v")) != -1) {
while ((opt = getopt(argc, argv, "o:N:s:e:vR")) != -1) {
switch (opt) {
case 'o':
output_prefix = optarg;
@@ -68,6 +69,9 @@ int main(int argc, char **argv) {
case 'v':
verbose = true;
break;
case 'R':
rotation_index = true;
break;
default:
print_usage(logger);
exit(EXIT_FAILURE);
@@ -110,7 +114,7 @@ int main(int argc, char **argv) {
// Configure Indexing
IndexingSettings indexing_settings;
indexing_settings.Algorithm(IndexingAlgorithmEnum::FFT);
indexing_settings.RotationIndexing(true);
indexing_settings.RotationIndexing(rotation_index);
indexing_settings.GeomRefinementAlgorithm(GeomRefinementAlgorithmEnum::BeamCenter);
experiment.ImportIndexingSettings(indexing_settings);
@@ -207,7 +211,15 @@ int main(int argc, char **argv) {
if (!img) continue;
DataMessage msg = img->ImageData();
DataMessage msg;
msg.image = img->ImageData().image;
msg.number = img->ImageData().number;
msg.error_pixel_count = img->ImageData().error_pixel_count;
msg.image_collection_efficiency = img->ImageData().image_collection_efficiency;
msg.storage_cell = img->ImageData().storage_cell;
msg.user_data = img->ImageData().user_data;
total_uncompressed_bytes += msg.image.GetUncompressedSize();
auto image_start_time = std::chrono::high_resolution_clock::now();