CN106989832A - 用于窄带滤光器不同视场波长漂移的测量装置和改正方法 - Google Patents
用于窄带滤光器不同视场波长漂移的测量装置和改正方法 Download PDFInfo
- Publication number
- CN106989832A CN106989832A CN201710176328.2A CN201710176328A CN106989832A CN 106989832 A CN106989832 A CN 106989832A CN 201710176328 A CN201710176328 A CN 201710176328A CN 106989832 A CN106989832 A CN 106989832A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- light source
- wave length
- scaling
- light filter
- 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
Links
- 230000000007 visual effect Effects 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005259 measurement Methods 0.000 title claims abstract description 24
- 238000001228 spectrum Methods 0.000 claims abstract description 117
- 238000009792 diffusion process Methods 0.000 claims abstract description 48
- 230000003595 spectral effect Effects 0.000 claims description 18
- 230000004313 glare Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 8
- 238000013480 data collection Methods 0.000 claims description 6
- 239000005338 frosted glass Substances 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 4
- 239000011022 opal Substances 0.000 claims description 4
- 238000000862 absorption spectrum Methods 0.000 claims description 3
- 238000001093 holography Methods 0.000 claims description 3
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 3
- 229910052753 mercury Inorganic materials 0.000 claims description 3
- 201000009310 astigmatism Diseases 0.000 claims description 2
- 238000000295 emission spectrum Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 235000007926 Craterellus fallax Nutrition 0.000 description 1
- 240000007175 Datura inoxia Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000571 coke Substances 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明实施例涉及一种用于窄带滤光器不同视场波长漂移的测量装置和改正方法,该方法包括:将定标光源和扩散片切入窄带滤光器所在的光路;通过所述窄带滤光器扫描所述定标光源的光谱;根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移;将所述定标光源和扩散片移出所述光路;通过所述窄带滤光器扫描观测目标的光谱;根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息。本发明实施例可用于天地基的小口径大视场的太阳观测,如全日面光球色球像,也可用于天地基的日冕仪观测,最终实现精确的波长和速度场本底标定,还可用于基于窄带滤光器的小口径大视场的对地遥感观测。
Description
技术领域
本发明涉及遥感成像领域,包括对地或天文目标的遥感成像,尤其涉及一种用于窄带滤光器不同视场波长漂移的测量装置和改正方法。
背景技术
在遥感观测领域,如空间对地观测或空间、地基的天文观测,经常用窄带滤光器在某条吸收线或发射线开展窄带成像观测,以得到观测目标更为丰富的信息。比如太阳观测中,可见光波段连续谱均形成在光球,若想实现太阳高层大气如色球的观测,只能在色球吸收线的线心附近观测。
绝大多数窄带滤光器如窄带的干涉滤光片、Lyot滤光器、法布里-珀罗等都是基于干涉原理实现窄带滤光。此类器件在使用时遇到的一个问题就是不同方向的光线在透过窄带滤光器后的仪器轮廓峰值不同。例如,在太阳观测中,滤光器要么放在准直光路,要么放在远心光路。在准直光路中,探测器接收的不同视场的信息经过窄带滤光器的光程并不相同,这会造成不同视场的滤光器透过峰值位置不同,一般称为不同视场的波长漂移。对于远心光路,尽管不同视场经过滤光器的角度一样,理论上不存在不同视场的波长漂移,实际上受光学表面加工精度的影响,整个光路内存在波长漂移。
窄带滤光器引起的波长漂移对观测数据的波长定标十分重要,因为波长定标不准会导致由波长信息推导的物理参数如速度场精度不高。对于小视场大口径的太阳观测,因观测视场较小,一般选择日面中心的宁静区,随机晃动望远镜,多帧叠加后得到均匀的不同视场波长一致的定标数据,用于推算窄带滤光器不同视场引起的波长漂移。
但是对于大口径小视场的太阳观测,如全日面光球或色球、日冕观测,无法推算窄带滤光器不同视场引起的波长漂移。
发明内容
本发明实施例提供一种用于窄带滤光器不同视场波长漂移的测量装置和改正方法,可以用于天地基的小口径大视场的太阳观测,如全日面光球色球像,也可以用于天地基的日冕仪观测,最终实现精确的波长和速度场本底标定,还可用于基于窄带滤光器的小口径大视场的对地遥感观测。
第一方面,本发明实施例提供了一种用于窄带滤光器不同视场波长漂移的改正方法,该方法包括:将定标光源和扩散片切入窄带滤光器所在的光路,其中,所述定标光源用于生成与所述窄带滤光器的观测谱线一致的光谱,所述扩散片用于将所述定标光源自身的空间特征平滑,得到具有发射角的均匀的面光源,以及减少所述定标光源自身速度引起的波长漂移;通过所述窄带滤光器扫描所述定标光源的光谱;根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移;将所述定标光源和扩散片移出所述光路;通过所述窄带滤光器扫描观测目标的光谱;根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息,其中,所述观测目标的真实光谱为改正所述扫描观测目标光谱的波长漂移后的光谱。
在一个可能的实施例中,所述根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移,包括:通过探测器接收所述定标光源的光谱;对所述探测器的每个探测点得到的所述定标光源的光谱通过高斯轮廓或内插方法进行拟合,得到所述每个探测点对应的所述定标光源光谱的线心位置;根据所述定标光源光谱的线心位置确定所述窄带滤光器在不同视场下的波长漂移。
在一个可能的实施例中,所述根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息,包括:通过所述探测器接收所述观测目标的光谱;对所述探测器的每个探测点得到的所述观测目标的光谱通过高斯轮廓或内插方法进行拟合,得到所述每个探测点对应的所述观测目标光谱的线心位置;在所述每个探测点得到的所述观测目标光谱的线心位置的基础上,减去所述窄带滤光器在不同视场下的波长漂移,得到所述观测目标的真实光谱和速度场信息。
在一个可能的实施例中,所述高斯轮廓拟合通过如下公式实现:
其中,y(λ)为高斯轮廓函数,a为高斯曲线的峰高,λc为高斯曲线的峰位置,b为高斯曲线的半宽高,c为常数项,所述峰位置λc代表所述定标光源光谱或所述观测目标光谱谱线线心的位置。
第二方面,本发明实施例提供了一种用于窄带滤光器不同视场波长漂移的测量装置,该包括:定标光源、扩散片、望远系统、窄带滤光器、探测器以及数据采集和处理系统。所述定标光源,用于生成与所述窄带滤光器的观测谱线一致的光谱;所述扩散片,用于将所述定标光源自身的空间特征平滑,得到具有发射角的均匀的面光源,以及减少所述定标光源自身速度引起的波长漂移;所述望远系统和所述窄带滤光器,用于扫描所述定标光源或观测目标光谱;所述探测器,用于接收所述定标光源或所述观测目标的光谱;所述数据采集和处理系统,用于采集所述探测器的每个探测点的数据,以及根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移,并根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息,其中,所述观测目标的真实光谱为改正所述扫描观测目标光谱的波长漂移后的光谱。
在一个可能的实施例中,所述定标光源为人造光源或天然光源,其中,所述人造光源为可调谐激光器、定标汞灯中的至少一种;所述天然光源为太阳上的多个吸收谱或发射谱中的至少一种。
在一个可能的实施例中,所述扩散片为透射式或反射式扩散片中的一种,所述扩散片具体为大角度的乳白玻璃、毛玻璃、小角度的全息或工程散射片中的至少一种。
在一个可能的实施例中,所述窄带滤光器为Lyot双折射滤光器、法布里-珀罗滤光器或可调谐的干涉滤光片中的一种。
在一个可能的实施例中,该装置还包括:消杂散光光阑,所述消杂散光光阑用于将所述扩散片引入的额外杂散光减少,提高所述窄带滤光器不同视场波长漂移的测量精度。
基于上述技术方案,本发明实施例提供的用于窄带滤光器不同视场波长漂移的测量装置和改正方法,针对小口径大视场的窄带成像遥感观测,通过扩散片将定标光源(天然或人工发射或吸收线光源)扩散,消除定标光源空间强度、空间速度场的光谱不均匀性,扩散光经望远系统和窄带滤光器系统后照亮整个探测器靶面,通过消杂散光光阑将引入扩散片后引入的额外杂散光消除到最小程度,通过窄带滤光器的波长调节装置扫描扩散光的谱线轮廓,然后通过内插或高斯拟合的方法拟合出探测器不同视场处的真实线心位置,最终得到滤光器在不同视场处的波长漂移。该发明涉及的扩散片可以为透射式的,也可以为反射式的;扩散片可以为乳白玻璃、毛玻璃、也可以为全息或工程散射片。
本发明实施例提供的测量和处理方法简便易行,可广泛用于法布里-珀罗、Lyot类型的窄带滤光器不同视场的波长漂移测量,对基于窄带成像观测设备的高精度波长定标和速度场测量具有十分重要的意义。本发明既可以用于地面窄带成像观测的常规定标装置,精心设计后也可以用于空间窄带成像观测的星上定标装置。不仅可以用于小口径大视场的太阳观测,如全日面光球色球像,也可以用于日冕仪窄带观测不同视场的波长漂移,最终实现精确的波长和速度场本底标定。此外,本发明也可用于改正扫描成像光谱仪(如光栅光谱仪)引起的波长漂移,实现精确的波长和速度场标定。
附图说明
图1为本发明实施例提供的用于窄带滤光器不同视场波长漂移的测量装置结构示意图;
图2为本发明实施例提供的用于窄带滤光器不同视场波长漂移的改正方法流程示意图;
图3为本发明实施例提供的测量装置扫描的两维光谱图;
图4为本发明实施例提供的不同视场窄带滤光器波长漂移的分布示意图;
图5为本发明实施例提供的原始测量的全日面速度场示意图;
图6为本发明实施例提供的改正波长漂移后的全日面速度场示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为便于对本发明实施例的理解,下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
图1为本发明实施例提供的一种用于窄带滤光器不同视场波长漂移的测量装置结构示意图。如图1所示,该实施例提供的装置包括:定标光源110、扩散片120、望远系统130、消杂散光光阑140、窄带滤光器150、探测器160、数据采集和处理系统170。
如图1所示,该测量装置又可称为定标装置,用于标定窄带滤光器在不同视场下的波长漂移。
其中,窄带滤光器150可以为窄带可调谐滤光器。
该装置各部分的作用介绍如下:
定标光源110的作用在于生成和窄带滤光器常规观测所有谱线一致的光谱。定标光源可以为人造光源,如可调谐激光器、定标汞灯等;也可以为天然光源,如太阳上存在许多吸收谱和发射谱。
扩散片120的作用在于将定标光源自身的空间特征平滑掉,得到具有一定发射角的较为均匀的面光源。扩散片的另外一个作用在于模糊天然定标光源自身速度引起的波长漂移,如太阳本身存在速度场,经过扩散片后,扩散光源所形成的视场内可以消除太阳速度场导致的波长漂移。为保证在观测视场内扩散后的定标光源亮度和观测目标的亮度尽可能一致,扩散片可选择的类型为大角度的乳白玻璃、毛玻璃;小角度的全息或工程散射片。扩散片的类型可以为透射式的,也可以为反射式的。
望远系统130的作用是将观测目标放大后成像在探测器上,放大率和望远系统的焦比有关。
消杂散光光阑140的作用在于将扩散片120引入的额外杂散光减少到最小程度,提高窄带滤光器150不同视场波长漂移的测量精度。
窄带滤光器150的作用是在极窄的光谱范围内观测观测目标的信息,滤光器透过带半宽一般小于观测目标的谱线半宽,可以为Lyot双折射滤光器、法布里-珀罗滤光器,也可以为可调谐的干涉滤光片。如太阳色球常用谱线半宽达几常用的窄带滤光器透过带半宽小于可调谐的原因是实现波长扫描,便于通过高斯拟合或内插的方法计算谱线线心,从而推算波长漂移。
探测器160的作用在于接收经望远系统130和窄带滤光器150后成的图像,可以为电荷耦合器件(Charge-Coupled Device,CCD)或互补金属氧化物半导体器件(Complementary Metal Oxide Semiconductor,CMOS)相机。
数据采集和处理系统170的作用在于控制窄带滤光器150的波长扫描、记录存储波长扫描数据并进行数据预处理。
具体地,望远系统130和窄带滤光器150用于扫描定标光源110或观测目标110′光谱;探测器160用于接收定标光源110或所述观测目标110′的光谱。其中,观测目标110′为待观测的目标,可先通过定标光源的光谱确定窄带滤光器在不同视场的波长漂移,再对扫描得到的观测目标110′的光谱进行校准,得到观测目标110′真实的光谱和速度场信息。
具体地,数据采集和处理系统170用于采集所述探测器的每个探测点的数据,以及根据定标光源的光谱确定窄带滤光器在不同视场下的波长漂移,并根据窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到观测目标的真实光谱和速度场信息,其中,观测目标的真实光谱为改正所述扫描观测目标光谱的波长漂移后的光谱。
本发明实施例提供的用于窄带滤光器不同视场波长漂移的测量装置,通过扩散片将定标光源(天然或人工发射或吸收线光源)扩散,消除定标光源空间强度、空间速度场的光谱不均匀性,扩散光经望远系统和窄带滤光器系统后照亮整个探测器靶面,通过消杂散光光阑将引入扩散片后引入的额外杂散光消除到最小程度,通过窄带滤光器的波长调节装置扫描扩散光的谱线轮廓,然后通过内插或高斯拟合的方法拟合出探测器不同视场处的真实线心位置,最终得到窄带滤光器在不同视场处的波长漂移。
图2为本发明实施例提供的一种用于窄带滤光器不同视场波长漂移的改正方法流程示意图。如图2所示,包括步骤S201-步骤S206。
在步骤S201,将定标光源和扩散片切入窄带滤光器所在的光路。
其中,所述定标光源用于生成与所述窄带滤光器的观测谱线一致的光谱,所述扩散片用于将所述定标光源自身的空间特征平滑,得到具有发射角的均匀的面光源,以及减少所述定标光源自身速度引起的波长漂移。
在步骤S202,通过所述窄带滤光器扫描所述定标光源的光谱。
在步骤S203,根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移。
具体地,通过探测器接收所述定标光源的光谱;对所述探测器的每个探测点得到的所述定标光源的光谱通过高斯轮廓或内插方法进行拟合,得到所述每个探测点对应的所述定标光源光谱的线心位置;根据所述定标光源光谱的线心位置确定所述窄带滤光器在不同视场下的波长漂移。
在步骤S204,将所述定标光源和扩散片移出所述光路。
在步骤S205,通过所述窄带滤光器扫描观测目标的光谱。
在步骤S206,根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息。
其中,所述观测目标的真实光谱为改正所述扫描观测目标光谱的波长漂移后的光谱。
具体地,通过探测器接收所述观测目标的光谱;对所述探测器的每个探测点得到的所述观测目标的光谱通过高斯轮廓或内插方法进行拟合,得到所述每个探测点对应的所述观测目标光谱的线心位置;在所述每个探测点得到的所述观测目标光谱的线心位置的基础上,减去所述窄带滤光器在不同视场下的波长漂移,得到所述观测目标的真实光谱和速度场信息。
具体地,上述高斯轮廓拟合通过如下公式实现:
其中,y(λ)为高斯轮廓函数,a为高斯曲线的峰高,λc为高斯曲线的峰位置,b为高斯曲线的半宽高,c为常数项,所述峰位置λc代表所述定标光源光谱或所述观测目标光谱谱线线心的位置。
在可选的实施例中,可将本发明实施例提供的用于窄带滤光器不同视场波长漂移的改正方法划分为如下两部分:
第一部分,定标观测测量不同视场窄带滤光器的波长漂移:
1)将定标光源和扩散片切入光路。
2)开启窄带滤光器的波长扫描模式,扫描定标光源的光谱。
3)对探测器每个点得到的扫描光谱用高斯轮廓进行拟合或用内插方法,计算扫描光谱的线心位置,最终得到不同视场窄带滤光器的波长漂移。用到的高斯拟合函数如下,对于吸收线,可以将其取符号后加上最小值转换成高斯轮廓,拟合后再逆变换回吸收线轮廓,如公式(1),其中公式(1)峰位置λc代表的就是(定标光源的)观测谱线线心的位置。
第二部分,改正常规观测光谱或速度场中滤光器引起的波长漂移,得到观测目标真实的光谱和速度场:
1)将定标光源和扩散片移出光路,开启常规观测模式。
2)开启窄带滤光器的波长扫描模式,扫描观测目标的光谱。
3)对探测器每个点得到的扫描光谱用高斯轮廓进行拟合或用内插方法,计算扫描光谱的线心位置,得到观测视场观测目标的线心位置。
4)减去第一部分中得到的不同视场窄带滤光器的波长漂移,得到观测视场内观测目标真实的光谱和速度场信息。
相应地,图3为本发明实施例提供的测量装置扫描的两维光谱图。图3中右下角的子图表示图中心位置的光谱轮廓(星号点),实线代表对观测光谱轮廓的高斯拟合,竖线给出拟合的线心位置。
在一种可能的实施例中,以中国科学院怀柔太阳观测基地运行的全日面磁场望远镜为例详细说明本装置的使用方法。
具体地,全日面磁场望远镜的基本参数为:10cm口径,有效焦距0.77m,配置半宽的双折射滤光器,可以在532.4nm处观测太阳窄带光球单色像和矢量磁场,探测器像元数992×992,观测视场33′,常规观测偏振参量Stokes I V,Q,U的谱线位置为偏离线心为避免引起不同视场的波长漂移,窄带滤光器放在远心光路中,从下面的分析可以看出,因各种原因,仍然存在波长漂移,且对观测影响很大。此处我们仅仅以Stoeks I为例,通过扫描的Stokes I轮廓我们可以得到全日面的速度场。
第一部分,通过定标观测测量不同视场窄带滤光器的波长漂移:
将定标光源和扩散片切入光路(望远镜前),此时定标光源为太阳,为保证常规观测和定标观测曝光时间相差不大,观测视场33′(≈0.5°),扩散片选择1°的全息扩散片,此时在观测视场内光强衰减约1/10。
开启窄带滤光器的波长扫描模式,扫描范围仪器设定的线心0到扫描间隔共16帧。图3是扫描的不同波长处的二维光强图,图3中可见越靠近线心,强度越弱。
对于二维图像上的一个像元,由扫描观测均得到一个光谱轮廓,光谱分辨率和扫描间隔对应,此时是因光谱分辨率有限,我们利用高斯函数(或内插)的方法对观测轮廓进行了拟合,根据拟合结果,可以得到线心的真实位置。上述方法的结果见图3右下方子图,图中光谱轮廓对应空间位置为二维图像的中心点,观测轮廓见星号点,实线代表高斯函数拟合的轮廓,实竖线代表有高斯拟合函数得到的谱线线心位置,此处的值为由上述结果可知,设定的线心并不准确,真实的线心位置在左右。
利用上述方法对二维图像上的每个像元均做高斯拟合,可以得到扫描谱线的线心分布图,结果见图4。图4为不同视场窄带滤光器波长漂移的分布图,其中,最小值最大值真实线心位置左右。图中黑色区域为光阑位置,强度为零,没有参与高斯拟合过程。图中可以清晰的看出不同视场的光束经窄带滤光器后引起的波长漂移,对于33′的观测视场,拟合的线心范围为0.1到不同视场的波长偏移最大相差约
第二部分,常规观测测量太阳速度场:
望远镜对准观测目标,此时观测目标为太阳,开启窄带滤光器的波长扫描模式,扫描范围仪器设定的线心0到扫描间隔利用和上述同样的计算方法,可以计算出观测目标太阳本身的波长漂移和仪器本身引起的波长漂移。因此时的计算重点在于计算速度场,所以我们将波长漂移换算到速度场,换算公式为:v=Δλ×2.9979×108/5324,此时的波长单位速度单位m/s。计算的速度场分布见图5,图5为原始测量的全日面速度场,速度场分布区间为-3km/s到3km/s,太阳本身的速度场在东西边缘最大为2km/s。这也说明常规观测测量的速度场中夹杂这不同视场窄带滤光器自身引起的波长漂移,对比图5和图4也可以发现你,二者的分布有类似之处,这些相似之处正是不同视场窄带滤光器自身引起的波长漂移。
第三部分,扣除不同视场窄带滤光器的波长漂移得到真实的太阳速度场:
因常规观测得到的速度场中包含不同视场窄带滤光器自身引起的波长漂移。如果我们将常规观测得到的图5结果减去定标观测得到的图4结果,就可以扣除不同视场窄带滤光器的波长漂移得到真实的太阳速度场。改正结果见图6。图6为改正后的全日面速度场,其中,色表最小值:-1.5km/s,最大值2km/s。由图6可见太阳沿东西方向(横轴方向)速度场的变化,在东边缘太阳自转引起的蓝移最大(左侧),速度场为负,接近1.7km/s,(部分东边缘没有在观测视场内,所以速度场偏小);在日面西边缘太阳自转引起的红移最大,速度场为正,对应约2km/s,和我们迄今认识的太阳速度场结果基本一致。证实了我们基于1°的高斯扩散片改正不同视场滤光器波长漂移的有效性。
本领域技术人员应该可以意识到,在上述一个或多个示例中,本申请所描述的功能可以用硬件、软件、固件或它们的任意组合来实现。当使用软件实现时,可以将这些功能存储在计算机可读介质中或者作为计算机可读介质上的一个或多个指令或代码进行传输。计算机可读介质包括计算机存储介质和通信介质,其中通信介质包括便于从一个地方向另一个地方传送计算机程序的任何介质。存储介质可以是通用或专用计算机能够存取的任何可用介质。
最后应说明的是,以上实施例仅用以说明本申请的技术方案而非限制。尽管参照实施例对本申请进行了详细说明,本领域的普通技术人员应当理解,对本申请的技术方案进行修改或者等同替换,都不脱离本申请技术方案的精神和范围,其均应涵盖在本申请的权利要求范围当中。
Claims (9)
1.一种用于窄带滤光器不同视场波长漂移的改正方法,其特征在于,所述方法包括:
将定标光源和扩散片切入窄带滤光器所在的光路,其中,所述定标光源用于生成与所述窄带滤光器的观测谱线一致的光谱,所述扩散片用于将所述定标光源自身的空间特征平滑,得到具有发射角的均匀的面光源,以及减少所述定标光源自身速度引起的波长漂移;
通过所述窄带滤光器扫描所述定标光源的光谱;
根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移;
将所述定标光源和扩散片移出所述光路;
通过所述窄带滤光器扫描观测目标的光谱;
根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息,其中,所述观测目标的真实光谱为改正所述扫描观测目标光谱的波长漂移后的光谱。
2.根据权利要求1所述的方法,其特征在于,所述根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移,包括:
通过探测器接收所述定标光源的光谱;
对所述探测器的每个探测点得到的所述定标光源的光谱通过高斯轮廓或内插方法进行拟合,得到所述每个探测点对应的所述定标光源光谱的线心位置;
根据所述定标光源光谱的线心位置确定所述窄带滤光器在不同视场下的波长漂移。
3.根据权利要求2所述的方法,其特征在于,所述根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息,包括:
通过所述探测器接收所述观测目标的光谱;
对所述探测器的每个探测点得到的所述观测目标的光谱通过高斯轮廓或内插方法进行拟合,得到所述每个探测点对应的所述观测目标光谱的线心位置;
在所述每个探测点得到的所述观测目标光谱的线心位置的基础上,减去所述窄带滤光器在不同视场下的波长漂移,得到所述观测目标的真实光谱和速度场信息。
4.根据权利要求2或3所述的方法,其特征在于,所述高斯轮廓拟合通过如下公式实现:
其中,y(λ)为高斯轮廓函数,a为高斯曲线的峰高,λc为高斯曲线的峰位置,b为高斯曲线的半宽高,c为常数项,所述峰位置λc代表所述定标光源光谱或所述观测目标光谱谱线线心的位置。
5.一种用于窄带滤光器不同视场波长漂移的测量装置,其特征在于,包括:定标光源、扩散片、望远系统、窄带滤光器、探测器以及数据采集和处理系统;
所述定标光源,用于生成与所述窄带滤光器的观测谱线一致的光谱;
所述扩散片,用于将所述定标光源自身的空间特征平滑,得到具有发射角的均匀的面光源,以及减少所述定标光源自身速度引起的波长漂移;
所述望远系统和所述窄带滤光器,用于扫描所述定标光源或观测目标光谱;
所述探测器,用于接收所述定标光源或所述观测目标的光谱;
所述数据采集和处理系统,用于采集所述探测器的每个探测点的数据,以及根据所述定标光源的光谱确定所述窄带滤光器在不同视场下的波长漂移,并根据所述窄带滤光器在不同视场下的波长漂移校准所述观测目标的光谱,得到所述观测目标的真实光谱和速度场信息,其中,所述观测目标的真实光谱为改正所述扫描观测目标光谱的波长漂移后的光谱。
6.根据权利要求5所述的装置,其特征在于,所述定标光源为人造光源或天然光源,其中,所述人造光源为可调谐激光器、定标汞灯中的至少一种;所述天然光源为太阳上的多个吸收谱或发射谱中的至少一种。
7.根据权利要求5所述的装置,其特征在于,所述扩散片为透射式或反射式扩散片中的一种,所述扩散片具体为大角度的乳白玻璃、毛玻璃、小角度的全息或工程散射片中的至少一种。
8.根据权利要求5所述的装置,其特征在于,所述窄带滤光器为Lyot双折射滤光器、法布里-珀罗滤光器或可调谐的干涉滤光片中的一种。
9.根据权利要求5所述的装置,其特征在于,所述装置还包括:消杂散光光阑,所述消杂散光光阑用于将所述扩散片引入的额外杂散光减少,提高所述窄带滤光器不同视场波长漂移的测量精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710176328.2A CN106989832B (zh) | 2017-03-23 | 2017-03-23 | 用于窄带滤光器不同视场波长漂移的测量装置和改正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710176328.2A CN106989832B (zh) | 2017-03-23 | 2017-03-23 | 用于窄带滤光器不同视场波长漂移的测量装置和改正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106989832A true CN106989832A (zh) | 2017-07-28 |
CN106989832B CN106989832B (zh) | 2019-03-12 |
Family
ID=59412916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710176328.2A Expired - Fee Related CN106989832B (zh) | 2017-03-23 | 2017-03-23 | 用于窄带滤光器不同视场波长漂移的测量装置和改正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106989832B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490214A (zh) * | 2018-03-14 | 2018-09-04 | 中国科学院光电技术研究所 | 一种基于可调谐滤光器成像观测的太阳大气视向速度场测量方法 |
CN109374550A (zh) * | 2018-11-16 | 2019-02-22 | 中国科学院合肥物质科学研究院 | 基于透过率光谱的超分辨光谱仪光谱定标方法 |
CN111060289A (zh) * | 2019-12-27 | 2020-04-24 | 中国科学院长春光学精密机械与物理研究所 | 一种高灵敏度日冕仪杂光检测装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1288165A (zh) * | 1999-09-09 | 2001-03-21 | 三星电子株式会社 | 长周期光纤光栅滤光器装置 |
CN1720650A (zh) * | 2002-10-03 | 2006-01-11 | 英特尔公司 | 背面波长锁定装置的调整和装配方法 |
CN1961235A (zh) * | 2002-11-06 | 2007-05-09 | 阿兹纳公司 | 用于色散补偿光纤系统的功率源 |
CN101023388A (zh) * | 2004-02-20 | 2007-08-22 | 加利福尼亚技术学院 | 用于布拉格光栅可调谐滤光器的方法和装置 |
CN104267451A (zh) * | 2014-09-11 | 2015-01-07 | 李�远 | 一种光学扩散片 |
CN104062010B (zh) * | 2014-06-27 | 2015-08-19 | 杭州彩谱科技有限公司 | 一种优化定标算法的分光光源颜色照度测量仪器 |
CN105068359A (zh) * | 2015-08-21 | 2015-11-18 | 中国工程物理研究院流体物理研究所 | 一种中心波长可调的滤色镜及滤色方法 |
CN106330328A (zh) * | 2015-07-10 | 2017-01-11 | 青岛海信宽带多媒体技术有限公司 | 光模块及其光信号输出控制方法 |
-
2017
- 2017-03-23 CN CN201710176328.2A patent/CN106989832B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1288165A (zh) * | 1999-09-09 | 2001-03-21 | 三星电子株式会社 | 长周期光纤光栅滤光器装置 |
CN1720650A (zh) * | 2002-10-03 | 2006-01-11 | 英特尔公司 | 背面波长锁定装置的调整和装配方法 |
CN1961235A (zh) * | 2002-11-06 | 2007-05-09 | 阿兹纳公司 | 用于色散补偿光纤系统的功率源 |
CN101023388A (zh) * | 2004-02-20 | 2007-08-22 | 加利福尼亚技术学院 | 用于布拉格光栅可调谐滤光器的方法和装置 |
CN104062010B (zh) * | 2014-06-27 | 2015-08-19 | 杭州彩谱科技有限公司 | 一种优化定标算法的分光光源颜色照度测量仪器 |
CN104267451A (zh) * | 2014-09-11 | 2015-01-07 | 李�远 | 一种光学扩散片 |
CN106330328A (zh) * | 2015-07-10 | 2017-01-11 | 青岛海信宽带多媒体技术有限公司 | 光模块及其光信号输出控制方法 |
CN105068359A (zh) * | 2015-08-21 | 2015-11-18 | 中国工程物理研究院流体物理研究所 | 一种中心波长可调的滤色镜及滤色方法 |
Non-Patent Citations (2)
Title |
---|
徐稚等: "抚仙湖一米红外太阳望远镜Hα 窄带滤光器", 《天文研究与技术》 * |
欧阳艳东等: "新型体全息窄带滤光器的光学特性研究", 《功能材料与器件学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490214A (zh) * | 2018-03-14 | 2018-09-04 | 中国科学院光电技术研究所 | 一种基于可调谐滤光器成像观测的太阳大气视向速度场测量方法 |
CN109374550A (zh) * | 2018-11-16 | 2019-02-22 | 中国科学院合肥物质科学研究院 | 基于透过率光谱的超分辨光谱仪光谱定标方法 |
CN111060289A (zh) * | 2019-12-27 | 2020-04-24 | 中国科学院长春光学精密机械与物理研究所 | 一种高灵敏度日冕仪杂光检测装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106989832B (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Blackman et al. | Performance verification of the extreme precision spectrograph | |
de la Cruz Rodriguez et al. | CRISPRED: A data pipeline for the CRISP imaging spectropolarimeter | |
Hook et al. | The Gemini–North Multi‐Object Spectrograph: Performance in imaging, long‐slit, and multi‐object spectroscopic modes | |
Egels et al. | Digital photogrammetry | |
Golimowski et al. | The Solar Neighborhood. IX. Hubble Space Telescope Detections of Companions to Five M and L Dwarfs within 10 Parsecs of the Sun | |
US20090022359A1 (en) | Vegetation index image generation methods and systems | |
Hora et al. | Photometry using the infrared array camera on the Spitzer space telescope | |
CN106989832A (zh) | 用于窄带滤光器不同视场波长漂移的测量装置和改正方法 | |
Biller et al. | A high-contrast search for variability in HR 8799bc with VLT-SPHERE | |
Moiseev | Reduction of CCD observations made with a scanning Fabry–Perot interferometer. III. Wavelength scale refinement | |
CN107727231A (zh) | 一种对光谱成像仪条带噪声进行修正的方法 | |
Arribas et al. | Differential atmospheric refraction in integral-field spectroscopy: Effects and correction-Atmospheric refraction in IFS | |
Sihler et al. | In-operation field-of-view retrieval (IFR) for satellite and ground-based DOAS-type instruments applying coincident high-resolution imager data | |
RU2531024C1 (ru) | Способ дистанционного зондирования земли (дзз) | |
CN105044113B (zh) | 一种二氧化硫气体成像仪 | |
CN105136294A (zh) | 地基可见高光谱分辨率月球观测系统 | |
Schlichenmaier et al. | Effects of solar evolution on finite acquisition time of Fabry–Perot interferometers in high resolution solar physics | |
Caron et al. | The CarbonSat candidate mission: radiometric and spectral performances over spatially heterogeneous scenes | |
Erskine et al. | High-resolution broadband spectroscopy using externally dispersed interferometry at the Hale telescope: part 1, data analysis and results | |
Korhonen et al. | Ellipsoidal primary of the RS CVn binary ζ Andromedae-Investigation using high-resolution spectroscopy and optical interferometry | |
Weber et al. | Polarization upgrade of specMACS: calibration and characterization of the 2D RGB polarization-resolving cameras | |
Johanneson et al. | Very high spatial resolution two-dimensional solar spectroscopy with video CCDs | |
Weber et al. | Polarization upgrade of specMACS: calibration and characterization of the 2D RGB polarization resolving cameras | |
AU2001286870C1 (en) | Passive ranging to a target reflecting solar radiation | |
Clermont et al. | Out-of-field stray light correction in optical instruments: the case of Metop-3MI |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190312 |