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-243 DRP 2-D Pipeline <p>Measure the optical aberrations of individual exposures of PFS data from out-of-focus images.</p> <p>The current plan is to follow the procedure by Joshua Mayers which was used for HSC images (e.g., <a href="https://jira.lsstcorp.org/browse/DM-8565" class="external-link" rel="nofollow">https://jira.lsstcorp.org/browse/DM-8565</a> and connected tickets). This means forward modeling of the PFS images taken at LAM with Galsim.</p> <p>I am documenting my progress at <a href="https://github.com/nevencaplar/PFS_Work_In_Progress/blob/master/Zernike.ipynb" class="external-link" rel="nofollow">https://github.com/nevencaplar/PFS_Work_In_Progress/blob/master/Zernike.ipynb.</a></p> <p>Update: July 5, 2018, report: <a href="https://www.overleaf.com/17550215gqrxkdfywqkp#/66627172/" class="external-link" rel="nofollow">https://www.overleaf.com/17550215gqrxkdfywqkp#/66627172/</a></p> PIPE2D-243

Characterize PSF using out-of-focus images

Story Normal Won't Fix Won't Fix ncaplar ncaplar Mon, 25 Sep 2017 16:45:21 +0000 Fri, 5 Apr 2019 15:04:41 +0000 Thu, 4 Apr 2019 20:56:40 +0000 0 2 <p>At the moment, I have been successful in running galsim and creating some mock images. I have also made a mock-up version of the pupil (see Pupil00.png). Jim G. should provide me at some point with illumination files which will give correct proportions of struts etc... At the moment I am struggling to understand if I am properly modeling the system, i.e., to set the all the parameters of the simulation (pupil size, pixel scale....) properly, so that the simulated images look somewhat like the LAM images.</p> <p>I have fully analyzed one of the test images taken at LAM (v3617).</p> <p>Below we can see one image which shows poststamp images of the defocused spots across the detector. The full black squares are the positions that were either too faint, or for which the simulation on the cluster did not finish (because I put too strict time limits to run my job). The numbers above are x and y positions on the chip. I am showing both the linear and log scaling. </p> <p><span class="image-wrap" style=""><a id="10816_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10816/10816_Originals+Log.png" title="Originals Log.png" file-preview-type="image" file-preview-id="10816" file-preview-title="Originals Log.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10819/_thumb_10819.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10812_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10812/10812_Models.png" title="Models.png" file-preview-type="image" file-preview-id="10812" file-preview-title="Models.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10816/_thumb_10816.png" style="border: 0px solid black" /></a></span> </p> <p>Then I model the post stamps. I have been modeling those with Galsim, with 16 parameters (z4 to z11 Zernike coefficients, size of the central obscuration, size of struts, x and y position of the central obscuration in the focal plane, x and y position on the chi, “size of the telescope”, size of the input fiber image). Results are quite good, capturing the main features of the actual data. There is only one obvious failure (at 3620,2347), and even data in corner of the detector is reproduced. Below I am showing results, with linear and log scaling.</p> <p><span class="image-wrap" style=""><a id="10811_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10811/10811_Models+Log.png" title="Models Log.png" file-preview-type="image" file-preview-id="10811" file-preview-title="Models Log.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10812/_thumb_10812.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10814_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10814/10814_Residuals%2C+no+rescale.png" title="Residuals, no rescale.png" file-preview-type="image" file-preview-id="10814" file-preview-title="Residuals, no rescale.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10811/_thumb_10811.png" style="border: 0px solid black" /></a></span> </p> <p>When we look at the residuals in linear scale, some of the images are almost complete erased. Log scale reveals that some structure is still left.</p> <p><span class="image-wrap" style=""><a id="10817_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10817/10817_Pupils.png" title="Pupils.png" file-preview-type="image" file-preview-id="10817" file-preview-title="Pupils.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10814/_thumb_10814.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10817_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10817/10817_Pupils.png" title="Pupils.png" file-preview-type="image" file-preview-id="10817" file-preview-title="Pupils.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10813/_thumb_10813.png" style="border: 0px solid black" /></a></span> </p> <p>One can also look at the “predicted” images of the pupils and their movement across the detector.</p> <p><span class="image-wrap" style=""><a id="10803_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10803/10803_2_z6.png" title="2_z6.png" file-preview-type="image" file-preview-id="10803" file-preview-title="2_z6.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10817/_thumb_10817.png" style="border: 0px solid black" /></a></span> </p> <p>Below I am showing the estimates for the behavior of different Zernike coefficients across the detector (from z5 to z11). Size of the spots is the same for all of the images, the information is only in the color. </p> <p><span class="image-wrap" style=""><a id="10804_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10804/10804_3_z7.png" title="3_z7.png" file-preview-type="image" file-preview-id="10804" file-preview-title="3_z7.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10802/_thumb_10802.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10804_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10804/10804_3_z7.png" title="3_z7.png" file-preview-type="image" file-preview-id="10804" file-preview-title="3_z7.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10803/_thumb_10803.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10805_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10805/10805_4_z8.png" title="4_z8.png" file-preview-type="image" file-preview-id="10805" file-preview-title="4_z8.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10804/_thumb_10804.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10806_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10806/10806_5_z9.png" title="5_z9.png" file-preview-type="image" file-preview-id="10806" file-preview-title="5_z9.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10805/_thumb_10805.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10807_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10807/10807_6_z10.png" title="6_z10.png" file-preview-type="image" file-preview-id="10807" file-preview-title="6_z10.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10806/_thumb_10806.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10808_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10808/10808_7_z11.png" title="7_z11.png" file-preview-type="image" file-preview-id="10808" file-preview-title="7_z11.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10807/_thumb_10807.png" style="border: 0px solid black" /></a></span> <br/> <span class="image-wrap" style=""><a id="10808_thumb" href="https://pfspipe.ipmu.jp/jira/secure/attachment/10808/10808_7_z11.png" title="7_z11.png" file-preview-type="image" file-preview-id="10808" file-preview-title="7_z11.png"><img src="https://pfspipe.ipmu.jp/jira/secure/thumbnail/10808/_thumb_10808.png" style="border: 0px solid black" /></a></span> </p> <p>Future improvements, roughly in the order of priority that I would suggest:<br/> 1. Analysis on better images with larger spots and stronger signal, also use images with different defocus to compare results<br/> 2. Check if the images are uniformly illuminated by comparing data taken from both sides of focus<br/> 3. Eliminate the pupil parameters from the fit and use the data which will be provided from France about the movement of the central obscuration and struts across the detector<br/> 4. I am using “size of the telescope” as a free parameter - obviously this is wrong and I am using it hide some of my lack of knowledge about how the spectrographs work and abusing the code which was done for other purpose. This needs to be better understood (and other parameters in galsim which I am not touching at the moment, e.g., wavelength of the light). <br/> 5. Do a bit more realistic modeling of input fiber (currently it is top hat, realistically is that some wings)<br/> 6. Analyzing with the images with corrected gains (currently variance in the scientific images is different than what would be inferred from variance image)</p> <p>This is an old overarching ticket summarizing the development of the PSF pipelne from defocused images. The work has been split to many other tickets. </p> Blocks PIPE2D-284 PIPE2D-269 Development Epic Link PIPE2D-1 Rank 0|02qpt6:508yjx Reviewers hassan Sprint 2DDRP-2019 D