Changeset 33 for mire/process.py


Ignore:
Timestamp:
10/28/08 17:50:04 (16 years ago)
Author:
meynadie
Message:

filtrage des données couronne

File:
1 edited

Legend:

Unmodified
Added
Removed
  • mire/process.py

    r31 r33  
    349349        if (x_rel < 0): 
    350350            tmp['theta']+=math.pi 
    351         if (tmp['theta'] < 0): 
     351        # Offset de pi pour se mettre en convention "astro" 
     352        # theta = 0 -> P droit quand nord vers le haut 
     353        tmp['theta']+=math.pi 
     354        while (tmp['theta'] < 0): 
    352355            tmp['theta'] += 2*math.pi 
    353         if (tmp['theta'] >= 2*math.pi): 
     356        while (tmp['theta'] >= 2*math.pi): 
    354357            tmp['theta'] -= 2*math.pi 
    355         couronne.append(tmp) 
     358        # On n'écrit pas les valeurs aberrantes 
     359        if ((tmp['ScFact'] > 1.05) and (tmp['ScFact'] < 1.07)): 
     360            couronne.append(tmp) 
    356361 
    357362output_couronne = file(options.ring,'w') 
Note: See TracChangeset for help on using the changeset viewer.