transform_spectral_filter() [1/2]

template<OperationReadyData DataType>

DataType MayaFlux::Yantra::transform_spectral_filter	(	DataType &	input,
		double	low_freq,
		double	high_freq,
		double	sample_rate,
		uint32_t	window_size,
		uint32_t	hop_size,
		std::vector< std::vector< double > > &	working_buffer
	)

Spectral filtering using existing FFT infrastructure with C++20 ranges (OUT-OF-PLACE)

Template Parameters

DataType

OperationReadyData type

Parameters

input	Input data - will NOT be modified
low_freq	Low cutoff frequency (Hz)
high_freq	High cutoff frequency (Hz)
sample_rate	Sample rate (Hz)
window_size	FFT window size
hop_size	Hop size for overlap-add
working_buffer	Buffer for operations (will be resized if needed)

Returns

Filtered data

Definition at line 227 of file SpectralHelper.hpp.

{
    auto [target_data, structure_info] = OperationHelper::setup_operation_buffer(input, working_buffer);
 
    auto processor = [low_freq, high_freq, sample_rate](Eigen::VectorXcd& spectrum, size_t) {
        auto bin_indices = std::views::iota(0, static_cast<int>(spectrum.size()));
        std::ranges::for_each(bin_indices, [&](int bin) {
            double freq = (bin * sample_rate) / (2.0 * (double)spectrum.size());
            if (freq < low_freq || freq > high_freq) {
                spectrum[bin] = 0.0;
            }
        });
    };
 
    working_buffer.resize(target_data.size());
    for (size_t i = 0; i < target_data.size(); ++i) {
        working_buffer[i] = process_spectral_windows(target_data[i], window_size, hop_size, processor);
    }
 
    return OperationHelper::reconstruct_from_double<DataType>(working_buffer, structure_info);
}

References process_spectral_windows(), and MayaFlux::Yantra::OperationHelper::setup_operation_buffer().

Here is the call graph for this function: