Geometrical visibility start and end times - DEPRECATED
[visi_dates] = CL_ev_visibility(cjd,mean_kep,stations,stations_masks, ... sim_period, [visi_min,prec,propag_model,er,mu,zonals,obla])
This function is deprecated.
Replacement function: CL_ev_visibilityExtr
Computes the periods of time (start and end times) when a satellite is visible from a given set of ground stations.
The satellite is visible from a ground station when its elevation is over a given threshold (stations_masks).
The satellite trajectory is computed using an analytical model (specified by propag_model).
The results are computed in the simulation period defined by sim_period.
The intervals returned in visi_dates define the visibility periods start and end times of the satellite by at least one ground station. It means that the visibility intervals for ground stations considered independently are concatenated.
Notes:
- Visibility intervals with length less than visi_min are not computed. This parameter is also used for the detection of visibility intervals. So that choosing a small value for visi_min will increase the computation time.
- Setting planet oblateness (obla) to 0 (so that computation are done assuming a spherical planet) results in faster computation.
- The inertial reference frame where the orbit parameters are defined is necessarily Gamma50 (Veis) as CL_mod_SidTimeG50 is used to convert positions to the rotating (planet fixed) reference frame.
- The orbit parameters are the "classical" orbital elements.
Modified julian days from 1950.0 (TUC) (1x1)
Satellite's keplerian mean orbital elements at time cjd in Gamma50 (Veis) frame [sma;ecc;inc;pom;raan;anm] (6x1)
Stations positions in the rotating (planet fixed) reference frame in elliptical (geodetic) coordinates [long,lat,alt] [rad,rad,m] (3xN)
Station minimum elevations (above which there can be visibility) [rad] (1xN or 1x1)
Simulation time interval ([cjd_start; cjd_end]) in modified julian days from 1950.0 (TUC) (2x1)
(optional) Minimum visibility duration (default is 60 seconds) [sec] (1x1)
(optional) Computation accuracy on start/end visibility times (default is 1 second) [sec] (1x1)
(optional) Propagation model: 'kep' for keplerian, 'j2' for secular J2, 'lyd' for lyddane or 'eh' for Eckstein Hechler (default is lyddane) (1x1)
(optional) Planet equatorial radius (default is %CL_eqRad) [m] (1x1)
(optional) Gravitational constant [m^3/s^2] (default value is %CL_mu)
(optional) Vector of zonals coefficients J1 to Jn (troncated to J5) to be used (default is %CL_j1jn(1:5)) (1 x N)
(optional) Planet oblateness (default is %CL_obla) (1x1)
Visibility start and end times: [cjd_visi_start ; cjd_visi_end] in modified julian days from 1950.0 (TUC) (2xM)
CNES - DCT/SB
t0 = 21915; mean_kep0 = [7070.e3 ; 0.001 ; CL_deg2rad([98;90;10;15])]; // Definition of ground stations sta1 = [CL_deg2rad(2);CL_deg2rad(70);200]; // high latitude sta2 = [CL_deg2rad(20);CL_deg2rad(0);400]; // equator stations = [sta1,sta2]; stations_masks = [ CL_deg2rad(10) , CL_deg2rad(2) ]; sim_period = [21915 ; 21918 ]; // 3 days // Visibility computation [visi_dates] = CL_ev_visibility(t0, mean_kep0, stations, .. stations_masks, sim_period); // Plot visibility duration (mn) as function of time scf(); plot2d3(visi_dates(1,:) - t0, .. (visi_dates(2,:) - visi_dates(1,:)) * 1440, style=2); // Same computation with obla=0 (faster) [visi_dates] = CL_ev_visibility(t0, mean_kep0, stations, .. stations_masks, sim_period, obla=0); // Plot visibility duration (mn) as function of time plot2d3(visi_dates(1,:) - t0, .. (visi_dates(2,:) - visi_dates(1,:)) * 1440, style=5); |