Introduction

This simulation should give an answer to the question how the quality of the reconstruction depends on the magnitude of a calibrator. The difficulties arising from a non constant sky background due to an object halo, or an imperfect AO are ignored. The setup used for this experiment is (for the common test setup see table 1):

This experiment uses monochromatic PSFs for the J- and K-Band with a strehl of 0.3. In contrast to all other simulations the calibrator and target uses the same PSF. Therefore only the brightness of the target and calibrator has an influence on the deconvolution. The setup for the PSF generation is described in table 1.

ParameterDescription
Atmosphere/AOThe AO delivers a strehl of about 30 percent.
ObservationIn order to test the sensitivity to the calibrator brightness, magnitudes from 14mag to 22 mag in 1mag steps are used.
Table 1: Setup for the simulation in J-Band and K-Band, .

The experiment uses a target which consists of a scaled down image of NGC4151 and an overlayed image of a simulated dust torus. A second object shows a small star cluster where different distances and brightness levels do exist. These objects are observed (we assume this) in J-Band and K-band five times at different position angles. The image contains exactly one star which can be used as a calibrator. It was additionally assumed, that the background around this star could be perfectly compensated (no inhomogenious sky background or target halo).

The raw data were generated according to the common scheme described in section Overview. In addition, the images for the target and calibrator are generated separately which means, that the psf-star image is not influenced by the target (no halo, etc.).

All simulated input data for the LN DRS pipeline are available as tar-files (ex1_j_input.tar.gz (16MB) and ex1_k_input.tar.gz (18MB)). The corresponding results are also available as tar files (ex1_j_results.tar.gz (4.8MB) and ex1_k_results.tar.gz (5.2MB)).

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.

Calibrator

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 (J-Band) and figure 2 (K-Band) the psf-stars for some stellar brightness levels are shown.

J-Band, calibrator, 14magJ-Band, calibrator, 15magJ-Band, calibrator, 16mag
J-Band, calibrator, 17magJ-Band, calibrator, 18magJ-Band, calibrator, 19mag
J-Band, calibrator, 20magJ-Band, calibrator, 21magJ-Band, calibrator, 22mag
Figure 1: Central 128x128 pixels of the generated J-Band images of a calibrator at different magnitudes (from top left to bottomright: 14, 15, 16, 17, 18, 19, 20, 21, and 22mag).
K-Band, calibrator, 14magK-Band, calibrator, 15magK-Band, calibrator, 16mag
K-Band, calibrator, 17magK-Band, calibrator, 18magK-Band, calibrator, 19mag
K-Band, calibrator, 20magK-Band, calibrator, 21magK-Band, calibrator, 22mag
Figure 2: Central 128x128 pixels of the generated K-Band images of a calibrator at different magnitudes (from top left to bottomright: 14, 15, 16, 17, 18, 19, 20, 21, and 22mag).

Raw Frames and Results for an AGN

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, NGC4151, position angle 108 degreeJ-Band, NGC4151, position angle 144 degreeJ-Band, NGC4151, position angle 180 degreeJ-Band, NGC4151, position angle 216 degreeJ-Band, NGC4151, position angle 252 degree
K-Band, NGC4151, position angle 108 degreeK-Band, NGC4151, position angle 144 degreeK-Band, NGC4151, position angle 180 degreeK-Band, NGC4151, position angle 216 degreeK-Band, NGC4151, position angle 252 degree
Figure 3: The simulated raw images for position angles of 108, 144, 180, 216, and 252 degree of NGC4151 at 10.1mag including sky background (top row J-Band, bottom row K-Band).

A comparison of the reconstructions depending on the calibrator brightness is presented in figure 4 (J-Band) and figure 5 (K-Band).

