CN102997994B - 一种基于人造光源的天空光光谱模拟方法 - Google Patents
一种基于人造光源的天空光光谱模拟方法 Download PDFInfo
- Publication number
- CN102997994B CN102997994B CN201210483543.4A CN201210483543A CN102997994B CN 102997994 B CN102997994 B CN 102997994B CN 201210483543 A CN201210483543 A CN 201210483543A CN 102997994 B CN102997994 B CN 102997994B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- light
- skylight
- metal halide
- sky
- 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.)
- Active
Links
Landscapes
- Non-Portable Lighting Devices Or Systems Thereof (AREA)
Abstract
本发明提供了一种基于人造光源的天空光光谱模拟方法，包含以下步骤：（1）利用MODTRAN辐射传输模型计算天空光的光谱特性；（2）选择金属卤化物光源作为400nm-2500nm范围内天空光光谱的主要模拟光源；（3）分段分析金属卤化物光源光谱与天空光光谱的异同；（4）使用大功率固体照明元件LED作为400nm-600nm范围内的光谱补偿光源；（5）选择相应滤光片，对600nm-2500nm范围内的金属卤化物光源光谱进行修正。该方法与同类系统相比，能够以较低的成本实现400nm-2500nm范围内的全光谱天空光模拟，适用于光学遥感物理或半物理仿真系统中天空光照明环境的研究。
Description
（一）技术领域
本发明涉及一种基于人造光源的天空光光谱模拟方法，属于光学遥感仿真技术领域，适用于光学遥感物理或半物理仿真系统中天空光照明环境的模拟研究。
（二）背景技术
遥感仿真技术是指通过数字或物理的方法对遥感过程的各个环节进行仿真实现。现有的遥感仿真技术主要有数字仿真、物理仿真和半物理仿真，其中数字仿真是通过数学的方法对遥感过程中涉及到的遥感器、大气传输等进行建模仿真；物理仿真是通过建立实物模型对遥感器、遥感平台、照明环境和观测目标等进行物理实现；半物理仿真技术是将数字仿真和物理仿真相结合，通过物理手段完成包含遥感器、遥感平台、照明环境和观测目标的观测场景模拟，通过数字手段完成大气传输等的模拟，可以在目标、环境、观测仪器均可控、可知的条件下，开展遥感系统性能指标和工作模式需求分析、遥感系统数据质量误差分析、系统工作模式和数据处理方法改进等工作。
照明环境作为遥感半物理仿真系统中物理仿真的重要组成部分，主要包括太阳光模拟和天空光模拟，二者相互匹配，共同完成遥感照明。太阳光模拟可用太阳模拟器实现，其模拟光源一般选用氙灯光源，国内以小型太阳模拟器为主，光谱为特制滤光片修正后的氙灯光谱；国外的大型太阳模拟器光谱为未修正的氙灯光谱。由于天空光是经大气散射后形成，其光谱特性与太阳光光谱特性有着密切的联系。关于天空光的光谱研究，国内主要使用非成像光谱仪在室外测量天空光光谱，分析某些光谱特征的异常，进而服务于大气污染监测等领域，仅有长光所在本世纪初为建立遥感成像模拟实验室提出了人造天空的构想，但目前尚无物理模拟实现方面的报道；国外只有美国Itek公司建立了可见光范围的天空光模拟设施，用于早期的多光谱成像仿真研究，但该设施未考虑天空光的光谱特性。
天空光作为遥感仿真系统中照明环境的重要组成部分，其光谱特性对获得的遥感数据质量有着非常重要的影响。因此，实现基于人造光源的天空光光谱模拟的最大挑战在于设计一种光源和其它光学元件的组合方式，实现400nm-2500nm范围内的天空光光谱特性模拟。
（三）发明内容
本发明的目的是提供一种基于人造光源的天空光光谱模拟方法，它可以实现400nm-2500nm范围内的天空光光谱特性模拟，服务于光学遥感物理或半物理仿真系统中天空光照明环境的模拟。
本发明的技术解决方案为：在利用大气辐射传输模型（moderatespectral resolution atmospheric transmittance code，以下简称MODTRAN）对天空光照明特性计算分析的基础上，选择合适的人造光源作为天空光的光谱模拟光源，通过对标准天空光光谱与所选光源光谱特征的比较，分段对天空光光谱进行模拟，在400nm-600nm，使用光谱补偿方法消除该范围内金属卤化物光源与天空光光谱的差别，在600nm-2500nm，使用合适的滤光片对该范围内金属卤化物光源与天空光光谱的差别进行修正，最终实现400nm-2500nm范围的全光谱天空光光谱特性模拟。
本发明一种基于人造光源的天空光光谱模拟方法，其步骤如下：
（1）利用MODTRAN辐射传输模型计算天空光的光谱特性；
（2）选择金属卤化物光源作为400nm-2500nm范围内天空光光谱的主要模拟光源；
（3）分段分析金属卤化物光源光谱与天空光光谱的异同；
（4）使用大功率固体照明元件LED作为400nm-600nm范围内的光谱补偿光源；
（5）选择相应滤光片，对600nm-2500nm范围内的金属卤化物光源光谱进行修正。
其中，步骤（1）中所述的“利用MODTRAN辐射传输模型计算天空光的光谱特性”，其含义说明如下：输入参数设置为中纬度夏季（北纬45°，7月）/乡村气溶胶模型VIS=23km/太阳高度角70°/天顶角45°/相对方位角45°，然后利用大气辐射传输方程计算大气中气体分子和气溶胶颗粒对光线的散射和吸收，从而获得天空光的光谱数据。
其中，步骤（2）中所述的“选择金属卤化物光源作为400nm-2500nm范围内天空光光谱的主要模拟光源”，其含义说明如下：金属卤化物光源的光谱范围为400nm-2500nm，光源的具体型号为飞利浦公司生产的CDM-R反射性陶瓷金属卤化物灯PAR30/35w/830。
其中，步骤（3）中所述的“分段分析金属卤化物光源光谱与天空光光谱的异同”，其含义说明如下：分段分析主要是针对400nm-600nm范围内金属卤化物光源的相对光谱能量所占的百分比较低，600nm-2500nm范围内金属卤化物光源的相对光谱能量所占的百分比较高，根据两者的不同，以600nm为分割点，将400nm-2500nm分为400nm-600nm和600nm-2500nm两个波段范围进行分段分析，并在下述步骤（4）和步骤（5）中分别进行补偿和修正。
其中，步骤（4）中所述的“使用大功率固体照明元件LED作为400nm-600nm范围内的光谱补偿光源”，其含义说明如下：LED光谱补偿光源的中心波长使用曲线分峰拟合的方法进行计算，其计算过程如下：由于LED的光谱曲线可近似认为符合高斯分布，因此分峰拟合的过程中，使用高斯函数作为LED的光谱分布曲线，其基本原理是假设目标光谱Y（v）（其中v是光谱频率）是由若干个单峰谱线相互叠加所形成，分峰拟合的任务是求解一组单峰谱线Xi（v）（i=1，2，...n），使得下式成立：
Y(v)=∑Xi(v)
但在实际情况下，上式两边不可能达到严格相等，在实际计算过程中，使目标光谱Y（v）与拟合光谱∑Xi(v)之间的误差在400nm-600nm范围内最小。
其中，步骤（5）中所述的“选择相应滤光片，对600nm-2500nm范围内的金属卤化物光源光谱进行修正”，其含义说明如下：滤光片的选择分两步完成，首先计算经步骤（4）补偿后的金属卤化物光源相对光谱能量分布与MODTRAN计算的标准天空光相对光谱能量分布在600nm-2500nm范围内的比值，再根据计算所得的滤光片透过率曲线选择相应的滤光片进行光谱修正。
本发明与现有技术相比的优点在于：补充了目前国内天空光物理仿真的空白，一定程度上可以克服传统的遥感试验只能在野外环境下才能获得自然光照明的局限，它具有以下优点：（1）利用金属卤化物光源，摒弃了常用的全光谱模拟光源——氙灯，与同类系统相比，可以很大程度地节约系统的体积和成本；（2）利用大功率LED对该金属卤化物光源400nm-600nm范围内的光谱进行补偿，操作安全且易于实现（3）使用光谱滤光片对该金属卤化物光源600nm-2500nm范围内的光谱进行修正，可使获得的模拟光谱更加接近真实天空光光谱（4）通过设计相应光谱补偿光源和光谱滤光片的参数要求，实现了400nm-2500nm的全光谱天空光模拟，打破了现有同类系统中不考虑光谱特性或只考虑局部光谱的功能障碍。
附图说明
图1为本发明所述方法流程框图。
图2为MODTRAN计算的天空光光谱。
图3为所选金属卤化物光源的光谱特征。
图4为所选滤光片的透过率曲线。
图4（a）为QB19滤光片透过率曲线；
图4（b）为GRB3滤光片透过率曲线。
图5为本发明中光谱补偿和光谱修正前后照明光谱对比。
具体实施方式
为了更好的说明本发明涉及的基于人造光源的天空光光谱模拟方法，利用Light tools光学设计软件对光源进行了建模仿真，模拟光源参数设置选用飞利浦公司生产的CDM-R反射性陶瓷金属卤化物，实现400nm-2500nm范围内的天空光光谱特性模拟。
本发明一种基于人造光源的天空光光谱模拟方法，见图1所示，具体实现步骤如下：
（1）利用MODTRAN辐射传输模型计算天空光的照明特性：MODTRAN计算天空光的参数设置为中纬度夏季（北纬45°，7月）/乡村气溶胶模型VLS=23km/太阳高度角70°/天顶角45°/相对方位角45°，利用大气辐射传输方程计算各层大气的消光、衰减以及大气中气体分子（水汽、CO2，O2等）和气溶胶颗粒对太阳光的吸收和散射，从而获得天空光的光谱数据，如图2所示，光谱范围为400nm-2500nm，光谱分辨率为10nm；
（2）选择金属卤化物光源作为400nm-2500nm范围内天空光光谱的主要模拟光源：为了保证所选光源在400nm-2500nm内的光谱能量分布接近标准天空光的光谱分布，对常见光源在400nm-2500nm内的光谱能量分布曲线与标准天空光光谱能量分布曲线进行了对比，通过分析最终选用飞利浦公司生产的CDM-R反射性陶瓷金属卤化物光源PAR30/35w/830作为天空光光谱的主要模拟光源，并使用SVC HR-1024光谱仪对该光源的光谱特性进行了测量，如图3所示。
（3）分段分析金属卤化物光源光谱与天空光光谱的异同：通过所选金属卤化物光源与MODTRAN辐射传输模型计算的天空光光谱对比分析可以看出：金属卤化物光源在400nm-600nm范围内的光谱能量所占的百分相对较低，600nm-2500nm范围内的光谱能量所占的百分比则比标准天空光高，为解决两者之间不同的变化趋势，将天空光光谱分为400nm-600nm、600nm-2500nm两段分别进行补偿和修正。
针对400nm-600nm金属卤化物光源光谱能量相对较低的部分，使用大功率固体照明元件LED光源进行光谱补偿；针对600nm-2500nm金属卤化物光源光谱能量相对较高的部分，则使用低透的滤光片进行光谱修正。通过光谱补偿和光谱滤光相结合的方法对天空光光谱进行分段模拟。
（4）使用大功率固体照明元件LED作为400nm-600nm范围内的光谱补偿光源：大功率固体照明元件LED的中心波长是使用Origin软件的曲线分峰拟合模块计算完成的，其基本原理是假设目标光谱Y（v）（其中v是光谱频率）是由若干个单峰谱线相互叠加所形成，分峰拟合的任务是找到一组单峰谱线Xi（v）（i=1，2，...n），使得下式成立：
Y(v)=∑Xi(v)
在实际情况下，上式两边不可能达到严格相等，因此在实际计算过程中，使目标光谱Y（v）与拟合光谱∑Xi(v)之间的误差最小。由于LED的光谱曲线可近似认为符合高斯分布，因此分峰拟合的过程中，使用高斯函数作为LED的光谱分布曲线，通过计算得到LED光谱补偿光源的中心波长为450nm、475nm、523nm和618nm，光谱半波宽为20nm。
（5）选择相应滤光片，对600nm-2500nm范围内的金属卤化物光源光谱进行修正：根据MODTRAN计算的标准天空光光谱曲线与补偿后的金属卤化物光源光谱曲线相比获得理想光谱滤光片的透过率要求，如表1所示，中心波长600nm处的透过率约为25％，半宽度约为200nm；400nm和900n m附近的透过率T≥80％；1000nm后的透过率约为70％。结合现有的较为成熟的滤光片加工工艺，选择图4（a）、（b）中所示江苏某玻片厂制作的QB19和GRB3组合滤光片对金属卤化物光源的光谱进行修正，光谱补偿和修正后的光谱特征如图5所示。
表1理想滤光片的透过率要求
通过图5中曲线的数据可以计算出，本发明设计的一种基于人造光源的天空光光谱特性模拟方法可实现400nm-2500nm范围内的天空光光谱模拟，其光谱匹配误差<35％,满足现有同类系统（如太阳模拟器等）对光谱匹配误差的基本要求。
Claims (6)
1.一种基于人造光源的天空光光谱模拟方法，其特征在于，它包含以下步骤：
(1)利用MODTRAN辐射传输模型计算天空光的光谱特性；
(2)选择金属卤化物光源作为400nm-2500nm范围内天空光光谱的主要模拟光源；
(3)分段分析金属卤化物光源光谱与天空光光谱的异同；
(4)使用大功率固体照明元件LED作为400nm-600nm范围内的光谱补偿光源；
(5)选择相应滤光片，对600nm-2500nm范围内的金属卤化物光源光谱进行修正。
2.根据权利要求1所述的一种基于人造光源的天空光光谱模拟方法，其特征在于：在步骤(1)中所述的“利用MODTRAN辐射传输模型计算天空光的光谱特性”，其含义说明如下：输入参数设置为中纬度夏季，即北纬45°，7月，/乡村气溶胶模型VIS＝23km/太阳高度角70°/天顶角45°/相对方位角45°，然后利用大气辐射传输方程计算大气中气体分子和气溶胶颗粒对光线的散射和吸收，从而获得天空光的光谱数据。
3.根据权利要求1所述的一种基于人造光源的天空光光谱模拟方法，其特征在于：在步骤(2)中所述的“选择金属卤化物光源作为400nm-2500nm范围内天空光光谱的主要模拟光源”，其含义说明如下：金属卤化物光源的光谱范围为400nm-2500nm。
4.根据权利要求1所述的一种基于人造光源的天空光光谱模拟方法，其特征在于：在步骤(3)中所述的“分段分析金属卤化物光源光谱与天空光光谱的异同”，其含义说明如下：分段分析主要是针对400nm-600nm范围内金属卤化物光源的相对光谱能量所占的百分比较低，600nm-2500nm范围内金属卤化物光源的相对光谱能量所占的百分比较高，根据两者的不同，以600nm为分割点，将400nm-2500nm分为400nm-600nm和600nm-2500nm两个波段范围进行分段分析。
5.根据权利要求1所述的一种基于人造光源的天空光光谱模拟方法，其特征在于：在步骤(4)中所述的“使用大功率固体照明元件LED作为400nm-600nm范围内的光谱补偿光源”，其含义说明如下：LED光谱补偿光源的中心波长使用曲线分峰拟合的方法进行计算，其计算过程如下：由于LED的光谱曲线近似认为符合高斯分布，因此分峰拟合的过程中，使用高斯函数作为LED的光谱分布曲线，假设v是光谱频率，则目标光谱Y(v)是由单峰谱线相互叠加所形成，分峰拟合的任务是求解一组单峰谱线Xi(v)(i＝1，2，…n)，使得下式成立：
Y(v)＝∑Xi(v)
但在实际情况下，上式两边不可能达到严格相等，在实际计算过程中，使目标光谱Y(v)与拟合光谱∑Xi(v)之间的误差在400nm-600nm范围内最小。
6.根据权利要求1所述的一种基于人造光源的天空光光谱模拟方法，其特征在于：在步骤(5)中所述的“选择相应滤光片，对600nm-2500nm范围内的金属卤化物光源光谱进行修正”，其含义说明如下：滤光片的选择分两步完成，首先计算经步骤(4)补偿后的金属卤化物光源相对光谱能量分布与MODTRAN计算的标准天空光相对光谱能量分布在600nm-2500nm范围内的比值，再根据计算所得的滤光片透过率曲线选择相应的滤光片进行光谱修正。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210483543.4A CN102997994B (zh) | 2012-11-23 | 2012-11-23 | 一种基于人造光源的天空光光谱模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210483543.4A CN102997994B (zh) | 2012-11-23 | 2012-11-23 | 一种基于人造光源的天空光光谱模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102997994A CN102997994A (zh) | 2013-03-27 |
| CN102997994B true CN102997994B (zh) | 2014-11-26 |
Family
ID=47926855
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210483543.4A Active CN102997994B (zh) | 2012-11-23 | 2012-11-23 | 一种基于人造光源的天空光光谱模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102997994B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104735848B (zh) * | 2013-12-20 | 2017-04-12 | 致茂电子股份有限公司 | 发光模块及其驱动方法 |
| CN104748853A (zh) * | 2015-04-15 | 2015-07-01 | 常州市武进区半导体照明应用技术研究院 | 日光光谱测量方法和装置 |
| CN105158811B (zh) * | 2015-09-24 | 2016-07-06 | 河北省科学院地理科学研究所 | 模拟现实场景的地物波谱采集装置和采集方法 |
| CN106369486B (zh) * | 2016-09-28 | 2018-11-16 | 河北工业大学 | 合成全波段匹配太阳模拟器光谱的方法 |
| CN108279221B (zh) * | 2017-12-07 | 2021-04-13 | 中国科学院国家天文台 | 一种局部天区大气透明度获取方法 |
| US20220112993A1 (en) * | 2018-08-17 | 2022-04-14 | Avalon St Sàrl | Solar simulator |
| CN110686776B (zh) * | 2019-10-09 | 2021-05-04 | 北京航空航天大学 | 一种室内高光谱遥感成像仿真装置与方法 |
| CN112305655B (zh) * | 2020-10-29 | 2023-02-17 | 沈阳仪表科学研究院有限公司 | 光伏用b级太阳模拟滤光片及其镀制方法 |
| CN117452541B (zh) * | 2023-11-16 | 2024-07-09 | 上海泊睿科学仪器有限公司 | 光伏用a+级太阳模拟滤光片及其镀制方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007183710A (ja) * | 2006-01-04 | 2007-07-19 | Hitachi Software Eng Co Ltd | 緑地の彩度補正システム |
| CN101598543A (zh) * | 2009-07-29 | 2009-12-09 | 中国科学院对地观测与数字地球科学中心 | 一种实用的遥感影像大气校正方法 |
| CN101694466A (zh) * | 2009-10-19 | 2010-04-14 | 浙江大学 | 背照式ccd阵列光谱仪标准具效应的校正方法 |
| CN102435318A (zh) * | 2011-09-16 | 2012-05-02 | 哈尔滨工业大学 | 考虑天空光对地表反射影响的偏振遥感数据仿真方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2547359C (en) * | 2003-11-26 | 2012-11-27 | Florida Environmental Research Institute, Inc. | Spectral imaging system |
-
2012
- 2012-11-23 CN CN201210483543.4A patent/CN102997994B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007183710A (ja) * | 2006-01-04 | 2007-07-19 | Hitachi Software Eng Co Ltd | 緑地の彩度補正システム |
| CN101598543A (zh) * | 2009-07-29 | 2009-12-09 | 中国科学院对地观测与数字地球科学中心 | 一种实用的遥感影像大气校正方法 |
| CN101694466A (zh) * | 2009-10-19 | 2010-04-14 | 浙江大学 | 背照式ccd阵列光谱仪标准具效应的校正方法 |
| CN102435318A (zh) * | 2011-09-16 | 2012-05-02 | 哈尔滨工业大学 | 考虑天空光对地表反射影响的偏振遥感数据仿真方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102997994A (zh) | 2013-03-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102997994B (zh) | 一种基于人造光源的天空光光谱模拟方法 | |
| Norton et al. | Comparison of solar spectral irradiance measurements using the average photon energy parameter | |
| CN109253976B (zh) | 基于光感模块的高光谱实时辐射定标方法 | |
| CN104251824B (zh) | 一种多光谱作物生长传感器温度补偿模型的构建方法 | |
| Igawa | Improving the All Sky Model for the luminance and radiance distributions of the sky | |
| CN104902153B (zh) | 一种多光谱相机色彩校正方法 | |
| CN111077093A (zh) | 一种基于多光谱技术的煤矸石快速检测的方法与装置 | |
| CN110866467A (zh) | 一种航空中红外高光谱数据温度和发射率反演方法 | |
| CN113837426A (zh) | 一种基于天气分型的光伏功率预测方法 | |
| Theristis et al. | Performance of a concentrating photovoltaic monomodule under real operating conditions: Part I–Outdoor characterisation | |
| CN115266610B (zh) | 一种植物生长监测方法及系统 | |
| CN102322957B (zh) | 一种干涉型高光谱成像仪的光谱漂移检测方法 | |
| Wittmer et al. | Analysis of PV grid installations performance, comparing measured data to simulation results to identify problems in operation and monitoring | |
| Yu et al. | Investigation of aerosol optical depth (AOD) and Ångström exponent over the desert region of northwestern China based on measurements from the China Aerosol Remote Sensing Network (CARSNET) | |
| CN105007676A (zh) | 基于led混光呈色模型的光谱功率分布提取方法及系统 | |
| CN104378861B (zh) | 一种色温可调黑体光源及其色温调节方法 | |
| Zhao et al. | Data mining techniques for novel local reference solar irradiance spectra calculations: A case study in Beijing | |
| CN105960046A (zh) | 基于led大平面光源的空间任意点照度匹配方法 | |
| Sevillano-Bendezú et al. | Average photon energy assessment based on modelled spectra from the National Solar Radiation Database for Lima, Peru | |
| CN103678793B (zh) | 一种对色散聚光型太阳能发电系统的仿真建模方法 | |
| CN113420253B (zh) | 一种多波段合成的黑碳气溶胶直接辐射效应计算方法 | |
| Tourasse et al. | Development of a system measuring the solar radiation spectrum in 5 planes for daylight and PV applications | |
| CN108872093A (zh) | 基于被动遥感的o4吸收校正及气溶胶消光廓线反演方法 | |
| CN109299431B (zh) | 基于希尔伯特边际谱特征的光伏组件性能一致性评价方法 | |
| CN110765641B (zh) | 一种定量计算光学观测光污染的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20130327 Assignee: Beijing Weishi Honghe Space Information Technology Co.,Ltd. Assignor: BEIHANG University Contract record no.: X2022990000751 Denomination of invention: A Simulation Method of Skylight Spectrum Based on Artificial Light Source Granted publication date: 20141126 License type: Common License Record date: 20221009 |