基于傅立叶变换的光谱数据分析

1绪论

本章介绍课题的研究背景，总结阐述光谱分析技术的发展应用，以及光谱测量仪

器的分类和各自特点，特别是傅里叶光谱仪及应用情况，简要介绍傅里叶变换光谱仪

的研究现状及成果;最后阐述本课题的研究目的、意义以及主要研究内容和技术指标

要求。

1.1选题的背景、目的和意义

在现代高技术战争中，激光武器及其对抗已显得日益重要，面对战场上激光战术

侦察、激光武器和激光制导武器等激光威胁，加速发展激光侦察告警技术己成为激光

对抗的首要任务。准确、可靠、迅速地掌握对方激光的属性己成为交战双方开战的重要前提，因此采用先进技术提高激光告警设备敌我识别的性能、抗干扰能力和反应速度是非常必要的。

激光告警技术是是光电对抗的重要组成部分。研究激光告警技术的目的是快速探测敌方激光威胁的存在，尽可能确定出其方位、波长、强度、脉冲特性(脉宽、重复频率等)等信息，以便我方能及时采取保护或反击措施。激光告警设备硬件通常由激光接收系统、光电传感器、信号处理器、显示与告警装置等部分组成。目前，告警设备在软件上基本都采用解方程组或者查表的方法，求解来袭激光的波长、角度和次数等基本信息。例如，典型相干识别法的迈克尔逊型、法布里一拍罗(F-P)型和光栅衍射型告警机，利用形成的干涉条纹间距确定入射激光的波长，利用干涉图的横向位移量确定入射激光方向等。当激光以一定波长和方向入射时，特定条纹在光电探测器上的位置的不同或者条纹阳间距的不同，制作波长和与入射方向对应的查找表，这样处理器只需计算目标条纹的成像位置和间距，便可通过软件查表实现波长和角度的测定。这种方法原理简单、编程容易;但是无法求出目标激光的光谱特征，从而无法得到威胁激光的时、空特性和类型[1]。

为了实时获取来袭激光的光谱分布和类型，提高告警系统的信噪比和探测率，需

要研究具有高速、准确、性能可靠的新型激光告警系统。

1.2激光光谱探测技术的国内外研究现状

目前，激光信号光谱的探测，主要通过光谱仪来实现，光谱仪从原理上可分为色

散型和干涉型两大类。

色散型光谱仪以棱镜或光栅为色散器件，利用色散元件将复色光色散成一系列谱线，然后再用探测器测量每一谱线的强度，具有结构简单、性能可靠等优点，但因该类光谱仪均含有入射狭缝，狭缝越窄，光谱分辨率越高，进入系统的光通量就变小，信噪比降低，即光谱分辨率和光通量是色散型光谱仪中相互制约的一对矛盾。由于这类光谱仪发展较早，技术相对成熟，在各领域已有诸多应用[2,3]。

干涉型光谱仪，是通过对入射光束引入不同光程差进而形成干涉图，并对干涉图进行处理，获得光谱信息。典型的干涉型光谱仪有傅里叶变换光谱仪，通过该光谱仪形成的干涉条纹和激光光谱间存在傅里叶变换关系，因此对条纹进行傅里叶变换，便可以得到光源光谱信息，具有高光谱分辨率和高光通量等优点，能探测微弱激光信号，在激光被动侦察中有重要意义。理论分析表明，在相同情况下，傅里叶变换光谱仪的能量通过率比典型色散型光谱仪高300倍左右，而光谱分辨率一般也高两个数量级以上。从实现的方式上，傅里叶变换光谱仪，又可分为时间调制干涉光谱仪和空间调制干涉光谱仪两大类[4,5]。

