transform_deconvolve() [2/2]

template<OperationReadyData DataType>

DataType MayaFlux::Yantra::transform_deconvolve	(	DataType &	input,
		const std::vector< double > &	impulse_to_remove,
		double	regularization = 1e-6
	)

Deconvolution using frequency domain division (IN-PLACE) Useful for removing known impulse responses.

Template Parameters

DataType

OperationReadyData type

Parameters

input	Input data - WILL BE MODIFIED
impulse_to_remove	Impulse response to remove
regularization	Regularization factor for numerical stability

Returns

Deconvolved data

Definition at line 444 of file ConvolutionHelper.hpp.

{
    auto [target_data, structure_info] = OperationHelper::extract_structured_double(input);
 
    auto deconvolution_op = [regularization](const Eigen::VectorXcd& input_fft,
                                const Eigen::VectorXcd& kernel_fft,
                                Eigen::VectorXcd& result_fft) {
        std::ranges::transform(input_fft, kernel_fft, result_fft.begin(),
            [regularization](const std::complex<double>& signal, const std::complex<double>& kernel) {
                double magnitude_sq = std::norm(kernel);
                if (magnitude_sq < regularization) {
                    return std::complex<double>(0.0, 0.0);
                }
                return signal * std::conj(kernel) / (magnitude_sq + regularization);
            });
    };
 
    for (auto& span : target_data) {
        auto result = fft_convolve_helper(span, impulse_to_remove, deconvolution_op);
        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().

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

Here is the call graph for this function:

Here is the caller graph for this function: