1 / 20

MIPS Enhancer

MIPS Enhancer. MIPS_Enhancer Overview MIPS Data Processing and Calibration Workshop April 12, 2004. How to Run MIPS_ENHANCER. Command line options mips_enhancer –f filename.fits (40+ options) Namelist options mips_enhancer –n namelist.in. What have we learned since launch ?.

blake
Download Presentation

MIPS Enhancer

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. MIPS Enhancer MIPS_Enhancer Overview MIPS Data Processing and Calibration Workshop April 12, 2004

  2. How to Run MIPS_ENHANCER • Command line options • mips_enhancer –f filename.fits (40+ options) • Namelist options • mips_enhancer –n namelist.in

  3. What have we learned since launch ? • Need to make MIPS_enhancer more efficient in terms of: • Memory • Speed • Improve method of combing surface bright nesses.

  4. Memory and Speed • Tiling • -x #, -y # • Only operating on tile sub-pixels • Do not hold mosaic in memory • -m • Memory limited to current tile • Write tiles out • Stitch mosaic from tiles • Necessary to change format of mosaic to limit memory needed. • Primary header (history, information from first DCE) • Extensions: surface brightness, uncertainty, pixel flag, overlap coverage)

  5. Conserving Surface Brightness • Improvements in: • Outlier rejection algorithm • Weighting overlapping pixels in mosaic

  6. Outlier Rejection Algorithm • For each sub-pixel there is a stack of overlapping BCD pixels.

  7. Sub-Pixel Outlier Rejection N BCD overlapping pixel values Start Iteration i = 1. • Find the median value: Medianorg • Find Standard Deviation: STDEVorg • Find N differences: |(Medianorg – fluxk)| • Sort differences • Withhold 10% largest differences (minimum 1) • Find new medianand standard deviation of remaining points: Mediani, STDEVi • Test Most deviant point: • (Fluxmost deviant point – Mediani)/STDEVi> Sigma Clipping (SC) • If reject – iterate

  8. Outlier Rejection cont. • Continue iterating until: • Remain points < Minimum Points (SP) • Iteration # > Maximum # iterations (SM) • ((STDEVorg – STDEVnew)/STDEVorg )*100 > Sigma Tolerance (ST)

  9. Variation on Outlier Rejection • Option: SK Test Most deviant point: (Fluxmost deviant point – Mediani)/STDEVi> Sigma Clipping (SC) • Option: SU Instead of using the standard deviation of stack: STDEVi Combine in quadrature: uncertainty of BCD pixel and standard deviation of Median

  10. Obsolete Outlier Rejection options • Quartiles • May retire thresholding (t) • (Fluxi/Median Flux) > t continue rejecting • (Fluxi/Median Flux) < t stop rejection

  11. Conserving Surface Brightness • For each sub-pixel we have N BCD pixels overlapping • Surface Brightness of BCD • Uncertainty associated with surface brightness • Overlap Area • How do we combine information to conserve surface brightness in mosaic ?

  12. Distortions

  13. Weighting Methods Most probable value is a weighted average of points: • = wixi/  wi • W1= overlap area, uncertainties • W2=overlap area • W3=overlap area, S/N (SB/uncertainty) • W4=overlap area, SQRT(S/N) • W5=overlap area, SQRT(S)

  14. Useful Options for MIPS_enhancer • -f filename (-i list of files) • -o output_name • -RR (remove rotation) • -s # (scale factor, recommended = 0.25) • dy (correct for distortion, default) • -x #, -y #, -m (memory & speed options) • -r # (70 microns, # BCDs to reject after stim, recommened 1) • -IB (ignore bias boost, 24 microns, default) • -W2 (W1,W3,W4,W5) Weight options • -SC #, ST #, SP #, SM # SK, SU (outlier rejection options Recommened:SC 4, ST 10, SP 3, SM 30, SK) • -FM (FD) sub-pixel flux weighted mean (median) • -D, -C, -FL (sets other appropriate options)

  15. Near Term Improvements • SED mode • Distortion (verify or tweak distortion coefficients produced from CODE-V) • Study on surface brightness conservation • Add weighting which only uses the number of reads that went into the slope determination. • Develop testing routine to see how well we have Use simulating routine (Emeric) to add sources to cal file. • Move on to 160 micron data • Start enhancing framework

  16. Suggestions ?

More Related