时间调制干涉光谱仪，是以迈克尔逊干涉仪为基础的，光的调制靠镜面的机械扫描运动来实现，这就决定了这种调制方式的扫描速度不可能很高，同时对机械扫描系统的加工、装配等的精度提出了较高的要求。在实际应用中，时间调制干涉光谱仪暴露出两大缺点:1)动镜驱动系统对倾斜、晃动非常敏感;2)实时性不好，干涉图的形成需要动镜运动一个周期，故不适合快速变化激光的光谱测量。它的优点是能达到较高的光谱分辨率[6]。

空间调制千涉光谱仪也称静态干涉仪，是九十年代以来随着探测器阵列的飞速发展而出现的，与时间调制型干涉成像光谱仪相比，空间调制型干涉光谱仪是在同一时刻探测器阵列的不同位置(空间)上，获得目标辐射光谱。这种静态的干涉光谱谱仪的类型很多，典型的有基于振幅分割的静态Michelson干涉仪（如图1.1）、三角共路Sagnac 干涉仪（如图1. 2)、Mach-Zehnder干涉仪(如图1. 3)、基于偏振干涉的双折射晶体的偏振干涉仪(如图1. 4 )和光楔光谱仪(如图1.5)

图1.1静态Michelson干涉仪

图1.2三角共路Sagnac干涉仪

图1.3 Mach-Zehnder干涉仪

图1.4双折射晶体偏振干涉仪

图1.5光楔光谱仪

这些静态干涉仪产生的干涉图，是按空间分布的强度信号，因此可实现对时变或脉冲辐射的实时监测：它们没有运动部件和扫描机构，使得仪器的结构变得简单、紧凑和稳定，抗干扰和震动能力强，减轻了设计和加工难度、体积和重量小，降低了成本。

鉴于以上诸多优点，空间调制干涉光谱技术一出现，即引起了美国、加拿大、日本、芬兰等国家的高度重视。目前国际上主要有:美国佛罗里达工学院和夏威夷大学基于Sagnac 干涉具的空间调制干涉光谱仪，其光谱范围为m m μμ0.5~0.1，分辨率为1100-cm ;美国华盛顿大学基于双折射元件的数字阵列扫描干涉仪，仪器的光谱范围为m m μμ2.2~4.0，分辨率为1300-cm ，空间分辨率为1. 5'，视场为︒5，体积只有40cm X 10cm X 15cm;美国宾夕法尼亚大学的研究人员研制的基于握拉斯顿棱镜阵列的干涉装置，可改善普通握拉斯顿棱镜光谱仪在分辨率上的不足，其结构紧凑，分辨率高;另外，美国Hughes Santa Barbara 和NASA Langley 研究中心分别研制成功光楔滤光片光谱仪WIS21 WIS22和基于地球静止卫星F-P 干涉光谱仪。英国圣安德鲁斯大学研制了基于Wollaston 棱镜的空间调制光谱仪;日本大阪大学基于Savat plate 的多通道红外空间调制光谱仪[7]。

在国内，中国科学院上海技术物理研究所、西安光机所和长春光机所及北京空间机电研究所等单位也开展了干涉光谱仪的研究工作。上海技术物理研究所主要研究光栅推扫式光谱仪和傅里叶变换光谱仪，该光谱仪采用时间调制形式，利用直线电机带动角镜

作一维直线运动，光谱范围m m μμ153-，分辨率11-cm ;西安光机所研制的大孔径静态空间调制干涉光谱仪采用Sagnac 结构，他们吸收了1996年美国科学应用国际公司((SAIC)的专利，光谱范围m m μμ75.05.0-;沈为民等人研制了一种紧凑型空间调制傅里叶变换光谱仪，它是由两块半五角形棱镜构成的环形共光路，具有无机械扫描、结构简单、体积小、性能稳定的特点。长春光机所应用光学国家重点实验室正在研制的是猫眼动镜驱动式傅里叶变换光谱仪的原理样机，该样机主要是用来测量紫外一真空紫外的光谱，光谱范围m m μμ3.017.0-，分辨率1025.0-cm 。浙江大学光电系在20世纪90年代初研制过弹簧片扫描结构的傅里叶变换光谱仪，但后来再没有相关报道。

各种静态干涉光谱仪性能比较如表1. 1所示。从比较结果可见，光楔型光谱仪的整体性能较好，工艺要求相对较低。因此，本项目中采用等效光楔(斜楔)干涉具做干涉 元件，获得干涉条纹[8]。

表1.1各类静态干涉仪光谱仪性能比较

1.3研究主要内容

为了实时获取来袭激光的光谱分布并识别其类型，提高激光告警系统的信噪比和探测率，本课题研究基于相干探测原理，以光纤束为光学天线、用等效斜楔干涉具做相干识别器的高速、高性能被动激光光谱探测技术。特定入射方向的入射光被光纤束光学天线接收，经准直后垂直照射干涉具进行相干叠加，聚焦透镜(柱面镜)将从干涉具出射的平行光束聚焦在光电传感器线阵上，然后由高速信号处理系统求得光源光谱信息[9]。

2光谱数据采集

2.1傅里叶变换光谱测量原理

Fourier 变换光谱仪是通过对双光束干涉仪产生的干涉图进行Fourier 变换来获取光谱的。在Fourier 变换光谱学的发展过程中，人们提出过多种形式的Fourier 变换光谱仪，尽管它们形式各异，但其物理原理和基本理论却是一致的。图2. 1是最简单而又能完全说明Fourier 变换光谱学基本原理的Michelson 干涉仪。下面通过对Michelson 干涉仪产生的干涉图的定量分析来阐述Fourier 变换光谱学的基本原理[10,11]。

图２.１ 迈克尔逊光谱仪原理

设有一振幅为A ，波数为v 的单色平行光以平行于静镜M1入射到反射系数为r ，透射系数为t 的理想分束器BS 。分束器将入射光分成振幅为ra 的反射光束和振幅为ta 的透射光束，分别经静镜M1和动镜M2反射后再次回到分束器，又分别被分束器分为透射和反射的两束。其中部分向着光源的方向传输，另部分则沿与入射方向垂直的方向传输，在探测器D 上进行干涉，形成干涉条纹。探测器接收到的信号为:

φi D e A t tA r A -⋅⋅+⋅= 式２.1

式φ是从动镜和静镜反射的两束相干光之间的位相差:

vx πφ2= 式2.2

式中x 为干涉仪两臂之间的光程差。探测器上信号的强度为

)2cos()(2)(2)

