Exemple #1
0
end_i = np.abs(km_utc_times - utc_timeT[-1]).argmin()

for section in xrange(section, section + 1):
    print section
    mean_x, mean_y, mean_alt, mean_pitch, mean_roll, mean_vel = ro.posav_section_info(
        m, posAV[km_idxs[section]:km_idxs[section + 1]])
    print 'Mean altitude:', mean_alt
    print 'Mean pitch:', mean_pitch
    print 'Mean roll:', mean_roll

    if (abs(mean_alt - 500) < 200) & (abs(mean_pitch) < 5) & (abs(mean_roll) <
                                                              5):
        poly_path, vertices, sides = ro.get_pos_poly(xp, yp, km_idxs[section],
                                                     km_idxs[section + 1])

        xatm_km, yatm_km, elevation_km = ro.get_atm_poly(
            xT, yT, elevationT, km_utc_times, utc_timeT, poly_path, section)
        if (size(xatm_km) > 0):
            lonDMS, latDMS = m(xatm_km[0], yatm_km[0], inverse=True)
            lonDMS_str = '%.2f' % lonDMS
            latDMS_str = '%.2f' % latDMS
            print lonDMS_str, latDMS_str
            xatm_km = xatm_km - np.amin(xatm_km)
            yatm_km = yatm_km - np.amin(yatm_km)
            spacing, minS, meanS, nearmaxS, maxS = ro.shot_spacing(
                xatm_km, yatm_km)

            mean_altSTR = ' %.0f' % mean_alt
            mean_pitchSTR = ' %.0f' % mean_pitch
            mean_rollSTR = ' %.0f' % mean_roll
            mean_velSTR = ' %.0f' % mean_vel
            meanSSTR = ' %.2f' % meanS
Exemple #2
0
xT, yT = m(lonT, latT)
#GET POSAV INDICES COINCIDING WITH START AND END OF ATM FILE. ADD PLUS/MINUS 1 FOR SOME LEEWAY.
start_i = np.abs(km_utc_times - utc_timeT[0]).argmin()
end_i = np.abs(km_utc_times - utc_timeT[-1]).argmin()

for section in xrange(section, section+1):
	print section
	mean_x, mean_y, mean_alt, mean_pitch, mean_roll, mean_vel = ro.posav_section_info(m, posAV[km_idxs[section]:km_idxs[section+1]]	)
	print 'Mean altitude:', mean_alt
	print 'Mean pitch:', mean_pitch
	print 'Mean roll:', mean_roll

	if (abs(mean_alt-500)<200) & (abs(mean_pitch)<5) & (abs(mean_roll)<5):
		poly_path, vertices, sides = ro.get_pos_poly(xp, yp, km_idxs[section], km_idxs[section+1])

		xatm_km, yatm_km, elevation_km = ro.get_atm_poly(xT, yT, elevationT, km_utc_times, utc_timeT, poly_path, section)
		if (size(xatm_km)>0):
			lonDMS, latDMS = m(xatm_km[0], yatm_km[0], inverse=True)
			lonDMS_str = '%.2f' %lonDMS
			latDMS_str = '%.2f' %latDMS
			print lonDMS_str, latDMS_str
			xatm_km = xatm_km - np.amin(xatm_km)
			yatm_km = yatm_km - np.amin(yatm_km)
			spacing, minS, meanS, nearmaxS, maxS = ro.shot_spacing(xatm_km, yatm_km)

			mean_altSTR=' %.0f' % mean_alt
			mean_pitchSTR=' %.0f' % mean_pitch
			mean_rollSTR=' %.0f' % mean_roll
			mean_velSTR=' %.0f' % mean_vel
			meanSSTR = ' %.2f' % meanS
			nearmaxSSTR = ' %.2f' % nearmaxS