,, ,,, and Rebeca de Sousa: Miguel R. Alarcon Javier Licandro Miquel Serra-Ricart Enrique Joven Vicens Gaitan
,, ,,, and Rebeca de Sousa: Miguel R. Alarcon Javier Licandro Miquel Serra-Ricart Enrique Joven Vicens Gaitan
3 Aplicaciones en Informática Avanzada (AIA), E-08172 Sant Cugat del Vallès, Catalonia, Spain
4 Departamento de Historia y Filosofı́a de la Ciencia, la Educación y el Lenguaje, Universidad de La Laguna (ULL), E-38206 La Laguna,
Canarias, Spain
ABSTRACT
Scientific complementary metal-oxide-semiconductor (CMOS) detectors have developed quickly in
recent years thanks to their low cost and high availability. They also have some advantages over charge-
coupled devices (CCDs), such as high frame rate or typically lower readout noise. These sensors started
to be used in astronomy following the development of the first back-illuminated models. Therefore,
it is worth studying their characteristics, advantages, and weaknesses. One of the most widespread
CMOS sensors are those from the Sony IMX series, which are included in large astronomical survey
projects based on small and fast telescopes because of their low cost, and capability for wide-field and
high-cadence surveys. In this paper, we aim to characterize the IMX455M and IMX411M sensors,
which are integrated into the QHY600 and QHY411 cameras, respectively, for use in astronomical
observations. These are large (36 × 24 and 54 × 40 mm) native 16 bit sensors with 3.76 µm pixels
and are sensitive in the optical range. We present the results of the laboratory characterization of
both cameras. They showed a very low dark current of 0.011 and 0.007 e− px−1 s−1 @–10◦ C for the
QHY600 and QHY411 cameras, respectively. They also show the presence of warm pixels, ∼0.024% in
the QHY600 and 0.005% in the QHY411. Warm pixels proved to be stable and linear with exposure
time, and are therefore easily corrected using dark frames. Pixels affected by the Salt & Pepper noise
are ∼2% of the total and a method to correct for this effect is presented. Both cameras were attached
to night telescopes and several on-sky tests were performed to prove their capabilities. On-sky tests
demonstrate that these CMOS behave as well as CCDs of similar characteristics and (for example)
they can attain photometric accuracies of a few milli-magnitudes.
of the planet Uranus at a wavelength of 8900 Å (Janesick out noise (1.6 e− ), good linearity (99.7%), and a stable
et al. 1987). The development and refinement of CCDs bias. This sensor family has a dual-amplifier structure
have increasingly made these sensors the most widely in which two 12 bit images are taken simultaneously and
used option for the manufacture of astronomical instru- are then merged into a single 16 bit high dynamic range
ments. image. However, this mechanism is susceptible to jumps
Over the last two decades, the rise of alternative and instabilities around the transition region between
technologies has undermined the prevalence of CCD high and low gain. In addition, edge glowing and charge
sensors. Complementary metal-oxide-semiconductor persistence over long periods of time is also present. The
(CMOS) image sensors started to be developed in the capabilities of these sensors decrease considerably with
1990s (Fossum 1997). However, they still had a num- exposure times longer than several seconds, which is a
ber of disadvantages when compared with CCDs, such consequence of these effects and the increase of dark
as lower dynamic range (DR), and poorer linearity and current; however, they are still suitable instruments for
sensitivity (Bigas et al. 2006). Although CMOS sen- high frame rate observations.
sors soon established themselves in the consumer mar- Many leading manufacturers are now developing in-
ket, their inherent constraints restricted their applica- strumentation based on next-generation sCMOS sen-
tion in certain fields, especially those related to sci- sors, but their suitability for general use in astronomy
ence. To overcome the typical limitations of CMOS, is still largely unexplored. One of the most widespread
the so-called scientific CMOS (sCMOS) were introduced sensors are those from the Sony IMX series, which are
in 2009 as a result of a collaboration between An- being included in large projects such as the Argus Op-
dor Technology, Fairchild Imaging (BAE Systems) and tical Array (Law et al. 2022), Large Array Survey Tele-
PCO Imaging (Coates et al. 2009). This new genera- scope (Ofek et al. 2023) or the next generation of tele-
tion of sensors combined high frame rates, reasonable scopes for the ATLAS project (Tonry et al. 2018) , such
pixel and sensor sizes, and quantum efficiencies com- as the one that will be installed at Teide Observatory,
parable to CCDs—especially back-illuminated (BI) sC- ATLAS-Teide (Licandro et al. 2023a), because of their
MOS (Princeton Instruments 2016)—and a considerable low cost and capability for wide-field and high-cadence
reduction in the noise levels that are traditionally asso- surveys.
ciated with CMOS. In this paper, we present the results of the laboratory
Scientific CMOS are beginning to be accepted in as- characterization of the sCMOS BI Sony IMX455M and
tronomy. The first analyses in this respect mainly IMX411M sensors when integrated into the QHY600M
involve two sensors: CIS2051 (later rebranded as and QHY411M cameras. In Section 2, the two de-
CIS2521), which was developed by Andor Technol- vices are tested, and the laboratory set-up and the tele-
ogy and integrated into the Neo series detectors; and scopes that we used are presented. Sections 3.1 and
GSENSE, which was developed by Gpixel and integrated 3.2 describe the tests performed under dark conditions:
by several companies, e.g., Finger Lakes Instrumenta- spatio-temporal variation of bias and dark and contam-
tion, QHYCCD, and Andor Technology. Experiments inating effects, such as random telegraph noise and the
performed with Andor Neo showed low readout noise presence of warm pixels. In Section 3.3, the main op-
(∼1 e− ) and linearity deviations up to the saturation erating features such as the gain fix pattern noise and
point within the expected range of ±1%. However, two linearity are verified using the photon transfer curve.
features were noticed: the transfer curve was observed The quantum efficiency measures are shown in Section
to bend, even at low signal levels, which was possibly 3.4 and the charge persistence effect reported in other
caused by the non-linearity of an amplifier stage; and sCMOS sensors is reviewed in Section 3.5. Finally, sev-
the signal shows certain irregularities with high vari- eral approaches to processing telescope data and on-sky
ance in the transition region between high and low gain results based on this analysis are discussed in Section 5.
modes (Schildknecht et al. 2013). This, along with a
limited QE due to being front-illuminated and a lower
fill factor (Qiu et al. 2013), means that this sensor is not
the most suitable for regular observations, and its use-
fulness is mostly constrained to bright objects and high
frame rates. The back-illuminated GSENSE2020 sensor,
integrated in the Andor Marana, was recently tested for
its application in astronomy by Qiu et al. (2021) and 2. METHODS
Karpov et al. (2021). This sensor showed a low read- 2.1. Instruments
3
Figure 1. Diagram of the optical test bench set-up, with the (1) Newport 68945 digital power supply, (2) Newport M-66881
QTH lamp, (3) Newport 76994 shutter, (4) Newport Oriel Cornerstone monochromator, (5) Hamamatsu S1336-5B1 photodiode,
(6) Labsphere SC6000, (7) Labsphere US-080-SF/SL integrating sphere, (9) Hamamatsu S2281 photodiode, and (10) QHY600M
Pro / QHY411M
Cooling system
(temperatures below ambient) Air cooling (-30C) Air cooling (-35C)
Water cooling (-45C)
4
The QHY600M1 camera is based on the back- integrating sphere, placed next to the monochromator,
illuminated IMX455 monochrome sensor that is man- had incorporated a Hamamatsu S1336-5B1 photodiode,
ufactured by SONY, which is a full-frame (35 mm for- which was in turn connected to a Keysight B2980A pi-
mat) sensor with 9576 × 6388, 3.76 µm square pixels. coammeter. Next, a dark box was located, followed by
The QHY411M2 camera is based on SONY’s IMX411 the test of the camera. A second photodiode (Hama-
monochrome sensor, which is also backilluminated but matsu S2281) was connected to the exit port of the dark
with a larger sensor size, 14304 × 10748 (equivalent to box. The conversion factors between the intensity mea-
medium-format cameras, 54 × 40 mm). Both of these sured by the picoammeter on the first photodiode and
sensors include an overscan region of 33 and 91 rows, the radiant power received on the second photodiode at
respectively, and they are native ADC sampled at 16 the exit of the dark box were characterized, in addition
bit, which is a significant change from previous genera- to the correction to the distance between the photodiode
tions of sCMOS sensors that were based on 12 bit im- and the sensor back focus. During the tests, the labora-
age merging (see Karpov et al. (2021)). The main fea- tory temperature was about 23◦ C and the humidity was
tures of the cameras are listed in Table 1. We used the 40–50%. Both cameras were air-cooled, with operating
QHY600M Pro version, which allows a faster 2×10 Gbps temperatures of −5◦ C for the QHY600M Pro and 0◦ C
fiber connection to a frame grabber with an additional 4 for the QHY411M.
GB of DDR3 memory—although in this paper all tests
were done via USB 3.0 connection—and triggering the 2.3. Telescopes
rolling shutter via an external GPS with an accuracy of The sky tests were performed with two of the robotic
more than a microsecond. It is worth mentioning that telescopes (Telescopios Abiertos Robticos, TAR) of
both sensors, especially the IMX455, are used in cam- Teide Observatory (Tenerife, Canary Islands, Spain).
eras from other manufactures, e.g., Atik Apx60, ZWO The QHY600M Pro was installed on the prime focus
ASI6200MM Pro. of TAR03, 0.46-m f /2.8 C18 reflector telescope on a
These cameras can operate in several modes and gain Planewave L500 altaz mount. An ROI of 4096×4096 was
settings, which essentially change their gain, readout used, giving a FOV of 410 .1 × 410 .1 with 000 .6 px−1 . The
noise (RON), and full-well capacity (FWC). In this pa- QHY411M was mounted on a Meade LX200-ACF 1600
per, we have focused on those that are considered most 0.406-m f /10, on an APM GE-300 Direct Drive equa-
appropriate for use in astronomy because they maintain torial mount, with UV/IR-Cut/L, SDSS g0 , r0 and i0 50
a good balance between these characteristics, which are mm filterswhich were manufactured by Baader Plane-
Mode#1 (High Gain Mode) and gain setting 0 (here- tarium GmbH. The sensor was trimmed to 40 × 40 mm,
after, #1@0) on the QHY600M Pro and Mode #4, and giving a FOV of 340 .2 × 340 .2 arcmin with 0.4000 px−1 .
gain setting 0 (hereafter, #4@0) on the QHY411M, al- Both cameras were connected via USB 3.0 and oper-
though some results for other modes are also shown. ated with air cooling at −10◦ C. After the first tests, the
2.2. Optical test bench QHY411M was installed in one of the 80 cm telescopes
of the Two-meter Twin Telescope (TTT) project, which
The evaluation of the sensors was performed using is a 0.80-m f /6.85 Ritchey–Chrétien altaz telescope that
existing experimental equipment, which are available was manufactured by ASA Astrosysteme GmbH, with a
at the Laboratory of Imaging and Sensors for Astron- total FOV of 330 .4×250 .0 and a plate scale of 000 .14 px−1 .
omy (LISA), at the Instituto de Astrofı́sica de Canarias
(Tenerife, Spain). The set-up is schematically shown in 3. RESULTS
Figure 1. A Newport M-66881 QTH lamp was connected
3.1. Bias stability and ”Salt & Pepper” effect
to A Newport 68945 digital power supply, with inten-
sity stabilization. In the Newport Oriel Cornerstone The bias frame stability in both cameras was tested
monochromator, the option to apply an order sorting in the laboratory under conditions of complete darkness.
filter was selected. This filter blocked higher diffraction In all of the tests that are presented here, including the
orders from interfering with the selected wavelengths. telescope runs, the offset setting was fixed at 10. First,
A Newport 76994 shutter was placed between the lamp 21 continuous unbinned full frames bias were obtained
and the monochromator to control the light beam before and stacked with a 3σ-clipping median to obtain the
taking each dark frame. A Labsphere US-080-SF/SL master bias. The result is shown in Figure 2. Both
cameras show dark signal non-uniformity (DSNU), with
a notable column-to-column pattern and pixel-to-pixel
1 https://www.qhyccd.com/scientific-camera-qhy600pro-imx455/ variations. Gradients can be observed in the QHY600M
2 https://www.qhyccd.com/scientific-camera-qhy411-qhy461/ Pro toward the upper region, with variations of less than
5
0.01
Frequency (%)
10 0.050 12.0 < ADU 16.5
0.025
0.2 7.5 < ADU 12.0
160 180 200 220 0.01 0.1 1 10 0.1
10 Frequency (%)
Frequency (%)
21 frames is retrieved. The distribution of each of these a column, given that the average of the signal would
groups is shown in the right-hand plot of Figure 5. For always lie in the central range. To show a kite-shaped
pixels showing a temporal deviation of the order of the dispersion, there must be some pixels that tend to
RON (σADU ≤7.5), the distribution of the measured sig- exhibit more deviations toward one of the Gaussians,
nal is Gaussian, similar to that of a conventional CCD. either the upper or the lower Gaussian, thus skewing
The large majority of the pixels are in these first two the average to either side and resulting in the widen-
ranges (45% of the pixels with σADU ≤3 and 98% with ing of the tails of the mean signal distribution, which
σADU ≤7.5). For the other ranges (with larger devia- is shown in yellow in the lower left-hand plot of Figure 3.
tions), the pixels begin to show a triple Gaussian dis-
tribution, with two smaller peaks that are centered on For a more detailed insight into the behavior of these
both sides of the central peak and with similar widths. pixels with high standard deviation, 1000 consecutive
For larger values of σADU the separation of the peaks is bias frames were obtained with the QHY600M Pro.
greater. In these anomalous pixels, the pixel distribu- The temporal average and standard deviation of all
tion is indicative of the following behavior: while most the frames has been taken in a small central region of
of the time these pixels return values that are around the 20 × 20, which is shown in the upper part of Figure 4.
mean value of the bias level with a normal noise equiva- Most of the pixels have an average signal around 171
lent to readout, their returned values occasionally jump ADU and a dispersion below the RON, as expected.
toward larger or smaller values with defined separations However, some of them show anomalous patterns in the
and different probabilities for each pixel. average value, the standard deviation, or both at the
The symmetry in the three Gaussian distributions same time. Several pixels have been selected as samples,
does not necessarily mean that all of the anomalous showing the time evolution of the signal in the lower
pixels jump between higher and lower levels with the plots. As a reference, a pixel with normal values of
same probability. In that case, the standard versus average signal and deviation has been taken, following
average plot would not show a kite-shaped scatter but a normal distribution over the 1000 frames, with mean
7
5 C 5 C
10 0 C 10 10 C
5 C 15 C
10 C 20 C
1 1
Frecuency (%)
Frecuency (%)
0.1 0.1
0.01 0.01
Figure 5. Dark current distributions for the QHY600M Pro (left-hand panel) and the QHY411M (right-hand panel). They
were obtained as the median stacked set of five dark frames with an exposure time of 1000 s.
Noise (ADU)
101 101
100 101 102 103 104 105 100 101 102 103 104 105
Signal (ADU) Signal (ADU)
Figure 8. Photon transfer curves obtained from sets of three bias-subtracted images, median stacked, and taken under uniform
illumination with increasing exposure times. The QHY600M Pro Mode #1@0 (left-hand panel) and the QHY411M Mode #4@0
(right-hand panel). The principal components of the noise have been distinguished as colored lines, whose characteristic values
are also included. The gain and FPN has been obtained from the fitting to the total noise function.
60000 Saturation level = 65532±1 ADU 60000 Saturation level = 65523±1 ADU
FWC = 50.0±0.2 ke FWC = 67.1±0.7 ke
50000 50000
Signal (ADU)
Signal (ADU)
40000 40000
30000 30000
20000 20000
10000 10000
0 0
100 100
Linearity (%)
Linearity (%)
99 99
98 98
0 20 40 60 80 100 120 140 0 20 40 60 80 100 120 140 160
Exposure time (ms) Exposure time (ms)
Figure 9. Average signal in sets of three images taken with increasing exposure time (top) for the QHY600M Pro Mode #1@0
(left-hand panel) and the QHY411M Mode #4@0 (right-hand panel). The saturation level is indicated by a horizontal-dashed
line. A straight line was fitted to points with values between 100 and 60,000 ADU. The relative rms of the differences between
the measurements and the fitted line is shown in the lower plots.
expression (1) on a logarithmic scale. This procedure component between 100 and 10,000 ADU. Thereafter,
was repeated in windows of 200 × 200 pixels along the the contribution of the FPN becomes significant and
ROI to obtaining a curve for each of them, 360 in to- a deviation from Poissonian behavior is observed. In
tal. The rms of the parameter distributions was used the QHY600M Pro, this source is the main contributor
to estimate their uncertainty. In both plots, three re- above 50,000 ADU, with a factor of 0.55% ± 0.02% ,
gions can be distinguished. At low illumination, below being lower, 0.31% ± 0.02%, in the QHY411M, where
10 ADU, RON, which does not depend on the signal, is the shot noise is dominant almost up to the saturation
the dominating noise source. From then on, shot noise point. A PTC has also been obtained for other operat-
starts to become important, and is essentially the main
11
100
Table 3. Main features of the different operating modes of QHY600M PRO
the QHY600M Pro (top) and QHY411M (bottom), obtained QHY411M
in laboratory tests. including the gain, readout noise (RON),
fixed pattern noise (FPN) and full-well capacity (FWC). 80
Figure 11. Etaloning pattern observed when illuminating the QHY600M Pro in a central ROI of 1000 × 1000 pixels 4 × 4
binned with ±1 nm bandwidth light at different wavelengths. Similar behaviour is observed in the QHY411M.
ing to crosstalk between adjacent pixels. This is the nm and above. This is a known effect in thinned back-
reason for the poor performance at longer wavelengths. illuminated devices, resulting from multiple reflections
It should be noted that the curves that are obtained produced inside the depletion region by a mismatch be-
here represent an overall reduction of 9% over the full tween its refractive indices and the adjacent layers. This
absolute QE curve published by QHYCCD, which has effect has not been observed in sky images taken with
a peak of 92% at 450 nm, 46% at 700 nm, and 15% TAR04 using a wider bandwidth, such as SDSS i0 (690–
at 900 nm. Nonetheless, the conditions and configura- 850 nm).
tion used by them, as well as the uncertainties in their
measurements, are unknown and the calibration method 3.5. Charge persistence
to obtain the QE curves is different from the one used Charge persistence is an effect that occurs when a por-
here, which has been used to calibrate multiple astro- tion of the signal remains in the detector element after
nomical instruments before. The same procedure and the sensor has been read out. It is a consequence of the
test bench have been used, for instance, to calculate the creation of traps at the interface between the photodi-
QE in sCMOS cameras such as Andor Marana, FLI- ode and the transfer gates, which capture free electrons
Kepler, and ORCA-Hamamatsu. They have also been and gradually release them, resulting in a decay of the
used in several deep depletion CCDs, such as the well- residual signal in the subsequent images after the illu-
known Teledyne e2v 4482 and 231-84 BI. In all cases, minating source has been removed. This effect has been
QE fits rather well with data supplied by the manufac- observed with tests performed in the laboratory with the
turers. Gill et al. (2022) have very recently presented FLI Kepler KL400 and Andor Marana cameras, both
a low-cost method to obtain, among other features, the with GSENSE400BSI sensors, showing a behavior simi-
absolute QE of detectors being applied to the IMX455 lar to that reported by Karpov et al. (2021). Although
sensor. Their results also show a lower performance than they did not observe any effect in sky images, previ-
that presented by QHYCCD, especially in the red part, ous tests performed with these cameras on the same
although they had a similar peak efficiency, with 93% at telescopes used in our work did show a smearing ef-
480 nm, 42% at 700 nm, and 10% at 900 nm, making fect in pixels that had previously been exposed close to
an overall deviation from our results of 5%. Finally, Be- the FWC, which remained visible up to tens of minutes
toule et al. (2023) have studied the quantum efficiency later.
of the IMX411 sensor in great detail and obtained a re- This effect has been tested in the laboratory for both
sult that is very similar to Figure 10, with a peak of 83% sensors. To do so, following the same set-up used in
at 490 nm, 37% at 700 nm and 6% at 900 nm. Section §3.3, i.e., removing the monochromator and ex-
During the tests, optical etaloning was observed at posing the cameras to the QTH lamp light, a pinhole
longer wavelengths, as shown in Figure 11 for the was placed inside the black box, just in front of the sen-
QHY600M Pro. Similar behavior was observed also for sor. Consecutive 1 s images were taken with the shutter
the QHY411M and reported by Betoule et al. (2023). closed. The shutter was then opened for about 10 s, sat-
The maximum variation over the average frame mean urating the illuminated pixels, and closed again. This
value is around 1% at 650 nm, and reaches 10% at 850 sequence is shown in Figure 12. On the top are the
13
16.85
16.90
16.95
Magnitude
17.00
17.05
17.10
17.15 Period = 3.60 ± 0.03 hr
Amp = 0.13 mag
Figure 12. Average signal (black dots) in 1 second frames 0.0 0.2 0.4 0.6 0.8 1.0
Phase
taken continuously during sensor flooding (central shaded
area). Images before (left-hand), during (center), and after
(right-hand) are shown at the top. The region in which the Figure 13. Phased light curve of the asteroid (3200)
average signal has been calculated is indicated by dashed Phaethon, obtained with the QHY600M mounted at the
rectangles. prime focus of a 0.46-m f /2.2 telescope, no filter installed.
0
tral region of the illuminated spot area, indicated by a
dashed square, the mean signal has been measured and 5
is shown in black dots below. The signal before and after 10
the flooding shows a stable trend, with very small oscil-
lations, not a sharp increase in the signal right after the 15
exposition followed by a slow decay to the original level, 20
as shown in Figure 10 of Karpov et al. (2021). Hence, 25 0.00
we conclude that the IMX455M and IMX411M sensors 0.05 0.10 0.15 0.20 0.25 0.30
JD - 2459936.3319
do not exhibit charge persistence.
Figure 14. Transit of the exoplanet TOI-1135 observed
4. ON-SKY TESTS with QHY411M mounted on a Nasmyth focus of the TTT-2
Both cameras have been extensively tested on images telescope at 0.80-m and f /6.85, SDSS g0 filter.
taken with telescopes, obtaining photometric accuracies
as expected for the characteristics described above. The low readout noise and negligible readout time, it is pos-
QHY600M Pro was installed at the prime focus of a 0.46- sible to take continuous short frames and combine them
m f /2.2 telescope. Figure 13 shows the phased light by aligning with the object or the stars, thus allowing
curve of asteroid (3200) Phaethon, which was observed fainter objects to be reached with very little time lost.
for 6 consecutive hours. At the time of observation, the This also allows the study of, for instance, very fast ro-
object had an apparent magnitude of V = 17 and was tating objects with good temporal sampling (Licandro
moving at a speed of 100 .6 minute−1. The photometric et al. 2023b) or using “shift-and-add” techniques, such
uncertainties that we obtained are in the order of hun- as synthetic tracking (Shao et al. 2014), to improve the
dredths of a magnitude and the rotation of the asteroid, detection performance of faint fast-moving objects.
with an amplitude of 0.13 mag, is clearly detectable.
This camera is a very suitable choice for very fast tele- Figure 14 shows a transit of the exoplanet TOI-1135
scopes with primary focus because its small size and that was observed with the QHY411M mounted on one
compact shape drastically reduce obscuration. With a of the Nasmyth foci of the TTT-2 telescope, 0.80-m
pixel size of 3.76 µm, it also allows us to obtain a plate f /6.85, SDSS g0 filter, with a plate scale of 000 .14 px−1 .
scale that is very suitable for sites with excellent seeing, Groups of four images of 5 s exposure were stacked to
such as the Teide Observatory. Furthermore, given its improve the SNR. The standard deviation of the differ-
14
ential aperture photometry measurements in pre-transit, tometric measurements. Figure 15(a) shows an image
during, and post-transit is 2.2, 3.3, and 2.0 mmag, re- of the galaxy M33 taken with the QHY411M on the
spectively (M. Mallorquin et al., in preparation). With TAR04 telescope, with a UV/IR-Cut/L filter, in a sin-
such a small plate scale, the psf was oversampled, with gle exposure of 60 s. When zooming in on a small area of
approximately 8 px of FWHM, spreading the flux of the 100×100 pixels (b), several pixels are observed with val-
star over a larger number of pixels. For very bright tar- ues that are clearly higher than those of their neighbors,
gets, such as this V = 9.6 mag star, this allows slightly as described in Section 3.1. This effect could influence
longer exposures to be taken without reaching the satu- the photometry of faint objects because the amplitude
ration point, thus reducing the random and flat-fielding of this random fluctuation may be comparable to the
errors of the telescope that, with a larger scale, would source signal in the pixel.
perhaps need some defocusing (Southworth et al. 2009). SSeveral strategies can be pursued to mitigate this
In addition, by taking short exposures, the time lost on problem, especially by using several consecutive frames.
a CCD could be equivalent to, or even longer than, the Figure 15(c) shows a sequence of 3 × 20s, where the
exposure time, which makes it very inefficient. With the S&P is clearly visible in each of them. If simple av-
QHY411M, the exposure time could be extended, thus eraging were performed, then the outliers would skew
improving photometric accuracy, and the readout time the signal obtained. Algorithms such as 3σ clipping
is almost non-existent. Although such small plate scales that use outlier-sensitive dispersion measures generally
are generally undesirable, they allow these cameras to do not work either because these metrics are biased by
be used in other scientific applications (e.g., fast-moving these fluctuations and define a very wide clipping range.
object astrometry or lucky imaging). This is seen in Figure 15(d) (top plot), where the three
frames have been combined with an average after 3σ-
5. DISCUSSION clipping. Some of the pixels that exhibited S&P also
show deviating values in the stacked frame.
The two instruments with sCMOS sensors that are
It is common to use spatial filters for this kind of
analyzed here present characteristics that are compat-
localized noise. The most typical for S&P is a median
ible with their use in astronomy: they are linear over
filter, where the value of each pixel is replaced by the
the whole dynamic range, have a high full-well capac-
median of the values of this and its close neighbors. This
ity, and are slightly affected by dark current, despite
is done in the second row of Figure 15(c). Note that
being able to work at higher temperatures than CCDs.
the image has been smoothed and the anomalous pixels
Regarding the quantum efficiency, although the curve
may have disappeared. However, this has been done at
obtained here is slightly lower than that reported by the
the expense of: (1) changing the signal and noise pat-
manufacturer, 80% at 500 nm is an acceptable perfor-
tern of the image and (2) correlating the nearby pixels.
mance for many scientific programmes and is in general
Although this is a very useful filter for improving the
similar or better than other CCD sensors in the same
cosmetics of the images, the photometry measurement
cost range. An improvement in efficiency toward redder
in the resulting image may be highly biased because the
wavelengths should be achieved in the next few years,
fluxes of each pixel have been altered by its neighbors.
so that sCMOS sensors can be used on a wider variety
of observational targets.
In this work, a solution based on convolutions and
The small pixel size means that these sensors are gen-
these two previous ideas is proposed to try to mitigate
erally not the best solution for slow focal length systems,
the S&P effect. First, it should be considered that the
except for dedicated programmes such as high spatial
S&P effect mostly impacts low noise areas, such as the
resolution or lucky imaging. Binning in sCMOS sensors
sky background or faint sources. For bright sources, the
is done after exposure, so it does not improve the read-
shot noise becomes higher than the random telegraph
out noise or frame rate. In general, in case it is needed,
noise and dominates all of the other fluctuations. The
it is better to do this by software after the exposure, so
median filter can be used to obtain a reference frame to
that the statistics can be preserved and the 16 bit limit is
identify outliers because, in well-sampled fields, they are
not reached. Nevertheless, they can be very valuable in
a good approximation to a smoothed frame. This can be
fast telescopes with larger fields and higher plate scales,
seen in Figure 15(c), where the third row shows the dif-
which allows better sampling of the PSF. In addition,
ference between the raw frame and the one filtered with
their manageability, and small size and weight are very
a 3 × 3 median kernel. The residual pattern is generally
interesting, e.g., for prime focus telescopes.
homogeneous, both in the sky area and in the vicinity of
The S&P effect is one of the main issues that affect
the sources. To identify the S&P, a threshold of 12 ADU
the use of these IMX455 and IMX411 sensors for pho-
15
Figure 15. (a) Image of M33 taken with the QHY411M at the TAR04 telescope with an exposure time of 60 s. (b) Zoomed
region 100 × 100 pixels where the Salt & Pepper effect is visible. (c) From top to bottom: 3 × 20 s sequence of the same
zoomed region; convolution with a 3 × 3 median filter; difference between the frames in the previous rows, showing the S&P
contaminated pixels in red; and raw frames with those pixels masked in red. (d) Result of stacking the three frames with an
average 3σ-sigma clipping (top) and with an average after masking the S&P affected pixels.
16
has been set because it is at this point that the distribu- 2. Their dark current is very low, as is the number of
tion of three Gaussians in Figure 3 (right-hand) starts warm pixels. They are stable for at least several
to be revealed. Hence, all of the pixels whose absolute weeks and their signal scales linearly with expo-
difference between the raw value and that resulting from sure time, so they may be quite easily removed
the convolution with the median filter is greater than 12 with dark subtraction. This is an improvement
ADU are masked. In the bottom row of Figure 15(c), from previous sCMOS sensors and makes it pos-
the raw frames are shown with the pixels masked in red. sible to take images with longer exposure times
By having a sequence of frames, the average of the un- without being affected by dark current.
masked values of each pixel can be taken to get a stacked
3. Its quantum efficiency peaks at 80% at 475 nm
image, which is shown on the bottom right-hand of the
and then drops rapidly at longer wavelengths, with
figure, where the S&P contamination has been highly
40% at 700 nm and 10% at 900 nm.
reduced. In cases where a pixel shows S&P in all of the
frames of the sequence and is therefore masked, the as- 4. They do not exhibit charge persistence or edge
tronomer has to decide, for instance, either not to take glow.
that pixel into account in the photometry or to replace
its value by an approximation, such as an interpolation 5. These sensors are affected by random telegraph
of the neighboring pixels or their median. In the exam- noise, which can introduce non-negligible de-
ple included here, this happens in only two pixels out viations in photometric measurements of low-
of 104 . It should be noted that this method may be brightness sources. Simple frame averaging, even
less accurate in fields with critically sampled sources. A with algorithms such as σ-clipping, is not enough
further review is currently underway for future work. to mitigate its effect.
Developing new algorithms or even using those already 6. They show promising performance on photomet-
available in common packages such as IRAF (Tody 1986) ric observations done with both fast and slow tele-
or Astropy (Astropy Collaboration et al. 2013) may not scopes.
be trivial when working with these cameras. The size of
Their low cost, power consumption, and replicability
the raw 16 bit images of the QHY600M Pro are about
make both cameras a very suitable solution for astro-
120 MB, while those of the QHY411M are 300 MB. If
nomical applications, especially with regard to their
a sequence is taken at a high frame rate, then the data
high frame range, near-zero readout time, and low read-
set may not be manageable with commonly available
out noise. These sensors are very good options for fast
CPU capacities. The advantage of convolution-based
small telescopes with large fields because of their small
approaches and simple arithmetical operations, such as
pixel size and large formats. The combination of such
median filtering or frame differencing, is that they are
telescopes and cameras permit very large field-of-view
easily deployable on GPUs, which allows the data to be
images with plate scales with reasonably good sampling
processed more efficiently and faster. The development
of the PSF. For instance, a 1100 f /2.2 telescope such
of GPU algorithms for astronomical image processing is
as the Celestron RASA11 with a QHY600 camera can
essential for further progress in the use of large sensors,
produce images covering a FOV of 7.5 deg2 with a plate
such as these sCMOS.
scale of 100 .27 px−1 , which are excellent options for sur-
6. CONCLUSION veys such as ATLAS-Teide (Licandro et al. 2023a). Even
so, having the necessary computational tools to process
In the previous sections, the key features of the
the data, especially GPU developments, is essential to
QHY600M Pro and QHY411M cameras as scientific in-
take advantage of the full performance of these cameras.
struments have been discussed in detail. For astronomy,
they have characteristics that make them very suitable
The authors declare no conflict of interest or rela-
for general use, although certain issues need to be ad-
tionship with the manufacturers of the cameras tested.
dressed. Our main conclusions are that:
M.R.A, M.S-R and J.L. acknowledge support from the
1. The built-in electronics in the pixels of sCMOS ACIISI, Consejerı́a de Economı́a, Conocimiento y Em-
sensors require that each pixel should be consid- pleo del Gobierno de Canarias and the European Re-
ered as an individual detector and this should gional Development Fund (ERDF) under grant with ref-
be taken into account when performing processes erence ProID2021010134 and from the Agencia Estatal
such as bias or dark subtraction. Spatial inhomo- de Investigacion del Ministerio de Ciencia e Ińnovacion
geneities in darkness are detectable all over the (AEI-MCINN) under grant ”Hydrated Minerals and Or-
frame. ganic Compounds in Primitive Asteroids” with reference
17
PID2020-120464GB-100. This research has been par- presented here. This article includes observations made
tially funded by Light Bridges, SL. which provided the in the Two meter Twin Telescope (TTT) at the IAC’s
QHY411M and Andor iKon-L 936 cameras for the tests Teide Observatory that Light Bridges, SL, operates on
the Island of Tenerife, Canary Islands (Spain).
REFERENCES
Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., Licandro, J., Popescu, M., Tatsumi, E., et al. 2023b,
et al. 2013, A&A, 558, A33, MNRAS, 521, 3784, doi: 10.1093/mnras/stad708
doi: 10.1051/0004-6361/201322068 Martin-Martinez, J., Rodriguez, R., & Nafria, M. 2020,
Betoule, M., Antier, S., Bertin, E., et al. 2023, A&A, 670, Advanced characterization and analysis of random
A119, doi: 10.1051/0004-6361/202244973 telegraph noise in CMOS devices (Springer), 467–493
Bigas, M., Cabruja, E., Forest, J., & Salvi, J. 2006, Ofek, E. O., Ben-Ami, S., Polishook, D., et al. 2023, arXiv
Microelectronics J., 37, 433, e-prints, arXiv:2304.04796,
doi: https://doi.org/10.1016/j.mejo.2005.07.002 doi: 10.48550/arXiv.2304.04796
Boyle, W. S., & Smith, G. E. 1970, Bell Syst. Tech. J., 49, Princeton Instruments. 2016, New Scientific CMOS
587, doi: 10.1002/j.1538-7305.1970.tb01790.x Cameras with Back-Illuminated Technology. Technical
Coates, C., Fowler, B., & Holts, G. 2009, White Paper: Note., Available at: https:
//www.princetoninstruments.com/wp-content/uploads/
Scientific CMOS technology, a high-performance imaging
2020/11/TechNote sCMOSBackIlluminatedTech.pdf
breakthrough, Available at:
(accessed July 8, 2022)
https://www.pco.de/fileadmin/user upload/db/
Qiu, P., Mao, Y.-N., Lu, X.-M., Xiang, E., & Jiang, X.-J.
download/scmos white paper 2mb.pdf (accessed Oct 3,
2013, Research in Astronomy and Astrophysics, 13, 615,
2022)
doi: 10.1088/1674-4527/13/5/012
Fossum, E. R. 1997, Nucl. Instrum. Methods Phys. Res. A,
Qiu, P., Zhao, Y., Zheng, J., Wang, J.-F., & Jiang, X.-J.
395, 291, doi: 10.1016/S0168-9002(97)00812-7
2021, RAA, 21, 268, doi: 10.1088/1674-4527/21/10/268
Gill, A. S., Shaaban, M. M., Tohuvavohu, A., et al. 2022, in Schildknecht, T., Hinze, A., Schlatter, P., et al. 2013, in
Society of Photo-Optical Instrumentation Engineers Proceedings of 6th European Conference on Space
(SPIE) Conference Series, Vol. 12191, X-Ray, Optical, Debris, Darmstadt, Germany
and Infrared Detectors for Astronomy X, ed. A. D. Shao, M., Nemati, B., Zhai, C., et al. 2014, ApJ, 782, 1,
Holland & J. Beletic, 1219114, doi: 10.1117/12.2627564 doi: 10.1088/0004-637X/782/1/1
Janesick, J. R. 2007, Photon Transfer DN → λ Southworth, J., Hinse, T. C., Jørgensen, U. G., et al. 2009,
(Bellingham, WA: SPIE Press) MNRAS, 396, 1023,
Janesick, J. R., Elliott, T., Collins, S., Blouke, M. M., & doi: 10.1111/j.1365-2966.2009.14767.x
Freeman, J. 1987, Optical Engineering, 26, 692, Tody, D. 1986, in Society of Photo-Optical Instrumentation
doi: 10.1117/12.7974139 Engineers (SPIE) Conference Series, Vol. 627,
Karpov, S., Bajat, A., Christov, A., Prouza, M., & Beskin, Instrumentation in astronomy VI, ed. D. L. Crawford,
G. 2021, arXiv e-prints, arXiv:2101.01517. 733, doi: 10.1117/12.968154
https://arxiv.org/abs/2101.01517 Tonry, J. L., Denneau, L., Heinze, A. N., et al. 2018, PASP,
Law, N. M., Corbett, H., Galliher, N. W., et al. 2022, 130, 064505, doi: 10.1088/1538-3873/aabadf
PASP, 134, 035003, doi: 10.1088/1538-3873/ac4811 Uren, M., Day, D., & Kirton, M. 1985, Applied physics
letters, 47, 1195
Licandro, J., Tonry, J., Alarcon, M. R., Serra-Ricart, M., &
Yokogawa, S., Oshiyama, I., Ikeda, H., et al. 2017, Scientific
Denneau, L. 2023a, arXiv e-prints, arXiv:2302.07954,
Reports, 7, 3832, doi: 10.1038/s41598-017-04200-y
doi: 10.48550/arXiv.2302.07954