2cos(22),(002222220vx v RTB v RTB vx A t r A t r I I A A v x I a d D D ππ+=+=+=⋅= 式2.3

其中22,t T r R ==表示分束器的反射率和透射率，)(0v B 是入射光束强度分布。可见，探测器产生的信号是由直流分量)(20v RTB 和交流分量组成)2cos()(20vx v RTB π由于只有交流分量反映输入光谱的形状，故去掉直流分量得到:

)2cos()(2),(0vx v RTB v x I a π= 式2.4

此交流分量代表了一定波数v 的入射光束，经Michelson 干涉仪后，在探测器D 上的强度I 。随光程差x 的变化规律。式2.4是在理想状态下，单色光的干涉图。实际上探测器得到的干涉条纹强度除与光源辐射强度Bo (v)有关外，还与分束板性质、探测器光谱响应，前置光学系统、电子线路等环节有关。因此对于实际情况，需要对式2.4进行修正，引入修正因子H(v) (H(v)

)2cos()()(2),(00vx v B v RTH v x I π= 式2.5 令)()(2)(0v B v RTH v B =表示修正好的光源强度，则：

⎰∞∞-=)2cos()(),(vx v B v x I z π 式2.6 这里B(v)表示光源的功率谱分布。不难看出，干涉图I(x)与光谱分布)(v B 是又是 Fourier 变换对，即

⎰∞∞

-=dv vx v B x I )2cos()()(π 式2.7 ⎰∞

∞-=dx

vx x I v B )2cos()()(π 式2.8 以上两式是Fourier 变换光谱学的基本关系式，也是傅里叶激光光谱探测技术的理论基础。以上是激光光谱探测技术的基础，实际应用中为了获得激光的功率谱分部，还必须附之以相应功能软件来完成频谱的复原和波长计算[12,13]。

2.2激光光谱探测系统方案

为了完成前面提出的研究任务，本文提出了基于等效斜楔干涉具的静态傅里叶变换光谱探测方案，在该系统中采用了光纤接收天线收集来袭激光并将其转换成方向-致的平行光:信号处理系统以具有较强信号处理功能的数字信号处理芯片为核心。以傅里叶变换光谱技术原理为基础，设计了适当的光路和干涉具，使入

射激光产生干涉，对干涉图进行傅里叶变换，从而得得入射激光的波长和光普信息[14]。

2.2.1整体结构

傅里叶变换激光光谱探测系统原理如图2.2所示，系统由光学接收天线、等效斜楔干涉具、聚焦透镜、线阵CCD 相机和信号处理系统组成。线阵CCD 置于聚焦透镜的焦平面上。光学天线用来收集激光并将其转换为平行光，经准直后照射到干涉具表面，在干涉具输出面处产生周期性的空间干涉条纹，经聚焦透镜，成像在线阵CCD 上，通过CCD 相机进行光电转换和模数转换后，运用TI 公司的TMS320VC5509定点型DSP 进行处理，得到入射激光的光谱信息和波长[15]。

图２.２傅里叶变换激光光谱探测系统

２.2.2可行性分析

图2. 3是该系统中等效斜楔干涉具原理示意图，干涉其由两个三角棱镜粘合成，棱镜I 是等腰直角棱镜，M1是其光入射面，M2是其全反射而;棱镜2是非直角棱镜,M3是其个反射面，Ml 是其光出射面，它的一个锐角是45O ，另一个锐角为α-O 45谱仪数据计算量较大，而系统对计算速度要求又较高，因此在数据采集与处理中使用高性能的DSP 处理器，以满足系统对速度的要求。整个系统集成化程度高、响应速度快、抗干扰能力强。下面对系统整体结构及可行性进行阐述。

(α为一小角);两个三角棱镜通过长边粘合在一起，粘合处镀有半透半反膜，形成一个分光面BS。光学天线收集的入射光，经准直后变成平行光，该平行光透过M1照射到BS上，经反射、透射分成光强近似相等的两束，再分别由M2和M3反射后，经过BS透射、反射后，从M4射出干涉具。由于角度α的原因，使得两束激光在干涉具中传播时产生了连续的光程差，从干涉具出射后，出射激光产生双光束干涉，形成干涉条纹[16]。

图２.３等效斜楔干涉具原理

本干涉具中，干涉条纹与其光谱之间满足式2.8，仍是一对傅里叶变换:

⎰∞∞-=dv

(

)

(π

cos(

)

vx

v

B

x

I)

2

式2.9

⎰∞∞-=dx

cos(

)

(

2

(π

)

vx

B)

x

I

v

式2.10 同样对干涉条纹进行傅里叶变换就可以得到光源的光谱信息。

干涉具产生的干涉条纹是空间分布的周期信号，使用电荷祸合器件(CCD)，将这种空间周期信号变换为时域周期信号，再将时域周期信号经过A/D转换后，运用离散傅里叶变换(DFT)，可求出入射激光的光谱信息和波长。将得到的数据存为1024.txt。该干涉具，具有非扫描、体积小、结构简单、成本低、光谱检测速度快等优点，可满足高速探测窄脉冲激光信号光谱的要求。

２.３等效斜楔干涉具的设计

由图２.４中(a)中可知，干涉具中角度α的存在是其产生干涉条纹的关键，它可等效为一个斜楔，而又不是一个真正的斜楔，因此我们称其为等效斜楔干涉

具。在该傅里叶变换光谱探测技术中，干涉具是核心部件，系统的光谱分辨率、波长计算误差大小等都与干涉具有直接的关系。

２.３. 1干涉具中的光程差分析

如图２.４中(a)所示，从干涉具出射的两束光传播方向不同，所以发生干涉的两条光线并不是由同一条入射光线分束形成，而是分别由两条入射光线经BS分束形成。

a

b

图２.４等效斜槲干涉具原理

光线1入射到BS 的A 点上，经过BS 反射到达M2上B 点，经M2反射后透过BS 到达柱面镜上C 点;光线2入射到BS 的D 点上，透过BS 照射到M3的E 点上，由于α角的原因，光线不能原路返回，而是会相对于原路线偏转大小为α2的角度，照射到BS 的F 点上，再经BS 反射到达C 点，与光线I 相交，发生干涉。做FG 土AC ，交点为G 点，AC 交DE 于H, EF 交于I ，延长DE 交M2’于J, M2’是M2相对于分光面BS 的镜像面，与M3相交于0点。由光线反射理论及相关三角关系可得

α2=∠=∠=∠EFG DCA DEF

45=∠=∠FAG AFG

α2cos FI FG AG == 式2.11 GC AG AC += 式2.12 α2cos EI HE = 式2.13 α2cos FC GC = 式2.14 EJ HE HJ AB +== 式2.15

从图中可以看出在A 点、H 点 之前两束光传播路经完全相同，光程差产生于此后的传播中。光线1在A 点以后传播的路径总长度为

AC AB I +=21 式2.16 光线2在F 点以后传播的路径总长度为

FC

IF EI HE FC EF HE I +++=++=2 式2.17

将式2.9 2.10 2.12 带入式2.14可得

α2cos )(21FC IF AB I ++= 式2.18

将式2.16减去式2.15，再结合式2.9、2.11、2.13，可以得到在J 点发生干涉的两条光线之间的路径差

)

12)(cos (2)12)(cos ()12(cos 2-++=-++-+=∆αααFC EF EJ FC IF EI EJ l 式2.19

一般α取一个很微小的值，否则干涉条纹太密，无法用CCD 采集。此时，12cos ≈α，因此EJ l 2≈∆，EJ 的值与J 点相对于0点的距离OJ 及α角的大小有关，此时干涉可等效于一个斜楔，如图2.4中(b)所示，并以0为原点建立如图直角坐标系。

αOJtg EJ =

光程差为光线路径差与干涉具材料折射率的乘积，可以求得由等效斜楔干涉具出射两束光之间光程差与角α和'2M 上坐标点之间的关系为

αtg nx l n x '2=∆= 式2.20

其中n 是干涉具材料的折射率，x 表示光程差[17]。 ２.３.2干涉条纹的傅里叶变换

当光程差为入射激光波长的整数倍时，即:

αλ

ctg n

x 2'=

∆

如图2.4中(a)所示，柱面镜相当于M2相对于BS 的镜像面，故这两个面之间的点存在一一对应关系。M2上光程差为几的两点间距等于柱面镜上干涉条纹间距d 。

αλ

ctg n

d 2=

式2.21

等效斜楔干涉具出射的两束激光光强为021I I I ==，则探测面上的光强分布为: 其中:

λ

1

=

f 式2.22

表示入射激光的波数，x 表示光程差。由傅里叶变换关系可以求得入射激光的波数光谱分布:

⎰∞

∞

--=dx e x I f E fx π2)()( 式2.23

上式说明只要对线阵CCD 采集到的由干涉条纹进行傅里叶变换就可以得到入射激光的光谱信息[18]。 ２.３.３α角的选择

α角是该等效斜楔干涉具最重要的参数，以下研究α角选择的理论依据。由

傅里叶变换理论可知傅里叶光谱仪的光谱分辨率是最大有效光程差的倒数。在设计等效斜楔干涉具时，要提高光谱分辨率，就要尽可能的增加最大有效光程差。由式3.19可知增加光程差的途径有两个，增大α角和增加'x 的值[20]。

对干涉具而言，'x 的最大值应该是干涉具的边长，但是只有照射到线阵CCD 光敏面上的光才有效，所以'x 的最大值应该是CCD 光敏面的有效长度。所以在选

定线阵CCD 的后，只能通过增加α角的值来增大最大有效光程差，但是由式3.20可看出，增大α角会导致干涉条纹间距的减小。由于线阵CCD 像元尺寸的限制，条纹间距不可能无限减小。根据香农抽样定理，在一个周期至少要采样两次，才可能准确地恢复信号，即最少需要两个像元才能检测一条干涉条纹，干涉条纹间距要大于等于两倍像元的宽度[21]。

研究中我们选用ATMEL 公司的CamraLink 相机采集干涉条纹，该相机共有2048个像元,像元尺寸为14 x 14m μ，有效探测面总长度为

cm m L 6872.21414=⨯=μ 式2.24

使用时只能使用奇数像元或者偶数像元，故最小条纹间距少大于等于4倍相元宽度，即m d μ56m in =等效斜楔干涉具材料为BK7光学玻璃，其折射率为1. 4581。代入式3. 20 可得α角最大为

min

min 2nd arctg

λ

α= 式3. 25

在最小干涉条纹间距和干涉据材料确定后，确定α的大小还与入射激光的波长有关，由式2.20可以看出，激光波长与条纹间距成正比，波长越小干涉条纹间距越小。

３干涉条纹的高速采集与处理

光电转换是傅里叶光谱测量的重要环节，数字信号处理方法是得到光谱信息得必需手段。由于激光告警接收机的特殊条件，必须采集随时到来的激光信号，并且后续处理系统也要在处理当前数据的同时接收下一帧数据，这就要求在系统设计时要考虑到高速信号的处理能力。由于DSP具有高速的运算能力，使系统实时高速成为可能。本章主要围绕干涉条纹的高速采集与处理，在理论结合试验的基础上，详细介绍干涉具条纹的采集存储、光谱快速获取等方法[22]。

3. 1光电探测器

空间调制型静态傅里叶变换光谱仪的光电信号转换通常由CCD阵列探测器完成。CCD探测器由于集成度高，光谱响应范围宽，光谱分辨率高而得到广泛的应用，是现代光电子学与测试技术中应用最广的光电传感器。

3. 1. 1 CCD相机的选取

1)阵列类型的选择

CCD可分为一维(线阵)的和二维(面阵)两种，其功能是把光学图像信号转变为电信号输出。由于本装置中干涉条纹是一系列分布方向固定的条纹，在垂直条纹方向上的某一行数据进行运算即可获得足够信息，因此选用线阵CCD作探测器，这样既保证了足够的信息量，又有有效降低了数据量，提高处理运算速度。

2)性能要求

根据要求光谱探测范围为nm

400-，故应选取光谱响应范围包含此频

nm1100

段的Si线阵CCD相机作为探测器。另外，系统中相机探测的激光信号具有光谱范围宽、激光能量低、持续时间短的特点，因此CCD相机应具有光谱响应范围宽、动态范围大、噪声小、暗电流小、灵敏度高等特点。

3. 1. 2 CCD相机性能

本系统中选用的法国Atmel公司的AV i i VA M2 CL型线阵CCD相机(见图3.

1 a)，其主要特性及参数如下:

1.光谱响应范围宽，可达nm

nm1100

250-;

2.像元数为1024或2048个(可通过相机设置来选择);

3.分辨率8bit, lObit, 12bit可选;

4.光敏感面积为m μ1414*;

5.线阵长30mm;

6.积分时m m μμ135-;

7.数据读出率可达60Mpixel/s

相机采用 Cameralink 接口，是目前工业相机普遍采用的接口格式。采用低压差分信号(LVDS )形式传输数据，数据稳定性高、传输速度快。图3.1b 为AViiVA M2 CL 型线阵CCD 相机的光谱响应曲线[23]。

a

b

图3. 1 AViiVA h12 CL 型相机及其光谱响应曲线 3.2干涉条纹的采集

光电探测器(CCD )将采集到的光信号转化为电信号后，需要送入处理系统进行处理。因此如何快速、实时地与处理系统进行通信是装置的关键技术之一。

数据采集是指从CCD相机输出的串行LVDS电平信号到单片机接收模块之间的过程。数据的采集、传输和处理要协调进行，并且要满足以下条件:(1)数据的采集速度不能低于数据的输出速度;(2)在相机的积分时间内必须要将上一帧的数据处理并输出完毕。如果上述条件不满足，就会导致数据丢失，失去实时性的要求[24]。

4. MATLAB与VC++实验结果对比

4.1用MATLAB进行实验

图4.1读入采集数据得到的源图像

运行程序

plot(A);

axis([0 1024 0 4096]);

xlabel('探测单元数');

ylabel('灰度值');

fA=fft(A,1024);

fAA=abs(fA);

fAA(1:19)=0;

fAA(1005:1024)=0;

n=1:241;

n=n+28;

此段程序将得到的光谱数据读入，然后转化为以探测单元数为横坐标，以灰度值为纵坐标的时域光谱图。

图4.2运行福利叶变换程序后得到的频谱图

运行程序

figure;plot(n,fAA(10:250));

M=max(fAA);find(fAA>=M)

xlabel('波长 \lambda(nm)');

ylabel('强度值');

fB=fAA(10:250);

fC=[];

for i=1:241

j=242-i;

fC(j)=fB(i);

end

m=269:-1:29;

