notebook.community

Edit and run

Flattening echelle spectra

Often-times, the model spectrum is just not good enough to get an adequate fit using the normal Fitters.RVFitter. When that is the case, the automated flattening just doesn't work. Let's give this a shot without using models. The thing we want to minimize is the offset between the flux at one order and the flux in the adjacent order.

$$L = \sum_{i=1}^{N_{orders}-1} \left(\frac{O_i(\lambda)}{M(i, \lambda| \Theta)} - \frac{O_{i+1}(\lambda)}{M(i+1, \lambda | \theta)}\right)^2 $$

Where $M(i, \lambda | \theta)$ is perhaps a 2d polynomial function of some order where $i$ is for the order number, $\lambda$ is the wavelength (or maybe pixel number?), and the $\theta$ are all the polynomial coefficients.

TRY:

Try using Robust Linear Models from statsmodels. That might help it to ignore the big Balmer line. I will need to put the Chebyshev factors into a feature vector to use this...

Those might help a bit. How about a hybrid approach where I divide by a stellar model (or empirical spectrum of similar type as Adam suggests), and then do a full robust 2D fit.

Robust linear models works well enough for me!

I was able to sufficiently flatten a few of my spectra to get RV fits. I need to calibrate the RV though, so pretty much need to flatten every spectrum and fit the RV. This will take a while...





In [1]:



    


import SpecFlattener
import glob
import StarData
from astropy.io import fits
import SpectralTypeRelations
import logging

logger = logging.getLogger()
logger.setLevel(logging.INFO)



    


















    






Module 'anfft' (FFTW Python bindings) could not be imported.
To install it, try running 'easy_install anfft' from the terminal.
Falling back on the slower 'fftpack' module for ND Fourier transforms.










    






:0: FutureWarning: IPython widgets are experimental and may change in the future.

















In [2]:



    


#hdf5_lib = '/media/ExtraSpace/Kurucz_FullGrid/CHIRON_grid_full.hdf5'
hdf5_lib = '/Volumes/DATADRIVE/Kurucz_Grid/TS23_grid_full.hdf5'
star_list = [f for f in glob.glob('../201*/H*.fits') if 'smoothed' not in f and 
             'uncorrected' not in f and 'flattened' not in f and 'B.fits' not in f]
print(len(star_list))
#star_list.index('../20131019/HIP_22913.fits')
for s in star_list:
    print s



    


















    






