transform_cross_correlate() [2/2]

template<OperationReadyData DataType>

DataType MayaFlux::Yantra::transform_cross_correlate	(	DataType &	input,
		const std::vector< double > &	template_signal,
		bool	normalize = true
	)

Cross-correlation using FFT (convolution with time-reversed impulse) (IN-PLACE)

Template Parameters

DataType

OperationReadyData type

Parameters

input	Input data - WILL BE MODIFIED
template_signal	Template signal for correlation
normalize	Whether to normalize the result

Returns

Cross-correlated data

Definition at line 235 of file ConvolutionHelper.hpp.

{
    auto [target_data, structure_info] = OperationHelper::extract_structured_double(input);
 
    std::vector<double> reversed_template(template_signal.rbegin(), template_signal.rend());
 
    auto correlation_op = [](const Eigen::VectorXcd& input_fft,
                              const Eigen::VectorXcd& kernel_fft,
                              Eigen::VectorXcd& result_fft) {
        std::ranges::transform(input_fft, kernel_fft, result_fft.begin(),
            [](const std::complex<double>& a, const std::complex<double>& b) {
                return a * std::conj(b);
            });
    };
 
    for (auto& span : target_data) {
        auto result = fft_convolve_helper(span, reversed_template, correlation_op);
 
        if (normalize && !result.empty()) {
            auto [min_it, max_it] = std::ranges::minmax_element(result);
            double max_abs = std::max(std::abs(*min_it), std::abs(*max_it));
            if (max_abs > 0.0) {
                std::ranges::transform(result, result.begin(),
                    [max_abs](double val) { return val / max_abs; });
            }
        }
 
        std::copy(result.begin(), result.begin() + std::min(result.size(), span.size()), span.begin());
    }
 
    auto reconstructed_data = target_data
        | std::views::transform([](const auto& span) {
              return std::vector<double>(span.begin(), span.end());
          })
        | std::ranges::to<std::vector>();
 
    return OperationHelper::reconstruct_from_double<DataType>(reconstructed_data, structure_info);
}

References MayaFlux::Yantra::OperationHelper::extract_structured_double(), fft_convolve_helper(), and MayaFlux::normalize().

Referenced by MayaFlux::Yantra::ConvolutionTransformer< InputType, OutputType >::transform_implementation(), transform_matched_filter(), and transform_matched_filter().

Here is the call graph for this function:

Here is the caller graph for this function: