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
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