J-Band reference imageJ-Band coadded raw image
J-Band, Richardson-Lucy reconstruction, 14 mag calibratorJ-Band, Richardson-Lucy reconstruction, 18 mag calibratorJ-Band, Richardson-Lucy reconstruction, 22 mag calibrator
J-Band, Building-Blocks reconstruction, 14 mag calibratorJ-Band, Building-Blocks reconstruction, 18 mag calibratorJ-Band, Building-Blocks reconstruction, 22 mag calibrator
Figure 4: AGN: The J-Band reconstructions depending ion the calibrator brightness are shown in this figure. The top left image shown the reference image of NGC4151. The top right image shown the five coadded raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The calibrator brightness was 14 mag, 18 mag, and 22 mag (from left to right). In the third row, the Building-Block algorithm was used.
K-Band reference imageK-Band coadded raw image
K-Band, Richardson-Lucy reconstruction, 14 mag calibratorK-Band, Richardson-Lucy reconstruction, 18 mag calibratorK-Band, Richardson-Lucy reconstruction, 22 mag calibrator
K-Band, Building-Blocks reconstruction, 14 mag calibratorK-Band, Building-Blocks reconstruction, 18 mag calibratorK-Band, Building-Blocks reconstruction, 22 mag calibrator
Figure 5: AGN: The K-Band reconstructions depending ion the calibrator brightness are shown in this figure. The top left image shown the reference image of NGC4151. The top right image shown the five coadded raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The calibrator brightness was 14 mag, 18 mag, and 22 mag (from left to right). In the third row, the Building-Block algorithm was used.

In table 2 (J-Band) and table 3 (K-Band) the errors for the AGN depending on the psf star brightness are presented. In figure 6 the image errors depending on the calibrator brightness are shown.

TestcaseRichardson-LucyBuilding-Block
Mag (C)Image ErrorImage Error
14.00.011 (0.011)0.004 (0.004)
15.00.024 (0.024)0.005 (0.005)
16.00.055 (0.052)0.006 (0.006)
17.00.093 (0.081)0.010 (0.010)
18.00.126 (0.111)0.025 (0.025)
19.00.160 (0.144)0.072 (0.072)
20.00.190 (0.179)0.153 (0.153)
21.00.233 (0.226)0.241 (0.241)
22.00.297 (0.286)0.461 (0.461)
Table 2: Reconstruction errors for the AGN depending on the psf star brightness (J-Band).
TestcaseRichardson-LucyBuilding-Block
Mag (C)Image ErrorImage Error
14.00.033 (0.019)0.004 (0.004)
15.00.056 (0.054)0.007 (0.007)
16.00.117 (0.107)0.015 (0.015)
17.00.173 (0.165)0.048 (0.048)
18.00.211 (0.207)0.115 (0.115)
19.00.246 (0.244)0.219 (0.219)
20.00.302 (0.290)0.305 (0.305)
21.00.361 (0.349)0.652 (0.652)
22.00.519 (0.502)0.825 (0.825)
Table 3: Reconstruction errors for the AGN depending on the psf star brightness (K-Band).
Image errors depending on the calibrator brightness
Figure 6: In this figure, the image errors depending on the calibrator brightness are shown.

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 calibratior brightness on the reconstruction error and the optimal iteration number is clearly visible.

Reconstruction error using Richardson-Lucy (J-Band)Reconstruction error using Richardson-Lucy (K-Band)
Reconstruction error using Building-Block (J-Band)Reconstruction error using Building-Block (K-Band)
Figure 7: Reconstruction error depending on the iteration number for the J-Band (left side) and the K-Band (right side). The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.

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.

Profile through the reconstructed galaxy using Richardson-Lucy (J-Band)Central part of the profile through the reconstructed galaxy using Richardson-Lucy (J-Band)
Profile through the reconstructed galaxy using Building-Block (J-Band)Central part of the profile through the reconstructed galaxy using Building-Block (J-Band)
Figure 8: Horizontal profile through the center of the reconstructed galaxy compared to the ideal image (J-Band). On the right side only the central part is shown. The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.
Profile through the reconstructed galaxy using Richardson-Lucy (K-Band)Central part of the profile through the reconstructed galaxy using Richardson-Lucy (K-Band)
Profile through the reconstructed galaxy using Building-Block (K-Band)Central part of the profile through the reconstructed galaxy using Building-Block (K-Band)
Figure 9: Horizontal profile through the center of the reconstructed galaxy compared to the ideal image (K-Band). On the right side only the central part is shown. The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.

Raw Frames and Results for a Star Cluster

The basis of this simulation is a small star cluster (see figure 2). The simulated raw images are shown in figure 10.