144
../20130106/HIP_28691.fits
../20130107/HIP_16147.fits
../20130107/HIP_18724.fits
../20130107/HIP_21402.fits
../20130107/HIP_25048.fits
../20130107/HIP_25281.fits
../20130107/HIP_28614.fits
../20130107/HIP_29038.fits
../20130107/HIP_3300.fits
../20130107/HIP_36850_A.fits
../20130107/HIP_45688.fits
../20130107/HIP_55266.fits
../20130107/HIP_55434.fits
../20130823/HIP_116971.fits
../20130824/HIP_116805.fits
../20131017/HIP_101716.fits
../20131017/HIP_14576.fits
../20131017/HIP_15444.fits
../20131018/HIP_111841.fits
../20131018/HIP_116631.fits
../20131018/HIP_118243.fits
../20131018/HIP_12719.fits
../20131018/HIP_14576.fits
../20131018/HIP_16244.fits
../20131018/HIP_2912.fits
../20131018/HIP_98194.fits
../20131018/HIP_99080.fits
../20131019/HBC_722.fits
../20131019/HIP_106711.fits
../20131019/HIP_12706.fits
../20131019/HIP_97496.fits
../20131020/HIP_105966.fits
../20131020/HIP_111169.fits
../20131020/HIP_1366.fits
../20131020/HIP_18805.fits
../20131020/HIP_3881.fits
../20131020/HIP_4436.fits
../20131020/HIP_95853.fits
../20131020/HIP_96288.fits
../20131020/HR_604.fits
../20131021/HIP_109521.fits
../20131021/HIP_113788.fits
../20131021/HIP_117371.fits
../20131021/HIP_20789.fits
../20131021/HIP_2505.fits
../20131021/HIP_9312.fits
../20131021/HIP_98055.fits
../20131021/HR_545.fits
../20131022/HIP_114520.fits
../20131022/HIP_5131.fits
../20131022/HIP_5132.fits
../20131022/HIP_93747.fits
../20131022/HIP_98325.fits
../20131122/HIP_10670.fits
../20131122/HIP_106711.fits
../20131122/HIP_13209.fits
../20131122/HIP_13879.fits
../20131122/HIP_813.fits
../20140110/HIP_22509.fits
../20140110/HIP_41307.fits
../20140111/HIP_19799.fits
../20140111/HIP_31278.fits
../20140111/HIP_36760.fits
../20140111/HIP_39567.fits
../20140111/HIP_70400.fits
../20140111/HR_545.fits
../20140111/HR_8799.fits
../20140112/HIP_25695.fits
../20140112/HIP_3478.fits
../20140112/HIP_38538.fits
../20140112/HIP_52638.fits
../20140112/HIP_56633.fits
../20140112/HIP_59819.fits
../20140112/HIP_65477.fits
../20140113/HIP_13327.fits
../20140113/HIP_25555.fits
../20140113/HIP_36393.fits
../20140113/HIP_43142.fits
../20140113/HIP_52911.fits
../20140113/HIP_54849.fits
../20140119/HIP_12803.fits
../20140119/HIP_16210.fits
../20140119/HIP_19968.fits
../20140119/HIP_26093.fits
../20140119/HIP_31434.fits
../20140119/HIP_36393.fits
../20140119/HIP_39847.fits
../20140119/HIP_47006.fits
../20140119/HIP_52422.fits
../20140119/HIP_52457.fits
../20140119/HIP_56034.fits
../20140119/HIP_60595.fits
../20140120/HIP_12332.fits
../20140120/HIP_21928.fits
../20140120/HIP_26563.fits
../20140120/HIP_34769.fits
../20140120/HIP_44127.fits
../20140120/HIP_50303.fits
../20140120/HIP_51685.fits
../20140120/HIP_65728.fits
../20140121/HIP_16340.fits
../20140121/HIP_18141.fits
../20140121/HIP_22840.fits
../20140121/HIP_33372.fits
../20140121/HIP_50860.fits
../20140121/HIP_60957.fits
../20140121/HIP_75178.fits
../20140122/HIP_14764.fits
../20140122/HIP_21819.fits
../20140122/HIP_22833.fits
../20140122/HIP_29997.fits
../20140122/HIP_39236.fits
../20140122/HIP_46225.fits
../20140123/HIP_35341.fits
../20140124/HIP_13879.fits
../20140124/HIP_42313.fits
../20140124/HIP_70327.fits
../20141101/HIP_106786.fits
../20141101/HIP_109521.fits
../20141101/HIP_116631.fits
../20141101/HIP_17489.fits
../20141101/HIP_20542.fits
../20141101/HIP_20648.fits
../20141101/HIP_20711.fits
../20141101/HIP_20901.fits
../20141101/HIP_21029.fits
../20141101/HIP_21589.fits
../20141101/HIP_21670.fits
../20141101/HIP_21683.fits
../20141101/HIP_23497.fits
../20141101/HIP_2505.fits
../20141101/HIP_93225.fits
../20141102/HIP_17401.fits
../20141102/HIP_17489.fits
../20141102/HIP_17527.fits
../20141102/HIP_17531.fits
../20141102/HIP_3478.fits
../20141103/HIP_17573.fits
../20141103/HIP_17588.fits
../20141103/HIP_17608.fits
../20141103/HIP_17702.fits
../20141103/HIP_17776.fits
../20141103/HIP_17851.fits
../20141103/HIP_3881.fits

















In [3]:



    


# Guess stellar properties
MS = SpectralTypeRelations.MainSequence()
def guess_teff_logg(fname):
    header = fits.getheader(fname)
    data = StarData.GetData(header['OBJECT'])
    spt = data.spectype
    teff = MS.Interpolate('Temperature', spt)
    logg = 3.5 if 'I' in spt else 4.0
    return teff, logg



    











In [80]:



    


teff, logg = guess_teff_logg('../20131021/HIP_109521.fits')
print(teff, logg)



    


















    






INFO:requests.packages.urllib3.connectionpool:Starting new HTTP connection (1): simbak.cfa.harvard.edu










    






(array(8079.999999999998), 4.0)

















In [104]:



    


reload(SpecFlattener)
reload(StellarModel)
output = SpecFlattener.flatten_spec(star_list[1], hdf5_lib, teff=teff, logg=logg, normalize_model=False,
                                       first_order=0, last_order=15)
final_orders, flattened, shifted_orders, mcf = output

%matplotlib notebook
for order in final_orders:
    plt.plot(order.x, order.y, 'k-', alpha=0.5)



    


















    






INFO:root:Initializing continuum for RV guess.
INFO:root:Estimating the RV and vsini by cross-correlation










    






Determine Chunk Log: Wl is 131072
Warning, 1996 of 1996 bins contained negative fluxes; they have been set to zero.









    






