Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

option for a stellar mask #3

Open
mzechmeister opened this issue Mar 20, 2018 · 1 comment
Open

option for a stellar mask #3

mzechmeister opened this issue Mar 20, 2018 · 1 comment
Assignees

Comments

@mzechmeister
Copy link
Owner

No description provided.

@mzechmeister
Copy link
Owner Author

A stellar mask seems now feasible, since template masking is now elaborated.
The badT flag is propagated (flagbroad).

For now one can modify the QMAP extension in the serval template fits files, which is then loaded in

return 1*hdu['WAVE'].data, hdu['SPEC'].data, (hdu['QMAP'].data if 'QMAP' in hdu else None), hdu[0].header

ww, ff, qq, head = read_template(tpl+(os.sep+os.path.basename(tpl.rstrip(os.sep))+'.fits' if os.path.isdir(tpl) else ''))

And finally masked in:
self.msk_fast = interp(self.wk, 1.*(BK>0))

More convenient would be an -mskstar option. Ups, there is already a -mskspec option; but it does nothing yet

argopt('-mskspec', help='Ascii 0-1 spectrum.'+default, default='')

The mask may need a shift to the frame of the template (RVabs, maybe using the keyword HIERARCH SERVAL TARG RV.

The mask could by applied to the bk (bad knot) map outside or inside the Tpl instance creation

def __init__(self, wk, fk, initfunc, evalfunc, bk=None, mask=False, berv=None, vsini=None, R=None, vrange=None):

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant