CN113281004A - 天文光学望远镜光电效率计算及实测验证方法 - Google Patents

天文光学望远镜光电效率计算及实测验证方法 Download PDF

Info

Publication number
CN113281004A
CN113281004A CN202110482419.5A CN202110482419A CN113281004A CN 113281004 A CN113281004 A CN 113281004A CN 202110482419 A CN202110482419 A CN 202110482419A CN 113281004 A CN113281004 A CN 113281004A
Authority
CN
China
Prior art keywords
telescope
star
transmittance
formula
wavelength
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110482419.5A
Other languages
English (en)
Other versions
CN113281004B (zh
Inventor
刘伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Purple Mountain Observatory of CAS
Original Assignee
Purple Mountain Observatory of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Purple Mountain Observatory of CAS filed Critical Purple Mountain Observatory of CAS
Priority to CN202110482419.5A priority Critical patent/CN113281004B/zh
Publication of CN113281004A publication Critical patent/CN113281004A/zh
Application granted granted Critical
Publication of CN113281004B publication Critical patent/CN113281004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • G01M11/02Testing optical properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Analytical Chemistry (AREA)
  • Computational Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本是发明公开了天文光学望远镜光电效率计算方法及实测验证方法,其中天文光学望远镜光电效率计算方法包括以下步骤:步骤一、统计天文光学望远镜各透光元器件在不同波段的透过率数据,根据数据,绘制成各透光元器件透过率‑波长曲线;步骤二、采用matlab进行编程,先加载这些透过率‑波长曲线的数据,再对各透过率‑波长曲线进行插值提取出各透光元器件在该插值处的透过率的值,将这些值进行相乘,得到插值处的总体效率,最后将总体效率拟合成透光元器件综合透过率曲线。本发明在完成望远镜的设计与镀膜之后对望远镜的整体的透过率进行计算,得到透光元器件综合透过率曲线,进而获得天文光学望远镜各波长处的整体光电效率,方法简单实用。

Description

天文光学望远镜光电效率计算及实测验证方法
技术领域
本发明涉及天文设备的技术领域,尤其涉及天文光学望远镜光电效率计算及实测验证方法。
背景技术
望远镜的光学设计可以分为折射式,反射式和折返式三种。但是无论采用哪种形式的光学结果,望远镜或多或少都存在像差和色差的问题。为了消去像差(球面像差,彗形像差)色差的影响,一些多透镜物镜组等被设计到望远镜光路中以改善望远镜的成像质量。这些附加的镜片都有一个光线透过率的问题,镜片对于入射光线有透射,反射以及被镜片吸收三种结果。反射和吸收都会造成光能量的损失,所以为增加光线的透过率,目前广泛采用的方法式对镜片进行镀膜,以降低反射率吸收率,增大透过率。设计多个改正镜,对望远镜各个光学元件镀膜是目前改善望远镜成像质量,提高透过率的普遍方法。
望远镜是由各种光学元件所组成,尤其是多波段同时成像的望远镜,其光学元器件的构成更多。这里以西藏阿里50cm HinOTORI三波段同时成像望远镜为例,介绍其光学结构并对其透过率进行分析。HinOTORI其结构采用卡塞格林反射式的结构。除了主镜和副镜为三个波段公用以外,望远镜每一个波段的所包含的光学元件不尽相同。三个波段的波段分离是通过两个分色镜实现。其中入射光被分色镜1反射,又通过一个平面镜反射到u'波段CCD相机成像;入射光透射过分色镜1,又被分色镜2反射,再到Rc波段CCD相机成像;入射光透射过分色镜1和分色镜2,形成Ic波段光束,并再Ic波段CCD相机成像。
在这些波段上,所包含的元器件各不相同,并且针对在每一个波段的光线透过波长,对各个元件进行了相应的光学镀膜处理,以提高其整体的透过率。
所以在完成望远镜的设计与镀膜之后对望远镜的整体的透过率进行计算是十分必要的。
发明内容
本发明的目的是为了解决背景技术中提及的问题,提供一种天文光学望远镜光电效率计算及实测验证方法。
为实现上述技术目的,本发明采取的技术方案为:
天文光学望远镜光电效率计算方法,包括以下步骤:
步骤一、统计天文光学望远镜各透光元器件在不同波段的透过率数据,根据数据,绘制成各透光元器件透过率-波长曲线;
步骤二、采用matlab进行编程,先加载这些透过率-波长曲线的数据,再对各透过率-波长曲线进行插值提取出各透光元器件在该插值处的透过率的值,将这些值进行相乘,得到插值处的总体效率,最后将总体效率拟合成透光元器件综合透过率曲线,进而获得天文光学望远镜各波长处的整体光电效率。
为优化上述技术方案,采取的具体措施还包括:
步骤一中,天文光学望远镜的透光元器件包括望远镜主镜,望远镜副镜,各分色镜、Ic波段镀膜以及CCD探测器,相应的,各透光元器件透过率-波长曲线分别包括望远镜主镜和副镜的透过率-波长曲线、分色镜的透过率-波长曲线、Ic波段镀膜的透过率-波长曲线以及CCD量子效率-波长曲线。
天文光学望远镜光电效率实测验证方法,包括以下步骤:
步骤a、定义采用权利要求1所述的透光元器件综合透过率曲线所得到的为望远镜整体效率ηA,用于观测实测验证的为望远镜效率ηB;选定目标星体作为观测对象;
步骤b、在大气外条件下,计算望远镜接收到目标星体的理想光子数;
步骤c、用望远镜在天文台上进行对目标星体进行实测,选取不同的天顶距进行成像,对图像进行处理,计算出各个波段上的大气消光因子K;
步骤d、利用大气消光因子K,将台址得到的目标星体的仪器星等归算到大气外,获得大气外的目标星体的仪器星等;
步骤e、计算出望远镜的CCD图像的Count值,再根据Count值计算CCD探测器所产生出的电子数,
步骤f、通过步骤b计算出的理想光子数以及步骤e得到的电子数,计算得到望远镜整体光电效率的值ηB,
步骤g、将望远镜的整体光电效率的值ηB与望远镜光电效率ηA进行比较,确认二者数值是否接近。
步骤b中,在大气外条件下,计算望远镜接收到目标星体的理想光子数的具体方法为:
Figure BDA0003048952300000021
其中NP为理想光子数,A为望远镜的集光面积,h为普朗克常数(6.626196×10- 27erg·s)mAB为目标天体的AB星等,星等采用大气外的查表星等,积分项是一个与波长相关的表达式,R(λ)表示在波长积分范围内的整体透过率;
为了计算积分项的数值,在matlab中进行编程计算,在Johnson-Cousins波长系统中的Ic波段,设定Ic波段积分的范围为700nm至900nm,采用sum函数得出R(λ)在700nm~900nm范围内积分项的值为0.089,代入上式(1),从而得到望远镜在单位时间内接收到的一个固定星等星的光子数公式为:
NP_Ic=8.426*108*10-0.4mAB 式:(2)
公式(2)对于Sloan波长系统同样适用,当望远镜的滤光片为Sloan中的u’波段时,设置u’积分波长范围为300nm~400nm;当望远镜的滤光片为Johnson-Cousins中的Rc波段时,设置Rc积分波长范围为550nm~800nm;当望远镜的滤光片为Johnson-Cousins中的Ic波段时,设置Rc积分波长范围为550nm~800nm,得到这三个波段的理想的光子NP_u',NP_Rc和NP_Ic如下式:
Figure BDA0003048952300000031
步骤c中,计算出各个波段上的大气消光因子K的具体方法为:
根据大气消光的理论可以得到大气消光的公式:
m(z)=m0+K*F(z) 式:(4)
其中m(z)为天体在天顶距z处的大气内星等,m0为天体的大气外星等,K为消光因子,F(z)为大气质量,它与天顶距Z的角度有直接的关系;
大气质量F(z)的计算公式为:
F(Z)=secz-0.0018167(secz-1)-0.002875(secz-1)2-0.0008083(secz-1)3 式:(5)
sec表示对天顶距Z取正割运算;
天顶距z处的大气内星等m(z)可以用仪器星等minst+加常数ZP表示,ZP为星等零点;图像中的仪器星等的计算公式为:
minst=-2.5*log10Count 式:(6)
Count为星在图像中的Count值;
因此大气外星等m0的整体计算公式可以表示为:
m0=minst+ZP-K*F(z) 式:(7)
为了获得各个波段的大气消光系数K,在实测过程中,改变天顶距Z,即望远镜的俯仰角,对多个标准星目标进行CCD成像,通过查表可得m0,并计算其仪器星等minst,再通过上式进行拟合,计算出各波段的大气消光系数K。
步骤d中,将台址得到的目标星体的仪器星等归算到大气外的具体方法为:
M0inst-band=Minst-band-F(Z)*Kband 式:(9)
M0inst-band为归算到大气层外后的天体仪器星等,Minst-band为台址上某一个大气质量下测量的天体仪器星等,F(Z)为当时的大气质量,其与天顶距有关,Kband为某一个波段的大气消光系数,通过实测得到。
步骤e、计算出望远镜的CCD图像的Count值,再根据Count值计算CCD探测器所产生出的电子数的具体方法为:
根据式(6)可以推导台址上各个波段所观测的含有大气消光影响下的天体目标其Count数Countband计算公式为
Figure BDA0003048952300000041
所以如果排除大气消光的影响,在得到归算后的大气外仪器星等之后,可以推导出在大气外,某一天体的理想Count值Count0band的计算公式为:
Figure BDA0003048952300000042
得到理想的Count数之后即可得到理想的电子数,计算公式为:
Figure BDA0003048952300000043
e-为电子数,ADU为Count数值。此参数在CCD相机出厂时经过测试;
已知CCD的增益以及理想的大气外Count0值之后,可以采用式(13)计算理想的大气外CCD转换生成的电子数e-0:
e-0=Count0*gain 式:(13)。
步骤f、通过步骤b计算出的理想光子数以及步骤e得到的电子数,计算得到望远镜整体光电效率的值η的具体方法为:
Figure BDA0003048952300000044
Ne-为整体接收的电子数,其计算公式为(11)和(13);Np为整体接收的光子数,其计算公式为式(3)。
本发明具有以下优点:
1、本发明在完成望远镜的设计与镀膜之后对望远镜的整体的透过率进行计算,得到透光元器件综合透过率曲线,进而获得天文光学望远镜各波长处的整体光电效率,方法简单实用。
2、本发明采用实测方法验证望远镜的光电效率,充分考虑到观测时无法避免的望远镜所处台址的大气透过率因素,从而能更准确的检测出望远镜的光电效率。
附图说明
图1是背景技术中多波段同时成像望远镜光学结构的结构示意图;
图2是Ic波段整体透过率理论仿真图;
图3是Rc波段整体透过率理论仿真图;
图4是u’波段整体透过率理论仿真图;
图5是大气消光影响示意图;
图6是天文光学望远镜光电效率计算方法的流程图;
图7是天文光学望远镜光电效率实测验证方法的流程图。
具体实施方式
以下结合附图对本发明的实施例作进一步详细描述。
这里以HinOTORI望远镜的Ic波段为例,详细介绍望远镜透过率以及光电效率的计算方法。
一:统计某一个波段各个元器件在不同波段的透过率曲线。如图1所示,
Ic光路上元器件包括主镜,副镜,分色镜1的透过率,分色镜2的透过率,Ic波段镀膜透过率曲线。获得这些透过率数据,并绘制成透过率-波长曲线。
二:综合各因素,绘制整体光电效率曲线。
由于Ic波段的中心波长约为798nm,为了过滤掉其他波段的光子,在CCD探测的前端设计了透过Ic波长的滤光片,光子经过滤光片之后直接进入CCD探测器,在CCD探测器光电作用之下,入射的光子被转化为电子,再经过数字化成像之后送到计算机进行图像处理分析,CCD将光子转换为电子的能力表现在其QE(量子效率)曲线上。在分析整体透过率时,滤光片的透过率曲线和CCD QE曲线也必须包含在内。获得这些元件的透过率和波长的之后,对其进行列表编号:其中第1组为Ic波段滤光片透过率;第2组为CCD量子效率;第3组为分色镜1的透过率;第4组为分色镜2的透过率;第5组为Ic波段镀膜的透过率;第6组为主镜和副镜的透过率,这些透过率及效率都是随波长而改变,与波长有着对应的关系。采用matlab进行编程,先加载这些透过率的数据,再用插值提取出这些透过率的各项之值,再在各插值处将各项因子相乘,得到插值处的总体效率,最后拟合成透过率曲线,见图2中黑色实线部分。
由于Rc波段与u’波段在光学元件和光路上与Ic波段有所不同,综合Rc和u’波段各透过率影响因素得到的整体透过率曲线如图3-4所示。
整体透过率的仿真处理光学透过率的因素,还包括了CCD QE量子效率这一项,QE是电子数与光子数的比值,整体效率包含QE,即可以理解为望远镜将收集光子转为电子的整体光电效率。另外必须指出,此总体透过率之与望远镜的仪器本身有关,并没有考虑观测时无法避免的望远镜所处台址的大气透过率因素。
此曲线所示的透过率曲线也是随各处波长的变化而有所不同。以上是通过编程获得望远镜整体光电效率的仿真图形。下面介绍望远镜采用实测方法验证望远镜的光电效率。
一:首先介绍对透过率曲线进行量化,并计算光子数的方法。
一架望远镜的集光力(收集光子的能力)直接与即口径(D)的大小有关。口径越大收集的光线越多。理论上收集的光子数Np可有表达式(1)计算。
Figure BDA0003048952300000061
其中NP为理想光子数,A为望远镜的集光面积,h为普朗克常数(6.626196×10- 27erg·s)mAB为目标天体的AB星等,积分项是一个与波长相关的表达式,R(λ)表示在波长积分范围内的整体透过率。
为了计算积分项的数值,在matlab中进行编程计算,设定Ic波段积分的范围为700nm至900nm,R(λ)为整体透过率,采用sum函数得出在700nm~900nm范围内积分项的值为0.089,代入上式(1),从而得到望远镜在单位时间内接收到的一个固定星等(AB系统)星的光子数公式为:
NP_Ic=8.426*108*10-0.4mAB 式:(2)
为了计算积分项的数值,在matlab中进行编程计算,以Johnson-Cousins波长系统中的Ic波段为例,设定Ic波段积分的范围为700nm至900nm,采用sum函数得出R(λ)在700nm~900nm范围内积分项的值为0.089,代入上式(1),从而得到望远镜在单位时间内接收到的一个固定星等星的光子数公式为:
NP_Ic=8.426*108*10-0.4mAB 式:(2)
此公式对于Sloan波长系统(u',g',r',i',and z')同样适用。例如当望远镜的滤光片为Sloan中的u’波段时,设置u’积分波长范围为300nm~400nm;当望远镜的滤光片为Johnson-Cousins中的Rc波段时,设置Rc积分波长范围为550nm~800nm;当望远镜的滤光片为Johnson-Cousins中的Ic波段时,设置Rc积分波长范围为550nm~800nm。得到这三个波段的理想的光子NP_u',NP_Rc和NP_Ic如下式:
Figure BDA0003048952300000071
二:计算大气消光因子。
不同的台址其大气的消光系数不同,大气消光取决于空气密度和辐射通过的大气层厚度,并且大气消光对不同波段的光线影响有很大不同,由于红外波段的穿透力较强,波长较长的波段其大气消光的影响较小,而紫外波段大气消光的影响较大。另外大气消光与观测的天顶距(Z)有很大关系,当望远镜指向天顶时,(天顶距Z=0°,天顶距Z=90°-俯仰角EL),大气消光的影响最小,而随着望远镜的俯仰角变小,望远镜逐渐指向水平,天顶距Z变大,大气消光的影响越加明显,即星等变大,亮度降低的现象。这也是很多望远镜都将台址选在大气稀薄,海拔较高地方的原因。
根据大气消光的理论可以得到大气消光的公式:
m(z)=m0+K*F(z) 式:(4)
其中m(z)为天体在天顶距z处的大气内星等,m0为天体的大气外星等,K为消光因子,F(z)为大气质量,它与天顶距Z的角度有直接的关系。
我们采用的大气质量F(z)的计算公式为:
F(Z)=secz-0.0018167(secz-1)-0.002875(secz-1)2-0.0008083(secz-1)3 式:(5)
sec表示对天顶距Z取正割运算。
天顶距z处的大气内星等m(z)可以用仪器星等minst+加ZP(星等零点)表示;图像中的仪器星等的计算公式为:
minst=-2.5*log10Count 式:(6)
Count为星在图像中的Count值;
因此整体计算公式可以表示为:
m0=minst+ZP-K*F(z) 式:(7)
为了获得各个波段的大气消光系数K,在实测过程中,改变天顶距Z,即望远镜B的俯仰角,对多个标准星目标进行CCD成像,查得其查表星等m0,并计算其仪器星等minst,再通过上式进行拟合,计算出各波段的大气消光系数K。
以这三波段中对大气消光最敏感的u'波段为例,在西藏阿里5100m海拔台址上实测的部分结果(3个天体目标)如下:
Figure BDA0003048952300000081
按式子进行拟合并且对不同的目标取均值得到u’波段的大消光系数:
K_u’=0.37mag/airmass。
按照同样的观测及计算方法,拟合得到Rc,Ic波段的大气消光系数K_Rc=0.061mag/airmass,K_Ic=0.026mag/airmass。可见紫外波段大气穿透力较弱,大气消光的影响最为明显。
由此可得各个波段大气消光星等表达式为,F(z)为与天顶距相关的大气质量。
Figure BDA0003048952300000082
三:将台址仪器星等归算到大气外星等,并计算其Count(ADU)值
在图像中,通过计算某个天体的流量count值(ADU)可以计算出某个天体的仪器星等,如式(6)所示。此计算的公式是利用了某一处台址的测量值,这其中的仪器星等包含了大气消光的影响。下面我们用台址实际测量的大气消光系数去排除大气消光对仪器星等的影响。
由式(5)可知,大气消光只和大气质量和消光系数有关,而大气消光的系数已按照上述的方法求得,所以可以通过下式将台址上测量得到得仪器星等归算到大气外F(z)=0(大气质量为零)时的仪器星等。
M0inst-band=Minst-band-F(Z)*Kband 式:(9)
M0inst-band为归算到大气层外后的天体仪器星等(大气外仪器星等),Minst-band为台址上某一个大气质量下测量的天体仪器星等,F(z)为当时的大气质量(只与天顶距有关),Kband为某一个波段的大气消光系数(需要通过实测得到)。
根据式(6)可以推导台址上各个波段所观测的含有大气消光影响下的天体目标其Count数Countband
Figure BDA0003048952300000091
所以如果排除大气消光的影响,在得到归算后的大气外仪器星等之后,可以推导出在大气外(airmass=0)外,某一天体的理想Count值Count0band的计算公式为:
Figure BDA0003048952300000092
以西藏HinOTORI望远镜的Johnson-Cousins Ic波段为例,其大气质量与仪器的关系,以及归算到大气外的仪器星等部分数据如下所示。
Figure BDA0003048952300000093
Figure BDA0003048952300000101
从中发现,将星等和Count数归算到大气之后,由于排除了大气消光的影响,大气外的星等变亮,Count数增多。
四:计算大气外电子数e-值:
得到理想的Count数之后即可得到理想的电子数。CCD的作用是将收集到的光子数转化为电子数并进行数字化处理。CCD将光子转换为电子的参数为CCD的增益(gain),其物理上的单位是(e-/ADU),即一个光子转换为多少个电子的能力。计算公式为:
Figure BDA0003048952300000102
e-为电子数,ADU为Count数值。此参数一般在CCD相机出厂时经过测试。
已知CCD的增益以及理想的大气外Count0值之后,可以采用式(13)计算理想的大气外CCD转换生成的电子数e-0:
e-0=Count0*gain 式:(13)
在大气质量airmass=1.67mag/airmass时,Ic波段的Count数和电子数部分统计如下:
Figure BDA0003048952300000103
Figure BDA0003048952300000111
五:计算望远镜整体光电效率η,并与仿真图形比对
到此为止,通过实际观测值可以得到某一个望远镜观测某一个天体时大气外的电子数,又由式(3)已经计算了大气外的理想光子数,所以望远镜的整体光电效率η的计算公式为:
Figure BDA0003048952300000112
Ne-为整体接收的电子数,其计算公式为(11)和(13);Np为整体接收的光子数,其计算公式为式(3)。
望远镜指向NGC188天区进行跟踪拍摄,得到三个波段的光电效率。部分数据结果如下:
Figure BDA0003048952300000113
以上是实际观测值推导出的结论,与图2,图3及图4拟合后的结果进行对比,观测所得出的结果基本与理论拟合出的结构相符合。这里计算η所用的Ne-和Np都是归算到了大气外的物理量,即光电效率是望远镜本身的元件性能决定的,排除了观测时大气透过率的影响,与观测时的台址条件没有关系,因此验证更加准确。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。

Claims (8)

1.天文光学望远镜光电效率计算方法,其特征是:包括以下步骤:
步骤一、统计天文光学望远镜各透光元器件在不同波段的透过率数据,根据数据,绘制成各透光元器件透过率-波长曲线;
步骤二、采用matlab进行编程,先加载这些透过率-波长曲线的数据,再对各透过率-波长曲线进行插值提取出各透光元器件在该插值处的透过率的值,将这些值进行相乘,得到插值处的总体效率,最后将总体效率拟合成透光元器件综合透过率曲线,进而获得天文光学望远镜各波长处的整体光电效率。
2.根据权利要求1所述的天文光学望远镜光电效率计算方法,其特征是:步骤一中,天文光学望远镜的透光元器件包括望远镜主镜,望远镜副镜,各分色镜、Ic波段镀膜以及CCD探测器,相应的,各透光元器件透过率-波长曲线分别包括望远镜主镜和副镜的透过率-波长曲线、分色镜的透过率-波长曲线、Ic波段镀膜的透过率-波长曲线以及CCD量子效率-波长曲线。
3.天文光学望远镜光电效率实测验证方法,其特征是:包括以下步骤:
步骤a、定义采用权利要求1所述的透光元器件综合透过率曲线所得到的为望远镜整体效率ηA,用于观测实测验证的为望远镜效率ηB;选定目标星体作为观测对象;
步骤b、在大气外条件下,计算望远镜接收到目标星体的理想光子数;
步骤c、用望远镜在天文台上进行对目标星体进行实测,选取不同的天顶距进行成像,对图像进行处理,计算出各个波段上的大气消光因子K;
步骤d、利用大气消光因子K,将台址得到的目标星体的仪器星等归算到大气外,获得大气外的目标星体的仪器星等;
步骤e、计算出望远镜的CCD图像的Count值,再根据Count值计算CCD探测器所产生出的电子数,
步骤f、通过步骤b计算出的理想光子数以及步骤e得到的电子数,计算得到望远镜整体光电效率的值ηB,
步骤g、将望远镜的整体光电效率的值ηB与望远镜光电效率ηA进行比较,确认二者数值是否接近。
4.根据权利要求3所述的天文光学望远镜光电效率实测验证方法,其特征是:步骤b中,在大气外条件下,计算望远镜接收到目标星体的理想光子数的具体方法为:
Figure FDA0003048952290000021
其中NP为理想光子数,A为望远镜的集光面积,h为普朗克常数(6.626196×10-27erg·s)mAB为目标天体的AB星等,星等采用大气外的查表星等,积分项是一个与波长相关的表达式,R(λ)表示在波长积分范围内的整体透过率;
为了计算积分项的数值,在matlab中进行编程计算,在Johnson-Cousins波长系统中的Ic波段,设定Ic波段积分的范围为700nm至900nm,采用sum函数得出R(λ)在700nm~900nm范围内积分项的值为0.089,代入上式(1),从而得到望远镜在单位时间内接收到的一个固定星等星的光子数公式为:
NP_Ic=8.426*108*10-0.4mAB 式:(2)
公式(2)对于Sloan波长系统同样适用,当望远镜的滤光片为Sloan中的u’波段时,设置u’积分波长范围为300nm~400nm;当望远镜的滤光片为Johnson-Cousins中的Rc波段时,设置Rc积分波长范围为550nm~800nm;当望远镜的滤光片为Johnson-Cousins中的Ic波段时,设置Rc积分波长范围为550nm~800nm,得到这三个波段的理想的光子NP_u',NP_Rc和NP_Ic如下式:
Figure FDA0003048952290000022
5.根据权利要求4所述的天文光学望远镜光电效率实测验证方法,其特征是:步骤c中,计算出各个波段上的大气消光因子K的具体方法为:
根据大气消光的理论可以得到大气消光的公式:
m(z)=m0+K*F(z) 式:(4)
其中m(z)为天体在天顶距z处的大气内星等,m0为天体的大气外星等,K为消光因子,F(z)为大气质量,它与天顶距Z的角度有直接的关系;
大气质量F(z)的计算公式为:
F(Z)=secz-0.0018167(secz-1)-0.002875(secz-1)2-0.0008083(secz-1)3式:(5)
sec表示对天顶距Z取正割运算;
天顶距z处的大气内星等m(z)可以用仪器星等minst+加常数ZP表示,ZP为星等零点;图像中的仪器星等的计算公式为:
minst=-2.5*log10Count 式:(6)
Count为星在图像中的Count值;
因此大气外星等m0的整体计算公式可以表示为:
m0=minst+ZP-K*F(z) 式:(7)
为了获得各个波段的大气消光系数K,在实测过程中,改变天顶距Z,即望远镜的俯仰角,对多个标准星目标进行CCD成像,通过查表可得m0,并计算其仪器星等minst,再通过上式进行拟合,计算出各波段的大气消光系数K。
6.根据权利要求5所述的天文光学望远镜光电效率实测验证方法,其特征是:步骤d中,将台址得到的目标星体的仪器星等归算到大气外的具体方法为:
M0inst-band=Minst-band-F(Z)*Kband 式:(9)
M0inst-band为归算到大气层外后的天体仪器星等,Minst-band为台址上某一个大气质量下测量的天体仪器星等,F(Z)为当时的大气质量,其与天顶距有关,Kband为某一个波段的大气消光系数,通过实测得到。
7.根据权利要求6所述的天文光学望远镜光电效率实测验证方法,其特征是:步骤e、计算出望远镜的CCD图像的Count值,再根据Count值计算CCD探测器所产生出的电子数的具体方法为:
根据式(6)可以推导台址上各个波段所观测的含有大气消光影响下的天体目标其Count数Countband计算公式为
Figure FDA0003048952290000031
所以如果排除大气消光的影响,在得到归算后的大气外仪器星等之后,可以推导出在大气外,某一天体的理想Count值Count0band的计算公式为:
Figure FDA0003048952290000032
得到理想的Count数之后即可得到理想的电子数,计算公式为:
Figure FDA0003048952290000033
e-为电子数,ADU为Count数值。此参数在CCD相机出厂时经过测试;
已知CCD的增益以及理想的大气外Count0值之后,可以采用式(13)计算理想的大气外CCD转换生成的电子数e-0:
e-0=Count0*gain 式:(13)。
8.根据权利要求7所述的天文光学望远镜光电效率实测验证方法,其特征是:步骤f、通过步骤b计算出的理想光子数以及步骤e得到的电子数,计算得到望远镜整体光电效率的值η的具体方法为:
Figure FDA0003048952290000041
Ne-为整体接收的电子数,其计算公式为(11)和(13);Np为整体接收的光子数,其计算公式为式(3)。
CN202110482419.5A 2021-04-30 2021-04-30 天文光学望远镜光电效率计算及实测验证方法 Active CN113281004B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110482419.5A CN113281004B (zh) 2021-04-30 2021-04-30 天文光学望远镜光电效率计算及实测验证方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110482419.5A CN113281004B (zh) 2021-04-30 2021-04-30 天文光学望远镜光电效率计算及实测验证方法

Publications (2)

Publication Number Publication Date
CN113281004A true CN113281004A (zh) 2021-08-20
CN113281004B CN113281004B (zh) 2022-08-02

Family

ID=77277998

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110482419.5A Active CN113281004B (zh) 2021-04-30 2021-04-30 天文光学望远镜光电效率计算及实测验证方法

Country Status (1)

Country Link
CN (1) CN113281004B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1999014578A1 (en) * 1997-09-16 1999-03-25 Labsphere Inc. Ultraviolet transmittance analyzing method and instrument
US20020084430A1 (en) * 2000-11-09 2002-07-04 Canesta, Inc. Methods for CMOS-compatible three-dimensional image sensing using quantum efficiency modulation
CN1982916A (zh) * 2006-04-18 2007-06-20 中国科学院安徽光学精密机械研究所 机载大气环境探测激光雷达及其探测方法
CN101251436A (zh) * 2008-03-28 2008-08-27 中国科学院上海技术物理研究所 卡塞格林二反射镜光学系统加工在线检验方法
CN102680209A (zh) * 2011-02-24 2012-09-19 泰勒斯公司 用于望远镜的简化检验平台及自检验望远镜
CN103424380A (zh) * 2013-08-05 2013-12-04 中国人民解放军63655部队 一种离轴式大气湍流强度廓线实时测量装置与方法
CN104215430A (zh) * 2014-08-26 2014-12-17 中国科学院长春光学精密机械与物理研究所 一种大口径光电望远镜主镜反射面光学效率检测方法
CN108957426A (zh) * 2018-07-17 2018-12-07 中科院合肥技术创新工程院 一种激光雷达光电探测系统探测性能测试方法及装置
CN109308297A (zh) * 2018-11-20 2019-02-05 中国科学院国家天文台 一种基于标定数据库和插值计算的望远镜反射面控制算法
CN110793754A (zh) * 2019-11-01 2020-02-14 中国科学院光电技术研究所 一种基于相移调制的拼接式望远镜系统偏心误差探测方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1999014578A1 (en) * 1997-09-16 1999-03-25 Labsphere Inc. Ultraviolet transmittance analyzing method and instrument
US20020084430A1 (en) * 2000-11-09 2002-07-04 Canesta, Inc. Methods for CMOS-compatible three-dimensional image sensing using quantum efficiency modulation
CN1982916A (zh) * 2006-04-18 2007-06-20 中国科学院安徽光学精密机械研究所 机载大气环境探测激光雷达及其探测方法
CN101251436A (zh) * 2008-03-28 2008-08-27 中国科学院上海技术物理研究所 卡塞格林二反射镜光学系统加工在线检验方法
CN102680209A (zh) * 2011-02-24 2012-09-19 泰勒斯公司 用于望远镜的简化检验平台及自检验望远镜
CN103424380A (zh) * 2013-08-05 2013-12-04 中国人民解放军63655部队 一种离轴式大气湍流强度廓线实时测量装置与方法
CN104215430A (zh) * 2014-08-26 2014-12-17 中国科学院长春光学精密机械与物理研究所 一种大口径光电望远镜主镜反射面光学效率检测方法
CN108957426A (zh) * 2018-07-17 2018-12-07 中科院合肥技术创新工程院 一种激光雷达光电探测系统探测性能测试方法及装置
CN109308297A (zh) * 2018-11-20 2019-02-05 中国科学院国家天文台 一种基于标定数据库和插值计算的望远镜反射面控制算法
CN110793754A (zh) * 2019-11-01 2020-02-14 中国科学院光电技术研究所 一种基于相移调制的拼接式望远镜系统偏心误差探测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
崔洪鲁 等: "基于HITRAN数据库的大气激光雷达信号仿真", 《空间科学学报》 *
崔洪鲁 等: "基于HITRAN数据库的大气激光雷达信号仿真", 《空间科学学报》, vol. 40, no. 6, 31 December 2020 (2020-12-31), pages 1046 - 1051 *

Also Published As

Publication number Publication date
CN113281004B (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
US9651763B2 (en) Co-aperture broadband infrared optical system
Mao et al. A multichannel visible spectroscopy system for the ITER-like W divertor on EAST
CN101963529A (zh) 极远紫外扫描成像光谱仪
CN103995344B (zh) 透射式非制冷被动消热差长波红外光学系统
CN109254384A (zh) 一种星敏感器小型化光学系统
CN104656235B (zh) 一种长焦距宽矩形视场的离轴三反光学系统
CN112683796A (zh) 一种基于地球同步轨道观测的差分吸收光谱仪光学系统
CN113538314B (zh) 四波段共光轴光电成像平台及其图像融合处理方法
CN109239898B (zh) 一种紧凑型同轴折反射式望远物镜
Borovicka et al. Data on 824 fireballs observed by the digital cameras of the European Fireball Network in 2017-2018. I. Description of the network, data reduction procedures, and the catalog
CN113281004B (zh) 天文光学望远镜光电效率计算及实测验证方法
Komiyama et al. A photometric and spectroscopic study of dwarf and giant galaxies in the Coma Cluster. I. Wide-area photometric survey: Observation and data analysis
Bowen Absolute radiometric calibration of the IKONOS sensor using radiometrically characterized stellar sources
CN113325578B (zh) 一种光电吊舱的光学系统
CN115683358A (zh) 基于制冷型红外成像系统的辐射定标装置及定标方法
CN210376857U (zh) 一种高精度小型化长焦距星敏感器光学系统
CN113834639A (zh) 基于像素级偏振片的消畸变超微纳光学系统
EP3799620B1 (en) System for the real-time high precision measurement of the atmospheric attenuation of electromagnetic radiation
CN208013534U (zh) 一种宽幅低杂光全天时星跟踪器光学结构
JP7170401B2 (ja) 光源角度測定装置、光源位置検出装置、並びに人工衛星
CN110703410A (zh) 一种无遮拦长焦距星敏感器光学系统
CN113390508B (zh) 光谱-空间分辨像质优化的短波红外成像方法与装置
CN117109735B (zh) 一种半球空间月球微光辐照度仪光学系统及设计方法
RU2798087C1 (ru) Зеркально-линзовый объектив
Scuderi et al. Feasibility study of a camera for the UV Italian Sky Surveyor (UVISS) on the International Space Station

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 210023 No. 10, Yuanhua Road, Qixia District, Nanjing, Jiangsu Province

Applicant after: PURPLE MOUNTAIN OBSERVATORY, CHINESE ACADEMY OF SCIENCES

Address before: 210008 No. 2 West Beijing Road, Gulou District, Jiangsu, Nanjing

Applicant before: PURPLE MOUNTAIN OBSERVATORY, CHINESE ACADEMY OF SCIENCES

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant