leonarski_f
21a8ea51ee
Build Packages / build:rpm (ubuntu2404_nocuda) (push) Failing after 8m56s
Build Packages / build:rpm (rocky8_nocuda) (push) Failing after 10m5s
Build Packages / build:rpm (rocky9_nocuda) (push) Failing after 11m41s
Build Packages / build:rpm (ubuntu2204_nocuda) (push) Failing after 11m39s
Build Packages / build:rpm (rocky9_sls9) (push) Failing after 11m42s
Build Packages / build:rpm (rocky8_sls9) (push) Failing after 11m45s
Build Packages / build:rpm (rocky8) (push) Failing after 11m47s
Build Packages / build:rpm (rocky9) (push) Failing after 10m21s
Build Packages / Generate python client (push) Successful in 14s
Build Packages / build:rpm (ubuntu2204) (push) Failing after 9m58s
Build Packages / Create release (push) Skipped
Build Packages / XDS test (neggia plugin) (push) Successful in 8m14s
Build Packages / Build documentation (push) Successful in 36s
Build Packages / XDS test (durin plugin) (push) Successful in 8m52s
Build Packages / build:rpm (ubuntu2404) (push) Failing after 9m4s
Build Packages / XDS test (JFJoch plugin) (push) Successful in 9m7s
Build Packages / DIALS test (push) Successful in 11m49s
Build Packages / Unit tests (push) Successful in 1h8m58s
82 lines
2.6 KiB
C++
82 lines
2.6 KiB
C++
// SPDX-FileCopyrightText: 2025 Filip Leonarski, Paul Scherrer Institute <filip.leonarski@psi.ch>
|
|
// SPDX-License-Identifier: GPL-3.0-only
|
|
|
|
#pragma once
|
|
|
|
#include <vector>
|
|
#include <unordered_set>
|
|
#include <map>
|
|
#include <mutex>
|
|
#include <optional>
|
|
|
|
#include "../common/TopPixels.h"
|
|
#include "JFJochReaderDataset.h"
|
|
#include "../common/CrystalLattice.h"
|
|
#include "../common/Histogram.h"
|
|
|
|
constexpr static int32_t GAP_PXL_VALUE = INT32_MIN + 1;
|
|
constexpr static int32_t ERROR_PXL_VALUE = INT32_MIN;
|
|
constexpr static int32_t SATURATED_PXL_VALUE = INT32_MAX;
|
|
|
|
struct JFJochReaderRawImage {
|
|
std::vector<uint8_t> image_buffer;
|
|
CompressedImage image;
|
|
};
|
|
|
|
class JFJochReaderImage {
|
|
std::shared_ptr<const JFJochReaderDataset> dataset;
|
|
|
|
std::vector<int32_t> image; // Image in the reader must be 32-bit signed, uncompressed
|
|
DataMessage message;
|
|
|
|
std::unordered_set<int64_t> saturated_pixel;
|
|
std::unordered_set<int64_t> error_pixel;
|
|
std::vector<std::pair<int32_t, int32_t>> valid_pixel;
|
|
|
|
// Fast stats without storing/sorting all valid pixels
|
|
int32_t valid_min = 0;
|
|
int32_t valid_max = 0;
|
|
size_t valid_count = 0;
|
|
bool has_valid = false;
|
|
|
|
// For overlay: track top pixels with a tiny O(K) structure; export to vector for UI
|
|
TopPixels top_pixels_acc{20};
|
|
std::vector<std::pair<int32_t, int32_t>> top_pixels;
|
|
|
|
// This histogram operates in square root of count from 0 to 10^20
|
|
Histogram count_histogram{100000};
|
|
|
|
constexpr static float auto_foreground_range = 99.0f;
|
|
int32_t auto_foreground;
|
|
void CalcAutoContrast();
|
|
|
|
template <class T>
|
|
void ProcessInputImage(const void* image, size_t npixel, int64_t sat_value, int64_t special_value);
|
|
void ProcessInputImage(const CompressedImage& image);
|
|
public:
|
|
JFJochReaderImage(const DataMessage &msg, const std::shared_ptr<const JFJochReaderDataset> &dataset);
|
|
JFJochReaderImage(const JFJochReaderImage &other);
|
|
|
|
const DataMessage &ImageData() const;
|
|
DataMessage &ImageData();
|
|
|
|
const std::vector<int32_t> &Image() const;
|
|
const std::unordered_set<int64_t> &SaturatedPixels() const;
|
|
const std::unordered_set<int64_t> &ErrorPixels() const;
|
|
const JFJochReaderDataset &Dataset() const;
|
|
|
|
std::optional<std::pair<int32_t, int32_t>> ValidMinMax() const;
|
|
const std::vector<std::pair<int32_t, int32_t>> &GetTopPixels() const;
|
|
|
|
void AddImage(const JFJochReaderImage& other);
|
|
std::vector<float> GetAzInt1D() const;
|
|
std::vector<float> GetAzInt1D_BinToQ() const;
|
|
std::shared_ptr<JFJochReaderDataset> CreateMutableDataset();
|
|
|
|
int32_t GetAutoContrastValue() const;
|
|
|
|
std::vector<float> GetHistogram() const;
|
|
};
|
|
|
|
|