walyssonBarbosa
commented
10 years ago Used planets_and_moons.txt and 'jup-', 'gany-', 'sat-', 'copernicus', 'gan-'] as parameters for planets or moons observations. Found 494 distinct observations in the filesystem using the following script:
#! /usr/bin/env python
from mtpipeline.get_settings import SETTINGS
import glob, os
from astropy.io import fits
locations = [('acs_input_path', 'flt.fits'), ('wfpc2_output_path', 'c0m.fits'), ('wfc3_input_path', 'flt.fits')]
files_list = []
for l, e in locations:
all_files_list = glob.glob(os.path.join(SETTINGS[l], '*_*/*{}'.format(e)))
files_list += [filename for filename
in all_files_list
if len(filename.split('/')[-1]) == 18
and filename.split('/')[-1].split('_')[-1] == '{}'.format(e)]
files_dict = {}
for f in files_list:
hdu = fits.open(f)
head = hdu[0].header['TARGNAME']
if head not in files_dict:
files_dict[head] = {'list': [f], 'category': ''}
else:
if f not in files_dict[head]['list']:
files_dict[head]['list'].append(f)
planets_moons_list = ['jup-', 'gany-', 'sat-', 'copernicus', 'gan-']
pm_file = open('mtpipeline/ephem/planets_and_moons.txt', 'r')
for pm in pm_file:
obj = pm.split(' ')
if len(obj) > 3:
planets_moons_list.append(obj[1])
for pla in planets_moons_list:
if pla == 'io':
pla = 'io-'
found = {}
for pm in planets_moons_list:
for obj in files_dict:
if pm.upper() in obj:
found[obj] = 0
print len(found)I looked over the list of files that are neither planet nor moon observation, they all seem to be asteroids, comets and other things. Next step is to create the list of files that are in this observations.
acviana
ktfhale
I'm not sure how to interpret the labels because I can't see any object inside them.
Create a script that goes over all
fitsandfltfiles in the filesystem (ACS, WFC3 and WFPC2). Create a dictionary and for each file check itsTARGNAME, categorize the file aspmornotpm(Planet or Moon, Not Planet nor Moon).