INFO:root:T = 8178.93960024
logg = 3.84080221678
RV = -100.0
INFO:root:LogL = -88140997.3266

INFO:root:T = 8179.03960024
logg = 3.84080221678
RV = -100.0
INFO:root:LogL = -88375191.6723

INFO:root:T = 8178.93960024
logg = 3.94080221678
RV = -100.0
INFO:root:LogL = -87605511.767

INFO:root:T = 8178.93960024
logg = 3.84080221678
RV = -99.9
INFO:root:LogL = -88363871.2827

INFO:root:T = 7000.0
logg = 4.5
RV = -100.0
INFO:root:LogL = -87256894.1132

INFO:root:T = 7000.1
logg = 4.5
RV = -100.0
INFO:root:LogL = -87250341.3278

INFO:root:T = 7000.0
logg = 4.6
RV = -100.0
INFO:root:LogL = -87634892.6582

INFO:root:T = 7000.0
logg = 4.5
RV = -99.9
INFO:root:LogL = -87232399.2596

INFO:root:T = 7780.5926732
logg = 4.0635357743
RV = -100.0
INFO:root:LogL = -76708174.915

INFO:root:T = 7780.6926732
logg = 4.0635357743
RV = -100.0
INFO:root:LogL = -76709224.4276

INFO:root:T = 7780.5926732
logg = 4.1635357743
RV = -100.0
INFO:root:LogL = -75243775.9745

INFO:root:T = 7780.5926732
logg = 4.0635357743
RV = -99.9
INFO:root:LogL = -76674354.7683

INFO:root:T = 7773.43018163
logg = 4.5
RV = -96.7030717746
INFO:root:LogL = -69202814.3337

INFO:root:T = 7773.53018163
logg = 4.5
RV = -96.7030717746
INFO:root:LogL = -69204643.2297

INFO:root:T = 7773.43018163
logg = 4.6
RV = -96.7030717746
INFO:root:LogL = -67699213.0429

INFO:root:T = 7773.43018163
logg = 4.5
RV = -96.6030717746
INFO:root:LogL = -69176479.8005

INFO:root:T = 7765.90029633
logg = 4.5
RV = -95.4285963207
INFO:root:LogL = -68718426.8809

INFO:root:T = 7766.00029633
logg = 4.5
RV = -95.4285963207
INFO:root:LogL = -68718879.6971

INFO:root:T = 7765.90029633
logg = 4.6
RV = -95.4285963207
INFO:root:LogL = -67241219.2858

INFO:root:T = 7765.90029633
logg = 4.5
RV = -95.3285963207
INFO:root:LogL = -68682452.4755

INFO:root:T = 7735.78075514
logg = 4.5
RV = -90.3306945052
INFO:root:LogL = -67161091.582

INFO:root:T = 7735.88075514
logg = 4.5
RV = -90.3306945052
INFO:root:LogL = -67162588.5135

INFO:root:T = 7735.78075514
logg = 4.6
RV = -90.3306945052
INFO:root:LogL = -65753218.6257

INFO:root:T = 7735.78075514
logg = 4.5
RV = -90.2306945052
INFO:root:LogL = -67147824.9646

INFO:root:T = 7679.26249257
logg = 4.5
RV = -80.9844104621
INFO:root:LogL = -63168914.963

INFO:root:T = 7679.36249257
logg = 4.5
RV = -80.9844104621
INFO:root:LogL = -63168207.8503

INFO:root:T = 7679.26249257
logg = 4.6
RV = -80.9844104621
INFO:root:LogL = -62567039.1323

INFO:root:T = 7679.26249257
logg = 4.5
RV = -80.8844104621
INFO:root:LogL = -63147666.7176

INFO:root:T = 7486.51752669
logg = 4.5
RV = -47.6339352452
INFO:root:LogL = -60645795.6113

INFO:root:T = 7486.61752669
logg = 4.5
RV = -47.6339352452
INFO:root:LogL = -60644516.3874

INFO:root:T = 7486.51752669
logg = 4.6
RV = -47.6339352452
INFO:root:LogL = -60579297.2887

INFO:root:T = 7486.51752669
logg = 4.5
RV = -47.5339352452
INFO:root:LogL = -60653908.4822

INFO:root:T = 7578.00586899
logg = 4.5
RV = -63.4640747167
INFO:root:LogL = -60777767.1557

INFO:root:T = 7578.10586899
logg = 4.5
RV = -63.4640747167
INFO:root:LogL = -60775823.8591

INFO:root:T = 7578.00586899
logg = 4.6
RV = -63.4640747167
INFO:root:LogL = -60374386.9721