J-Band, Star Cluster, position angle 108 degreeJ-Band, Star Cluster, position angle 144 degreeJ-Band, Star Cluster, position angle 180 degreeJ-Band, Star Cluster, position angle 216 degreeJ-Band, Star Cluster, position angle 252 degree
K-Band, Star Cluster, position angle 108 degreeK-Band, Star Cluster, position angle 144 degreeK-Band, Star Cluster, position angle 180 degreeK-Band, Star Cluster, position angle 216 degreeK-Band, Star Cluster, position angle 252 degree
Figure 10: The simulated raw images for position angles of 108, 144, 180, 216, and 252 degree of a star cluster including sky background (top row J-Band, bottom row K-Band).

A comparison of the reconstructions depending on the calibrator brightness is presented in figure 4 (J-Band) and figure 5 (K-Band).

J-Band reference imageJ-Band coadded raw image
J-Band, Richardson-Lucy reconstruction, 14 mag calibratorJ-Band, Richardson-Lucy reconstruction, 18 mag calibratorJ-Band, Richardson-Lucy reconstruction, 22 mag calibrator
J-Band, Building-Blocks reconstruction, 14 mag calibratorJ-Band, Building-Blocks reconstruction, 18 mag calibratorJ-Band, Building-Blocks reconstruction, 22 mag calibrator
Figure 11: Star cluster: The J-Band reconstructions depending ion the calibrator brightness are shown in this figure. The top left image shown the reference image of the star cluster. The top right image shown the five coadded raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The calibrator brightness was 14 mag, 18 mag, and 22 mag (from left to right). In the third row, the Building-Block algorithm was used.
K-Band reference imageK-Band coadded raw image
K-Band, Richardson-Lucy reconstruction, 14 mag calibratorK-Band, Richardson-Lucy reconstruction, 18 mag calibratorK-Band, Richardson-Lucy reconstruction, 22 mag calibrator
K-Band, Building-Blocks reconstruction, 14 mag calibratorK-Band, Building-Blocks reconstruction, 18 mag calibratorK-Band, Building-Blocks reconstruction, 22 mag calibrator
Figure 12: Star cluster: The K-Band reconstructions depending ion the calibrator brightness are shown in this figure. The top left image shown the reference image of the star cluster. The top right image shown the five coadded raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The calibrator brightness was 14 mag, 18 mag, and 22 mag (from left to right). In the third row, the Building-Block algorithm was used.

In table 4 (J-Band) and table 5 (K-Band) the errors for the star cluster depending on the psf star brightness are presented. The image errors are also summarized in figure 13 and the photometric errors in figure 14.

