Python yanny Exemples

Exemple #1

Fichier : run_rm.py Projet : joyce2015/redmonster

                            star_flags)
    # Write flags file
    prihdu = fits.PrimaryHDU(zpick.zwarning)
    thdulist = fits.HDUList([prihdu])
    thdulist.writeto(
        '/uufs/astro.utah.edu/common/home/u0814744/scratch/screens/flags-%s.fits'
        % plate,
        overwrite=True)
    output = io.WriteRedmonster(
        zpick,
        dest='/uufs/astro.utah.edu/common/home/u0814744/scratch/screens',
        overwrite=True)


# Read yanny file
x = y.yanny(filename='spInspect_alltest_bolton.par.txt', np=True)

# Get fibers, zpipe, zperson for each plate
args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0]
fibers3686 = []
zpipe3686 = []
zperson3686 = []
for i in args:
    fibers3686.append(x['BOSSOBJECT'][i][2] - 1)
    zpipe3686.append(x['BOSSOBJECT'][i][5])
    zperson3686.append(x['BOSSOBJECT'][i][6])

args = n.where(x['BOSSOBJECT']['plate'] == 3687)[0]
fibers3687 = []
zpipe3687 = []
zperson3687 = []

Exemple #2

Fichier : dchi2_optimize.py Projet : julienguy/redfit

import matplotlib.pyplot as p
p.interactive(True)
from redmonster.sandbox import yanny as y
from astropy.io import fits
from redmonster.datamgr import spec, io
from redmonster.physics import zfinder, zfitter, zpicker
from redmonster.physics import misc
from time import gmtime, strftime
from os.path import join
from os import environ
from scipy.optimize import curve_fit
from os.path import join
from os import environ

# Read yanny file
x = y.yanny(filename='spInspect_alltest_bolton.par.txt', np=True)

# Get fibers, zpipe, zperson for each plate
args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0]
fibers3686 = []
zpipe3686 = []
zperson3686 = []
classperson3686 = []
comments3686 = []
for i in args:
    fibers3686.append( x['BOSSOBJECT'][i][2]-1)
    zpipe3686.append( x['BOSSOBJECT'][i][5])
    zperson3686.append( x['BOSSOBJECT'][i][6])
    classperson3686.append( x['BOSSOBJECT'][i][7])
    comments3686.append( x['BOSSOBJECT'][i][8])

Exemple #3

Fichier : yanny_to_arrays.py Projet : desihub/redmonster

import numpy as n

from redmonster.sandbox import yanny as y

# Read yanny file
x = y.yanny(filename='/uufs/astro.utah.edu/common/home/u0814744/boss/\
            spInspect_alltest_bolton.par.txt', np=True)

# Get fibers, zpipe, zperson for each plate
args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0]
fibers3686 = []
zpipe3686 = []
zperson3686 = []
for i in args:
    fibers3686.append( x['BOSSOBJECT'][i][2])
    zpipe3686.append( x['BOSSOBJECT'][i][5])
    zperson3686.append( x['BOSSOBJECT'][i][6])

args = n.where(x['BOSSOBJECT']['plate'] == 3687)[0]
fibers3687 = []
zpipe3687 = []
zperson3687 = []
for i in args:
    fibers3687.append( x['BOSSOBJECT'][i][2])
    zpipe3687.append( x['BOSSOBJECT'][i][5])
    zperson3687.append( x['BOSSOBJECT'][i][6])

args = n.where(x['BOSSOBJECT']['plate'] == 3804)[0]
fibers3804 = []
zpipe3804 = []
zperson3804 = []

Exemple #4

Fichier : chi2_optimize.py Projet : julienguy/redfit

import matplotlib.pyplot as p
p.interactive(True)
from redmonster.sandbox import yanny as y
from astropy.io import fits
from redmonster.datamgr import spec, io
from redmonster.physics import zfinder, zfitter, zpicker
from redmonster.physics import misc
from time import gmtime, strftime
import multiprocessing as mp


startfib = 100
nfib = 1

# Read yanny file
x = y.yanny(filename='/Users/boltonlab3/Downloads/spInspect_alltest_bolton.par.txt', np=True)

# Get fibers, zpipe, zperson for each plate
args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0]
fibers3686 = []
zpipe3686 = []
zperson3686 = []
for i in args:
    fibers3686.append( x['BOSSOBJECT'][i][2])
    zpipe3686.append( x['BOSSOBJECT'][i][5])
    zperson3686.append( x['BOSSOBJECT'][i][6])

# Only keep 50 fibers
fiberid = fibers3686[startfib:startfib+nfib]
zpipe = zpipe3686[startfib:startfib+nfib]
zperson = zperson3686[startfib:startfib+nfib]

Exemple #5

import matplotlib.pyplot as p
p.interactive(True)
from redmonster.sandbox import yanny as y
from astropy.io import fits
from redmonster.datamgr import spec, io
from redmonster.physics import zfinder, zfitter, zpicker
from redmonster.physics import misc
from time import gmtime, strftime
import multiprocessing as mp


startfib = 100
nfib = 1

# Read yanny file
x = y.yanny(filename='/Users/boltonlab3/Downloads/spInspect_alltest_bolton.par.txt', np=True)

# Get fibers, zpipe, zperson for each plate
args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0]
fibers3686 = []
zpipe3686 = []
zperson3686 = []
for i in args:
    fibers3686.append( x['BOSSOBJECT'][i][2])
    zpipe3686.append( x['BOSSOBJECT'][i][5])
    zperson3686.append( x['BOSSOBJECT'][i][6])

# Only keep 50 fibers
fiberid = fibers3686[startfib:startfib+nfib]
zpipe = zpipe3686[startfib:startfib+nfib]
zperson = zperson3686[startfib:startfib+nfib]

Exemple #6

import numpy as n

from redmonster.sandbox import yanny as y

# Read yanny file
x = y.yanny(filename='/uufs/astro.utah.edu/common/home/u0814744/boss/\
            spInspect_alltest_bolton.par.txt',
            np=True)

# Get fibers, zpipe, zperson for each plate
args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0]
fibers3686 = []
zpipe3686 = []
zperson3686 = []
for i in args:
    fibers3686.append(x['BOSSOBJECT'][i][2])
    zpipe3686.append(x['BOSSOBJECT'][i][5])
    zperson3686.append(x['BOSSOBJECT'][i][6])

args = n.where(x['BOSSOBJECT']['plate'] == 3687)[0]
fibers3687 = []
zpipe3687 = []
zperson3687 = []
for i in args:
    fibers3687.append(x['BOSSOBJECT'][i][2])
    zpipe3687.append(x['BOSSOBJECT'][i][5])
    zperson3687.append(x['BOSSOBJECT'][i][6])

args = n.where(x['BOSSOBJECT']['plate'] == 3804)[0]
fibers3804 = []
zpipe3804 = []