INFO:root:T = 7578.00586899
logg = 4.5
RV = -63.3640747167
INFO:root:LogL = -60764070.9611

INFO:root:T = 7562.29113351
logg = 4.5
RV = -60.3173829948
INFO:root:LogL = -60491567.5896

INFO:root:T = 7562.39113351
logg = 4.5
RV = -60.3173829948
INFO:root:LogL = -60490200.5592

INFO:root:T = 7562.29113351
logg = 4.6
RV = -60.3173829948
INFO:root:LogL = -60128947.2036

INFO:root:T = 7562.29113351
logg = 4.5
RV = -60.2173829948
INFO:root:LogL = -60479722.2645

INFO:root:T = 7499.43219159
logg = 4.5
RV = -47.7306161072
INFO:root:LogL = -60590235.1223

INFO:root:T = 7499.53219159
logg = 4.5
RV = -47.7306161072
INFO:root:LogL = -60589907.4703

INFO:root:T = 7499.43219159
logg = 4.6
RV = -47.7306161072
INFO:root:LogL = -60557092.0764

INFO:root:T = 7499.43219159
logg = 4.5
RV = -47.6306161072
INFO:root:LogL = -60607993.6955

INFO:root:T = 7530.24688017
logg = 4.5
RV = -53.9008965791
INFO:root:LogL = -60190763.4419

INFO:root:T = 7530.34688017
logg = 4.5
RV = -53.9008965791
INFO:root:LogL = -60190003.9342

INFO:root:T = 7530.24688017
logg = 4.6
RV = -53.9008965791
INFO:root:LogL = -59931784.2839

INFO:root:T = 7530.24688017
logg = 4.5
RV = -53.8008965791
INFO:root:LogL = -60186602.0328

INFO:root:T = 7528.25650825
logg = 4.5
RV = -53.2842961561
INFO:root:LogL = -60176897.4756

INFO:root:T = 7528.35650825
logg = 4.5
RV = -53.2842961561
INFO:root:LogL = -60176236.5899

INFO:root:T = 7528.25650825
logg = 4.6
RV = -53.2842961561
INFO:root:LogL = -59926519.3998

INFO:root:T = 7528.25650825
logg = 4.5
RV = -53.1842961561
INFO:root:LogL = -60172416.2418

INFO:root:T = 7520.29502058
logg = 4.5
RV = -50.8178944643
INFO:root:LogL = -60189362.853

INFO:root:T = 7520.39502058
logg = 4.5
RV = -50.8178944643
INFO:root:LogL = -60188433.1293

INFO:root:T = 7520.29502058
logg = 4.6
RV = -50.8178944643
INFO:root:LogL = -60000237.0647

INFO:root:T = 7520.29502058
logg = 4.5
RV = -50.7178944643
INFO:root:LogL = -60191904.9417

INFO:root:T = 7524.60786669
logg = 4.5
RV = -52.1539777925
INFO:root:LogL = -60170129.2731

INFO:root:T = 7524.70786669
logg = 4.5
RV = -52.1539777925
INFO:root:LogL = -60169184.9752

INFO:root:T = 7524.60786669
logg = 4.6
RV = -52.1539777925
INFO:root:LogL = -59937176.0507

INFO:root:T = 7524.60786669
logg = 4.5
RV = -52.0539777925
INFO:root:LogL = -60163580.8987

INFO:root:T = 7523.2014771
logg = 4.5
RV = -51.7182901661
INFO:root:LogL = -60168704.4598

INFO:root:T = 7523.3014771
logg = 4.5
RV = -51.7182901661
INFO:root:LogL = -60168044.6642

INFO:root:T = 7523.2014771
logg = 4.6
RV = -51.7182901661
INFO:root:LogL = -59948917.2311

INFO:root:T = 7523.2014771
logg = 4.5
RV = -51.6182901661
INFO:root:LogL = -60165398.4391

INFO:root:T = 7516.42565031
logg = 4.5
RV = -47.9272074152
INFO:root:LogL = -60575305.2365

INFO:root:T = 7516.52565031
logg = 4.5
RV = -47.9272074152
INFO:root:LogL = -60575544.7178

INFO:root:T = 7516.42565031
logg = 4.6
RV = -47.9272074152
INFO:root:LogL = -60440474.5682

INFO:root:T = 7516.42565031
logg = 4.5
RV = -47.8272074152
INFO:root:LogL = -60590655.7956

INFO:root:T = 7522.85580477
logg = 4.5
RV = -51.5248861094
INFO:root:LogL = -60165205.8806

INFO:root:T = 7522.95580477
logg = 4.5
RV = -51.5248861094
INFO:root:LogL = -60163675.0357

INFO:root:T = 7522.85580477
logg = 4.6
RV = -51.5248861094
INFO:root:LogL = -59952126.2478

INFO:root:T = 7522.85580477
logg = 4.5
RV = -51.4248861094
INFO:root:LogL = -60163562.4209

INFO:root:T = 7541.31539271
logg = 4.5
RV = -41.9806411689
INFO:root:LogL = -60175542.7901

INFO:root:T = 7541.41539271
logg = 4.5
RV = -41.9806411689
INFO:root:LogL = -60175399.2196

INFO:root:T = 7541.31539271
logg = 4.6
RV = -41.9806411689
INFO:root:LogL = -59917947.8069

INFO:root:T = 7541.31539271
logg = 4.5
RV = -41.8806411689
INFO:root:LogL = -60173508.7621

INFO:root:T = 7527.49783514
logg = 4.5
RV = -49.124795976
INFO:root:LogL = -60399583.0937

INFO:root:T = 7527.59783514
logg = 4.5
RV = -49.124795976
INFO:root:LogL = -60399808.7528

INFO:root:T = 7527.49783514
logg = 4.6
RV = -49.124795976
INFO:root:LogL = -60266772.6624

INFO:root:T = 7527.49783514
logg = 4.5
RV = -49.024795976
INFO:root:LogL = -60413712.7863

INFO:root:T = 7523.33144592
logg = 4.5
RV = -51.2789631972
INFO:root:LogL = -60156454.4426

INFO:root:T = 7523.43144592
logg = 4.5
RV = -51.2789631972
INFO:root:LogL = -60155358.3842

INFO:root:T = 7523.33144592
logg = 4.6
RV = -51.2789631972
INFO:root:LogL = -59951028.7429

INFO:root:T = 7523.33144592
logg = 4.5
RV = -51.1789631972
INFO:root:LogL = -60156406.8893

INFO:root:T = 7524.14163117
logg = 4.5
RV = -51.1656938068
INFO:root:LogL = -60147996.6176

INFO:root:T = 7524.24163117
logg = 4.5
RV = -51.1656938068
INFO:root:LogL = -60147037.5045

INFO:root:T = 7524.14163117
logg = 4.6
RV = -51.1656938068
INFO:root:LogL = -59946893.1288

INFO:root:T = 7524.14163117
logg = 4.5
RV = -51.0656938068
INFO:root:LogL = -60147990.3743

INFO:root:T = 7529.60019658
logg = 4.5
RV = -50.4336399664
INFO:root:LogL = -60148464.0085

INFO:root:T = 7529.70019658
logg = 4.5
RV = -50.4336399664
INFO:root:LogL = -60148277.4811

INFO:root:T = 7529.60019658
logg = 4.6
RV = -50.4336399664
INFO:root:LogL = -60023389.0752

INFO:root:T = 7529.60019658
logg = 4.5
RV = -50.3336399664
INFO:root:LogL = -60162317.5996

INFO:root:T = 7527.02226276
logg = 4.5
RV = -50.7793693315
INFO:root:LogL = -60131941.8783

INFO:root:T = 7527.12226276
logg = 4.5
RV = -50.7793693315
INFO:root:LogL = -60131256.5521

INFO:root:T = 7527.02226276
logg = 4.6
RV = -50.7793693315
INFO:root:LogL = -59961954.412

INFO:root:T = 7527.02226276
logg = 4.5
RV = -50.6793693315
INFO:root:LogL = -60137679.4285

INFO:root:T = 7530.10572753
logg = 4.5
RV = -50.968546714
INFO:root:LogL = -60101706.1156

INFO:root:T = 7530.20572753
logg = 4.5
RV = -50.968546714
INFO:root:LogL = -60101575.121

INFO:root:T = 7530.10572753
logg = 4.6
RV = -50.968546714
INFO:root:LogL = -59931180.46

INFO:root:T = 7530.10572753
logg = 4.5
RV = -50.868546714
INFO:root:LogL = -60108248.1658

INFO:root:T = 7534.60515117
logg = 4.5
RV = -52.0283946788
INFO:root:LogL = -60049800.7336

INFO:root:T = 7534.70515117
logg = 4.5
RV = -52.0283946788
INFO:root:LogL = -60048852.0693

INFO:root:T = 7534.60515117
logg = 4.6
RV = -52.0283946788
INFO:root:LogL = -59830006.9243

INFO:root:T = 7534.60515117
logg = 4.5
RV = -51.9283946788
INFO:root:LogL = -60050191.5269

