# List RTL-SDR devices
rtl_eeprom -d 0

# Test reception
rtl_fm -d 0 -f 137620000 -s 48000 -g 40 - | play -r 48000 -t raw -e s16 -c 1 -V1 -

# Record NOAA pass
noaa-apt -d 0 -f 137620000 -g 40 -o recording.wav

# /etc/systemd/system/noaa-recorder.service
[Unit]
Description=NOAA Satellite Recorder
After=network.target

[Service]
Type=simple
User=wez
WorkingDirectory=/home/wez/scripts
ExecStart=/usr/bin/python3 /home/wez/scripts/noaa_recorder.py
Restart=on-failure
RestartSec=30

[Install]
WantedBy=multi-user.target

#!/usr/bin/env python3
"""
NOAA Automatic Recorder
Predicts passes and records NOAA APT signals
"""

import subprocess
from datetime import datetime, timedelta
from orbit_predictor.sources import EtcTLESource
from orbit_predictor.locations import Location

# Ground station location
GROUND_STATION = Location(
    latitude=38.0,  # Configure for your location
    longitude=-97.0,
    elevation_m=300
)

# Active NOAA satellites
NOAA_SATELLITES = {
    'NOAA-15': {'norad_id': 25338, 'frequency': 137620000},
    'NOAA-18': {'norad_id': 28654, 'frequency': 137912500},
    'NOAA-19': {'norad_id': 33591, 'frequency': 137100000},
}

def predict_next_pass(satellite_name, hours=24):
    """Predict next satellite pass"""
    source = EtcTLESource.from_url('https://celestrak.org/NORAD/elements/weather.txt')
    predictor = source.get_predictor(satellite_name)
    
    passes = predictor.passes_over(GROUND_STATION, limit=hours*60)
    return passes

def record_pass(satellite_name, frequency, duration_seconds, output_file):
    """Record NOAA APT signal"""
    cmd = [
        'rtl_fm',
        '-d', '0',           # Device index
        '-f', str(frequency), # Frequency in Hz
        '-s', '48000',        # Sample rate
        '-g', '40',           # Gain
        '-p', '1',            # PPM correction
        output_file
    ]
    
    subprocess.run(cmd, timeout=duration_seconds + 30)

def decode_apt(wav_file, output_image):
    """Decode APT signal to image"""
    subprocess.run(['wxtoimg', '-t', 'auto', wav_file, output_image])

class Satellite(db.Model):
    """Satellite TLE data and metadata"""
    __tablename__ = 'satellites'
    
    id = db.Column(db.Integer, primary_key=True)
    norad_id = db.Column(db.Integer, unique=True, nullable=False, index=True)
    name = db.Column(db.String(100), nullable=False)
    tle_line1 = db.Column(db.Text, nullable=False)
    tle_line2 = db.Column(db.Text, nullable=False)
    satellite_type = db.Column(db.String(50))  # 'weather', 'ham', 'gps'
    frequency = db.Column(db.Float)  # MHz
    active = db.Column(db.Boolean, default=True)
    tle_updated_at = db.Column(db.DateTime, default=datetime.utcnow)

class SatellitePass(db.Model):
    """Predicted satellite pass over ground station"""
    __tablename__ = 'satellite_passes'
    
    id = db.Column(db.Integer, primary_key=True)
    satellite_id = db.Column(db.Integer, db.ForeignKey('satellites.id'))
    aos_time = db.Column(db.DateTime, nullable=False)  # Acquisition of signal
    los_time = db.Column(db.DateTime, nullable=False)  # Loss of signal
    max_elevation = db.Column(db.Float, nullable=False)  # Degrees
    aos_azimuth = db.Column(db.Float)  # Rise direction
    los_azimuth = db.Column(db.Float)  # Set direction
    duration_minutes = db.Column(db.Float)
    pass_type = db.Column(db.String(20))  # 'ascending', 'descending'

class SatelliteReception(db.Model):
    """Recorded satellite reception"""
    __tablename__ = 'satellite_receptions'
    
    id = db.Column(db.Integer, primary_key=True)
    satellite_id = db.Column(db.Integer, db.ForeignKey('satellites.id'))
    pass_id = db.Column(db.Integer, db.ForeignKey('satellite_passes.id'))
    user_id = db.Column(db.Integer, db.ForeignKey('users.id'))
    
    # Reception details
    start_time = db.Column(db.DateTime, nullable=False)
    end_time = db.Column(db.DateTime)
    frequency = db.Column(db.Float)  # MHz
    sample_rate = db.Column(db.Integer)  # Hz
    gain = db.Column(db.Float)  # dB
    
    # File paths
    raw_file = db.Column(db.String(255))  # .wav or .raw
    image_file = db.Column(db.String(255))  # .png or .jpg
    
    # Processing status
    decoded = db.Column(db.Boolean, default=False)
    uploaded = db.Column(db.Boolean, default=False)
    
    # Quality metrics
    snr = db.Column(db.Float)  # Signal-to-noise ratio
    max_signal = db.Column(db.Float)  # Peak signal strength

# Record NOAA-19 pass
import requests

# Schedule recording
response = requests.post(
    'https://satellite.stsgym.com/api/satellites/receptions',
    json={
        'satellite_id': 1,  # NOAA-19
        'pass_id': 123,
        'auto_record': True
    }
)

print(response.json())
# {
#   "reception_id": 456,
#   "status": "scheduled",
#   "start_time": "2024-03-15T08:30:00Z"
# }

# Update TLE data daily at 3 AM UTC
0 3 * * * /usr/bin/python3 /home/wez/scripts/update_tle.py

# Predict passes every 5 minutes
*/5 * * * * /usr/bin/python3 /home/wez/scripts/predict_passes.py

# Record temperature every 5 minutes
*/5 * * * * /usr/bin/python3 /home/wez/scripts/temp_monitor.py --record --alert

# Sync temperature data every 10 minutes
*/10 * * * * /usr/bin/python3 /home/wez/scripts/temp_monitor.py --sync --hours 1

# Daily temperature report at 5 AM
0 5 * * * /usr/bin/python3 /home/wez/scripts/temp_monitor.py --report

-- Speed up pass queries
CREATE INDEX idx_passes_time ON satellite_passes (aos_time, los_time);

-- Speed up reception queries
CREATE INDEX idx_receptions_satellite ON satellite_receptions (satellite_id, start_time);

# Cache pass predictions
from functools import lru_cache
from datetime import datetime, timedelta

@lru_cache(maxsize=128)
def get_cached_passes(satellite_id: int, date: str):
    # Cache passes for 1 hour
    return SatellitePass.query.filter(
        SatellitePass.satellite_id == satellite_id,
        SatellitePass.aos_time > datetime.utcnow()
    ).all()