Introduction

In some observation conditions the psf of the target is not equal to the psf of the calibrator. Several reasons for the discrepancies do exist. For this test case different AO-performances for the target and calibrator are assumed which result in different strehl values. Unequal spectra of the target and the calibrator or other sources for discrepancies are not coverd by this test case.

The setup of the simulation is described in table 1. The input data for the test uses a raw target image with a strehl of 30 percent and the strehl of the psf-star covers a range from 20 percent up to 40 percent. In addition an artificial perfect psf (strehl 100 percent) is used for the reconstruction. In this experiment a bright psf star is chosen (14 mag), so that the reconstructions are not influenced by the photon noise of the calibrator.

ParameterDescription
Atmosphere/AOThe AO delivers a strehl of about 30 percent for the target but for the psf-star strehl values of 20, 25, 27, 29, 30, 31, 33, 35, 40, and 100 percent are used.
ObservationThe brightness of the psf-star is set to 14mag.
Table 1: Setup for the simulation in J-Band and K-Band, for the common test setup see table 1.

The raw data were generated according to the common scheme described in section Overview. In addition the images for the target and psf-star 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 a tar-file (ex2_j_input.tar.gz (18MB) and ex2_k_input.tar.gz (21MB)). The corresponding results are also available as tar files (ex2_j_results.tar.gz (5.5MB) and ex2_k_results.tar.gz (5.7MB)).

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 calibrator for some strehl values are shown.

J-Band, calibrator, strehl 0.20J-Band, calibrator, strehl 0.25J-Band, calibrator, strehl 0.27
J-Band, calibrator, strehl 0.29J-Band, calibrator, strehl 0.30J-Band, calibrator, strehl 0.31
J-Band, calibrator, strehl 0.30J-Band, calibrator, strehl 0.35J-Band, calibrator, strehl 0.40
Figure 1: Central 128x128 pixels of the generated J-Band images of a calibrator at 14 mag and with different strehl (from top left to bottom right: 0.20, 0.25, 0.37, 0.29, 0.30, 0.31, 0.33, 0.35, and 0.40).
K-Band, calibrator, strehl 0.20K-Band, calibrator, strehl 0.25K-Band, calibrator, strehl 0.27
K-Band, calibrator, strehl 0.29K-Band, calibrator, strehl 0.30K-Band, calibrator, strehl 0.31
K-Band, calibrator, strehl 0.30K-Band, calibrator, strehl 0.35K-Band, calibrator, strehl 0.40
Figure 2: Central 128x128 pixels of the generated K-Band images of a calibrator at 14 mag and with different strehl (from top left to bottom right: 0.20, 0.25, 0.37, 0.29, 0.30, 0.31, 0.33, 0.35, and 0.40).

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

