#!/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)")
#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()
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:
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)