figure;plot(m,fC)

%M=max(fB);find(fB>=M)

xlabel('波长 \lambda(nm)');

ylabel('强度值');

axis([29 269 0 140000]);

此段程序是将上面的时域光谱图通过福利叶变换将其转化为图。以波长为横坐标，以强度为纵坐标。

4.2 用ＶＣ＋＋进行实验

运行下面的程序使窗口初始化

BOOL CReadWord1Dlg::OnInitDialog()

{

CDialog::OnInitDialog();

// Add "About..." menu item to system menu.

// IDM_ABOUTBOX must be in the system command range.

ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);

ASSERT(IDM_ABOUTBOX < 0xF000);

CMenu* pSysMenu = GetSystemMenu(FALSE);

if (pSysMenu != NULL)

{

CString strAboutMenu;

strAboutMenu.LoadString(IDS_ABOUTBOX);

if (!strAboutMenu.IsEmpty())

{

pSysMenu->AppendMenu(MF_SEPARATOR);

pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);

}

}

// Set the icon for this dialog. The framework does this automatically

// when the application's main window is not a dialog

SetIcon(m_hIcon, TRUE); // Set big icon

SetIcon(m_hIcon, FALSE); // Set small icon

// TODO: Add extra initialization here

return TRUE; // return TRUE unless you set the focus to a control

}

图4.３窗口初始化

运行下面的程序实现读入得到的光谱数据，光谱数据储存在１０２４.ｔｘｔ文件夹内。void CReadWord1Dlg::OnInput() //读数

{

// TODO: Add your control notification handler code here

CFileDialog dlg(TRUE, NULL,NULL,OFN_ALLOWMULTISELECT);

if (dlg.DoModal()!=IDOK)

{

EndDialog(IDCANCEL); //取消则退出程序

}

CString strLine, strEdit;

int i=0;

m_File.Open(dlg.GetPathName(),CFile::modeRead);

for(int m=0; m<=1023; m++)

{

m_File.ReadString( strLine );

m_Data = (double)atof((char *)(LPTSTR)(LPCTSTR)strLine);

m_flt[m] = m_Data;

m_Edit1.GetWindowText(strEdit);

m_Edit1.SetSel(strEdit.GetLength(), strEdit.GetLength());

strLine += "\r\n";

m_Edit1.ReplaceSel(strLine);

}