• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    The Impact of Bias Row Noise to Photometric Accuracy:Case Study Based on a Scientific CMOS Detector

    2024-03-22 04:11:48LiShaoHuZhanChaoLiuHaonanChiQiuyanLuoHuaipuMuandWenzhongShi

    Li Shao , Hu Zhan,2 , Chao Liu,3,4, Haonan Chi, Qiuyan Luo, Huaipu Mu, and Wenzhong Shi

    1 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China; shaoli@nao.cas.cn

    2 Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China

    3 University of Chinese Academy of Sciences, Beijing 100049, China

    4 Institute for Frontier in Astronomy and Astrophysics, Beijing Normal University, Beijing 100875, China

    Abstract We tested a new model of CMOS detector manufactured by the Gpixel Inc, for potential space astronomical application.In laboratory, we obtain some bias images under the typical application environment.In these bias images, clear random row noise pattern is observed.The row noise also contains some characteristic spatial frequencies.We quantitatively estimated the impact of this feature to photometric measurements, by making simulated images.We compared different bias noise types under strict parameter control.The result shows the row noise will significantly deteriorate the photometric accuracy.It effectively increases the readout noise by a factor of 2–10.However, if it is properly removed, the image quality and photometric accuracy will be significantly improved.

    Key words: instrumentation: detectors – methods: statistical – techniques: image processing

    1.Introduction

    Since the time of digitalization, CCD (charge coupled device) detectors have played an essential role in astronomy(Mackay 1986;Janesick 2001;Howell 2006;Lesser 2015).But after the invention of CMOS (complementary metal oxide semiconductor) image detector (Fossum 1997), CMOS has attracted more and more attentions in astronomy due to its advantages(Qiu et al.2013;Wang et al.2014;Mendikoa et al.2016; Jorden et al.2017; Shugarov 2020; Ardilanov et al.2021; Greffe et al.2022; Liu et al.2022), especially in time domain area (Pratlong et al.2016; Niu et al.2022; Song et al.2022).It also shows the potential to replace CCD (see e.g.,Gach et al.2022).

    During the development of the Chinese Space Station Telescope(CSST),both CCD and CMOS were considered as possible options for the CSST Survey Camera (CSC).The CSC is designed for large area optical survey.It uses 30 detectors to form a giant 2.3 billion pixel array to cover a field of view of more than one square degree.The optical survey will cover wavelengths ranging from 250 to 1000 nm.The detectors are required to have high quantum efficiency,low readout noise, low dark current and large dynamic range to fit for the proposed wide field blind survey (Zhan 2011,2018, 2021).In order to test whether the state-of-art CMOS detectors can satisfy the scientific requirement, we perform series of tests in the laboratory at the National Astronomical Observatories,Chinese Academy of Sciences,in Beijing.The detector we tested is an HR9090BSI scientific CMOS(sCMOS) detector5The detector we tested is a customized one.The standard version of this model is officially named as “GSENSE1081BSI.” A Chinese exoplanet searching project called “Earth 2.0” (Song et al.2022) uses the standard version.manufactured by the Gpixel Inc.6official website: https://www.gpixel.com/.This model is under active development in last few years.The Gpixel Inc.also produces some other science level CMOS detectors for (potential) astronomical applications (Shugarov 2020; Gill et al.2022).

    The readout noise is one of the most important parameters to evaluate the performance of a detector.The bias noise properties do have subtle consequences to final photometry or spectroscopy quality(Jiménez Mu?oz et al.2021).Only a few studies were done on CMOS detectors in astronomy field (Basden 2015; Basden &Morris 2016; Greffe et al.2022), probably due to limited application cases.More can be found in non-astronomy scientific applications(Diekmann et al.2017;Matsiaka&Plakhotnik 2021).In this paper,we will study the bias noise properties of this detector and its photometric performance with simulated images.Since this model is still under development, the content of this paper only reflects the status before the end of 2022.

    2.The Noise Properties of CMOS Bias Image

    2.1.The CMOS Detector and the Bias Images

    Figure 1.The setup of the detector test bench.

    The bias images were taken from an HR9090BSI sCMOS detector (serial number HR9090BSI-A10LV-SV-22001JC-0027).This one is designed to optimize the performance in visible bands.The pixel physical size is 10 μm×10 μm.The detector supports on-chip 16 bit column-parallel analog-todigital converter (ADC), with up to four low voltage differential signalling (LVDS) channels for image output.The detector uses rolling shutter.The supported minimal exposure time is about 0.6 ms, which is used for taking bias images.The readout time is about 3 s when using two LVDS channels at frequency of 250 MHz.The raw image has 8976×9222 pixels.On both left and right sides, there are 34 columns at the edge of array designed as “electrical black”(EB) region.The readings of these pixels only come from electronics.In principle,they can be used to calculate bias level as well (similar to overscan region in CCD case).However in practice we found the values in EB region were exceptionally large, and therefore not useful for our analysis.In this paper,we only use the central 8900×9120 pixels (“active” region).

    We prepared a test bench (Figure 1) for the test of CSC detectors in an ISO 7 clean room.It can satisfy multiple test requirements (e.g., pixel response at different wavelengths).The optical system was kept in a light tight box, which was installed on a vibration controlled optical table.All devices emitting light(e.g.,running indicator)were carefully separated.Light intensity monitor system was also used when necessary.It makes sure the bias and dark images are not contaminated by any external light source.The detector was installed in a vacuum dewar,cooled with liquid nitrogen.The temperature of the detector was kept at about 188 K, with typical temperature fluctuation no more than 1 K, similar to the designed on-board environment.

    The bias pattern is not very sensitive to the temperature.The variation of temperature has negligible contribution to the bias noise.The backend test board was also provided by the Gpixel Inc.(serial number 008).The test was performed in full dark environment to avoid external light contamination to the bias images.Given the typical dark current of ≯0.02 e-/s/pixel and high readout frequency, the dark current contribution during readout is negligible.In this paper, we use 30 bias images that were continuously taken on 2022 September 5.We note that other data sets show similar properties.

    In this paper we will use ADU for simplicity.The gain value(in astronomical definition) is about 1.3 e-/ADU (measured from photon transfer curve7In the CMOS detector, each pixel has its own gain value.The traditional photon transfer curve technique only gives a typical value.).Readers can convert all the results with this gain if they prefer unit of electrons.

    Figure 2.Left: master bias.Middle: individual bias image.Right: noise residual image (middle–left).

    2.2.Noise Residual Images and Row Noise

    We first stacked the 30 bias images to make a master bias.The stacking was performed pixel by pixel.In each pixel, we sorted the 30 values, excluded the top two and bottom two values to avoid any possible cosmic ray contamination, and then took the average of the rest.We note if we do not exclude any data or exclude more data,the result does not significantly change.On the master bias,we identified a pixel as“bad”if its value exceeded the 7σ threshold of typical(median)bias value.There are totally 52565 pixels marked as “bad,” including two rows.Most defected pixels are excluded in this way.The master bias has a strong feature of horizontal and vertical lines as well as square blocks.A similar feature was also observed in the other CMOS model with large arrays (e.g., Wang et al.2014; Karpov et al.2020).

    The master bias was then subtracted from each individual bias image to remove the fixed pattern,leaving only the noise fluctuation.Hereafter we call such a fixed-pattern-corrected bias image as “noise residual image.” On the right panel of Figure 2 there shows one example.The horizontal stripes can be clearly observed.We calculated the average value and the standard deviation in each row or column on the noise residual image.The bad pixels were excluded from the calculation.We merge all the row and column values from the 30 noise residual images together to make statistical plots(Figures 3 and 4).

    The result shows that the row-to-row variation(σ=3.59 ADU) is much larger than column-to-column variation(σ=0.05 ADU).It means the bias images contain spatially correlated noise, which is visually represented as horizontal stripes in noise residual image.The standard deviation in each column is comparable with the pixel value noise(which is based on Gaussian fitting of the pixel histogram,σ=4.46 ADU)and nominal“readout noise”(which is from the standard deviation of noise residual image, σ=4.77 ADU).Therefore the row-to-row fluctuation is the dominant component of the total readout noise.The standard deviation in each row is actually quite low (typically σ=2.84 ADU),implying that the detector may have the potential to decrease the overall readout noise to less than 3 ADU.

    We notice that the noise distribution is always non-Gaussian.The pixel-wise noise distribution is known to be somehow non-Gaussian.An internal technical report from the Gpixel Inc.(private communication) shows that the histogram of single pixel noise values has a long tail at the high end.A similar behavior was also reported in other CMOS models(e.g.,Ishida et al.2018).

    We calculated the pixel-to-pixel cross-correlation coefficient as a function of spatial shift(both x and y direction).The result is displayed in Figure 5.The result clearly shows the correlation in each row (i.e., x direction) is strong (>0.5).There are also many horizontal lines with weak but statistically significant correlations, positively or negatively.Hereafter in this paper we call such a feature as row noise (RN).We note that the pattern also shows some periodicity in the y direction,implying that the noise also contains characteristic frequency components.

    In order to confirm the existence of the characteristic spatial frequencies in the y direction, we calculated the spatial power spectra of the bias noise in Figure 6.The power spectrum is defined as the squared absolute amplitude of the Fourier transform component at given spatial frequency.We calculated the power spectrum in each row (or column) and derived the stacked spectrum with 20% trimmed mean value for all rows(or columns) in each individual bias image.Then we obtained 30 spectra(illustrated in gray lines)and derived the average of them(red line).In this way,we can reject any potential outliers(due to bad pixels) to get a robust estimation of the intrinsic power spectrum.

    Figure 3.Top right:pixel value histogram of all 30 noise residual images(fixed-pattern-corrected bias images).Bottom left:average value versus standard deviation in each row (the Pearson correlation coefficient is displayed).Top left: histogram of average (in each row)values.Bottom right: histogram of standard deviation (in each row) values.The red dashed line is Gaussian fitting.

    We find that the noise in each row(top panel of Figure 6)has a significant low frequency component, consistent with “1/f noise” (see e.g., Tian & El Gamal 2000).The noise in each column(bottom panel of Figure 6)has flat spectrum except that there are five characteristic frequencies (0.048, 0.133, 0,181,0.229, 0.410, in units of pixel-1).It may be caused by the power supply fluctuation, which is quite common in sCMOS detectors using column parallel readout architecture (i.e., each row is read in sequence, see Mikkonen 2016; Gilroy 2020).The readout could be affected by internal or external electrical circuits which has characteristic frequencies.However,we note that the white noise component is still the dominant noise source (the observed white noise component is equivalent to a noise source with rms of ≮4 ADU).

    2.3.Image Correction to Remove Row Noise

    The row noise, or more generally stripe noise, is commonly observed in images in remote sensing and/or taken by CMOS detectors (e.g., Simpson et al.1995; Rakwatin et al.2007;Gilroy 2020; Wei et al.2020).There are many developed algorithms to remove the stripe noise, for example, histogram matching with facet filter (Rakwatin et al.2007), midway histogram equalization (Tendero et al.2010), low-rank-based single-image decomposition (Chang et al.2016), group sparsity based regularization model (Chen et al.2017), wavelet deep neural network (Guan et al.2019), etc.It is also observed in some astronomical data (e.g., Schlawin et al.2020, for JWST NIRCam in time-series).The problem is not new and there are many possible solutions.Therefore it will be useful to investigate the properties of row noise corrected images as well.

    Figure 4.Same as Figure 3 but in columns.

    The origin of row noise in CMOS is complicated(Gilroy 2020).In this paper, only additive (or external signal independent) row noise is discussed (as only bias images are involved).8Actually, in another imaging test, we find there is row direction crosstalk,i.e.,signal dependent row noise,when there are many saturated pixels in a row.But it only appears when there is significantly saturated,i.e.,very bright,target.It does not affect the main conclusion of this paper,where we care more about faint objects with low signal-to-noise ratios.If we add the crosstalk in our simulation (presented in next section), the result will never be better.Our focus in this paper is not the correction method itself either.Therefore we use a very simple and straightforward method: subtracting the averaged empty pixel value in each row.

    To start, it is important to define the empty pixels, i.e., the image area not exposed to external light9Image area exposed to homogeneously distributed background light but not affected by astronomical objects may also be used for this purpose.However,in real application, such exercise has limitations: e.g., it is difficult for large extended objects (e.g., Magellanic Clouds) or over-crowded region (e.g., star cluster)where hardly any pixel is not contaminated by targets.In this paper,we do not discuss the algorithms to overcome such problems.Nevertheless, if used,it should be tested in real applications.On the other hand,our solution is the simplest way to avoid such problem..In case of HR9090BSI, the EB region (34 columns on each side) could be used for this purpose, after proper development in future(suggested by the Gpixel Inc.).Or alternatively,it is possible to mask some active pixels with opaque material to make an optically blind area at the edge of image.This option was proposed by the Gpixel Inc.(private communication) as well.In this paper,we assume the former option will be available.So we test the feasibility of this algorithm here.We use the leftmost 34 and rightmost 34 pixels of the active region (to mimic the size of EB)in each row to calculate the row average.We subtracted the average value of the 68 pixels in each row of the noise residual image.A similar algorithm is often used in CCD image processing, in which the reference empty area is from overscan instead of physical pixels (e.g., Wolf et al.2018).In Figure 7 the noise residual image before and after the correction are compared.Hereafter, we call these corrected images as “RN-corr” images.

    The row noise correction significantly reduces the readout noise (from 4.77 ADU to 2.89 ADU).The horizontal stripes become much less prominent visually.Similar power spectrum and correlation coefficient analysis were performed on RN-corr images for comparison (see Figures 8 and 9).

    Although there is still some small pixel-to-pixel correlation residual left in each row, the amplitude is greatly suppressed(Figure 8).The strength of the power spectrum in vertical direction(~6×104,Figure 9)is about one third of the original one (~2×105, Figure 6), becoming comparable with the horizontal one.Though as side effect,the power spectrum is no longer flat.Some characteristic frequencies are still there (e.g.,at 0.048 pixel-1),but their amplitudes are also greatly reduced.All these results demonstrate that the correction is very helpful.In next section, we will include the corrected images into our analysis as well.

    3.The Impact of Row Noise to Photometry

    3.1.The Simulated Noise Residual Images

    Here we use simulated images with galaxies and stars to better quantify the effect of row noise in astronomical application.We always simulate master bias corrected images(i.e.,object added to the noise residual image).

    For RN images, we use cutouts from original CMOS noise residual images.We randomly pick sub-regions of 196×196 pixels.Each random block must contain no more than 10 bad pixels and all bad pixels were set to zero.Random blocks should never overlap with each other so that the same noise pattern would not appear twice.There are totally 30446 random cutouts taken from 30 images.We also made cutouts for RN-corr images at the same positions for controlled comparison.

    For comparison, we also simulated pure Gaussian noise images.We made two sets of noise residual images with σ=4.77 ADU (Gauss-high) and σ=2.89 ADU (Gauss-low)respectively.They were used to compare with RN and RN-corr images respectively, in sense of nominal readout noise.We do not only observe the effect of overall readout noise reduction,but also the effect of RN itself.Figure 10 compares the different noise residual types.

    Figure 5.Pixel value correlation coefficient as a function of spatial shift.

    3.2.Real Galaxy Models

    Real galaxies were added to the simulated noise residual images.The image cutouts of real galaxies were directly taken from the website of the Legacy Survey10https://www.legacysurvey.org/(Dey et al.2019).15 galaxies were randomly picked.For simplicity and better signal-to-noise ratio, all of the three bands were merged together.Source detection and segmentation were performed to remove all nearby neighboring objects (except one galaxy α=187.5740, δ=73.0354 where a possible foreground star was kept to mimic the effect of possible confusion in real pipeline photometry).All pixels out of the target galaxy segment were set to zero.The cutout images were then normalized to unity total flux.Figure 11 shows all the galaxies we used.

    When adding the galaxy model to noise residual cutout, we randomly picked one model from the 15 models,multiplied by a random flux scale between 0.1 and 10.We use the same set of models for each type of noise residual image,to make sure that the difference only comes from noise residual.The left column of Figure 10 shows one example of such an image set.

    3.3.Sérsic Models

    Figure 6.Spatial power spectra of bias noise in rows (top)and columns (bottom).Grey lines are from individual bias images.Red lines are the average of them.In bottom panel, the values of five characteristic frequencies are noted.

    Sérsic model is a good way to parameterize the light distribution of elliptical galaxies, often used in simplified simulations (e.g., Holzschuh et al.2022).We use python package astropy to make 2D Sérsic model images,with random parameters listed in Table 1.Some examples are shown in Figure 12.Since the Sérsic model extends to infinite radius,we made a segment (consistent with an elliptical aperture) which included 90% of the total flux for further photometric measurement.Different from real galaxy models,Sérsic models have larger variability in galaxy size, probably more representative for field galaxies in a blind survey.The middle column of Figure 10 shows one example for different noise residual background.

    3.4.Star Models

    To simulate stars, we use 2D Gaussian model with σ=1.13 pixels (FWHM 2.66 pixels), which is the designed optical resolution (dynamic) of CSC.The model values were first assigned to a 0.05 pixel resolution sub-grid and then summed together to recover the true pixel value.Random amplitudes (between 10 and 103.5) and x/y center offsets(between -0.5 and 0.5) were used for each simulated star.Some examples were displayed in Figure 13.A circular aperture with radius of 4 pixels was used as segment for further photometric measurement.The right column of Figure 10 shows one example for different noise residual background.

    3.5.The Photometry and Basic Parameters

    We used the python package photutils v1.5 to measure several basic photometric parameters: position, elongation,position angle of the major axis, segment flux, Kron flux and circular aperture flux.

    The position measurement is very important in astrometry,i.e., determining the actual position of the object.Many other more complicated photometric parameters are sensitive to position measurement accuracy.The elongation is defined as the ratio between the length of major axis and the length of minor axis.It is often used as an indicator of the inclination angle of spiral galaxy.The position angle of the major axis is defined as the angle between the major axis and y-axis direction.It is also an important parameter to describe how the galaxy is positioned.Both the elongation and position angle are important morphological parameters,especially useful in weak lensing shear measurement (although the actual algorithm is more complicated).

    Figure 7.Before (left) and after (right) the row noise correction, displayed in the same scale.

    The segment flux is the sum of pixel values within the segment.It is a model-independent flux measurement,which is believed to be a robust estimation of the total flux.The Kron flux(Kron 1980)is another parameter to calculate the total flux.It sums up the flux within an elliptical aperture (i.e., assuming the galaxy shape can be represented by an elliptical).The parameters of the ellipse are determined from the first-order and second-order moment values.We use the default parameter in photutils.Comparing with segment flux, Kron flux is more sensitive to position and morphological measurements, but more reasonable in sense of measuring the total flux of a(generally symmetric)galaxy.It is widely used in many optical surveys (e.g., Guo et al.2013; Chambers et al.2016).The circular aperture flux was also calculated, which was preferred for stars.We always used the circular aperture with radius of 1.5 FWHM (3.99 pixels), which was sufficiently large to include most flux of a point source (more than 99.9% in our case) and was not too large to include too many noisy empty pixels.Most astrophysical studies require accurate flux measurement.It is important to check the reliability of the flux measurement.

    We used the fixed segment, i.e., the segment of the object(either galaxy or star) was fixed to the model image segment mentioned in previous subsections.If we use free segment or use resampled image,the conclusion does not change.We also note that Poisson resampling of the input flux was not used in our simulated images.In this way we make sure the difference only comes from the bias noise fluctuation within the same area.We will make a direct comparison between different noise types at the same nominal“readout noise”level and“signal-tonoise ratio.”

    3.6.Results

    We define the measurement error as the measurement difference between simulated image (model + noise) and pure model image, i.e., we treat the pure model result as “true value.”If the measurement is flux or elongation,the difference will be normalized by the “true value” and presented in percentage.

    Figure 8.Similar to Figure 5 but for row noise corrected images.

    Figure 9.Similar to Figure 6 but for row noise corrected images.

    Table 1 Random Parameters for Sérsic Models

    Figure 14 shows the measurement error as a function of average pixel flux in segment,which is defined as the total flux(in ADU) within the detection segment divided by the number of pixels of the segment.Here the position offset is defined as the absolute position offset with the sign of Y offset (i.e.,Each data point is binned from 1000 measurements.Error bars show the variation within the bin.Statistically, the mean value is accurate even for RN.But at any input brightness, the uncertainty of RN measurements are clearly worse than the Gaussian counterparts (note they have the same “readout noise”).Real galaxies and Sérsic models have similar results.

    Figure 15 shows the distribution of measurement error of all real galaxy images for six parameters.It shows the performance of RN is the worst for all parameters.The measurement error is especially large for Y position, segment flux and Kron flux.After correction, the performance of RN-corr is generally between Gauss-high and Gauss-low cases.It suggests that the row noise correction is not perfect but still very useful.

    The distribution is sometimes not symmetric.It may be related to the underlying algorithm and properties of the input galaxy templates.Such asymmetry does not affect our conclusion.Interestingly, the Y position uncertainty in RN case is significantly worse than X position.It is somehow expected, given the anisotropic nature of row noise.

    To better quantify the difference between different noise types, we plot the measurement uncertainty (the standard deviation of measurement error within the bin)as a function of average pixel flux in segment.Figures 16–18 show the results of all four noise types combined with three different models.

    The result is similar.Except the position angle of Sérsic images,RN is significantly worse than all other cases.RN-corr is much improved, comparable with Gaussian (between the high and low noise cases).The effect is much prominent in galaxies(including real galaxies and Sérsic models) than stars.It is probably due to the fact that the flux measurement apertures of stars are usually much smaller than that of galaxies.More pixels mean higher noise contribution from bias noise, and hence more vulnerable to row noise.It also clearly demonstrates that the correction is useful but the corrected result is still not as good as pure white noise.

    Figure 10.The same object with different noise residual as background (displayed in the same color scale).Left column: real galaxy model; middle column: Sérsic model; right column: star model.The images are zoomed in for better visibility.

    Figure 11.Real galaxy models used in simulation.Right ascension and decl.of the image center are shown for each galaxy (some galaxies are not centered).

    Figure 12.Example Sérsic models used in simulation.Blue line shows the segment which contains 90% of total flux.

    By comparing the flux uncertainties with the same input brightness,it is found that the uncertainty of RN case is about 3–10 times larger than Gauss-high for galaxies,and about 2 times larger for stars.In this paper,our test strictly controls the variables,leaving only the bias noise as the only parameter.So the effect of row noise we measured here can be directly translated into the increase of effective readout noise,by a factor of 2–10,depending on the target type.If we compare RN with Gauss-low, the difference is larger.

    Figure 13.Example Gaussian models used in simulation.All panels are in the same scale.The images are zoomed in.

    4.Discussion and Summary

    In this paper, we studied the noise properties of bias images taken from an HR9090BSI sCMOS detector.We found a strong pixel cross-correlation in each row, i.e., the row noise,and some minor fluctuations with fixed vertical spatial frequencies in the column direction.A simple row noise correction was applied to the real bias image,leading to a much lower readout noise and much smaller pixel-to-pixel correlation.We made a set of simulated galaxies and stars with different noise residual images as background.A series of photometry tests were performed.The results suggest that the row noise feature can significantly deteriorate the overall photometry performance, especially for galaxies.It effectively increases the readout noise by a factor of 2–10, depending on the type of object.

    Figure 14.Measurement error as a function of average pixel flux in detection segment(see text)for simulated stars.The blue data points are slightly shifted along Xaxis for better visibility.

    We demonstrate that the existence of row noise pattern, i.e.,pixel-to-pixel correlation, could greatly increase the photometric uncertainty.Such an effect is not simply due to the increase of nominal readout noise.In fact,at the same nominal readout noise level,the row noise case performs much worse than its Gaussian noise counterpart.It again proves the importance of pixel-to-pixel independence in astronomical high-precision photometry.It also makes a clear warning to anybody working on observation planning or camera design.The nominal“readout noise”may be misleading when the detector shows strong pixel-to-pixel correlation (no matter what the origin of such a correlation is).It may underestimate the required observation time for a given photometric accuracy, leading to an ambiguous observation result, or overestimate the performance of camera, leading to a failure of achieving the original design goal.

    We also note that, given the fact that the row noise is very common in sCMOS detectors, it is necessary to seriously consider the destriping method, either implemented in hardware or adopted at the data processing stage with efficient and robust algorithm.Our result suggests even the simplest destriping algorithm can make prominent improvement.Any future astronomical project using sCMOS detector should be aware of these facts.

    Our analysis also has limitations,although we note that some of them are out of the scope of this paper.First,our analysis of row noise is totally phenomenal.We treat the camera as a black box (due to some reasons).In order to fully understand the origin of row noise in this specific CMOS model,we will need more knowledge of the detector system, including the manufacturing details of the chip and the readout electronics.Second, we only consider the additive signal-independent row noise.Other multiplicative and/or signal-dependent noise may also affect the photometry.All these additional noise components will further worsen the photometry result if not properly corrected.In order identify these noise components,we will need more data, especially images from real astronomical observations.We note that in near future the“Earth 2.0” project (Song et al.2022) may give us a good chance.Third, this paper is focused on bias noise only.Our result is a good quantitative evaluator of the image readout performance.However, as a side effect, it makes the result looks too “exaggerated” when comparing to real situation.In practice,the bias noise usually only contributes a small fraction of the total photometric uncertainty.Many noise sources, such as Poisson noise,dark current,sky background and stray light,will make major contributions.Readers should be cautious when interpreting our results.To make a more realistic analysis of the photometric uncertainties,data from a more sophisticated simulation (e.g., Fang et al., in preparation, specifically designed for CSST11CSST simulation v2.0 release: https://csst-tb.bao.ac.cn/code/csst_sim/csst-simulation/-/tree/release_v2.0.) should be used.Finally, we do not discuss too much about the algorithm for row noise correction.When more and more CMOS data appears in the field of optical astronomy, it may be interesting to systematically study the effect of existing methods to high-precision photometry based on real observation images.

    Figure 15.The distribution of measurement errors for real galaxy models.Different bias noise types are displayed in different colors.

    Figure 16.The measurement uncertainty as a function of average pixel flux in detection segment for real galaxy models.Different bias noise types are displayed in different colors.

    Acknowledgments

    We thank the anonymous referee for her/his comments which help to improve this paper.

    This work is support by the National Key R&D Program of China No.2022YFF0503400.

    Figure 17.Same as Figure 16 but for simulated Sérsic images.

    The Legacy Surveys consist of three individual and complementary projects: the Dark Energy Camera Legacy Survey (DECaLS; Proposal ID #2014B-0404; PIs: David Schlegel and Arjun Dey), the Beijing-Arizona Sky Survey(BASS; NOAO Prop.ID #2015A-0801; PIs: Zhou Xu and Xiaohui Fan), and the Mayall z-band Legacy Survey (MzLS;Prop.ID #2016A-0453; PI: Arjun Dey).DECaLS, BASS and MzLS together include data obtained, respectively, at the Blanco telescope, Cerro Tololo Inter-American Observatory,NSF?s NOIRLab; the Bok telescope, Steward Observatory,University of Arizona; and the Mayall telescope, Kitt Peak National Observatory, NOIRLab.Pipeline processing and analyses of the data were supported by NOIRLab and the Lawrence Berkeley National Laboratory (LBNL).The Legacy Surveys project is honored to be permitted to conduct astronomical research on Iolkam Du’ag(Kitt Peak),a mountain with particular significance to the Tohono O’odham Nation.

    Figure 18.Same as Figure 16 but for simulated stars.

    NOIRLab is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.LBNL is managed by the Regents of the University of California under contract to the U.S.Department of Energy.

    This project used data obtained with the Dark Energy Camera (DECam), which was constructed by the Dark Energy Survey(DES)collaboration.Funding for the DES Projects has been provided by the U.S.Department of Energy, the U.S.National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom,the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign,the Kavli Institute of Cosmological Physics at the University of Chicago, Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University,Financiadora de Estudos e Projetos, Fundacao Carlos Chagas Filho de Amparo, Financiadora de Estudos e Projetos,Fundacao Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cientifico e Tecnologico and the Ministerio da Ciencia,Tecnologia e Inovacao, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey.The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenossische Technische Hochschule (ETH) Zurich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciencies de l’Espai (IEEC/CSIC),the Institut de Fisica d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig Maximilians Universitat Munchen and the associated Excellence Cluster Universe, the University of Michigan, NSF?s NOIRLab, the University of Nottingham, the Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University.

    BASS is a key project of the Telescope Access Program(TAP), which has been funded by the National Astronomical Observatories,the Chinese Academy of Sciences(the Strategic Priority Research Program “The Emergence of Cosmological Structures” grant # XDB09000000), and the Special Fund for Astronomy from the Ministry of Finance.The BASS is also supported by the External Cooperation Program of Chinese Academy of Sciences (grant # 114A11KYSB20160057), and the National Natural Science Foundation of China (grants #12120101003 and # 11433005).

    The Legacy Survey team makes use of data products from the Near-Earth Object Wide-field Infrared Survey Explorer(NEOWISE), which is a project of the Jet Propulsion Laboratory/California Institute of Technology.NEOWISE is funded by the National Aeronautics and Space Administration.

    The Legacy Surveys imaging of the DESI footprint is supported by the Director, Office of Science, Office of High Energy Physics of the U.S.Department of Energy under Contract No.DE-AC02-05CH1123, by the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility under the same contract; and by the U.S.National Science Foundation, Division of Astronomical Sciences under Contract No.AST-0950945 to NOAO.

    ORCID iDs

    Li Shao https://orcid.org/0000-0003-2015-777X

    Hu Zhan https://orcid.org/0000-0003-1718-6481

    一区二区三区免费毛片| 人人妻,人人澡人人爽秒播| 国产免费男女视频| 中文字幕精品亚洲无线码一区| 男人的好看免费观看在线视频| 亚洲精华国产精华精| 欧美黑人巨大hd| 亚洲精华国产精华精| 十八禁人妻一区二区| 嫁个100分男人电影在线观看| 国产成人影院久久av| 成人av在线播放网站| 少妇人妻精品综合一区二区 | 老司机在亚洲福利影院| 日本成人三级电影网站| 国内精品美女久久久久久| 亚洲自拍偷在线| 亚洲18禁久久av| 久久亚洲真实| 久9热在线精品视频| 最新在线观看一区二区三区| 午夜福利成人在线免费观看| 夜夜躁狠狠躁天天躁| h日本视频在线播放| 无遮挡黄片免费观看| 亚洲精品国产精品久久久不卡| 日韩高清综合在线| 99久国产av精品| 最新美女视频免费是黄的| 国产激情偷乱视频一区二区| 国产精品一区二区三区四区免费观看 | 成人av一区二区三区在线看| 国产精品久久久久久精品电影| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清专用| 在线播放国产精品三级| 小蜜桃在线观看免费完整版高清| 日韩国内少妇激情av| 两个人的视频大全免费| 日韩欧美一区二区三区在线观看| 搡老妇女老女人老熟妇| 三级毛片av免费| 亚洲在线观看片| 色老头精品视频在线观看| 悠悠久久av| 男女下面进入的视频免费午夜| 中文字幕熟女人妻在线| 一级作爱视频免费观看| 观看免费一级毛片| h日本视频在线播放| 五月伊人婷婷丁香| 久久精品影院6| 亚洲av第一区精品v没综合| 毛片女人毛片| 欧美中文日本在线观看视频| xxx96com| 国产探花极品一区二区| 老汉色∧v一级毛片| 国产高清激情床上av| 国产精品精品国产色婷婷| bbb黄色大片| 亚洲在线观看片| 老熟妇乱子伦视频在线观看| 亚洲真实伦在线观看| 亚洲国产欧洲综合997久久,| 亚洲人成伊人成综合网2020| 色老头精品视频在线观看| av欧美777| 亚洲国产精品成人综合色| 人妻久久中文字幕网| 国产精品一区二区三区四区免费观看 | 欧美日韩福利视频一区二区| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 乱人视频在线观看| 色av中文字幕| 成年女人看的毛片在线观看| 国产亚洲欧美在线一区二区| 国产成人欧美在线观看| 在线观看66精品国产| 高清日韩中文字幕在线| 禁无遮挡网站| 日韩大尺度精品在线看网址| 最好的美女福利视频网| 夜夜爽天天搞| 国产 一区 欧美 日韩| 在线天堂最新版资源| 宅男免费午夜| 麻豆国产av国片精品| 免费一级毛片在线播放高清视频| 亚洲自拍偷在线| 亚洲av熟女| 精品99又大又爽又粗少妇毛片 | 国产伦在线观看视频一区| 亚洲自拍偷在线| 两个人视频免费观看高清| 精品99又大又爽又粗少妇毛片 | 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 亚洲av免费高清在线观看| 美女被艹到高潮喷水动态| 国内少妇人妻偷人精品xxx网站| 琪琪午夜伦伦电影理论片6080| 国产中年淑女户外野战色| 韩国av一区二区三区四区| 日韩欧美精品v在线| 欧美zozozo另类| 国产精品 欧美亚洲| 全区人妻精品视频| 午夜福利免费观看在线| 亚洲avbb在线观看| 在线观看免费视频日本深夜| 亚洲无线在线观看| 久久伊人香网站| 国产视频一区二区在线看| 五月玫瑰六月丁香| 国产97色在线日韩免费| 麻豆成人av在线观看| 99精品欧美一区二区三区四区| 麻豆成人午夜福利视频| 偷拍熟女少妇极品色| 成年女人毛片免费观看观看9| 激情在线观看视频在线高清| 国产欧美日韩一区二区三| 亚洲最大成人手机在线| 看黄色毛片网站| 村上凉子中文字幕在线| 69av精品久久久久久| 亚洲第一电影网av| 97碰自拍视频| 欧美在线一区亚洲| 国产激情偷乱视频一区二区| 日韩成人在线观看一区二区三区| 最近最新免费中文字幕在线| 麻豆一二三区av精品| 国产三级在线视频| 香蕉丝袜av| 国产97色在线日韩免费| 变态另类成人亚洲欧美熟女| 亚洲狠狠婷婷综合久久图片| 一卡2卡三卡四卡精品乱码亚洲| 老司机在亚洲福利影院| 色综合站精品国产| 日本在线视频免费播放| 大型黄色视频在线免费观看| 一区二区三区免费毛片| 亚洲中文日韩欧美视频| 欧美午夜高清在线| 成人性生交大片免费视频hd| 欧美日韩精品网址| 九九热线精品视视频播放| www国产在线视频色| 国语自产精品视频在线第100页| 欧美黑人欧美精品刺激| 久久久国产成人精品二区| 一个人看的www免费观看视频| 国产成年人精品一区二区| 亚洲人与动物交配视频| 国产乱人伦免费视频| 国产成年人精品一区二区| 一级黄片播放器| 日韩精品中文字幕看吧| 中文字幕人妻丝袜一区二区| 99久久综合精品五月天人人| 精品国产超薄肉色丝袜足j| 久久精品综合一区二区三区| 99久久成人亚洲精品观看| 悠悠久久av| 18禁美女被吸乳视频| 天堂影院成人在线观看| 男人舔奶头视频| 69人妻影院| 亚洲天堂国产精品一区在线| 成人三级黄色视频| av女优亚洲男人天堂| 国产午夜精品论理片| 九九久久精品国产亚洲av麻豆| 午夜精品久久久久久毛片777| 搡老熟女国产l中国老女人| 国产高清激情床上av| av女优亚洲男人天堂| 伊人久久精品亚洲午夜| 女生性感内裤真人,穿戴方法视频| 久久精品国产99精品国产亚洲性色| 三级男女做爰猛烈吃奶摸视频| 99在线人妻在线中文字幕| 狂野欧美激情性xxxx| 国产精品亚洲一级av第二区| 久久久国产精品麻豆| 亚洲人成网站在线播| 欧美日韩黄片免| 观看美女的网站| 制服人妻中文乱码| 欧美在线黄色| 51国产日韩欧美| 成人高潮视频无遮挡免费网站| 97超视频在线观看视频| 少妇的丰满在线观看| 夜夜看夜夜爽夜夜摸| 亚洲狠狠婷婷综合久久图片| 最近最新中文字幕大全电影3| 男人舔奶头视频| 中亚洲国语对白在线视频| 成人一区二区视频在线观看| 他把我摸到了高潮在线观看| 成年免费大片在线观看| 老熟妇仑乱视频hdxx| 久久6这里有精品| 欧美一级a爱片免费观看看| 亚洲乱码一区二区免费版| 日韩高清综合在线| 色老头精品视频在线观看| 九色成人免费人妻av| 国产精品久久久久久久久免 | 欧美日本亚洲视频在线播放| 国产真实乱freesex| 国产精品乱码一区二三区的特点| 国产精品三级大全| xxx96com| 国产免费av片在线观看野外av| 精品无人区乱码1区二区| 黄色成人免费大全| 午夜亚洲福利在线播放| 一级毛片高清免费大全| 99精品欧美一区二区三区四区| 免费看日本二区| 999久久久精品免费观看国产| 一个人免费在线观看电影| 亚洲欧美日韩高清在线视频| 国产精品av视频在线免费观看| 国产美女午夜福利| 午夜福利在线观看免费完整高清在 | 最后的刺客免费高清国语| 18+在线观看网站| 国产免费av片在线观看野外av| 非洲黑人性xxxx精品又粗又长| 国产精品亚洲一级av第二区| 免费搜索国产男女视频| 久久精品国产亚洲av涩爱 | 亚洲成人久久爱视频| 波野结衣二区三区在线 | 人妻丰满熟妇av一区二区三区| 亚洲乱码一区二区免费版| 禁无遮挡网站| 日本成人三级电影网站| 欧美在线一区亚洲| 国产免费男女视频| 最近视频中文字幕2019在线8| 天堂影院成人在线观看| 亚洲片人在线观看| 欧美一级毛片孕妇| 欧美+亚洲+日韩+国产| 最后的刺客免费高清国语| 特级一级黄色大片| 乱人视频在线观看| 一区二区三区国产精品乱码| 国产av不卡久久| 亚洲国产精品久久男人天堂| 欧美又色又爽又黄视频| 青草久久国产| 免费av不卡在线播放| 深爱激情五月婷婷| 一区二区三区国产精品乱码| 变态另类丝袜制服| 高潮久久久久久久久久久不卡| eeuss影院久久| 中文字幕精品亚洲无线码一区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久人妻蜜臀av| 国产91精品成人一区二区三区| 久久精品国产99精品国产亚洲性色| ponron亚洲| 午夜老司机福利剧场| 国产亚洲av嫩草精品影院| 国内精品一区二区在线观看| 深爱激情五月婷婷| 91在线观看av| 欧美一区二区国产精品久久精品| 免费在线观看日本一区| 99精品欧美一区二区三区四区| 精品久久久久久成人av| 国产成人系列免费观看| 女人被狂操c到高潮| 国内少妇人妻偷人精品xxx网站| 俺也久久电影网| 岛国在线观看网站| 天美传媒精品一区二区| 国产成+人综合+亚洲专区| 有码 亚洲区| 天堂影院成人在线观看| 国产蜜桃级精品一区二区三区| 狂野欧美激情性xxxx| 日韩精品青青久久久久久| 免费观看的影片在线观看| 好看av亚洲va欧美ⅴa在| 国产黄a三级三级三级人| 在线播放国产精品三级| 九九久久精品国产亚洲av麻豆| 老司机午夜福利在线观看视频| 波多野结衣高清无吗| 亚洲精品久久国产高清桃花| 国产精品98久久久久久宅男小说| а√天堂www在线а√下载| 国产亚洲精品综合一区在线观看| 综合色av麻豆| 成人精品一区二区免费| 欧美日韩中文字幕国产精品一区二区三区| 精品国内亚洲2022精品成人| 亚洲av成人av| 亚洲欧美日韩高清在线视频| 十八禁网站免费在线| 久久久久精品国产欧美久久久| 精华霜和精华液先用哪个| 国产日本99.免费观看| 99久久久亚洲精品蜜臀av| 男插女下体视频免费在线播放| 亚洲av五月六月丁香网| 啪啪无遮挡十八禁网站| 欧美日韩精品网址| 男女视频在线观看网站免费| 欧美不卡视频在线免费观看| 亚洲欧美日韩卡通动漫| 国内精品一区二区在线观看| 欧美激情在线99| 欧美日韩国产亚洲二区| 制服丝袜大香蕉在线| 在线播放无遮挡| 亚洲精品国产精品久久久不卡| 天堂网av新在线| 丰满人妻一区二区三区视频av | 国产成人系列免费观看| 蜜桃亚洲精品一区二区三区| 亚洲国产精品999在线| 老司机午夜十八禁免费视频| 男人舔奶头视频| 女警被强在线播放| 又粗又爽又猛毛片免费看| 在线观看日韩欧美| 成人国产一区最新在线观看| 99热6这里只有精品| 成人av一区二区三区在线看| 国产精品美女特级片免费视频播放器| 好男人电影高清在线观看| 亚洲国产色片| ponron亚洲| 精品一区二区三区视频在线 | 国产69精品久久久久777片| 久久久久亚洲av毛片大全| 国产精品美女特级片免费视频播放器| 嫩草影视91久久| bbb黄色大片| 日韩人妻高清精品专区| 国产淫片久久久久久久久 | 91av网一区二区| e午夜精品久久久久久久| 又粗又爽又猛毛片免费看| 日韩有码中文字幕| 神马国产精品三级电影在线观看| 97超级碰碰碰精品色视频在线观看| 国产乱人视频| 一级黄片播放器| 久久久久亚洲av毛片大全| 国产中年淑女户外野战色| 此物有八面人人有两片| 久久久久性生活片| 国产午夜精品久久久久久一区二区三区 | 亚洲午夜理论影院| 免费观看人在逋| 国产成+人综合+亚洲专区| 超碰av人人做人人爽久久 | 国产三级黄色录像| 窝窝影院91人妻| 亚洲内射少妇av| 欧美最黄视频在线播放免费| 久久精品综合一区二区三区| 国产亚洲精品av在线| 精品久久久久久久毛片微露脸| 看片在线看免费视频| 久久香蕉国产精品| 9191精品国产免费久久| 18禁美女被吸乳视频| www日本在线高清视频| 99国产精品一区二区蜜桃av| 又紧又爽又黄一区二区| 国产精华一区二区三区| 一区二区三区高清视频在线| 狂野欧美白嫩少妇大欣赏| 亚洲专区中文字幕在线| 夜夜夜夜夜久久久久| 国产午夜精品久久久久久一区二区三区 | 又黄又粗又硬又大视频| 亚洲aⅴ乱码一区二区在线播放| 757午夜福利合集在线观看| 丝袜美腿在线中文| 午夜两性在线视频| 1024手机看黄色片| 伊人久久大香线蕉亚洲五| 午夜激情福利司机影院| 国产成人系列免费观看| 久久伊人香网站| 国产91精品成人一区二区三区| 国产精品一区二区免费欧美| 男女床上黄色一级片免费看| av专区在线播放| 亚洲av电影在线进入| 亚洲国产日韩欧美精品在线观看 | 女同久久另类99精品国产91| 麻豆国产av国片精品| 精品电影一区二区在线| 在线看三级毛片| 欧美日本亚洲视频在线播放| 亚洲人成网站在线播| 久久久成人免费电影| 午夜久久久久精精品| 99热只有精品国产| 午夜视频国产福利| 在线观看美女被高潮喷水网站 | 亚洲中文字幕一区二区三区有码在线看| 日韩亚洲欧美综合| 国产av在哪里看| 女人被狂操c到高潮| 久久精品国产自在天天线| 亚洲第一电影网av| 天堂√8在线中文| 丁香六月欧美| 天美传媒精品一区二区| 变态另类成人亚洲欧美熟女| 久9热在线精品视频| 日韩欧美国产一区二区入口| 特大巨黑吊av在线直播| www.色视频.com| 久久久精品大字幕| 露出奶头的视频| 亚洲一区二区三区不卡视频| 亚洲国产精品成人综合色| 女同久久另类99精品国产91| 欧美日韩福利视频一区二区| 精华霜和精华液先用哪个| 色av中文字幕| 美女高潮喷水抽搐中文字幕| 午夜日韩欧美国产| 99国产精品一区二区蜜桃av| 精品国产超薄肉色丝袜足j| 操出白浆在线播放| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 国产91精品成人一区二区三区| 免费在线观看日本一区| 欧美一区二区亚洲| 国语自产精品视频在线第100页| 免费观看的影片在线观看| 亚洲专区中文字幕在线| 99久国产av精品| 色综合亚洲欧美另类图片| 69人妻影院| av天堂中文字幕网| 久久久国产成人免费| 人人妻,人人澡人人爽秒播| 婷婷丁香在线五月| 深夜精品福利| 欧美高清成人免费视频www| 日本a在线网址| 一级黄色大片毛片| 国产成人aa在线观看| 欧美av亚洲av综合av国产av| 成人18禁在线播放| 国产精品三级大全| 欧美精品啪啪一区二区三区| 俺也久久电影网| 国产乱人视频| 人妻久久中文字幕网| 丰满人妻一区二区三区视频av | 日本精品一区二区三区蜜桃| 亚洲国产日韩欧美精品在线观看 | 国产高潮美女av| 欧美3d第一页| 久久中文看片网| 一级毛片高清免费大全| 日韩av在线大香蕉| 1000部很黄的大片| 久久草成人影院| 国产精品久久久人人做人人爽| 免费人成在线观看视频色| 超碰av人人做人人爽久久 | 日日干狠狠操夜夜爽| 国产精品久久久久久久电影 | 男人舔女人下体高潮全视频| 日韩大尺度精品在线看网址| 天美传媒精品一区二区| 亚洲电影在线观看av| 免费在线观看亚洲国产| 久久99热这里只有精品18| 黑人欧美特级aaaaaa片| 偷拍熟女少妇极品色| 在线十欧美十亚洲十日本专区| 美女cb高潮喷水在线观看| 久久久久亚洲av毛片大全| 亚洲人成网站在线播| 国产欧美日韩精品亚洲av| 久久精品亚洲精品国产色婷小说| 一本久久中文字幕| 婷婷精品国产亚洲av在线| 日本三级黄在线观看| 19禁男女啪啪无遮挡网站| 国产成+人综合+亚洲专区| 99riav亚洲国产免费| 神马国产精品三级电影在线观看| 亚洲精品乱码久久久v下载方式 | 国产精品一区二区三区四区免费观看 | 欧美激情久久久久久爽电影| 男人和女人高潮做爰伦理| 婷婷精品国产亚洲av在线| 色尼玛亚洲综合影院| 色综合站精品国产| 日韩精品中文字幕看吧| 在线观看日韩欧美| 午夜a级毛片| 啦啦啦韩国在线观看视频| 亚洲欧美一区二区三区黑人| 国内精品美女久久久久久| 性欧美人与动物交配| 欧美黄色片欧美黄色片| 特级一级黄色大片| 一本一本综合久久| 亚洲成人中文字幕在线播放| 人妻丰满熟妇av一区二区三区| 性色av乱码一区二区三区2| 久久香蕉精品热| 亚洲成a人片在线一区二区| 久久精品国产清高在天天线| 91九色精品人成在线观看| 日韩欧美精品免费久久 | 超碰av人人做人人爽久久 | 嫩草影院入口| 91字幕亚洲| 成人18禁在线播放| 最近视频中文字幕2019在线8| 欧美日韩综合久久久久久 | 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 狂野欧美激情性xxxx| 亚洲国产欧洲综合997久久,| 757午夜福利合集在线观看| av女优亚洲男人天堂| 一a级毛片在线观看| 免费一级毛片在线播放高清视频| 一个人观看的视频www高清免费观看| 午夜福利高清视频| 亚洲久久久久久中文字幕| 欧美一区二区国产精品久久精品| 亚洲专区国产一区二区| 在线免费观看不下载黄p国产 | 国产一区在线观看成人免费| 可以在线观看毛片的网站| 久久亚洲真实| av天堂在线播放| xxx96com| 两个人视频免费观看高清| 老熟妇仑乱视频hdxx| 亚洲精品乱码久久久v下载方式 | 黄色成人免费大全| av片东京热男人的天堂| 丰满人妻一区二区三区视频av | 成人特级av手机在线观看| 在线a可以看的网站| 成年免费大片在线观看| 亚洲国产欧洲综合997久久,| 国产精品国产高清国产av| 51国产日韩欧美| 久久精品国产清高在天天线| 国产免费男女视频| 国产真实乱freesex| 久久精品国产亚洲av香蕉五月| 久久精品国产综合久久久| a在线观看视频网站| 亚洲国产精品999在线| 欧美最黄视频在线播放免费| 日本 欧美在线| 人人妻人人澡欧美一区二区| 天堂网av新在线| 99久久无色码亚洲精品果冻| 亚洲精品在线观看二区| 亚洲在线自拍视频| 久久6这里有精品| 波多野结衣高清无吗| 真实男女啪啪啪动态图| 91在线精品国自产拍蜜月 | 熟女少妇亚洲综合色aaa.| 国产一区二区亚洲精品在线观看| 波野结衣二区三区在线 | 亚洲 国产 在线| 制服人妻中文乱码| 日日干狠狠操夜夜爽| 中文字幕av成人在线电影| 99久久成人亚洲精品观看| 制服丝袜大香蕉在线| 国内少妇人妻偷人精品xxx网站| 级片在线观看| 亚洲欧美一区二区三区黑人| 又粗又爽又猛毛片免费看| 久久国产精品人妻蜜桃| 99国产精品一区二区蜜桃av| 精品人妻1区二区| 免费在线观看日本一区| 欧美成人一区二区免费高清观看| 日韩成人在线观看一区二区三区| 国产av麻豆久久久久久久| a在线观看视频网站| 3wmmmm亚洲av在线观看| 欧美高清成人免费视频www| 男人和女人高潮做爰伦理| 国产伦精品一区二区三区视频9 | 女生性感内裤真人,穿戴方法视频| 女同久久另类99精品国产91|