Files
Jungfraujoch/common/ROIAzimuthal.cpp
leonarski_f 75e401f0e5
Build Packages / Unit tests (push) Successful in 1h31m59s
Build Packages / build:rpm (rocky8_nocuda) (push) Successful in 8m43s
Build Packages / build:rpm (rocky9_nocuda) (push) Successful in 10m5s
Build Packages / build:rpm (ubuntu2204_nocuda) (push) Successful in 9m27s
Build Packages / build:rpm (ubuntu2404_nocuda) (push) Successful in 8m56s
Build Packages / build:rpm (rocky8_sls9) (push) Successful in 9m24s
Build Packages / build:rpm (rocky9_sls9) (push) Successful in 10m27s
Build Packages / build:rpm (rocky8) (push) Successful in 9m20s
Build Packages / build:rpm (rocky9) (push) Successful in 10m50s
Build Packages / build:rpm (ubuntu2204) (push) Successful in 9m54s
Build Packages / build:rpm (ubuntu2404) (push) Successful in 8m38s
Build Packages / DIALS test (push) Successful in 12m13s
Build Packages / XDS test (durin plugin) (push) Successful in 7m8s
Build Packages / XDS test (JFJoch plugin) (push) Successful in 7m8s
Build Packages / XDS test (neggia plugin) (push) Successful in 7m50s
Build Packages / Generate python client (push) Successful in 16s
Build Packages / Build documentation (push) Successful in 50s
Build Packages / Create release (push) Skipped
v1.0.0-rc.153 (#63)
This is an UNSTABLE release. It includes many experimental features, as well as many AI generated fixes. We recommend using rc.152 for production use.

* jfjoch_broker: Add EXPERIMENTAL pixelrefine mode for image processing
* jfjoch_broker: Allow to load user mask from 8-bit and 16-bit TIFF files
* jfjoch_broker: Add ROI calculation in non-FPGA workflow
* jfjoch_broker: Fixes to TCP image pusher
* jfjoch_broker: Remove NUMA bindings
* jfjoch_broker: Improvements to indexing
* jfjoch_broker: For PSI EIGER, trimming energies are taken from the detector configuration (now compulsory) instead of hardcoded values
* jfjoch_writer: Save ROI definitions and the per-pixel ROI bitmap in the master file; azimuthal ROIs support phi (angular) sectors
* jfjoch_viewer: Major redesign with dockable panels and saved layouts, plus on-canvas creation/move/resize of box, circle and azimuthal ROIs
* jfjoch_viewer: Run jfjoch_process reprocessing jobs from inside the GUI and overlay per-run results

Reviewed-on: #63
2026-06-23 20:29:49 +02:00

81 lines
2.3 KiB
C++

// SPDX-FileCopyrightText: 2024 Filip Leonarski, Paul Scherrer Institute <filip.leonarski@psi.ch>
// SPDX-License-Identifier: GPL-3.0-only
#include <cmath>
#include "JFJochMath.h"
#include "ROIAzimuthal.h"
#include "JFJochException.h"
static float NormalizePhi_deg(float phi) {
phi = std::fmod(phi, 360.0f);
if (phi < 0)
phi += 360.0f;
return phi;
}
ROIAzimuthal::ROIAzimuthal(const std::string &in_name, float in_d_min_A, float in_d_max_A,
float in_phi_min_deg, float in_phi_max_deg)
: ROIElement(in_name),
phi_min_deg(NormalizePhi_deg(in_phi_min_deg)),
phi_max_deg(NormalizePhi_deg(in_phi_max_deg)) {
if ((in_d_min_A <= 0) || (in_d_max_A <= 0))
throw JFJochException(JFJochExceptionCategory::InputParameterInvalid,
"Resolution cannot be zero or negative");
if (in_d_min_A > in_d_max_A) {
d_max_A = in_d_min_A;
d_min_A = in_d_max_A;
} else {
d_max_A = in_d_max_A;
d_min_A = in_d_min_A;
}
}
float ROIAzimuthal::GetDMin_A() const {
return d_min_A;
}
float ROIAzimuthal::GetDMax_A() const {
return d_max_A;
}
bool ROIAzimuthal::HasPhi() const {
return phi_min_deg != phi_max_deg;
}
float ROIAzimuthal::GetPhiMin_deg() const {
return phi_min_deg;
}
float ROIAzimuthal::GetPhiMax_deg() const {
return phi_max_deg;
}
bool ROIAzimuthal::CheckROI(int64_t x, int64_t y, float resolution, float phi_deg) const {
if (resolution < d_min_A || resolution > d_max_A)
return false;
if (phi_min_deg == phi_max_deg)
return true; // full ring: all angles
if (phi_min_deg <= phi_max_deg)
return (phi_deg >= phi_min_deg) && (phi_deg <= phi_max_deg);
return (phi_deg >= phi_min_deg) || (phi_deg <= phi_max_deg); // sector wraps across 0
}
float ROIAzimuthal::GetQMax_recipA() const {
return 2.0f * PI / d_min_A;
}
float ROIAzimuthal::GetQMin_recipA() const {
return 2.0f * PI / d_max_A;
}
ROIConfig ROIAzimuthal::ExportMetadata() const {
double qmin = GetQMin_recipA();
double qmax = GetQMax_recipA();
return ROIConfig{
.type = ROIConfig::ROIType::Azim,
.name = name,
.azim = ROIConfigAzim{.qmin = qmin, .qmax = qmax,
.phi_min = phi_min_deg, .phi_max = phi_max_deg}
};
}