INFO:root:T = 7564.88705667
logg = 4.5
RV = -56.8149321308
INFO:root:LogL = -60059159.6707

INFO:root:T = 7564.98705667
logg = 4.5
RV = -56.8149321308
INFO:root:LogL = -60058395.7954

INFO:root:T = 7564.88705667
logg = 4.6
RV = -56.8149321308
INFO:root:LogL = -59728355.57

INFO:root:T = 7564.88705667
logg = 4.5
RV = -56.7149321308
INFO:root:LogL = -60049503.4685

INFO:root:T = 7547.95106694
logg = 4.5
RV = -54.1379292219
INFO:root:LogL = -60006189.3302

INFO:root:T = 7548.05106694
logg = 4.5
RV = -54.1379292219
INFO:root:LogL = -60004173.2466

INFO:root:T = 7547.95106694
logg = 4.6
RV = -54.1379292219
INFO:root:LogL = -59733172.4675

INFO:root:T = 7547.95106694
logg = 4.5
RV = -54.0379292219
INFO:root:LogL = -59999611.1801

INFO:root:T = 7559.12882016
logg = 4.5
RV = -55.9047511418
INFO:root:LogL = -60021899.825

INFO:root:T = 7559.22882016
logg = 4.5
RV = -55.9047511418
INFO:root:LogL = -60021252.6979

INFO:root:T = 7559.12882016
logg = 4.6
RV = -55.9047511418
INFO:root:LogL = -59713865.6659

INFO:root:T = 7559.12882016
logg = 4.5
RV = -55.8047511418
INFO:root:LogL = -60013193.048

INFO:root:T = 7550.31474518
logg = 4.5
RV = -54.23146165
INFO:root:LogL = -59984636.4346

INFO:root:T = 7550.41474518
logg = 4.5
RV = -54.23146165
INFO:root:LogL = -59983421.6061

INFO:root:T = 7550.31474518
logg = 4.6
RV = -54.23146165
INFO:root:LogL = -59710954.0284

INFO:root:T = 7550.31474518
logg = 4.5
RV = -54.13146165
INFO:root:LogL = -59979894.6405

INFO:root:T = 7553.98720852
logg = 4.5
RV = -54.4050646801
INFO:root:LogL = -59958967.8679

INFO:root:T = 7554.08720852
logg = 4.5
RV = -54.4050646801
INFO:root:LogL = -59958152.2749

INFO:root:T = 7553.98720852
logg = 4.6
RV = -54.4050646801
INFO:root:LogL = -59678687.148

INFO:root:T = 7553.98720852
logg = 4.5
RV = -54.3050646801
INFO:root:LogL = -59950317.2266

INFO:root:T = 7569.16439097
logg = 4.5
RV = -53.2369929996
INFO:root:LogL = -59932141.9311

INFO:root:T = 7569.26439097
logg = 4.5
RV = -53.2369929996
INFO:root:LogL = -59932982.6318

INFO:root:T = 7569.16439097
logg = 4.6
RV = -53.2369929996
INFO:root:LogL = -59753239.9398

INFO:root:T = 7569.16439097
logg = 4.5
RV = -53.1369929996
INFO:root:LogL = -59947875.2437

INFO:root:T = 7560.3510221
logg = 4.5
RV = -53.9152906098
INFO:root:LogL = -59882990.8784

INFO:root:T = 7560.4510221
logg = 4.5
RV = -53.9152906098
INFO:root:LogL = -59881848.0782

INFO:root:T = 7560.3510221
logg = 4.6
RV = -53.9152906098
INFO:root:LogL = -59621643.9906

INFO:root:T = 7560.3510221
logg = 4.5
RV = -53.8152906098
INFO:root:LogL = -59880477.4545

INFO:root:T = 7563.82723826
logg = 4.5
RV = -53.6477528332
INFO:root:LogL = -59869642.8327

INFO:root:T = 7563.92723826
logg = 4.5
RV = -53.6477528332
INFO:root:LogL = -59868951.2918

INFO:root:T = 7563.82723826
logg = 4.6
RV = -53.6477528332
INFO:root:LogL = -59638709.8781

INFO:root:T = 7563.82723826
logg = 4.5
RV = -53.5477528332
INFO:root:LogL = -59871723.9025

INFO:root:T = 7578.62416185
logg = 4.5
RV = -53.3315705002
INFO:root:LogL = -60011956.2611

INFO:root:T = 7578.72416185
logg = 4.5
RV = -53.3315705002
INFO:root:LogL = -60012582.4689

INFO:root:T = 7578.62416185
logg = 4.6
RV = -53.3315705002
INFO:root:LogL = -59806564.6962