A comparison of the reconstructions depending on the calibrator strehl 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 with 0.20 strehl calibratorJ-Band, Richardson-Lucy reconstruction with 0.30 strehl calibratorJ-Band, Richardson-Lucy reconstruction with 0.40 strehl calibratorJ-Band, Richardson-Lucy reconstruction with 1.00 strehl calibrator
J-Band, Building-Block reconstruction with 0.20 strehl calibratorJ-Band, Building-Block reconstruction with 0.30 strehl calibratorJ-Band, Building-Block reconstruction with 0.40 strehl calibratorJ-Band, Building-Block reconstruction with 1.00 strehl calibrator
Figure 4: AGN: The J-Band reconstructions depending on the calibrator strehl 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 strehl was 0.20, 0.30 (same as for the target), 0.40, and 1.0 (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 with 0.20 strehl calibratorK-Band, Richardson-Lucy reconstruction with 0.30 strehl calibratorK-Band, Richardson-Lucy reconstruction with 0.40 strehl calibratorK-Band, Richardson-Lucy reconstruction with 1.00 strehl calibrator
K-Band, Building-Block reconstruction with 0.20 strehl calibratorK-Band, Building-Block reconstruction with 0.30 strehl calibratorK-Band, Building-Block reconstruction with 0.40 strehl calibratorK-Band, Building-Block reconstruction with 1.00 strehl calibrator
Figure 5: AGN: The K-Band reconstructions depending on the calibrator strehl 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 strehl was 0.20, 0.30 (same as for the target), 0.40, and 1.0 (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 star cluster depending on the calibrator strehl are presented. In figure 6 the image errors depending on the calibrator strehl ratio are shown.

TestcaseRichardson-LucyBuilding-Block
Strehl (C)Image ErrorImage Error
0.200.268 (0.248)0.392 (0.392)
0.250.120 (0.114)0.171 (0.171)
0.270.073 (0.071)0.098 (0.098)
0.290.027 (0.027)0.033 (0.033)
0.300.013 (0.013)0.009 (0.009)
0.310.031 (0.031)0.028 (0.028)
0.330.083 (0.083)0.081 (0.081)
0.350.129 (0.129)0.127 (0.127)
0.400.226 (0.226)0.225 (0.225)
0.500.000 (0.349)0.348 (0.348)
0.700.511 (0.511)0.511 (0.511)
1.000.642 (0.642)0.643 (0.643)
Table 2: AGN: Errors depending on the calibrator strehl (J-Band).
TestcaseRichardson-LucyBuilding-Block
Strehl (C)Image ErrorImage Error
0.200.319 (0.284)0.417 (0.417)
0.250.151 (0.142)0.191 (0.191)
0.270.082 (0.081)0.103 (0.103)
0.290.034 (0.034)0.037 (0.037)
0.300.020 (0.020)0.015 (0.015)
0.310.035 (0.035)0.028 (0.028)
0.330.081 (0.081)0.076 (0.076)
0.350.120 (0.120)0.116 (0.116)
0.400.214 (0.214)0.212 (0.212)
0.500.342 (0.342)0.341 (0.341)
0.700.493 (0.493)0.493 (0.493)
1.000.611 (0.611)0.613 (0.613)
Table 3: AGN: Errors depending on the calibrator strehl (K-Band).
Image errors depending on the calibrator strehl ratio
Figure 6: In this figure, the image errors depending on the calibrator strehl ratio 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 calibrator strehl on the reconstruction error and the optimal iteration number is clearly visible.

Richardson-Lucy reconstruction error (J-Band)Richardson-Lucy reconstruction error (K-Band)
Building-Block reconstruction error (J-Band)Building-Block reconstruction error (K-Band)
Figure 7: Reconstruction error depending on the iteration number for the J-Band (left side) and K-Band (right side). In the first row, Richardson-Lucy was used for the reconstruction. 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 NHC4151 compared with the ideal reference image is shown.

Richardson-Lucy: Profile through the reconstructed galaxy (J-Band)CRichardson-Lucy: entral part of the profile through the reconstructed galaxy (J-Band)
Building-Block: Profile through the reconstructed galaxy (J-Band)Building-Block: Central part of the profile through the reconstructed galaxy (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. In the first row, Richardson-Lucy was used for the reconstruction. In the second row, the Building-Block method was used.
Richardson-Lucy: Profile through the reconstructed galaxy (K-Band)Richardson-Lucy: Central part of the profile through the reconstructed galaxy (K-Band)
Building-Block: Profile through the reconstructed galaxy (K-Band)Building-Block: Central part of the profile through the reconstructed galaxy (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. In the first row, Richardson-Lucy was used for the reconstruction. 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 strehl is presented in figure 11 (J-Band) and figure 12 (K-Band).

J-Band reference imageJ-Band coadded raw image
J-Band, Richardson-Lucy reconstruction with 0.20 strehl calibratorJ-Band, Richardson-Lucy reconstruction with 0.30 strehl calibratorJ-Band, Richardson-Lucy reconstruction with 0.40 strehl calibratorJ-Band, Richardson-Lucy reconstruction with 1.00 strehl calibrator
J-Band, Building-Block reconstruction with 0.20 strehl calibratorJ-Band, Building-Block reconstruction with 0.30 strehl calibratorJ-Band, Building-Block reconstruction with 0.40 strehl calibratorJ-Band, Building-Block reconstruction with 1.00 strehl calibrator
Figure 11: Star cluster: The J-Band reconstructions depending on the calibrator strehl 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 strehl was 0.20, 0.30 (same as for the target), 0.40, and 1.0 (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 with 0.20 strehl calibratorK-Band, Richardson-Lucy reconstruction with 0.30 strehl calibratorK-Band, Richardson-Lucy reconstruction with 0.40 strehl calibratorK-Band, Richardson-Lucy reconstruction with 1.00 strehl calibrator
K-Band, Building-Block reconstruction with 0.20 strehl calibratorK-Band, Building-Block reconstruction with 0.30 strehl calibratorK-Band, Building-Block reconstruction with 0.40 strehl calibratorK-Band, Building-Block reconstruction with 1.00 strehl calibrator
Figure 12: Star cluster: The K-Band reconstructions depending on the calibrator strehl 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 strehl was 0.20, 0.30 (same as for the target), 0.40, and 1.0 (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 calibrator strehl are presented.. The image errors are also summarized in figure 13 and the photometric errors in figure 14.

TestcaseRichardson-LucyBuilding-Block
Strehl (C)Image ErrorMagnitude ErrorImage ErrorMagnitude Error
0.200.103 (0.103)0.08, -0.07, 1.87, 4.81, 2.670.184 (0.182)0.53, -0.08, 1.75, 4.27, 2.30
0.250.062 (0.062)0.02, -0.04, 1.44, 1.67, 2.690.097 (0.091)0.13, -0.06, 1.78, 4.29, 2.35
0.270.036 (0.036)0.01, -0.03, 0.57, 0.70, 2.700.073 (0.068)0.04, -0.05, 1.79, 4.29, 2.36
0.290.014 (0.014)0.01, -0.01, 0.14, 0.24, 0.520.029 (0.029)0.02, -0.02, 0.56, 0.86, 2.39
0.300.014 (0.014)0.01, 0.00, 0.02, 0.11, 0.060.013 (0.013)0.01, -0.00, 0.12, 0.25, 0.42
0.310.038 (0.038)0.04, 0.03, -0.00, 0.10, -0.060.034 (0.034)0.04, 0.03, 0.01, 0.10, -0.02
0.330.094 (0.094)0.10, 0.10, -0.08, 0.12, -0.220.092 (0.092)0.10, 0.09, -0.07, 0.13, -0.21
0.350.146 (0.146)0.16, 0.16, -0.10, 0.13, -0.340.144 (0.144)0.16, 0.15, -0.09, 0.13, -0.35
0.400.252 (0.252)0.31, 0.30, -0.20, 0.19, -0.510.251 (0.251)0.30, 0.29, -0.20, 0.19, -0.51
0.500.383 (0.382)0.51, 0.50, -0.31, 0.31, -0.680.376 (0.376)0.50, 0.49, -0.31, 0.31, -0.67
0.700.556 (0.555)0.85, 0.84, -0.38, 0.40, -0.880.552 (0.552)0.84, 0.83, -0.38, 0.40, -0.88
1.000.698 (0.698)1.25, 1.23, -0.38, 0.44, -1.010.699 (0.699)1.25, 1.24, -0.39, 0.44, -1.01
Table 4: Star cluster: Errors depending on the calibrator strehl (J-Band).
TestcaseRichardson-LucyBuilding-Block
Strehl (C)Image ErrorMagnitude ErrorImage ErrorMagnitude Error
0.200.191 (0.191)0.37, -0.09, 1.50, 4.04, 2.020.210 (0.206)0.51, -0.08, 1.44, 3.80, 1.84
0.250.111 (0.111)0.10, -0.08, 1.51, 4.03, 2.040.113 (0.107)0.15, -0.07, 1.45, 3.79, 1.87
0.270.077 (0.077)0.03, -0.06, 1.52, 4.02, 2.060.079 (0.075)0.05, -0.06, 1.46, 3.79, 1.88
0.290.050 (0.050)0.00, -0.04, 0.47, 1.53, 2.080.044 (0.043)0.01, -0.04, 0.58, 2.10, 1.90
0.300.047 (0.047)0.00, -0.02, 0.11, 0.75, 2.080.035 (0.035)0.01, -0.01, 0.16, 0.93, 1.92
0.310.055 (0.055)0.02, 0.01, -0.02, 0.53, 0.580.042 (0.042)0.03, 0.02, 0.01, 0.60, 0.71
0.330.098 (0.098)0.07, 0.07, -0.09, 0.47, 0.120.089 (0.089)0.08, 0.08, -0.08, 0.51, 0.19
0.350.139 (0.139)0.12, 0.12, -0.16, 0.41, -0.120.133 (0.133)0.12, 0.13, -0.14, 0.46, -0.08
0.400.244 (0.244)0.25, 0.26, -0.28, 0.42, -0.430.241 (0.241)0.27, 0.27, -0.26, 0.44, -0.40
0.500.386 (0.386)0.47, 0.47, -0.36, 0.45, -0.680.387 (0.387)0.49, 0.48, -0.35, 0.46, -0.65
0.700.553 (0.553)0.81, 0.80, -0.44, 0.51, -0.850.554 (0.554)0.82, 0.81, -0.43, 0.51, -0.83
1.000.680 (0.680)1.16, 1.15, -0.46, 0.56, -0.970.682 (0.682)1.17, 1.16, -0.44, 0.57, -0.97
Table 5: Star cluster: Errors depending on the calibrator strehl (K-Band).
Image errors depending on the calibrator strehl ratio
Figure 13: In this figure, the image errors depending on the calibrator strehl ratio are shown.
Photometric errors depending on the calibrator strehl ratio using Richardson-Lucy (J-Band)Photometric errors depending on the calibrator strehl ratio using Richardson-Lucy (K-Band)
Photometric errors depending on the calibrator strehl ratio using Building-Block (J-Band)Photometric errors depending on the calibrator strehl ratio using Building-Block (K-Band)
Figure 14: Photometric errors depending on the calibrator strehl ratio 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 calibrator strehl on the reconstruction error and the optimal iteration number is clearly visible.

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

In the figures figure , figure 16, figure 17, figure 18, and figure 19 the image and photometric errors depending on the calibrator strehl 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 (J-Band, 0.20 strehl), Richardson-LucyImage and photometric errors depending on the iteration number (K-Band, 0.20 strehl), Richardson-Lucy
Image and photometric errors depending on the iteration number (J-Band, 0.20 strehl), Building-BlockImage and photometric errors depending on the iteration number (K-Band, 0.20 strehl), Building-Block
Figure 16: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The calibrator has a strehl of 0.20. In the first row, Richardson-Lucy was used for the reconstruction. In the second row, the Building-Block method was used.
Image and photometric errors depending on the iteration number (J-Band, 0.30 strehl), Richardson-LucyImage and photometric errors depending on the iteration number (K-Band, 0.30 strehl), Richardson-Lucy
Image and photometric errors depending on the iteration number (J-Band, 0.30 strehl), Building-BlockImage and photometric errors depending on the iteration number (K-Band, 0.30 strehl), Building-Block
Figure 17: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The calibrator has a strehl of 0.30. In the first row, Richardson-Lucy was used for the reconstruction. In the second row, the Building-Block method was used.
Image and photometric errors depending on the iteration number (J-Band, 0.40 strehl), Richardson-LucyImage and photometric errors depending on the iteration number (K-Band, 0.40 strehl), Richardson-Lucy
Image and photometric errors depending on the iteration number (J-Band, 0.40 strehl), Building-BlockImage and photometric errors depending on the iteration number (K-Band, 0.40 strehl), Building-Block
Figure 18: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The calibrator has a strehl of 0.40. In the first row, Richardson-Lucy was used for the reconstruction. In the second row, the Building-Block method was used.
Image and photometric errors depending on the iteration number (J-Band, 1.00 strehl), Richardson-LucyImage and photometric errors depending on the iteration number (K-Band, 1.00 strehl), Richardson-Lucy
Image and photometric errors depending on the iteration number (J-Band, 1.00 strehl), Building-BlockImage and photometric errors depending on the iteration number (K-Band, 1.00 strehl), Building-Block
Figure 19: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The calibrator has a strehl of 1.00. In the first row, Richardson-Lucy was used for the reconstruction. In the second row, the Building-Block method was used.