https://pfspipe.ipmu.jp/jira This file is an XML representation of an issue en-us 8.3.4 803005 13-09-2019 https://pfspipe.ipmu.jp/jira/browse/PIPE2D-648 DRP 2-D Pipeline <p>In order to do spectral fitting of uncalibrated spectra, in advance, the spectral whitening applies to spectra. We try a method in which Fourier transformation to a spectrum is performed and then a bandpass filter is applied to it. We use AMBRE templates and observed SDSS stellar spectra for this study. </p> <p>This ticket includes the application method to those spectra. Evaluation of the results and the fine tuning of a width of a bandpass filter will be covered in a ticket for spectral fitting.</p> PIPE2D-648

Spectral whitening of AMBRE templates and SDSS spectra

Story Normal Done Done Takuji Yamashita Takuji Yamashita flux-calibration Wed, 4 Nov 2020 05:33:54 +0000 Mon, 4 Jan 2021 20:23:22 +0000 Tue, 15 Dec 2020 06:44:54 +0000 0 3 <p>In order to subtract the continuum from a spectrum and to get a whitened spectrum, I did the spectral whitening for a test spectrum. The spectral whitening is performed in the following procedure:</p> <ul> <li>Fourier transformation of a spectrum</li> <li>Apply a bandpass filter in the wavenumber domain to extract only the large wavenumber (small-scale) components</li> <li>Inverse Fourier transform of the filtered spectrum, which should have been whitened.</li> </ul> <p>I used a simulated spectrum similar to the PFS outputs made from the observed SDSS spectrum of an F-star. I used two high-pass filters, the boxcar and Tukey. The results were check, particularly, at four absorption lines which are strongly depending on the stellar parameters. In any cases with the boxcar filter, the small undulations (sinc functions) are left on the continuum, while those undulations subside around the monitored lines in the case with Tukey filter. It is difficult that both Balmer lines and other narrow lines survive at the same time, because the large widths of Balmer lines are close to the middle-scale fluctuations of the continuum which we want to subtract. I found the resultant spectra of the spectral whitening is very sensitive to a parameter of an adopted filter (e.g., the cut-out wavenumber and the shape parameter). The details are as follows.</p> <h3><a name="Testspectrum"></a><b>Test spectrum</b></h3> <p>I am using similar spectra to the PFS outputs. To make simulated spectra, I multiplied the observed SDSS spectra of F-stars in Lee et al. (2011) by the atmospheric transmission and then divided it by the flat spectrum (from the quartz lamp). </p> <h3><a name="Parameterdependenceofabsorptionlines"></a><b>Parameter dependence of absorption lines</b></h3> <p>Before doing the spectral whitening, I checked the absorption lines that are strongly depending on the stellar parameters (Teff, log g, Z, and <span class="error">&#91;alpha/Fe&#93;</span>). Because these lines are important in stellar classification, I monitor these lines and the continua by them during the whitening process. Four figures below are the standard deviations (blue) and differences between max and min (red) of each wavelength among a subset of AMBRE templates. The left figure is made among a subset in which three parameters (log, Z, and <span class="error">&#91;alpha/Fe&#93;</span>) are fixed but Teff is free (from 5500 K to 8000 K). Similarly, other figures are free in each parameters (logg, Z, and <span class="error">&#91;alpha/Fe&#93;</span>). Balmer series are sensitive to Teff. Na D/He I and Paschen series are sensitive to Teff, Logg, and Z. Mg I triplets are to Logg, Z, and <span class="error">&#91;alpha/Fe&#93;</span>. Ca II triplets are sensitive to all the parameters.</p> <p><span class="image-wrap" style=""><a id="13234_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13234/13234_highpass_boxcar200.png" title="highpass_boxcar200.png" file-preview-type="image" file-preview-id="13234" file-preview-title="highpass_boxcar200.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13228/_thumb_13228.png" style="border: 0px solid black" /></a></span> <span class="image-wrap" style=""><a id="13233_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13233/13233_highpass_boxcar300.png" title="highpass_boxcar300.png" file-preview-type="image" file-preview-id="13233" file-preview-title="highpass_boxcar300.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13229/_thumb_13229.png" style="border: 0px solid black" /></a></span> <span class="image-wrap" style=""><a id="13232_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13232/13232_highpass_boxcar400.png" title="highpass_boxcar400.png" file-preview-type="image" file-preview-id="13232" file-preview-title="highpass_boxcar400.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13227/_thumb_13227.png" style="border: 0px solid black" /></a></span> <span class="image-wrap" style=""><a id="13232_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13232/13232_highpass_boxcar400.png" title="highpass_boxcar400.png" file-preview-type="image" file-preview-id="13232" file-preview-title="highpass_boxcar400.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13230/_thumb_13230.png" style="border: 0px solid black" /></a></span></p> <h3><a name="Spectralwhitening"></a><b>Spectral whitening</b></h3> <p>I tested two filters of the boxcar filter and Tukey filter (<a href="https://scipy.github.io/devdocs/generated/scipy.signal.windows.tukey.html" class="external-link" rel="nofollow">https://scipy.github.io/devdocs/generated/scipy.signal.windows.tukey.html</a>). Each figure shows the input spectrum (top panel), the power spectrum and filtered one (second row from the top), the whitened spectrum from the inverse Fourier transform of the filtered spectrum (third row), and four absorption lines I am monitoring (bottom). Different figures are results with different filters and different parameters of a filter. In the bottom panel, the spectra whose local continua were subtracted using a sigma-clipped median value.</p> <h4><a name="Boxcarfilter"></a><b>Boxcar filter</b></h4> <p>In these four cases with the boxcar filter, the different cut-out wavenumber are adopted: 1/log10(wavelength/AA) = 200, 300, 400, 500, which are close values to a wavenumber roughly comparable to or slightly less than the scale of Balmer lines.</p> <p><span class="image-wrap" style=""><a id="13235_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13235/13235_highpass_boxcar20.png" title="highpass_boxcar20.png" file-preview-type="image" file-preview-id="13235" file-preview-title="highpass_boxcar20.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13234/_thumb_13234.png" style="border: 0px solid black" /></a></span> <span class="image-wrap" style=""><a id="13236_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13236/13236_highpass_tukey_a0.02_cut0.0.png" title="highpass_tukey_a0.02_cut0.0.png" file-preview-type="image" file-preview-id="13236" file-preview-title="highpass_tukey_a0.02_cut0.0.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13233/_thumb_13233.png" style="border: 0px solid black" /></a></span> <span class="image-wrap" style=""><a id="13237_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13237/13237_highpass_tukey_a0.1_cut0.0.png" title="highpass_tukey_a0.1_cut0.0.png" file-preview-type="image" file-preview-id="13237" file-preview-title="highpass_tukey_a0.1_cut0.0.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13232/_thumb_13232.png" style="border: 0px solid black" /></a></span> <span class="image-wrap" style=""><a id="13236_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13236/13236_highpass_tukey_a0.02_cut0.0.png" title="highpass_tukey_a0.02_cut0.0.png" file-preview-type="image" file-preview-id="13236" file-preview-title="highpass_tukey_a0.02_cut0.0.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13231/_thumb_13231.png" style="border: 0px solid black" /></a></span></p> <p>I also show the case that the cut-out wavenumber of 20 is used to make the Balmer lines survive.</p> <p>  <span class="image-wrap" style=""><a id="13237_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13237/13237_highpass_tukey_a0.1_cut0.0.png" title="highpass_tukey_a0.1_cut0.0.png" file-preview-type="image" file-preview-id="13237" file-preview-title="highpass_tukey_a0.1_cut0.0.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13235/_thumb_13235.png" style="border: 0px solid black" /></a></span></p> <h3><a name="Tukeyfilter"></a><b>Tukey filter</b></h3> <p>I tried to delete the small undulations of sinc functions in the boxcar cases and I used a Tukey filter for that. We can control the width and the sharpness of the filter. The undulations subsided around the monitored lines if we choose a proper parameters of the filter. The high sharpness of the filter makes the spectrum flat, but that makes the lines narrow and removes Balmer lines. On the other hand, in the case of the low sharpness, Balmer lines survive but the continuum still remains. <br/>  <br/>    <span class="image-wrap" style=""><a id="13400_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13400/13400_ccf_t7000_4_0_0_v200.png" title="ccf_t7000_4_0_0_v200.png" file-preview-type="image" file-preview-id="13400" file-preview-title="ccf_t7000_4_0_0_v200.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13236/_thumb_13236.png" style="border: 0px solid black" /></a></span><br/>    <span class="image-wrap" style=""><a id="13531_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/13531/13531_deltav_450s_v200.png" title="deltav_450s_v200.png" file-preview-type="image" file-preview-id="13531" file-preview-title="deltav_450s_v200.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/13237/_thumb_13237.png" style="border: 0px solid black" /></a></span></p> <p> </p> <p>We also try the median filter to whiten a spectrum, because in this spectral whitening it is hard to make both Balmer lines (broad lines) and other narrow lines survive at the same time. </p> <p>We can close this ticket following the 2DDRP technical telecon. I will try the median filter as a alternative. </p> Blocks PIPE2D-651 PIPE2D-647 Development Epic Link PIPE2D-303 Rank 0|02qps2:s Reviewers hassan Sprint 2DDRP-2021 A