载荷探测器像元满阱参数测试

载荷探测器像元满阱参数测试 Measurement of pixel full-well parameters of detector used on spaceborne spectrometer 1 2 first-author 常　振（1988-），男， 安徽合肥人，博士，副研究员，2012年于安徽师范大学获得学士学位，2017年于中国科学技术大学获得博士学位，主要从事CCD成像技术方面的研究。E-mail：zhchang@aiofm.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965955&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965963&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965960&type= 22.01333237 32.00400543 常　振 （1988-），男， 安徽合肥人，博士，副研究员，2012年于安徽师范大学获得学士学位，2017年于中国科学技术大学获得博士学位，主要从事CCD成像技术方面的研究。E-mail： zhchang@aiofm.ac.cn zhchang@aiofm.ac.cn 2 corresp E-mail： yuwang@aiofm.ac.cn E-mail： yuwang@aiofm.ac.cn 王　煜（1971-），男，北京人，博士，研究员， 1993年于哈尔滨工业大学获得学士学位，2005年于英国斯特莱斯克莱德大学获得博士学位，现为安徽大学物质科学与信息技术研究院环境物理所副所长，主要从事航空、航天遥感设备的研发工作。E-mail：yuwang@aiofm.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965966&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965972&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965969&type= 22.01332855 32.00399780 王　煜 （1971-），男，北京人，博士，研究员， 1993年于哈尔滨工业大学获得学士学位，2005年于英国斯特莱斯克莱德大学获得博士学位，现为安徽大学物质科学与信息技术研究院环境物理所副所长，主要从事航空、航天遥感设备的研发工作。E-mail： yuwang@aiofm.ac.cn yuwang@aiofm.ac.cn 1 1 1 1 中国科学院 合肥物质科学研究院 安徽光学精密机械研究所 环境光学与技术重点实验室，安徽 合肥230031 中国科学院 合肥物质科学研究院 安徽光学精密机械研究所 环境光学与技术重点实验室，安徽 合肥230031 Key Laboratory of Environmental Optical and Technology， Anhui Institute of Optics and Fine Mechanics， Hefei Institutes of Physical Science， Chinese Academy of Sciences， Hefei 230031， China Key Laboratory of Environmental Optical and Technology， Anhui Institute of Optics and Fine Mechanics， Hefei Institutes of Physical Science， Chinese Academy of Sciences， Hefei 230031， China 安徽大学 物质科学与信息技术研究院 信息材料与智能感知安徽省实验室， 安徽 合肥 230039 安徽大学 物质科学与信息技术研究院 信息材料与智能感知安徽省实验室， 安徽 合肥 230039 Information Materials and Intelligent Sensing Laboratory of Anhui Province， Institutes of Physical Science and Information Technology， Anhui University， Hefei 230039， China Information Materials and Intelligent Sensing Laboratory of Anhui Province， Institutes of Physical Science and Information Technology， Anhui University， Hefei 230039， China 在CCD成像电路的研发过程中，为了确保CCD达到最优的满阱性能，需要不断调整CCD驱动信号参数并进行满阱测试，该过程通常需要重复几十甚至上百次。通用的光子转移曲线法需要积分球等设备搭建平场光源，系统复杂，满阱测试效率较低。提出了一种满阱测试方法——LED点测试法，该方法通过光子转移曲线法获取系统增益后，仅通过搭建LED点光源配合普通民品镜头即可实现CCD满阱测试。基于载荷CCD成像电路，使用LED点测试法进行测试。结果表明，载荷CCD的满阱电子数可达817.013ke - ，误差不大于0.643%。针对LED点测试的同组数据采用传统满阱测试方法进行测试，对比结果表明，LED点测试法和传统测试方法相对光子转移曲线法的误差分别为0.039 7%和1.9%。LED点测试法可用作简易条件下快速测量CCD满阱的通用方法，能够提升CCD成像电路的研发效率。 During the development of a charge-coupled device （CCD） imaging circuit， it is necessary to adjust the parameters of the CCD’s driving signal and repeatedly test the full-well parameters dozens or even hundreds of times to ensure the optimal full-well performance of the CCD. To achieve this， the general method employed， known as the photon transfer curve method， requires integrating a sphere and other equipment to create a flat-field light source system. However， the resulting system is complex， and hence， it is unsuitable for application in a full-well testing method during the development of an imaging system. In this paper， a novel full-well testing method， a light-emitting diode （LED） point test method， is proposed. Accordingly， CCD full-well testing can be easily realized by developing an LED point light source and common civilian lens. The results reveal that the number of full-well electrons of a spaceborne CCD can reach 817.013ke - ， and the error can be less than 0.643%. The results also indicate that the errors resulting from the LED point test method and the traditional method are， respectively， 0.0397% and 1.9% compared with the errors resulting from the photon transfer curve method. In summary， the LED point test method can be used as a general method to rapidly measure the full-well parameters of a CCD under simple conditions， and it can enhance the efficiency of CCD imaging circuit development. 电荷耦合器件 满阱容量 光子转移曲线 LED点测试法 数据拟合 charge-coupled device（CCD） full well capacity photon transfer curve LED-point test method data fitting 1 引 言 1 星载大气痕量气体差分吸收光谱仪（Environment Monitoring Instrument， EMI）是我国首台用于大气污染监测的星载高光谱载荷，数据产品精度优于97%，在空间分辨率和反演噪声上都显著高于OMI（Ozone Monitoring Instrument）卫星产品，相关性系数为0.96。该载荷填补了国产卫星探测全球及区域大气污染物分布的空白，可为大气污染防治、全球气候变化管理等提供有力的数据支撑 ［ 1 1 - 2 2 1-2 ］ 。为了进一步提升国产卫星探测全球大气污染分布的能力，实现更高空间分辨率的对地探测，新型号载荷使用了全新探测器，对地观测空间分辨率达到7公里，临边对地空间分辨率为2公里，这一分辨率精度对污染源位置辨识、排放通量计算等有着重要作用 ［ 3 3 ］ 。 新型号载荷使用深势阱的科学级电荷耦合器件（Charge-Coupled Device， CCD），可提供相对于EMI更高信噪比数据，从而提供更精确的污染物浓度分布数据。目前，国内尚无使用该探测的案例，作为国内首次使用该探测器的载荷，在成像电路的开发过程中需要对其电子学关键性能——满阱容量（Full Well Capacity， FWC），即满阱电子数进行测试及验证，以确认探测器电子学系统处于最佳状态。虽然CCD在出厂时标注了满阱参数，但是由于CCD驱动电路的设计参数如水平及垂直转移时钟的频率、驱动电流等均会影响CCD的满阱性能，因此在成像电路的研发过程中，需要不断调整参数并测试CCD的满阱，以寻找最优满阱性能对应的参数组合。参数的调整过程中有时需要进行几十次甚至上百次满阱测试，因此需要简单可行的满阱测试方法以提升CCD成像电路的研发效率。 CCD成像系统满阱性能的测试通常采用光子转移曲线（Photon Transfer Curve， PTC）法，该方法可通过简单的方法获取CCD的读出噪声、系统增益、满阱电子数和像元响应非一致性（Photo Response Non Uniformity， PRNU）等重要参数。PTC是一种起源于国外的较为公认的CCD多项参数测试方法。国内的有关科学级CCD成像技术的研究多基于卫星载荷项目等 ［ 4 4 - 5 5 4-5 ］ ，民用及实验室大都采用市场现有的CCD相机搭建成像系统 ［ 6 6 - 8 8 6-8 ］ ，因此，PTC法是一种应用广泛且可靠的CCD成像系统满阱测试方法 ［ 9 9 - 10 10 9-10 ］ 。当CCD的满阱输出未达到成像系统模数转换（ADC）的量程时，可采用EMVA1288标准定义计算满阱电子数，即根据DN值方差和辐射量关系曲线的3%或5%线性拟合偏差作为饱和点，但是这一人为选择因素很可能给测试结果带来较大偏差。 由于PTC方法测试系统复杂，在CCD成像电路的研发过程中，每次测试使用PTC方法进行满阱测试的效率较低。本文构建了一种简易的测试方法——LED点测试法。该方法使用LED光源作为目标光源获取数据，结合PTC法获取的增益参数，计算得到CCD的满阱容量。由于CCD驱动信号时钟参数的改变不会影响成像系统增益，因此研发过程中只需进行一次PTC测试获取系统增益，后续驱动信号的调整仅使用LED点测试法即可获得探测器满阱数值。通过实验验证了LED点测试法的可行性和准确性，与传统的CCD成像系统满阱测试法相比，LED点测试法基于对像元溢出物理过程的分析，更符合像元满阱的物理学定义，更接近像元满阱的真实情况，可信度更高。该方法可作为简易条件下快速测量CCD满阱的通用性方法，大幅提升CCD成像电路的研发效率。探测器满阱性能的测试结果可为后期载荷数据反演提供有力的数据支撑。 2 原 理 1 2.1　PTC方法 2 CCD的基本结构是由许多MOS管组合而成，MOS管内部有电势垒可保存电荷，保存电荷的多少与晶体掺杂浓度、静态工作的供电电压等有关 ［ 11 11 ］ 。光线照射到CCD的感光区域可产生光生电子，在时序电压的驱动下，产生的电子通过MOS管内部的载流子耦合可实现转移。同一光强CCD产生的光生电子数并非常量，而是呈现一定的随机分布，单位时间内入射光产生的光生电子数可认为服从泊松分布。根据泊松分布理论，CCD收集的电子数目的平均值 μ 和噪声 σ 的关系如下： <math display="block" id="M1"><mi>σ</mi><mo>=</mo><mroot><mrow><mi>μ</mi></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965978&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965975&type= 11.85333347 4.91066647 . 对于深势阱的科学级CCD，由于散粒噪声随光强的增大而增加，因此在一定光强下，可将它作为系统的主要噪声，忽略暗电流噪声、读出噪声、电源噪声等不随光强变化的噪声 ［ 12 12 ］ 。根据CCD的工作原理，CCD像元搜集到的电荷可转换为电压信号，电压信号经读出电路、偏置电路和采样电路转换为数字信号。CCD成像系统将电子数转换为DN值的过程包括两个重要参数：电子数转换成电压信号的增益 η 和电压经读出电路及采样电路转换为DN值的增益 k 。假定整个系统电子数到DN值的变换是线性的，则此处可定义系统增益 <math id="M2"><mi>g</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965984&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965981&type= 1.94733346 3.64066648 ，单位是DN/e - ，表示像元搜集的平均电子数与DN值的关系。于是有： <math display="block" id="M3"><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">N</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">g</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub><mo>=</mo><mi>g</mi><mo>·</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">E</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">g</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub><mo>+</mo><mi>B</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965990&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965987&type= 35.22133255 4.57200003 ， 其中： <math id="M4"><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">N</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">g</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965996&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32965993&type= 11.00666618 4.57200003 表示像元DN值的平均值， <math id="M5"><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">E</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">g</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966004&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966002&type= 9.31333351 4.57200003 表示像元搜集到的电子数的平均值， B 表示系统偏置。同理，像元搜集到的电子数的噪声部分 <math id="M6"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">E</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966010&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966008&type= 13.20800018 3.72533321 也符合系统增益的定义： <math display="block" id="M7"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">N</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi></mrow></msub><mo>=</mo><mi>g</mi><mo>·</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">E</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966017&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966012&type= 35.64466476 4.57200003 ， 其中： <math id="M8"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">N</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966023&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966020&type= 14.81666565 3.72533321 表示像元搜集到的电子数在单位时间内带来的DN值波动，即散粒噪声，与 式（2） 式（2） 不同，噪声部分是不包括偏置的。获得的DN值噪声中还包括电路的读出噪声、电源噪声等，这些噪声与CCD获取的光强无关，可统计为读出噪声 <math id="M9"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966029&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966026&type= 10.66800022 3.72533321 ，根据噪声的叠加原理，系统总噪声 <math id="M10"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966035&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966032&type= 5.50333309 3.72533321 可表示为： <math display="block" id="M11"><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>=</mo><msup><mrow><mi>g</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">E</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">h</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966041&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966038&type= 42.50266647 4.82600021 . 结合式（1）~ 式（4） 式（4） ，有： <math display="block" id="M12"><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>=</mo><mi>g</mi><mo>·</mo><mfenced separators="|"><mrow><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">N</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">g</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub><mo>-</mo><mi>B</mi></mrow></mfenced><mo>+</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi><mo>-</mo><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966047&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966044&type= 50.71533585 4.82600021 . 式（5） 式（5） 即PTC的原理公式。 根据上述公式可得到PTC曲线，如标准差-平均值对数曲线、方差-平均值曲线等。 国内外基于PTC技术对CCD的成像过程进行了大量研究。李洪博、董岩、刘云清等提出了一种可提高相机动态范围的图像融合方法，使用PTC技术分别测得了高、低增益下CCD成像系统的PTC，并根据曲线的线性部分得到了图像融合阈值，实现了相机动态范围的提升 ［ 13 13 ］ 。冯婕、李豫东、文林等使用PTC技术对辐照前后探测器的性能进行了测试，发现辐照导致转换增益发生了7.82%的退化 ［ 14 14 ］ 。Qiang Wen， Siqi Zhu，Shichang Liu等使用神经网络算法对PTC的非线性区进行校正，提升了图像传感器的线性范围和图像对比度 ［ 15 15 ］ 。Deyan Levski， Martin Wany， Bhaskar Choubey提出了一种在使用PTC技术进行CMOS传感器性能测试时补偿信号相关读出噪声的方法，优化了利用PTC求取读出噪声的算法 ［ 16 16 ］ 。PTC方法也存在一些缺陷，通常需要搭建平场面光源（如使用积分球等设备），测试系统结构复杂，事实上由于CCD对光的响应极为敏感，使用一般的积分球最低亮度照射CCD也会发生饱和，为此定制了低照度积分球进行测试。PTC方法的前提是光源为绝对的平场，这是因为该方法事实上是以空间维噪声作为时间维噪声，因此需要所有像元接收相同亮度的光，对测试环境的杂散光要求较高。使用PTC方法有时并不能得到理想的曲线，尤其对于科学级CCD，像元相关性等因素使PTC在线性区会发生非正常弯曲现象 ［ 17 17 - 18 18 17-18 ］ 。另外，由于PTC方法采用统计方法计算CCD满阱，无法得到单个像元的满阱数值。 2.2　LED点测试法 2 根据CCD原理，CCD像元的主要结构是MOS管结构，光线照射到像元的耗尽区诱发电子能级跃迁产生新的电子-空穴对，可理解为光生电荷。光生电荷在像元势阱的作用下聚集于像元电极下，当光线增强，光生电荷量增多接近CCD满阱电子数，像元势阱不足以束缚所有电荷时，该像元位置处的电荷将部分被周围像元所捕获，即向周围溢出。电荷首先会向相邻像元溢出，导致其电荷数量增多，继续增加光强，像元势阱被电荷填满，将不具备继续收集电荷的能力，电荷量不再上升，附近像元也会随着光强的继续增强达到饱和。根据这一原理，可使用点光源照射CCD成像系统，使点光源在CCD上聚焦于一个像元之内，通过研究该像元的成像情况即可测得该像元的满阱电子数。 构建绝对的点光源系统，需要平行光管、星点板、数值孔径匹配的镜头等设备，成本较高。为了简化测试过程，本文使用LED点测试法，即使用LED作为光强可调的点光源，结合普通的成像镜头使LED点光源在CCD图像中覆盖3×3个像元区域，面积为78 μm×78 μm。通常CCD会首先在横向或纵向相邻的像元溢出，则LED点光源图像中的横向或纵向3个像元的DN值如 图1 图1 所示。 图1 图1 （a）为LED点在水平或垂直3个像元中的成像示意图， 图1 图1 （b）为随着光强的增大3个像元DN值的变化规律。 图1 图1 （b）中，区域（1）是线性增长区，此时3个像元均未达到满阱，随着光强的增大，3个像元的DN值同步增强，具体规律应符合高斯曲线的增长规律；当 b 像元接近满阱时，其搜集电荷的能力减弱，光生电荷被临近像元捕获的可能性增大，此时 b 的增长曲线逐渐放缓，临近像元的增长曲线开始加快，通常读出方向一侧的像元 a 稍早进入区域（2），即满阱区；继续增大光强，此时 b 像元达到饱和，数值不再上升，随后 a ， c 像元也达到饱和，即区域（3）。在满阱区，对 图1 图1 （b）中率先受到溢出影响的像元DN值变化趋势进行多项式拟合，根据多项式拟合曲线结果可求取 图1 图1 （b）中该像元曲线的拐点，根据满阱区的拐点对应的DN值，结合系统增益，即可得到CCD在该点像元的满阱电子数。 需要注意的是，CCD可能存在坏像元（一般高等级的CCD，坏像元数量极少），测试时应避开（通过不同亮度的平场图像可直观看到坏像元）。上述方法获取的是CCD单个像元的满阱数值，由于CCD的PRNU参数表征了像元之间对同一光强的响应差异性，通常不大于3%，文中所述载荷CCD的PRNU仅为0.643%，因此使用单个像元的满阱数值作为CCD的满阱参数时，其偏差不大于0.643%。为获取不同像元的满阱数值，调整LED光源在物平面的位置进行重复测试即可。为了获取更深的势阱和更大的像元有效感光面积，科学级CCD通常不做抗溢出结构，文中所述探测器即不具有抗溢出结构。对于采用了抗溢出设计的CCD，其LED点测试法的数据会有所不同，不会出现相邻像元DN值的迅速增大现象，因此使用这一方法可以看出CCD是否具有抗溢出设计。对于采用了抗溢出设计的CCD，LED点测试法的测试数据同样可用于成像电路研发过程中的满阱计算，可参考传统满阱测试法，找到中心像元线性偏差为3%的像素值作为满阱DN值，再结合增益即可获取满阱电子数。 3 实验及结果分析 1 本文分别采用PTC法和LED点测试法对载荷探测器满阱性能进行了实验，并采用传统饱和灰度值方法计算满阱容量作为对比。为了获取完整的实验数据，实验采用光源可调的测试方法，分别搭建了 图2 图2 所示的两种测试系统。 图2 图2 （a）为PTC测试系统， 图2 图2 （b）为LED点测试系统。为了降低光源不稳定性带来的误差，测试中使用了高稳定性电源且数据采集时将多幅平均结果作为测试数据。 3.1　数据采集 2 在PTC测试实验中，为了快速获取更多的测试数据，搭建了 图2 图2 （a）所示的测试系统，并获取差分图像，具体步骤如下： （1）调整积分球光源强度，获取两幅平场图像，分别记为Flat_A，Flat_B； （2）将两幅平场图像相减得到差分平场图像，记为Diff_AB； （3）选取图像中相对平坦的中间区域的200×200个像元，求取Flat_A，Flat_B中在该区域所有像元的平均值，记为 μ 1 ； （4）求取差分平场图像，Diff_AB在步骤（3）所述区域中所有像元的方差，记为 <math id="M13"><msubsup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966071&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966068&type= 3.21733332 3.97933316 ； （5）更改积分球光源强度，重复步骤（1）~（4）获取不同光强下的 μ i 和 <math id="M14"><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966077&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966074&type= 3.21733332 3.97933316 ，（ i= 1~ N，N 表示所有光强等级数量）； （6）将载荷调整工作模式为0 s，获取2幅图像，求取（3）所述区域中所有像元的平均值作为系统偏置，记为 B ； （7）绘制所有数据点，以 <math id="M15"><msub><mrow><mi>μ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><mi>B</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966083&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966080&type= 10.07533360 4.57200003 作为横轴， <math id="M16"><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966089&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966086&type= 6.51933336 7.87400007 作为纵轴，绘制PTC曲线。 根据概率分布理论，当随机变量 X 与随机变量 Y 相互独立时，其差的方差 D （ X - Y ）= D （ X ）+ D （ Y ）。上述步骤中，获取的两幅平场图像实际上是服从同一分布的随机变量的两次采样，两次采样相互独立。使用样本方差作为随机变量方差的估计时，也符合方差公式。因此，差分图像的方差 <math id="M17"><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966095&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966092&type= 3.21733332 3.97933316 应为平场图像方差的2倍，即步骤（7）中计算得到的方差需要乘以1/2。 搭建 图2 图2 （b）所示的LED点测试法测试系统，以测量CCD某个像元的满阱电子数，具体步骤如下： （1）调整LED的电压为 U 1 ，点亮LED，调整镜头焦距及LED与探测器的距离，使LED在CCD图像上聚焦为一个亮点，占据3×3的像元面积； （2）观察LED照亮像元对应的DN值分布，调整其位置，使DN值分布符合 图3 图3 （a）中的规律，即左右及上下像元的DN值相近，获取16幅图像，计算9个像素点在16幅图像（或更多）中的DN值的平均值； （3）依次改变LED电压为 U i ，获取不同的光强等级，选取光强等级的原则为在饱和点附件选取更多的点，以减小数据拟合误差，重复步骤（2），获取9个像素点在不同光强下的数据，分别记为 <math id="M18"><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">1</mn><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966111&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966108&type= 4.23333359 3.72533321 ， <math id="M19"><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966117&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966114&type= 4.23333359 3.72533321 ，……， <math id="M20"><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">9</mn><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966123&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966120&type= 4.23333359 3.72533321 ，（ i =1~ M ， M 为所有光强等级数量）； （4）以电压 U i 为横轴、像元DN值 <math id="M21"><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">1</mn><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966127&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966125&type= 4.23333359 3.72533321 ， <math id="M22"><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966134&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966130&type= 4.23333359 3.72533321 ，……， <math id="M23"><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">9</mn><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966140&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32966136&type= 4.23333359 3.72533321 为纵轴绘制3条曲线，选取合适范围进行多项式拟合得到关键曲线拐点及其对应的中心像元DN值，该值即为中心像元达到满阱时对应的DN值。 3.2　结果及分析 2 3.2.1　PTC测试结果 3 按照上述方法，得到载荷探测器在全幅工作模式下的标准差-平均值对数PTC如 图3 图3 所示，方差-平均值PTC如 图4 图4 所示。 1,2; 对 图3 图3 （d）中的拟合曲线求导，得到其斜率为0的点为（1.97，0.676 34），对应的标准差为4.762DN，即电路的读出噪声（与光强无关的噪声）为4.762DN。对 图3 图3 （a）中数据的饱和值点（即最高点）之前的数据进行线性拟合，结果如 图3 图3 （b）所示，斜率为0.386，这与理论值0.5不符，但是可以看出，随着Log（Mean-B）值的增大曲线斜率呈增长趋势，这说明散粒噪声在总噪声中的比重不足，因此不适合用该曲线来计算成像系统增益。对数PTC曲线在起始阶段呈现明显的水平趋势，所以可用来测算成像系统的读出噪声。CCD的读出数据中包含Pre-scan和Over-scan的读出像元，选取Pre-scan和Over-scan区域的像元进行平均值和方差的计算，得到此区域像元的平均值为522.4DN，方差为21.22DN 2 。此结果与PTC曲线得到的读出噪声偏差不大，实际中可使用Pre-scan和Over-scan快速获取载荷的读出噪声。 对 图4 图4 （a）数据进行线性拟合，斜率为0.011 05，即成像系统增益 g =0.011 05；数据最高点对应的Mean-B数值为9 036DN，其对应的电子数即为CCD像元满阱电子数为817.738 ke - 。由于方差-平均值曲线中的截距拟合结果误差较大，此处取标准差-平均值对数曲线中的拟合结果，即4.762DN作为读出噪声，具体载荷探测器的满阱参数测试结果如 表1 表1 所示。 Parameters Value Readout Noise（DN） 4.762 System Gain（DN/e - ） 0.011 05 FWC（ke - ） 817.7 3.2.2　LED点测试结果 3 根据2.2所述方法获取LED点测试数据，得到LED点在CCD中的成像结果，如 图5 图5 所示。由 图5 图5 （a）可见，在LED光强增加的过程中，中间像元DN值迅速增大，上下两侧像元率先于水平方向两侧像元开始增大。 图5 图5 （b）可直观看出LED点光源在CCD垂直方向（帧转移型面阵CCD像元的读出过程一般包括帧转移和水平转移过程，帧转移有时也称垂直转移，沿垂直方向溢出，即是沿帧转移方向溢出）率先溢出。 根据CCD的像元搜集电荷的原理，CCD对光生电荷的约束是依靠像元MOS管结构的电势垒和电势阱，随着CCD像元搜集的电荷数的增加，接近像元饱和时，该像元结构搜集电荷的能力发生变化 ［ 19 19 ］ ，此时临近像元会收集到更多的电子。根据帧转移型CCD的工作原理，CCD帧转移方向为垂直方向，像元读出方向为水平方向 ［ 20 20 ］ ，因此饱和后的像元在水平或垂直方向上率先发生溢出 ［ 21 21 ］ ，而且是临近的某一个像元受到溢出影响最大，即DN值会迅速上升。 图6 图6 所示为LED点亮的9个像元在不同光强下DN值的增加趋势。 图6 图6 （a）为像元DN值走势图，可看出LED点光源投影到CCD像元上后的电荷向周围溢出过程：在中心像元未饱和之前，其DN值的增长与LED电压的增加呈线性比例（由于低电压处LED亮度稳定性差，因此在低DN值区域成像部分波动，但对满阱容量的计算并无影响，因为像元不可能在低值处饱和）；当中心像元接近饱和时，临近的下方像元首先收集到中心像元区域的电荷，其DN值迅速增长；光源继续增强，上方临近像元开始收集到中心像元区域电荷；当中心像元饱和溢出时，周围像元开始受到影响。 按照像元饱和过程，可选取首先受到中心像元饱和影响的像元（down_center pixel）DN值变化做分析。如 图6 图6 （b）所示，计算得到down_center像元在转折区域附近拟合曲线的拐点为（2.09，1 517），该点对应的中心像元DN拟合曲线的拟合值为9 562，即该像元发生饱和时对应的DN值。PTC方法得到的去除偏置的饱和DN值为9 036，LED点测试法扣除偏置值522.4后对应的DN值为9 039.6，其对应的满阱电子数即为818 063e - 。两者的误差仅为0.039 7%，证明LED点测试法的准确性。进一步地，通过调整LED点光源的位置，可获取CCD成像面上不同像元的满阱参数。4个不同像元采用LED点测试法的测试结果如 图7 图7 所示。 1;2; 综上，采用LED点测试法获取的5个像元测试结果如 表2 表2 所示。可见，每个像元饱和DN值相对平均值的偏差也符合出厂参数中的PRNU参数0.643%。 Pixel no. Saturated values/DN Mean value/DN Deviation/% FWC/ke - 1 9 562 5 9 550.4（减去偏置后对应满阱817.013 ke - ） 0.121 5 818.063 2 9 587 0.383 2 820.326 3 9 607 0.592 6 822.136 4 9 501 0.517 812.543 5 9 495 0.58 812 3.2.3　传统方法的测试结果 3 根据EMVA1288标准的CCD成像系统满阱测试法，即传统测试方法，首先需要求取成像系统的饱和DN值，然后通过线性拟合得到饱和值。随着光强的增大，成像系统获取的DN值在像元饱和前呈线性增长，接近饱和及发生溢出后DN值随光强的变化会变缓并最终停止增长。对线性区进行拟合，计算像元DN值与该线性拟合结果的误差比例，使得该误差比例达到5%或3%时的像元DN值为饱和DN值。为了对比LED点测试法与传统满阱测试方法，对中心像元的DN值采用传统方法进行处理，结果如 图8 图8 所示（彩图见期刊电子版）。 图8 图8 中，三角形符号为中心像元随光强的变化趋势，在其上升区进行线性拟合，拟合直线如红线所示，方形符号曲线为线性拟合误差比例。初始阶段由于电压较低，LED光源波动较大，因此拟合误差也会偏大，但并不影响计算结果。当误差比例在转折区域上升为3%和5%时对应的中心像元DN值分别为9 728和9 985，其对应的满阱电子数即为880 362e - 和903 620e - 。该结果与PTC方法的测试结果偏差分别为1.9%和4.7%，电子数偏差则为62.624ke - 和85.882ke - 。可见，这一方法误差较大，且选取不同拟合误差比例作为满阱定义存在较大区别。同型号芯片在不同实验组测试中可能选取不同定义，得到无法统一的数据结果。 4 结 论 1 LED点测试法基于像元饱和过程的物理现象和原理，是一种获得CCD像元满阱电子数的简易方法。使用LED点测试法对载荷探测器的多个像元进行多次测试，得到满阱电子数平均为817.013ke - ，与PTC法的测试结果817.738ke - ，偏差仅为0.088 7%；针对单个像元的LED点测试数据，分别采用LED点测试法和传统的满阱测试方法进行测试，其结果与PTC测试结果的误差分别为0.039 7%和1.9%，从而验证了LED点测试法的准确性。使用单个像元的满阱电子数作为CCD满阱参数时，其偏差不大于0.643%，即文中所述载荷探测器的PRNU出厂参数。综上，LED点测试法作为简易条件下CCD满阱测试方法，具有较高的准确性和可靠性，可大幅提升CCD成像电路的研发效率。 参考文献 黄珊 ， 司福祺 ， 赵敏杰 ， 等 . 星载大气痕量气体差分吸收光谱仪狭缝函数研究 ［J］. 光谱学与光谱分析 ， 2019 ， 39 （ 7 ）： 2008 - 2012 . 星载大气痕量气体差分吸收光谱仪狭缝函数研究 HUANG SH ， SI F Q ， ZHAO M J ， et al . Study on the slit function of atmospheric trace gas differential optical absorption spectrometer ［J］. Spectroscopy and Spectral Analysis ， 2019 ， 39 （ 7 ）： 2008 - 2012 . （in Chinese） Study on the slit function of atmospheric trace gas differential optical absorption spectrometer ZHAO M J ， SI F Q ， WANG Y ， et al . First year on-orbit calibration of the Chinese environmental trace gas monitoring instrument onboard GaoFen-5 ［J］. IEEE Transactions on Geoscience and Remote Sensing ， 2020 ， 58 （ 12 ）： 8531 - 8540 . doi: 10.1109/tgrs.2020.2988573 http://dx.doi.org/10.1109/tgrs.2020.2988573 First year on-orbit calibration of the Chinese environmental trace gas monitoring instrument onboard GaoFen-5 鲁月林 ， 王煜 ， 司福祺 ， 等 . 星载差分吸收光谱仪摆镜控制系统设计 ［J］. 光学 精密工程 ， 2019 ， 27 （ 3 ）： 630 - 636 . doi: 10.3788/ope.20192703.0630 http://dx.doi.org/10.3788/ope.20192703.0630 星载差分吸收光谱仪摆镜控制系统设计 LU Y L ， WANG Y ， SI F Q ， et al . Design of scanning mirror control system for satellite-borne DOAS spectrometer ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 3 ）： 630 - 636 . （in Chinese） . doi: 10.3788/ope.20192703.0630 http://dx.doi.org/10.3788/ope.20192703.0630 Design of scanning mirror control system for satellite-borne DOAS spectrometer 常振 ， 王煜 ， 司福祺 ， 等 . 基于科学级CCD的紫外成像系统设计与实现 ［J］. 中国激光 ， 2017 ， 44 （ 8 ）： 0804002 . doi: 10.3788/cjl201744.0804002 http://dx.doi.org/10.3788/cjl201744.0804002 基于科学级CCD的紫外成像系统设计与实现 CHANG ZH ， WANG Y ， SI F Q ， et al . Design and implementation of ultraviolet imaging system based on scientific grade CCD ［J］. Chinese Journal of Lasers ， 2017 ， 44 （ 8 ）： 0804002 . （in Chinese） . doi: 10.3788/cjl201744.0804002 http://dx.doi.org/10.3788/cjl201744.0804002 Design and implementation of ultraviolet imaging system based on scientific grade CCD 邱晓晗 ， 王煜 ， 常振 ， 等 . 机载紫外DOAS成像光谱仪CCD成像电路的设计及实施 ［J］. 红外与激光工程 ， 2017 ， 46 （ 5 ）： 0538002 . doi: 10.3788/irla201746.0538002 http://dx.doi.org/10.3788/irla201746.0538002 机载紫外DOAS成像光谱仪CCD成像电路的设计及实施 QIU X H ， WANG Y ， CHANG ZH ， et al . Design and implementation of CCD imaging circuit for airborne ultraviolet DOAS imaging spectrometer ［J］. Infrared and Laser Engineering ， 2017 ， 46 （ 5 ）： 0538002 . （in Chinese） . doi: 10.3788/irla201746.0538002 http://dx.doi.org/10.3788/irla201746.0538002 Design and implementation of CCD imaging circuit for airborne ultraviolet DOAS imaging spectrometer 刘丹丹 ， 黄印博 ， 孙宇松 ， 等 . 地表气压对温室气体柱总量反演的敏感性分析及敦煌地区温室气体的柱总量观测 ［J］. 光学 精密工程 ， 2020 ， 28 （ 3 ）： 573 - 583 . doi: 10.3788/ope.20202803.0573 http://dx.doi.org/10.3788/ope.20202803.0573 地表气压对温室气体柱总量反演的敏感性分析及敦煌地区温室气体的柱总量观测 LIU D D ， HUANG Y B ， SUN Y S ， et al . Analysis of earth surface pressure sensitivity on inversion of greenhouse gases columns and observation of greenhouse gases columns in Dunhuang in China ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 3 ）： 573 - 583 . （in Chinese） . doi: 10.3788/ope.20202803.0573 http://dx.doi.org/10.3788/ope.20202803.0573 Analysis of earth surface pressure sensitivity on inversion of greenhouse gases columns and observation of greenhouse gases columns in Dunhuang in China 闫欢欢 ， 李晓静 ， 张兴赢 ， 等 . 大气SO 2 柱总量遥感反演算法比较分析及验证 ［J］. 物理学报 ， 2016 ， 65 （ 8 ）： 084204 . doi: 10.7498/aps.65.084204 http://dx.doi.org/10.7498/aps.65.084204 大气SO2柱总量遥感反演算法比较分析及验证 YAN H H ， LI X J ， ZHANG X Y ， et al . Comparison and validation of band residual difference algorithm and principal component analysis algorithm for retrievals of atmospheric SO 2 columns from satellite observations ［J］. Acta Physica Sinica ， 2016 ， 65 （ 8 ）： 084204 . （in Chinese） . doi: 10.7498/aps.65.084204 http://dx.doi.org/10.7498/aps.65.084204 Comparison and validation of band residual difference algorithm and principal component analysis algorithm for retrievals of atmospheric SO2 columns from satellite observations 牟福生 ， 雒静 ， 李素文 . MAX-DOAS仪器观测O 4 吸收对气溶胶参量的敏感性研究 ［J］. 光子学报 ， 2020 ， 49 （ 2 ）： 0201002 . doi: 10.3788/gzxb20204902.0201002 http://dx.doi.org/10.3788/gzxb20204902.0201002 MAX-DOAS仪器观测O4吸收对气溶胶参量的敏感性研究 MOU F SH ， LUO J ， LI S W . Sensitivity study of aerosol parameters from observations of O 4 absorption by MAX-DOAS ［J］. Acta Photonica Sinica ， 2020 ， 49 （ 2 ）： 0201002 . （in Chinese） . doi: 10.3788/gzxb20204902.0201002 http://dx.doi.org/10.3788/gzxb20204902.0201002 Sensitivity study of aerosol parameters from observations of O4 absorption by MAX-DOAS 王伶俐 ， 李守贝 ， 张惠鸽 ， 等 . 科学级光学CCD信噪比特性分析 ［J］. 核聚变与等离子体物理 ， 2014 ， 34 （ 4 ）： 361 - 366 . doi: 10.3969/j.issn.0254-6086.2014.04.013 http://dx.doi.org/10.3969/j.issn.0254-6086.2014.04.013 科学级光学CCD信噪比特性分析 WANG L L ， LI S B ， ZHANG H G ， et al . Performance analysis on SNR of scientific grade optical CCD ［J］. Nuclear Fusion and Plasma Physics ， 2014 ， 34 （ 4 ）： 361 - 366 . （in Chinese） . doi: 10.3969/j.issn.0254-6086.2014.04.013 http://dx.doi.org/10.3969/j.issn.0254-6086.2014.04.013 Performance analysis on SNR of scientific grade optical CCD 夏果 ， 黄禅 ， 吴骕 ， 等 . 微型光谱仪增益及读出噪声的测量方法 ［J］. 光学学报 ， 2017 ， 37 （ 1 ）： 0112001 . doi: 10.3788/aos201737.0112001 http://dx.doi.org/10.3788/aos201737.0112001 微型光谱仪增益及读出噪声的测量方法 XIA G ， HUANG SH ， WU S ， et al . Method of measuring gain and readout noise of miniature spectrometer ［J］. Acta Optica Sinica ， 2017 ， 37 （ 1 ）： 0112001 . （in Chinese） . doi: 10.3788/aos201737.0112001 http://dx.doi.org/10.3788/aos201737.0112001 Method of measuring gain and readout noise of miniature spectrometer XIONG X ， KIM C ， SHIMODA H . Earth observing missions and sensors： development， implementation， and characterization ［J］. SPIE. Asia-pacific Remote Sensing IV ， 2012 . Earth observing missions and sensors： development， implementation， and characterization 周鹏骥 ， 王晓东 ， 董吉洪 ， 等 . 天问一号高分相机成像噪声分析与抑制 ［J］. 光学 精密工程 ， 2022 ， 30 （ 2 ）： 217 - 226 . 天问一号高分相机成像噪声分析与抑制 ZHOU P J ， WANG X D ， DONG J H ， et al . Imaging noise analyzing and suppressing for Tianwen-1 high-resolution camera ［J］. Opt. Precision Eng. ， 2022 ， 30 （ 2 ）： 217 - 226 . （in Chinese） Imaging noise analyzing and suppressing for Tianwen-1 high-resolution camera 李洪博 ， 董岩 ， 刘云清 ， 等 . 可提高相机动态范围的图像融合方法研究 ［J］. 红外技术 ， 2018 ， 40 （ 9 ）： 887 - 893 . 可提高相机动态范围的图像融合方法研究 LI H B ， DONG Y ， LIU Y Q ， et al . Image fusion method for improving dynamic range of camera ［J］. Infrared Technology ， 2018 ， 40 （ 9 ）： 887 - 893 . （in Chinese） Image fusion method for improving dynamic range of camera 冯婕 ， 李豫东 ， 文林 ， 等 . CMOS图像传感器光子转移曲线辐照后的退化机理 ［J］. 光学 精密工程 ， 2017 ， 25 （ 10 ）： 2676 - 2681 . doi: 10.3788/ope.20172510.2676 http://dx.doi.org/10.3788/ope.20172510.2676 CMOS图像传感器光子转移曲线辐照后的退化机理 FENG J ， LI Y D ， WEN L ， et al . Degradation mechanism for photon transfer curve of CMOS image sensor after irradiation ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 10 ）： 2676 - 2681 . （in Chinese） . doi: 10.3788/ope.20172510.2676 http://dx.doi.org/10.3788/ope.20172510.2676 Degradation mechanism for photon transfer curve of CMOS image sensor after irradiation WEN Q ， ZHU S Q ， LIU S C ， et al . Image sensor correction algorithm for photon transfer curve based on neural network ［J］. Optical Engineering ， 2020 ， 59 ： 067102 . doi: 10.1117/1.oe.59.6.067102 http://dx.doi.org/10.1117/1.oe.59.6.067102 Image sensor correction algorithm for photon transfer curve based on neural network LEVSKI D ， WÄNY M ， CHOUBEY B . Compensation of signal-dependent readout noise in photon transfer curve characterisation of CMOS image sensors ［J］. IEEE Transactions on Circuits and Systems II： Express Briefs ， 2021 ， 68 （ 1 ）： 102 - 105 . doi: 10.1109/tcsii.2020.3010366 http://dx.doi.org/10.1109/tcsii.2020.3010366 Compensation of signal-dependent readout noise in photon transfer curve characterisation of CMOS image sensors ASTIER P ， ANTILOGUS P ， JURAMY C ， et al . The shape of the photon transfer curve of CCD sensors ［J］. Astronomy & Astrophysics ， 2019 ， 629 ： A36 . doi: 10.1051/0004-6361/201935508 http://dx.doi.org/10.1051/0004-6361/201935508 The shape of the photon transfer curve of CCD sensors ANTILOGUS P ， ASTIER P ， DOHERTY P ， et al . The brighter-fatter effect and pixel correlations in CCD sensors ［J］. Journal of Instrumentation ， 2014 ， 9 （ 3 ）： C03048 . doi: 10.1088/1748-0221/9/03/c03048 http://dx.doi.org/10.1088/1748-0221/9/03/c03048 The brighter-fatter effect and pixel correlations in CCD sensors MA B ， SHANG Z H ， HU Y ， et al . A new method of CCD dark current correction via extracting the dark Information from scientific images ［C］. SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 9154， High Energy， Optical， and Infrared Detectors for Astronomy VI ， Montréal， Quebec， Canada . 2014 ， 9154 ： 593 - 600 . doi: 10.1117/12.2055416 http://dx.doi.org/10.1117/12.2055416 A new method of CCD dark current correction via extracting the dark Information from scientific images JORDEN P R ， JERRAM P A ， FRYER M ， et al . e2v CMOS and CCD sensors and systems for astronomy ［J］. Journal of Instrumentation ， 2017 ， 12 （ 7 ）： C07008 . doi: 10.1088/1748-0221/12/07/c07008 http://dx.doi.org/10.1088/1748-0221/12/07/c07008 e2v CMOS and CCD sensors and systems for astronomy 常振 ， 王煜 ， 司福祺 ， 等 . CCD数字像元的合并方法及其应用 ［J］. 光学 精密工程 ， 2017 ， 25 （ 8 ）： 2204 - 2211 . doi: 10.3788/ope.20172508.2204 http://dx.doi.org/10.3788/ope.20172508.2204 CCD数字像元的合并方法及其应用 CHANG ZH ， WANG Y ， SI F Q ， et al . Digital binning method for CCD pixels and its applications ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 8 ）： 2204 - 2211 . （in Chinese） . doi: 10.3788/ope.20172508.2204 http://dx.doi.org/10.3788/ope.20172508.2204 Digital binning method for CCD pixels and its applications 10.37188/OPE.20223013.1542 TP73 TN386.5 A 1004-924X（2022）13-1542-13 常振,王煜,林方等.载荷探测器像元满阱参数测试[J].光学精密工程,2022,30(13):1542-1554. 常振,王煜,林方等.载荷探测器像元满阱参数测试[J].光学精密工程,2022,30(13):1542-1554. CHANG Zhen,WANG Yu,LIN Fang,et al.Measurement of pixel full-well parameters of detector used on spaceborne spectrometer[J].Optics and Precision Engineering,2022,30(13):1542-1554. CHANG Zhen,WANG Yu,LIN Fang,et al.Measurement of pixel full-well parameters of detector used on spaceborne spectrometer[J].Optics and Precision Engineering,2022,30(13):1542-1554. 国家自然科学基金青年科学基金资助项目（No.61905256 ） 2022-03-24 2022-04-12 2022-07-10 2022-08-07 光学精密工程 13 30