INFO:root:T = 7578.62416185
logg = 4.5
RV = -53.2315705002
INFO:root:LogL = -60021813.7947

INFO:root:T = 7565.38743842
logg = 4.5
RV = -53.6144142999
INFO:root:LogL = -59867563.2426

INFO:root:T = 7565.48743842
logg = 4.5
RV = -53.6144142999
INFO:root:LogL = -59867548.3791

INFO:root:T = 7565.38743842
logg = 4.6
RV = -53.6144142999
INFO:root:LogL = -59649981.4894

INFO:root:T = 7565.38743842
logg = 4.5
RV = -53.5144142999
INFO:root:LogL = -59874247.022

INFO:root:T = 7566.75145014
logg = 4.5
RV = -53.8583213765
INFO:root:LogL = -59858469.4125

INFO:root:T = 7566.85145014
logg = 4.5
RV = -53.8583213765
INFO:root:LogL = -59858154.0226

INFO:root:T = 7566.75145014
logg = 4.6
RV = -53.8583213765
INFO:root:LogL = -59625227.2126

INFO:root:T = 7566.75145014
logg = 4.5
RV = -53.7583213765
INFO:root:LogL = -59862542.7034

INFO:root:T = 7591.22827052
logg = 4.5
RV = -57.2215351707
INFO:root:LogL = -59863291.7885

INFO:root:T = 7591.32827052
logg = 4.5
RV = -57.2215351707
INFO:root:LogL = -59862814.8304

INFO:root:T = 7591.22827052
logg = 4.6
RV = -57.2215351707
INFO:root:LogL = -59532419.6345

INFO:root:T = 7591.22827052
logg = 4.5
RV = -57.1215351707
INFO:root:LogL = -59858185.1308

INFO:root:T = 7575.59795586
logg = 4.5
RV = -55.0738669259
INFO:root:LogL = -59839482.2882

INFO:root:T = 7575.69795586
logg = 4.5
RV = -55.0738669259
INFO:root:LogL = -59839348.2871

INFO:root:T = 7575.59795586
logg = 4.6
RV = -55.0738669259
INFO:root:LogL = -59548976.1159

INFO:root:T = 7575.59795586
logg = 4.5
RV = -54.9738669259
INFO:root:LogL = -59837745.3597

INFO:root:T = 7574.81274924
logg = 4.5
RV = -54.9390275962
INFO:root:LogL = -59840885.6573

INFO:root:T = 7574.91274924
logg = 4.5
RV = -54.9390275962
INFO:root:LogL = -59839421.6541

INFO:root:T = 7574.81274924
logg = 4.6
RV = -54.9390275962
INFO:root:LogL = -59551749.5141

INFO:root:T = 7574.81274924
logg = 4.5
RV = -54.8390275962
INFO:root:LogL = -59842548.3666

INFO:root:T = 7575.28406037
logg = 4.5
RV = -55.0199633331
INFO:root:LogL = -59840191.3304

INFO:root:T = 7575.38406037
logg = 4.5
RV = -55.0199633331
INFO:root:LogL = -59839848.5849

INFO:root:T = 7575.28406037
logg = 4.6
RV = -55.0199633331
INFO:root:LogL = -59550357.5584

INFO:root:T = 7575.28406037
logg = 4.5
RV = -54.9199633331
INFO:root:LogL = -59838701.3517

INFO:root:T = 7575.56733134
logg = 4.5
RV = -55.0686079403
INFO:root:LogL = -59839354.4213

INFO:root:T = 7575.66733134
logg = 4.5
RV = -55.0686079403
INFO:root:LogL = -59839389.6891

INFO:root:T = 7575.56733134
logg = 4.6
RV = -55.0686079403
INFO:root:LogL = -59549502.9739

INFO:root:T = 7575.56733134
logg = 4.5
RV = -54.9686079403
INFO:root:LogL = -59837885.9464

INFO:root:T = 7575.53440174
logg = 4.5
RV = -55.0629531164
INFO:root:LogL = -59839636.4447

INFO:root:T = 7575.63440174
logg = 4.5
RV = -55.0629531164
INFO:root:LogL = -59839306.9812

INFO:root:T = 7575.53440174
logg = 4.6
RV = -55.0629531164
INFO:root:LogL = -59549657.6695

INFO:root:T = 7575.53440174
logg = 4.5
RV = -54.9629531164
INFO:root:LogL = -59837772.5162

INFO:root:T = 7575.56574922
logg = 4.5
RV = -55.0683362519
INFO:root:LogL = -59839346.4466

