This test investigates the dependency of the reconstruction error on the atmospheric dispersion (refraction). In addition, a very red target (0.5 mag brighter in the red) and a slightly red calibrator (0.2 mag brighter in the red) were used. For this test, the PSFs are generated using several (41) monochromatic PSFs. This simulation only tests the difficult cases where the target has a declination between -30 and +20 degree. For each target, three parallactic angles were used which corresponds to the largest angle which can be achieved going down to 75 degree zenith angle.
The setup of the simulation is described in table 1. The input data for the test uses a strehl value of 0.30 for target and calibrator. The psf star has a magnitude of 14. In addition, the target and psf star PSF are generated by using 11 (no atmospheric dispersion) or 41 (with atmospheric dispersion) monochromatic PSFs.
The multi-monochromatic PSFs are created by generating a monochromatic OPD screen for the shortest wavelength λ0. The OPD values are then rescaled for a given wavelength λ (we used 11 or 41 discrete wavelength to sample a spectral band) by λ0/λ. The result was put into an array which was enlarged by λ/λ0. Then the atmospheric dispersion was simulated by moving the center of the psf according to the assumed zenith angle and a model atmosphere. A PSF for a specific wavelength was then calculated by |FFT-1|2, weighted by the spectral intensity and all summed up to the final multi-monochromatic PSF.
Parameter | Description |
---|---|
Atmosphere/AO | The AO delivers a strehl of about 30 percent for the target and the psf star. The psfs are separately simulated and therefore not equal!. The model of the atmosphere is from Allan's Astrophysical Quantities, 1999, section 11.20 Atmospheric refraction and air path, page 262. |
Spectra | The target and calibrator show a spectra which was sampled at 41 equally spaced wavelength over a part of the J-Band (1.156 - 1.250 μm: J shorter half) and K-Band (2.0 - 2.4 μm). The spectrum of the target is 0.5 mag brighter on the red side and the spectrum of the calibrator is 0.2 mag brighter at the red side. |
Observation | The observation target is NGC4151 with an overlayed dust tori (simulated) at 10.1mag integral brightness and a small star cluster. The brightness of the psf-star is set to 14mag. For each position angle an image equivalent to 60s exposure time was created. The central wavelength was 1.203 micron (shorter side of the J-Band) and 2.2 micron (whole K-band). This target was placed at -30, -20, -10, 0, 10, and 20 degree declination. |
In table 2 all combinations of declination, culmination, smallest altitude, and parallactic angles are shown. In addition, each test case was also created without atmospheric dispersion.
Declination [deg] | Culmination [deg] | Smallest altitude [deg] | PA [deg] |
---|---|---|---|
-30 | 27.3 | 15.6 | 144 180 216 |
-20 | 37.3 | 16.5 | 135 180 225 |
-10 | 47.3 | 15.0 | 128 180 232 |
0 | 57.3 | 15.1 | 124 180 236 |
10 | 67.3 | 28.9 | 122 180 238 |
20 | 77.3 | 55.7 | 120 180 240 |
The raw data were generated according to the common scheme described in section Overview. In addition, the target and psf star PSF are generated by using 41 monochromatic PSFs.
All simulated input data for the LN DRS pipeline are available as a tar-file (ex8_j_input.tar.gz (3.9MB) and ex8_k_input.tar.gz (4.3MB)). The corresponding results are also available as tar files (ex8_j_results.tar.gz (552KB) and ex8_k_results.tar.gz (546KB)).
The generated calibrator raw images are shown in section Calibrator. The raw images and deconvolution results are shown in section Raw Frames and Results for an AGN for the AGN and section Raw Frames and Results for a Star Cluster for the star cluster.
The simulated raw LBT interferograms used for the deconvolution are ideal images, they are not influenced by detector effects like different pixel gain or bad pixels. In figure 1 and figure 2 the raw calibrator for different declination is shown.
The basis of this simulation is an image of NGC4151 with an overlayed dust torus (see figure 1). The simulated raw images are shown in figure 3 (J-Band) and in figure 4 (K-Band).
A comparison of the reconstructions depending on the atmospheric dispersion is presented in figure 5 and figure 6 for the J-Band and K-Band.
In figure 7 the image error describing the difference between the reconstruction and the reference image is shown for each iteration (in steps of 100 iterations). The effect of the atmospheric dispersion on the reconstruction error and the optimal iteration number is clearly visible.
In figure 8 (J-Band) and figure 9 (K-Band) a horizontal cut slightly above the intensity maximum of the reconstructions of NGC4151 compared with the ideal reference image is shown.
In table 3 (J-Band) and table 4 (K-Band) the errors for the AGN depending on the atmospheric dispersion are presented.
Testcase | Richardson-Lucy | Building-Block | ||
---|---|---|---|---|
Dec | PA | Altitudes | Image Error | Image Error |
-30 | 144 180 216 | 15.6 27.3 15.6 | 0.156 (0.155) | 0.120 (0.120) |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.143 (0.142) | 0.103 (0.103) |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.124 (0.123) | 0.088 (0.088) |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.091 (0.091) | 0.059 (0.059) |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.068 (0.068) | 0.047 (0.047) |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.037 (0.037) | 0.021 (0.021) |
Testcase | Richardson-Lucy | Building-Block | ||
---|---|---|---|---|
Dec | PA | Altitudes | Image Error | Image Error |
-30 | 144 180 216 | 15.6 27.3 15.6 | 0.157 (0.157) | 0.127 (0.127) |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.118 (0.117) | 0.079 (0.079) |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.095 (0.094) | 0.046 (0.046) |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.075 (0.075) | 0.041 (0.041) |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.044 (0.044) | 0.022 (0.022) |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.038 (0.038) | 0.017 (0.017) |
The basis of this simulation is a small star cluster (see figure 2). The simulated raw images are shown in figure 10 (J-Band) and in figure 11 (K-Band).
A comparison of the reconstructions depending on the atmospheric dispersion is presented in figure 12 (J-Band) and figure 13 (K-Band).
In figure 14 the image error describing the difference between the reconstruction and the reference image is shown for each iteration (in steps of 100 iterations). The effect of the atmospheric dispersion on the reconstruction error and the optimal iteration number is clearly visible.
In the figures figure 15, figure 16, figure 17, and figure 18 the image and photometric errors depending on the atmospheric dispersion and the iteration number is shown. The overall shape is the same, but the optimum number of iterations is difficult to select, because not all photometric errors have their minimum at the same iteration number.
In table 5 (J-Band) and table 6 (K-Band) the errors for the AGN depending on the atmospheric dispersion are presented.
Testcase | Richardson-Lucy | Building-Block | ||||||
---|---|---|---|---|---|---|---|---|
Dec | PA | Altitudes | Image Error | Magnitude Error | Position Error | Image Error | Magnitude Error | Position Error |
-30 | 144 180 216 | 15.6 27.3 15.6 | 0.100 (0.100) | 0.07, 0.04, 0.56, 0.10, 2.68 | 0.05, 0.04, 0.25, 0.57, 0.59 | 0.110 (0.110) | 0.07, 0.03, 0.85, 0.24, 2.43 | 0.05, 0.05, 0.49, 0.54, 0.60 |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.104 (0.104) | 0.04, 0.06, 0.53, 0.40, 0.98 | 0.06, 0.05, 0.18, 0.07, 0.57 | 0.108 (0.108) | 0.04, 0.05, 0.74, 0.54, 2.20 | 0.05, 0.06, 0.38, 0.01, 0.84 |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.070 (0.070) | 0.04, 0.04, 0.43, 0.40, 0.36 | 0.03, 0.02, 0.62, 0.05, 1.23 | 0.070 (0.070) | 0.03, 0.03, 0.70, 0.52, 0.96 | 0.03, 0.03, 0.36, 0.01, 1.17 |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.049 (0.049) | 0.02, 0.03, 0.36, 0.26, 0.86 | 0.00, 0.02, 0.75, 0.60, 0.07 | 0.049 (0.049) | 0.02, 0.01, 0.59, 0.40, 2.21 | 0.00, 0.03, 0.32, 0.42, 0.50 |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.072 (0.072) | 0.04, 0.03, 0.05, -0.02, 0.39 | 0.04, 0.04, 0.67, 0.36, 1.11 | 0.068 (0.068) | 0.04, 0.02, 0.15, 0.01, 0.68 | 0.04, 0.03, 0.60, 0.25, 0.92 |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.030 (0.030) | 0.02, 0.01, 0.10, 0.24, 0.58 | 0.03, 0.03, 0.06, 0.14, 0.51 | 0.026 (0.026) | 0.02, 0.00, 0.19, 0.34, 1.15 | 0.01, 0.03, 0.07, 0.01, 0.47 |
Testcase | Richardson-Lucy | Building-Block | ||||||
---|---|---|---|---|---|---|---|---|
Dec | PA | Altitudes | Image Error | Magnitude Error | Position Error | Image Error | Magnitude Error | Position Error |
-30 | 144 180 216 | 15.6 27.3 15.6 | 0.167 (0.167) | 0.12, 0.03, 0.90, 0.55, 2.08 | 0.29, 0.15, 1.82, 1.19, 0.36 | 0.150 (0.150) | 0.13, 0.04, 0.81, 0.63, 1.94 | 0.27, 0.14, 1.62, 1.30, 0.36 |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.093 (0.093) | 0.03, 0.03, 0.88, 0.65, 2.07 | 0.08, 0.06, 0.75, 0.59, 0.39 | 0.065 (0.065) | 0.04, 0.03, 0.80, 0.76, 1.92 | 0.04, 0.05, 0.71, 0.48, 0.44 |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.070 (0.070) | 0.05, 0.03, 1.17, 0.23, 1.90 | 0.10, 0.05, 0.59, 1.39, 0.94 | 0.052 (0.052) | 0.05, 0.04, 1.06, 0.33, 1.51 | 0.09, 0.05, 0.24, 1.27, 1.72 |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.075 (0.075) | 0.06, 0.03, 0.17, 0.09, 2.01 | 0.08, 0.03, 0.24, 1.29, 0.42 | 0.061 (0.061) | 0.07, 0.05, 0.19, 0.13, 1.88 | 0.08, 0.01, 0.19, 1.36, 0.45 |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.067 (0.067) | 0.03, 0.02, 0.75, 0.70, 1.69 | 0.07, 0.05, 1.57, 2.26, 1.79 | 0.053 (0.053) | 0.03, 0.03, 0.68, 0.75, 1.75 | 0.06, 0.06, 1.89, 2.32, 1.07 |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.057 (0.057) | 0.05, 0.03, 0.91, 0.98, 1.92 | 0.05, 0.02, 0.71, 0.60, 0.65 | 0.046 (0.046) | 0.05, 0.03, 0.87, 1.07, 1.90 | 0.03, 0.02, 0.61, 0.78, 0.45 |