VeggieVampire/RTL-SDR_Signal_Detection_Amateur_Bands.py

## RTL-SDR_Signal_Detection_Amateur_Bands.py
import rtlsdr
import numpy as np
from scipy import signal

# Connect to the SDR and configure it
sdr = rtlsdr.RtlSdr()
sdr.sample_rate = 2.4e6  # 2.4 MHz sample rate
sdr.gain = 'auto'

# Define the frequency ranges to scan
ranges = [
    {'lower_freq': 29.7e6, 'upper_freq': 49.99e6, 'mode': 'NFM', 'step': 10e3, 'band': 'VHF Low Band'},
    {'lower_freq': 50e6, 'upper_freq': 53.98e6, 'mode': 'NFM', 'step': 20e3, 'band': '6 Meter Amateur Band'},
    {'lower_freq': 54e6, 'upper_freq': 71.95e6, 'mode': 'WFM', 'step': 50e3, 'band': 'VHF TV Broadcast 2 – 4'},
    {'lower_freq': 144e6, 'upper_freq': 147.995e6, 'mode': 'NFM', 'step': 5e3, 'band': '2 Meter Amateur Band'},
    {'lower_freq': 225e6, 'upper_freq': 379.975e6, 'mode': 'AM', 'step': 25e3, 'band': 'UHF Aircraft Band'},
    {'lower_freq': 1240e6, 'upper_freq': 1300e6, 'mode': 'NFM', 'step': 25e3, 'band': '25 cm Amateur Band'},
]

# Initialize the array to hold the results
results = []

# Calculate the total number of steps to take
num_steps = sum((r['upper_freq'] - r['lower_freq']) // r['step'] for r in ranges)

# Initialize the step counter
steps_taken = 0

# Scan each frequency range
for r in ranges:
    # Set the lower and upper bounds of the frequency range
    lower_freq = r['lower_freq']
    upper_freq = r['upper_freq']
    mode = r['mode']
    band = r['band']
    step_size = r['step']

    # Scan the frequency range in steps
    freq = lower_freq
    while freq <= upper_freq:
        # Set the center frequency of the SDR
        sdr.center_freq = freq

        # Read in samples from the SDR
        samples = sdr.read_samples(256 * 1024)

        # Calculate the power spectrum of the signal
        f, Pxx_den = signal.welch(np.real(samples), fs=sdr.sample_rate)

        # Find the indices of the power spectrum corresponding to the frequencies of interest
        indices = np.where(np.logical_and(f >= lower_freq, f <= upper_freq))[0]

        # Get the frequencies and power values at those indices
        frequencies = f[indices]
        power = Pxx_den[indices]

        # Check if there are any frequencies in the specified range
        if len(frequencies) > 0:
            # Find the frequencies with the highest power
            max_power_indices = np.argmax(power)
            max_frequencies = frequencies[max_power_indices]

            # Append the results to the array and save them to a file
            for freq, pwr in zip(max_frequencies, power[max_power_indices]):
                results.append((freq, pwr, mode, band))
                print(f"Signal found at frequency {freq / 1e6:.2f} MHz in mode {mode} in the {band} band")
                with open('frequencies.txt', 'a') as f:
                    f.write(f"{freq}: {pwr}: {mode}: {band}\n")

        # Increment the frequency counter by the step size
        lower_freq += step_size
        steps_taken += 1

        # Print the progress and current frequency
        total_progress = 100 * steps_taken / num_steps
        progress = 100 * (lower_freq - r['lower_freq']) / (upper_freq - r['lower_freq'])

        print(f"Frequency: {lower_freq / 1e6:.2f} MHz |  mode: {mode} | {band} progress: {progress:.2f}% | Total progress: {total_progress:.2f}%")

# Clean up
sdr.close()
print ("Scan completed. Check the frequencies.txt file for the found frequencies if any.")
del sdr
	import rtlsdr
	import numpy as np
	from scipy import signal

	# Connect to the SDR and configure it
	sdr = rtlsdr.RtlSdr()
	sdr.sample_rate = 2.4e6 # 2.4 MHz sample rate
	sdr.gain = 'auto'

	# Define the frequency ranges to scan
	ranges = [
	{'lower_freq': 29.7e6, 'upper_freq': 49.99e6, 'mode': 'NFM', 'step': 10e3, 'band': 'VHF Low Band'},
	{'lower_freq': 50e6, 'upper_freq': 53.98e6, 'mode': 'NFM', 'step': 20e3, 'band': '6 Meter Amateur Band'},
	{'lower_freq': 54e6, 'upper_freq': 71.95e6, 'mode': 'WFM', 'step': 50e3, 'band': 'VHF TV Broadcast 2 – 4'},
	{'lower_freq': 144e6, 'upper_freq': 147.995e6, 'mode': 'NFM', 'step': 5e3, 'band': '2 Meter Amateur Band'},
	{'lower_freq': 225e6, 'upper_freq': 379.975e6, 'mode': 'AM', 'step': 25e3, 'band': 'UHF Aircraft Band'},
	{'lower_freq': 1240e6, 'upper_freq': 1300e6, 'mode': 'NFM', 'step': 25e3, 'band': '25 cm Amateur Band'},
	]

	# Initialize the array to hold the results
	results = []

	# Calculate the total number of steps to take
	num_steps = sum((r['upper_freq'] - r['lower_freq']) // r['step'] for r in ranges)

	# Initialize the step counter
	steps_taken = 0

	# Scan each frequency range
	for r in ranges:
	# Set the lower and upper bounds of the frequency range
	lower_freq = r['lower_freq']
	upper_freq = r['upper_freq']
	mode = r['mode']
	band = r['band']
	step_size = r['step']

	# Scan the frequency range in steps
	freq = lower_freq
	while freq <= upper_freq:
	# Set the center frequency of the SDR
	sdr.center_freq = freq

	# Read in samples from the SDR
	samples = sdr.read_samples(256 * 1024)

	# Calculate the power spectrum of the signal
	f, Pxx_den = signal.welch(np.real(samples), fs=sdr.sample_rate)

	# Find the indices of the power spectrum corresponding to the frequencies of interest
	indices = np.where(np.logical_and(f >= lower_freq, f <= upper_freq))[0]

	# Get the frequencies and power values at those indices
	frequencies = f[indices]
	power = Pxx_den[indices]

	# Check if there are any frequencies in the specified range
	if len(frequencies) > 0:
	# Find the frequencies with the highest power
	max_power_indices = np.argmax(power)
	max_frequencies = frequencies[max_power_indices]

	# Append the results to the array and save them to a file
	for freq, pwr in zip(max_frequencies, power[max_power_indices]):
	results.append((freq, pwr, mode, band))
	print(f"Signal found at frequency {freq / 1e6:.2f} MHz in mode {mode} in the {band} band")
	with open('frequencies.txt', 'a') as f:
	f.write(f"{freq}: {pwr}: {mode}: {band}\n")

	# Increment the frequency counter by the step size
	lower_freq += step_size
	steps_taken += 1

	# Print the progress and current frequency
	total_progress = 100 * steps_taken / num_steps
	progress = 100 * (lower_freq - r['lower_freq']) / (upper_freq - r['lower_freq'])

	print(f"Frequency: {lower_freq / 1e6:.2f} MHz \| mode: {mode} \| {band} progress: {progress:.2f}% \| Total progress: {total_progress:.2f}%")

	# Clean up
	sdr.close()
	print ("Scan completed. Check the frequencies.txt file for the found frequencies if any.")
	del sdr