INFO:root:T = 7575.66574922
logg = 4.5
RV = -55.0683362519
INFO:root:LogL = -59839393.8583

INFO:root:T = 7575.56574922
logg = 4.6
RV = -55.0683362519
INFO:root:LogL = -59549498.4987

INFO:root:T = 7575.56574922
logg = 4.5
RV = -54.9683362519
INFO:root:LogL = -59837886.557

INFO:root:T = 7575.56429512
logg = 4.5
RV = -55.0680865473
INFO:root:LogL = -59839337.0881

INFO:root:T = 7575.66429512
logg = 4.5
RV = -55.0680865473
INFO:root:LogL = -59839386.7387

INFO:root:T = 7575.56429512
logg = 4.6
RV = -55.0680865473
INFO:root:LogL = -59549493.0285

INFO:root:T = 7575.56429512
logg = 4.5
RV = -54.9680865473
INFO:root:LogL = -59837885.1895

INFO:root:T = 7575.54934843
logg = 4.5
RV = -55.0655198318
INFO:root:LogL = -59839280.1845

INFO:root:T = 7575.64934843
logg = 4.5
RV = -55.0655198318
INFO:root:LogL = -59839305.3753

INFO:root:T = 7575.54934843
logg = 4.6
RV = -55.0655198318
INFO:root:LogL = -59549533.8936

INFO:root:T = 7575.54934843
logg = 4.5
RV = -54.9655198318
INFO:root:LogL = -59837828.5287

INFO:root:T = 7575.54412536
logg = 4.5
RV = -55.0646229015
INFO:root:LogL = -59839267.0647

INFO:root:T = 7575.64412536
logg = 4.5
RV = -55.0646229015
INFO:root:LogL = -59839298.8919

INFO:root:T = 7575.54412536
logg = 4.6
RV = -55.0646229015
INFO:root:LogL = -59549590.1767

INFO:root:T = 7575.54412536
logg = 4.5
RV = -54.9646229015
INFO:root:LogL = -59837808.909

INFO:root:T = 7575.54400917
logg = 4.5
RV = -55.0646029504
INFO:root:LogL = -59839266.8689

INFO:root:T = 7575.64400917
logg = 4.5
RV = -55.0646029504
INFO:root:LogL = -59839299.1763

INFO:root:T = 7575.54400917
logg = 4.6
RV = -55.0646029504
INFO:root:LogL = -59549591.7532

INFO:root:T = 7575.54400917
logg = 4.5
RV = -54.9646029504
INFO:root:LogL = -59837808.3658

INFO:root:T = 7575.54389555
logg = 4.5
RV = -55.0645834391
INFO:root:LogL = -59839266.2155

INFO:root:T = 7575.64389555
logg = 4.5
RV = -55.0645834391
INFO:root:LogL = -59839299.3967

INFO:root:T = 7575.54389555
logg = 4.6
RV = -55.0645834391
INFO:root:LogL = -59549593.4346

INFO:root:T = 7575.54389555
logg = 4.5
RV = -54.9645834391
INFO:root:LogL = -59837808.142

INFO:root:T = 7575.53914865
logg = 4.5
RV = -55.0637682777
INFO:root:LogL = -59839573.7425

INFO:root:T = 7575.63914865
logg = 4.5
RV = -55.0637682777
INFO:root:LogL = -59839296.1946

INFO:root:T = 7575.53914865
logg = 4.6
RV = -55.0637682777
INFO:root:LogL = -59549637.9484

INFO:root:T = 7575.53914865
logg = 4.5
RV = -54.9637682777
INFO:root:LogL = -59837792.6192

INFO:root:T = 7575.54389555
logg = 4.5
RV = -55.0645834391
INFO:root:LogL = -59839266.2155

INFO:root:T = 7575.64389555
logg = 4.5
RV = -55.0645834391
INFO:root:LogL = -59839299.3967

INFO:root:T = 7575.54389555
logg = 4.6
RV = -55.0645834391
INFO:root:LogL = -59549593.4346

INFO:root:T = 7575.54389555
logg = 4.5
RV = -54.9645834391
INFO:root:LogL = -59837808.142

INFO:root:Flattening the spectrum using the best-fit values
INFO:root:Outputting flattened spectrum to file ../20130107/HIP_16147_flattened.fits










    















    




















    































In [102]:



    


%matplotlib notebook
for order in flattened:
    plt.plot(order.x, order.y, 'k-', alpha=0.5)



    


















    




















    




























    


    

        

            

            

                Content source: kgullikson88/TS23-Scripts
            

            
Similar notebooks:

            

            
                notebook.community |
                gallery |
                about