This test investigates the dependency of the reconstruction error on the atmospheric dispersion (refraction). For this test, the PSFs are generated using several (11 or 41) monochromatic PSFs. For the target and calbrator a constant spectrum was assumed. 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 constant spectra which was sampled at 11 or 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). |
Observation | The observation target is NGC4151 with an overlayed dust tori (simulated) at 10.1mag integral brightness. 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 10 monochromatic PSFs.
All simulated input data for the LN DRS pipeline are available as a tar-file (ex7_j_input.tar.gz (3.9MB) and ex7_k_input.tar.gz (4.3MB)). The corresponding results are also available as tar files (ex7_j_results.tar.gz (552KB) and ex7_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 with dispersion is shown.
The basis of the simulation is an image of NGC4151 with an overlayed dust torus (see figure 1. The simulated raw images for the smallest altitude with atmospheric dispersion and without atmospheric dispersion are shown in figure 3 for the J-Band and in figure 4 for the K-Band.
A comparison of the reconstructions depending on the atmospheric dispersion is presented in figure 5 for the J-Band. For the K-Band the results are presented in figure 6.
In table 3 (J-Band) and table 4 (K-Band) the errors for the star cluster depending on the atmospheric dispersion are presented. In figure 7 the image errors depending on the atmospheric dispersion are shown.
Testcase | Richardson-Lucy | Building-Block | ||
---|---|---|---|---|
Dec | PA | Altitudes | Image Error | Image Error |
-30 | 144 180 216 | 15.6 27.3 15.6 | 0.141 (0.140) | 0.076 (0.076) |
-30 | 144 180 216 | 90.0 90.0 90.0 | 0.047 (0.046) | 0.034 (0.034) |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.122 (0.122) | 0.050 (0.050) |
-20 | 135 180 225 | 90.0 90.0 90.0 | 0.038 (0.038) | 0.024 (0.024) |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.109 (0.108) | 0.040 (0.040) |
-10 | 128 180 232 | 90.0 90.0 90.0 | 0.028 (0.028) | 0.015 (0.015) |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.078 (0.078) | 0.033 (0.033) |
0 | 124 180 236 | 90.0 90.0 90.0 | 0.027 (0.027) | 0.012 (0.012) |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.058 (0.058) | 0.023 (0.023) |
10 | 122 180 238 | 90.0 90.0 90.0 | 0.029 (0.029) | 0.010 (0.010) |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.034 (0.034) | 0.015 (0.015) |
20 | 120 180 240 | 90.0 90.0 90.0 | 0.029 (0.029) | 0.012 (0.012) |
Testcase | Richardson-Lucy | Building-Block | ||
---|---|---|---|---|
Dec | PA | Altitudes | Image Error | Image Error |
-30 | 144 180 216 | 15.6 27.3 15.6 | 0.143 (0.139) | 0.092 (0.092) |
-30 | 144 180 216 | 90.0 90.0 90.0 | 0.095 (0.093) | 0.064 (0.064) |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.110 (0.107) | 0.069 (0.069) |
-20 | 135 180 225 | 90.0 90.0 90.0 | 0.071 (0.070) | 0.036 (0.036) |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.092 (0.090) | 0.040 (0.040) |
-10 | 128 180 232 | 90.0 90.0 90.0 | 0.051 (0.051) | 0.017 (0.017) |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.073 (0.073) | 0.026 (0.026) |
0 | 124 180 236 | 90.0 90.0 90.0 | 0.039 (0.039) | 0.017 (0.017) |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.042 (0.042) | 0.016 (0.016) |
10 | 122 180 238 | 90.0 90.0 90.0 | 0.037 (0.037) | 0.017 (0.017) |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.032 (0.032) | 0.017 (0.017) |
20 | 120 180 240 | 90.0 90.0 90.0 | 0.034 (0.034) | 0.018 (0.018) |
In figure 8 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 9 (J-Band) and figure 10 (K-Band) a horizontal cut slightly above the intensity maximum of the reconstructions of NHC4151 compared with the ideal reference image is shown.
The basis of the simulation is an image of a starcluster (see figure 2. The simulated raw images for the smallest altitude with atmospheric dispersion and without atmospheric dispersion are shown in figure 11 for the J-Band and in figure 12 for the K-Band.
A comparison of the reconstructions depending on the atmospheric dispersion is presented in figure 13 for the J-Band. For the K-Band the results are presented in figure 14.
In figure 15 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 16, figure 17, figure 18, and figure 19 the image and photometric errors depending on the atmospheric dispersion and the iteration number is shown.
In table 5 (J-Band) and table 6 (K-Band) the errors for the star cluster 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.059 (0.059) | 0.02, 0.01, 0.42, -0.01, 0.23 | 0.02, 0.02, 0.14, 1.05, 0.09 | 0.049 (0.048) | 0.02, 0.01, 0.50, 0.05, 0.42 | 0.02, 0.02, 0.42, 1.00, 0.04 |
-30 | 144 180 216 | 90.0 90.0 90.0 | 0.018 (0.018) | 0.02, 0.01, -0.08, 0.11, 0.53 | 0.01, 0.01, 0.47, 0.09, 0.08 | 0.016 (0.016) | 0.02, -0.00, 0.05, 0.23, 1.01 | 0.01, 0.01, 0.26, 0.02, 0.12 |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.047 (0.047) | 0.04, 0.03, 0.00, 0.25, 0.69 | 0.00, 0.02, 0.41, 0.59, 0.42 | 0.043 (0.043) | 0.03, 0.02, 0.05, 0.28, 0.91 | 0.00, 0.02, 0.41, 0.63, 0.20 |
-20 | 135 180 225 | 90.0 90.0 90.0 | 0.030 (0.030) | 0.04, 0.01, 0.16, 0.13, 0.89 | 0.01, 0.01, 0.12, 0.03, 0.05 | 0.021 (0.021) | 0.02, 0.01, 0.18, 0.14, 1.07 | 0.01, 0.01, 0.08, 0.01, 0.12 |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.026 (0.026) | 0.02, 0.00, 0.10, 0.62, 1.72 | 0.00, 0.01, 0.38, 0.27, 1.43 | 0.025 (0.025) | 0.02, 0.00, 0.15, 0.67, 2.42 | 0.00, 0.01, 0.31, 0.23, 0.68 |
-10 | 128 180 232 | 90.0 90.0 90.0 | 0.020 (0.020) | 0.02, -0.00, 0.12, 0.01, 0.61 | 0.00, 0.00, 0.18, 0.50, 0.27 | 0.018 (0.018) | 0.02, -0.01, 0.22, 0.13, 1.15 | 0.00, 0.00, 0.11, 0.38, 0.14 |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.030 (0.030) | -0.01, 0.01, 0.12, 0.15, 0.73 | 0.03, 0.02, 0.56, 0.08, 0.34 | 0.028 (0.028) | -0.00, 0.01, 0.21, 0.21, 1.33 | 0.03, 0.02, 0.41, 0.03, 0.17 |
0 | 124 180 236 | 90.0 90.0 90.0 | 0.031 (0.031) | 0.03, 0.02, 0.07, 0.10, 0.24 | 0.01, 0.00, 0.19, 0.25, 1.75 | 0.025 (0.025) | 0.02, 0.01, 0.09, 0.15, 0.50 | 0.01, 0.00, 0.17, 0.25, 1.85 |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.030 (0.030) | 0.03, 0.02, -0.05, 0.10, 0.47 | 0.00, 0.02, 0.56, 0.06, 0.44 | 0.026 (0.026) | 0.02, 0.01, -0.02, 0.15, 0.88 | 0.00, 0.02, 0.48, 0.02, 0.10 |
10 | 122 180 238 | 90.0 90.0 90.0 | 0.023 (0.023) | 0.02, 0.02, 0.10, 0.18, 0.18 | 0.00, 0.00, 0.04, 0.15, 0.11 | 0.019 (0.019) | 0.02, 0.01, 0.15, 0.26, 0.37 | 0.00, 0.00, 0.07, 0.03, 0.04 |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.019 (0.019) | 0.02, 0.01, 0.09, 0.24, 0.36 | 0.01, 0.00, 0.12, 0.13, 0.31 | 0.020 (0.020) | 0.02, 0.00, 0.18, 0.37, 0.87 | 0.00, 0.00, 0.08, 0.04, 0.14 |
20 | 120 180 240 | 90.0 90.0 90.0 | 0.026 (0.026) | 0.03, 0.02, 0.12, 0.17, 0.21 | 0.01, 0.00, 0.25, 0.47, 0.53 | 0.020 (0.020) | 0.01, 0.01, 0.16, 0.22, 0.42 | 0.00, 0.00, 0.16, 0.40, 0.09 |
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.126 (0.126) | 0.05, 0.03, 0.90, 1.41, 0.28 | 0.19, 0.04, 1.13, 2.23, 1.72 | 0.075 (0.075) | 0.05, 0.03, 0.73, 1.60, 0.44 | 0.17, 0.04, 0.59, 2.40, 1.12 |
-30 | 144 180 216 | 90.0 90.0 90.0 | 0.068 (0.068) | 0.05, 0.02, 0.29, 0.58, 2.07 | 0.08, 0.03, 1.77, 0.80, 0.44 | 0.038 (0.038) | 0.05, 0.02, 0.30, 0.65, 1.92 | 0.07, 0.03, 1.90, 0.77, 0.48 |
-20 | 135 180 225 | 16.5 37.3 16.5 | 0.107 (0.107) | 0.06, 0.03, 0.65, 1.05, 1.64 | 0.09, 0.08, 1.38, 2.30, 0.97 | 0.076 (0.076) | 0.06, 0.04, 0.72, 1.09, 1.64 | 0.07, 0.07, 1.00, 2.62, 1.61 |
-20 | 135 180 225 | 90.0 90.0 90.0 | 0.067 (0.067) | 0.03, 0.00, 0.41, 1.21, 2.04 | 0.08, 0.03, 1.13, 1.57, 0.53 | 0.045 (0.045) | 0.04, 0.01, 0.47, 1.33, 1.90 | 0.05, 0.02, 0.98, 1.23, 0.51 |
-10 | 128 180 232 | 15.0 47.3 15.0 | 0.087 (0.087) | 0.06, 0.03, 1.07, 0.31, 1.65 | 0.06, 0.05, 1.86, 0.90, 1.47 | 0.062 (0.062) | 0.07, 0.03, 0.87, 0.42, 1.87 | 0.05, 0.04, 1.97, 1.14, 0.67 |
-10 | 128 180 232 | 90.0 90.0 90.0 | 0.055 (0.055) | 0.03, 0.04, 0.45, 0.46, 0.96 | 0.07, 0.02, 1.01, 0.95, 3.46 | 0.041 (0.041) | 0.03, 0.04, 0.52, 0.54, 0.98 | 0.06, 0.03, 0.82, 0.95, 3.20 |
0 | 124 180 236 | 15.1 57.3 15.1 | 0.074 (0.074) | 0.04, 0.03, 0.43, 0.99, 0.35 | 0.08, 0.03, 1.17, 1.01, 3.83 | 0.054 (0.054) | 0.04, 0.04, 0.39, 1.11, 0.39 | 0.05, 0.01, 1.68, 1.41, 4.33 |
0 | 124 180 236 | 90.0 90.0 90.0 | 0.060 (0.060) | 0.06, -0.00, 0.44, 0.83, -0.23 | 0.04, 0.05, 0.79, 1.20, 1.22 | 0.050 (0.050) | 0.07, -0.00, 0.49, 0.90, -0.09 | 0.03, 0.05, 0.74, 0.81, 1.29 |
10 | 122 180 238 | 28.9 67.3 28.9 | 0.062 (0.062) | 0.04, 0.02, 0.78, 0.71, 2.09 | 0.07, 0.03, 0.30, 1.91, 0.40 | 0.052 (0.052) | 0.04, 0.03, 0.86, 0.83, 1.93 | 0.06, 0.01, 0.36, 1.87, 0.42 |
10 | 122 180 238 | 90.0 90.0 90.0 | 0.056 (0.056) | 0.05, 0.01, 1.07, 0.50, 1.66 | 0.09, 0.02, 1.11, 0.27, 2.00 | 0.046 (0.046) | 0.05, 0.01, 1.02, 0.58, 1.58 | 0.08, 0.02, 0.93, 0.24, 1.92 |
20 | 120 180 240 | 55.7 77.3 55.7 | 0.059 (0.059) | 0.04, 0.03, 0.81, 0.20, 2.10 | 0.05, 0.03, 1.33, 0.41, 0.43 | 0.051 (0.051) | 0.04, 0.03, 0.82, 0.26, 1.94 | 0.03, 0.02, 1.35, 0.65, 0.46 |
20 | 120 180 240 | 90.0 90.0 90.0 | 0.047 (0.047) | 0.00, 0.00, 0.21, 0.94, 2.09 | 0.06, 0.03, 0.25, 0.37, 0.42 | 0.038 (0.038) | 0.02, 0.01, 0.27, 1.10, 1.93 | 0.05, 0.02, 0.44, 0.30, 0.44 |