TestcaseRichardson-LucyBuilding-Block
Mag (C)Image ErrorMagnitude ErrorImage ErrorMagnitude Error
14.00.011 (0.011)0.01, 0.01, -0.01, 0.16, 0.030.012 (0.012)0.01, 0.00, 0.09, 0.29, 0.41
15.00.010 (0.010)0.01, 0.00, -0.00, 0.17, 0.040.012 (0.012)0.01, 0.00, 0.09, 0.30, 0.41
16.00.009 (0.009)0.00, -0.00, 0.00, 0.20, 0.100.012 (0.012)0.01, 0.00, 0.08, 0.28, 0.36
17.00.015 (0.015)0.01, -0.02, 0.10, 0.38, 0.490.013 (0.013)0.01, 0.01, 0.07, 0.22, 0.24
18.00.044 (0.043)0.02, -0.04, 0.38, 1.73, 2.690.019 (0.019)0.01, 0.01, 0.01, 0.23, 0.15
19.00.096 (0.095)0.10, -0.08, 1.86, 4.81, 2.670.078 (0.078)0.01, 0.01, 0.13, 0.54, 0.52
20.00.177 (0.176)0.35, -0.08, 1.84, 4.83, 2.640.319 (0.308)0.07, 0.05, 0.38, 0.46, 0.41
21.00.308 (0.308)1.16, -0.10, 1.81, 4.88, 2.600.626 (0.621)0.24, 0.14, 0.82, 4.65, 2.66
22.00.464 (0.463)8.28, 0.07, 1.95, 4.99, 2.570.826 (0.826)0.76, 0.88, -0.34, 4.65, 2.40
Table 4: Reconstruction errors for the star cluster depending on the psf star brightness (J-Band).
TestcaseRichardson-LucyBuilding-Block
Mag (C)Image ErrorMagnitude ErrorImage ErrorMagnitude Error
14.00.050 (0.050)0.05, 0.01, 0.18, 0.38, 0.630.038 (0.038)0.06, 0.01, 0.23, 0.44, 0.76
15.00.048 (0.048)0.04, 0.00, 0.17, 0.41, 0.790.039 (0.039)0.06, 0.02, 0.21, 0.46, 0.78
16.00.050 (0.050)0.04, -0.02, 0.33, 1.04, 2.080.038 (0.038)0.06, 0.01, 0.23, 0.51, 0.93
17.00.076 (0.076)0.09, -0.06, 1.49, 4.03, 2.050.057 (0.057)0.09, 0.03, 0.15, 0.56, 1.24
18.00.143 (0.143)0.46, -0.08, 1.47, 4.04, 2.010.160 (0.158)0.15, 0.07, 0.66, 0.80, 1.72
19.00.275 (0.275)2.17, -0.10, 1.40, 4.08, 1.960.459 (0.450)0.42, 0.43, 0.07, 0.95, 2.04
20.00.358 (0.358)7.47, 0.13, 1.55, 4.15, 1.920.701 (0.700)1.16, 0.89, 0.68, 3.94, 2.02
21.00.575 (0.575)7.20, 0.50, 1.71, 4.25, 1.900.905 (0.905)2.47, 2.21, -0.52, 2.81, 1.42
22.00.688 (0.687)7.22, 1.42, 1.72, 4.30, 1.870.986 (0.945)5.38, 4.13, 1.10, 4.58, 0.83
Table 5: Reconstruction errors for the star cluster depending on the psf star brightness (K-Band).
Image errors depending on the calibrator brightness
Figure 13: In this figure, the image errors depending on the calibrator brightness are shown.
Photometric errors depending on the calibrator brightness using Richardson-Lucy (J-Band)Photometric errors depending on the calibrator brightness using Richardson-Lucy (K-Band)
Photometric errors depending on the calibrator brightness using Building-Block (J-Band)Photometric errors depending on the calibrator brightness using Building-Block (K-Band)
Figure 14: Photometric errors depending on the calibrator bightness in the J-Band (left side) and K-Band (right side). The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.

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 calibratior brightness on the reconstruction error and the optimal iteration number is clearly visible.

Reconstruction error using Richardson-Lucy (J-Band)Reconstruction error using Richardson-Lucy (K-Band)
Reconstruction error using Building-Block (J-Band)Reconstruction error using Building-Block (K-Band)
Figure 15: Reconstruction error depending on the iteration number for the J-Band (left side) and the K-Band (right side). The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.

In figure 16, figure 17, and figure 18 the image and photometric errors depending on the calibrator brightness 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.

Image and photometric errors depending on the iteration number using Richardson-Lucy (J-Band, 14 mag calibrator)Image and photometric errors depending on the iteration number using Richardson-Lucy (K-Band, 14 mag calibrator)
Image and photometric errors depending on the iteration number using Building-Block (J-Band, 14 mag calibrator)Image and photometric errors depending on the iteration number using Building-Block (K-Band, 14 mag calibrator)
Figure 16: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The calibrator brightness is 14 mag. The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.
Image and photometric errors depending on the iteration number using Richardson-Lucy (J-Band, 18 mag calibrator)Image and photometric errors depending on the iteration number using Richardson-Lucy (K-Band, 18 mag calibrator)
Image and photometric errors depending on the iteration number using Building-Block (J-Band, 18 mag calibrator)Image and photometric errors depending on the iteration number using Building-Block (K-Band, 18 mag calibrator)
Figure 17: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The calibrator brightness is 18 mag. The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.
Image and photometric errors depending on the iteration number using Richardson-Lucy (J-Band, 22 mag calibrator)Image and photometric errors depending on the iteration number using Richardson-Lucy (K-Band, 22 mag calibrator)
Image and photometric errors depending on the iteration number using Building-Block (J-Band, 22 mag calibrator)Image and photometric errors depending on the iteration number using Building-Block (K-Band, 22 mag calibrator)
Figure 18: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The calibrator brightness is 22 mag. The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.