Python eiscat_3d_module Exemples

Exemple #1

#!/usr/bin/env python

import numpy as n
import matplotlib.pyplot as plt
import debris
import radar_library as rl

radars=[]
radars.append(rl.eiscat_uhf())
radars.append(rl.tromso_space_radar())
radars.append(rl.eiscat_svalbard())
radars.append(rl.eiscat_3d_module())
radars.append(rl.eiscat_3d())

range_m=10**n.linspace(5,8,num=100)

for r in radars:
    min_diam=n.zeros(len(range_m))
    gain_tx=r._tx[0].beam.I_0
    gain_rx=r._rx[0].beam.I_0
    rx_noise=r._tx[0].rx_noise
    B=r._tx[0].coh_int_bandwidth
    wavelength=r._tx[0].wavelength
    tx_power=r._tx[0].tx_power
    print(r._tx[0])

    for ri,rng_m in enumerate(range_m):
        min_diam[ri]=debris.target_diameter(gain_tx, gain_rx, wavelength, tx_power, rng_m, rng_m, enr=1.0, bandwidth=B, rx_noise_temp=rx_noise)
    plt.loglog(range_m/1e3,min_diam,label=r.name)
    
plt.xlabel("Range (km)")

Exemple #2

Fichier : DEV_demo2.py Projet : zzwei1/SORTSpp

#Project 2: Lets put an EISCAT 3D module in my home island!

import radar_library
import radar_config

import matplotlib.pyplot as plt

radar = radar_library.eiscat_3d_module()

radar._tx[0].lat = 59.993693
radar._tx[0].lon = 20.166953

radar._rx[0].lat = 59.993693
radar._rx[0].lon = 20.166953

radar._rx[1].lat = 60.008134
radar._rx[1].lon = 18.566308
radar._rx[2].lat = 60.555240
radar._rx[2].lon = 21.428962

#Now we dont want to plt
#radar_config.plot_radar(radar)
#plt.show()

Exemple #3

import copy

import numpy as np
import matplotlib.pyplot as plt

sys.path.insert(0, "/home/danielk/IRF/IRF_GITLAB/SORTSpp")

import population_library as plib
import radar_library as rlib
from propagator_orekit import PropagatorOrekit
from sorts_config import p as default_propagator
import dpt_tools as dpt
radars = [
    rlib.eiscat_3d(),
    #rlib.eiscat_3d(stage = 2),
    rlib.eiscat_3d_module(),
]
plot_out = '/home/danielk/IRF/E3D_PA/FP_sims/plots/population_filter/'
master_in = "./master/celn_20090501_00.sim"

base_pop = plib.master_catalog_factor(
    input_file=master_in,
    treshhold=0.01,
    seed=65487945,
    propagator=default_propagator,
    propagator_options={},
)

SNRs = [1.0, 10.0]

for SNR_lim in SNRs:

Exemple #4

Fichier : DEV_radar_library.py Projet : zzwei1/SORTSpp

import numpy as np
import scipy.constants as c

import matplotlib.pyplot as plt

from radar_config import plot_radar
import radar_library as rl

radars = [
    rl.eiscat_3d(beam='gauss'),
    #rl.eiscat_3d(beam = 'array'),
    rl.eiscat_3d_module(beam='gauss'),
    #rl.eiscat_3d_module(beam = 'array'),
    rl.eiscat_svalbard(),
    rl.eiscat_uhf(),
]

for radar in radars:
    plot_radar(radar)