CN105606342B - 高轨卫星星地激光远场分布测试方法 - Google Patents

高轨卫星星地激光远场分布测试方法 Download PDF

Info

Publication number
CN105606342B
CN105606342B CN201511021877.XA CN201511021877A CN105606342B CN 105606342 B CN105606342 B CN 105606342B CN 201511021877 A CN201511021877 A CN 201511021877A CN 105606342 B CN105606342 B CN 105606342B
Authority
CN
China
Prior art keywords
light spot
computer
satellite
orbit satellite
spot information
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
Application number
CN201511021877.XA
Other languages
English (en)
Other versions
CN105606342A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201511021877.XA priority Critical patent/CN105606342B/zh
Publication of CN105606342A publication Critical patent/CN105606342A/zh
Application granted granted Critical
Publication of CN105606342B publication Critical patent/CN105606342B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01M11/02Testing optical properties
    • G01M11/0242Testing optical properties by measuring geometrical properties or aberrations

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Image Processing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

高轨卫星星地激光远场分布测试方法,涉及星地激光通信领域。完成了对孔径受限系统极远光场分布特性测试和微畸变光场的极远距离传输及分布特性测试。计算机控制转台使转台上的望远镜对高轨卫星进行预瞄准,然后望远镜上的相机实时记录高轨卫星星地激光远场分布的光斑信息,计算机实时读取相机中的光斑信息,并根据光斑信息计算过域像素数和光斑信息的灰度质心坐标,计算机采用螺旋扫描方定位高轨卫星的位置,然后计算机根据过阈像素数判断转台是否捕获到高轨卫星,若转台捕获到高轨卫星,则通过计算机控制转台移动,实现准确捕获,完成高轨卫星的闭环实时跟踪,即获得高轨卫星星地激光远场分布测试。它适用于星地激光分布测试。

Description

高轨卫星星地激光远场分布测试方法
技术领域
本发明涉及星地激光通信领域。
背景技术
空间光通信系统中,光束的远场传输特性是一个重要的问题。尤其当光束经过大气等介质信道时,受到了空间条件的约束,对其远场分布将会产生严重影响。对于星地激光通信等链路,光束将在自由空间(真空)与约束空间之间进行超长距离的传输,其远场分布特性的研究非常重要。对高轨卫星星地激光远场分布测试十分重要,目前尚没有十分有效的测试手段,本方法主要是解决星地激光通信光信号传输四万两千公里后的光场分布测试难题。
发明内容
本发明针对空间激光通信光场分布特点,目的是完成对孔径受限系统极远光场分布特性测试和微畸变光场的极远距离传输及分布特性测试,提出了高轨卫星星地激光远场分布测试方法。
该方法主要解决以下两方面问题:
1)孔径受限系统极远光场分布特性测试。实际的空间光通信中的发射系统都是孔径受限系统。经过孔径受限系统的发射光束传输到极远距离后的光场分布特性,是卫星光通信以及深空光通信中必须要深入掌握的数据。
2)微畸变光场的极远距离传输及分布特性测试。空间光通信系统中,光学系统和光源等的缺陷都会导致发射光场产生微小畸变。微畸变的光场在极远距离传输过程中的传输特性和分布特性,对空间光通信系统优化设计具有重要意义。
高轨卫星星地激光远场分布测试方法,该方法是基于高轨卫星星地激光远场分布测试装置实现的,该装置包括望远镜1、相机2、转台3和计算机4;
望远镜1安装在转台3上,且随转台3移动;相机2安装在望远镜1的物镜上;转台3的转动信号输入输出端连接计算机4的转动信号输出输入端;相机2的拍照信号输入输出端连接计算机4的拍照信号输出输入端;该装置包括望远镜1、相机2、转台3和计算机4;
该方法包括下述步骤:
步骤一、计算机4向转台3发送控制指令,使控制转台3移动到指定方位、指定俯仰角度,对高轨卫星进行提前预瞄准;
步骤二、计算机4持续读取相机2中的光斑信息,并根据光斑信息计算过阈像素数和光斑信息的灰度质心坐标;
步骤二中计算机4根据光斑信息计算过阈像素数的具体过程如下;
设定光斑信息的光斑像素指定阈值d,将光斑信息的光斑像素灰度值与光斑信息的光斑像素指定阈值d进行比较;
若光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d,则计数一次;
若光斑信息的光斑像素灰度值≤光斑信息的光斑像素指定阈值d,则不计数;
对光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d的计数数目进行求和统计,获得的总的计数数目,即为过域像素数;
步骤二中计算机4根据光斑信息计算光斑信息的灰度质心坐标的具体过程如下;
设定光斑信息的大小为M1×M2维图像,第m1行m2列像元坐标记为(xm1,ym2);M1×M2维图像由m1行m2列像元构成;m1和m2均为整数;
其中,I(xm1,ym2)表示像元的质量密度,即图像上x、y处像素点的灰度值;x、y值均为整数;
则光斑总能量Wm的表达式为:
光斑信息的质心坐标的表达式(xc,yc)为:
步骤三、计算机4读取转台3的当前方位和俯仰角度,并且计算机4控制转台3采用螺旋扫描方式扫描高轨卫星,并计算机4根据过阈像素数判断转台3是否捕获到高轨卫星,若是,则执行步骤四;若否,则返回步骤一;
步骤三中的计算机4控制转台3采用螺旋扫描方式扫描高轨卫星,即采用等线速的方式从初始瞄准点开始螺旋向外侧进行扫描,具体过程如下:
步骤三一、以初始瞄准点作为坐标原点(0,0)建立极坐标系;螺旋扫描曲线可表示为:
式中,r——极坐标半径变量;θ——极坐标角变量;Iθ——扫描步长;
步骤三二、从初始瞄准点(0,0)到任意点(r,θ)的扫描时间近似表达为:
式中,△t——每步停留时间,其大小取决于扫描模式,任意点(r,θ)即为高轨卫星的坐标;
步骤三中计算机4根据过阈像素数判断转台3是否捕获到高轨卫星,若是,则执行步骤四;若否,则返回步骤一;
若过阈像素数﹥4,表示转台3捕获到高轨卫星;
若过阈像素数≤3,表示转台3未捕获到高轨卫星;
步骤四、计算机4根据步骤二获得的光斑信息的灰度质心坐标持续反馈控制转台3移动,完成高轨卫星闭环跟踪,即获得高轨卫星星地激光远场分布测试。
本发明所述的方法及所使用的装置,完成高轨卫星的捕捉,实时积累实测数据,为空间光通信系统的设计提供依据。
有益效果:高轨卫星星地激光远场分布测试方法基于高轨卫星星地激光远场分布测试装置实现,该装置的特点是光端机并不需要信号发射单元,只需要信号接收单元即可。为使结构简便零巧,接收单元仅有信标组件,这将大大降低空间激光通信光端机系统的性能要求,且装置具有轻小型化、便携性,可研制多套装置在不同地点进行测试。
高轨卫星星地激光远场分布测试方法,首先计算机控制转台使转台上的望远镜对高轨卫星进行预瞄准,然后望远镜上的相机实时记录高轨卫星星地激光远场分布的光斑信息,计算机实时读取相机中的光斑信息,并根据光斑信息计算过域像素数和光斑信息的灰度质心坐标,计算机采用螺旋扫描方定位高轨卫星的位置,然后计算机根据过阈像素数判断转台是否捕获到高轨卫星,若转台捕获到高轨卫星,则通过计算机控制转台移动,实现准确捕获,完成高轨卫星的闭环实时跟踪。
本发明适用于星地激光分布测试。
附图说明
图1为高轨卫星星地激光远场分布测试方法的流程图;
图2为高轨卫星星地激光远场分布测试装置的结构示意图;
图3为螺旋扫描方式示意图;
图4为高轨卫星星地激光远场分布测试装置结合高轨星地激光链路系统进行激光远场分布的数据采集试验分析原理图。
具体实施方式
具体实施方式一、参照图1和图2具体说明本实施方式,本实施方式所述的高轨卫星星地激光远场分布测试方法,该方法是基于高轨卫星星地激光远场分布测试装置实现的,高轨卫星星地激光远场分布测试装置,该装置包括望远镜1、相机2、转台3和计算机4;
望远镜1安装在转台3上,且随转台3移动;相机2安装在望远镜1的物镜上;转台3的转动信号输入输出端连接计算机4的转动信号输出输入端;相机2的拍照信号输入输出端连接计算机4的拍照信号输出输入端;
该方法包括下述步骤:
步骤一、计算机4向转台3发送控制指令,使控制转台3移动到指定方位、指定俯仰角度,对高轨卫星进行提前预瞄准;
步骤二、计算机4持续读取相机2中的光斑信息,并根据光斑信息计算过阈像素数和光斑信息的灰度质心坐标;
步骤三、计算机4读取转台3的当前方位和俯仰角度,并且计算机4控制转台3采用螺旋扫描方式捕获高轨卫星;同时,计算机4根据过阈像素数判断转台3是否捕获到高轨卫星,若是,则执行步骤四;若否,则返回步骤一;
步骤四、计算机4根据步骤二获得的光斑信息的灰度质心坐标持续反馈控制转台3移动,完成高轨卫星闭环跟踪,即获得高轨卫星星地激光远场分布测试。
本实施方式中,地面设置多个高轨卫星星地激光远场分布测试装置,并按照设定的位置排列,获取高轨卫星的具体位置。高轨卫星在完成星地激光传输的过程中,受到大气条件、气候因素等多个因素的限制,使得高轨卫星的跟踪非常困难,不利于对高轨卫星星地激光远场分布的测试。本发明中的计算机4控制相机2记录望远镜1拍摄图片,计算机控制转台3的移动,转台3转动的角度,仰俯角度等。
计算机控制望远镜1转动,使望远镜1的位置为观察星地激光信号的最佳位置,即预瞄准。然后望远镜1及相机2获得高轨卫星星地激光的远场分布图,计算机4控制相机2每隔0.001s拍摄一次远场分布图,然后将拍摄的图片实时发送至计算机4,然后计算机4进行分析处理。
工作过程:计算机控制转台使转台上的望远镜对高轨卫星进行预瞄准,然后望远镜上的相机实时记录高轨卫星星地激光远场分布的光斑信息,计算机实时读取相机中的光斑信息,并根据光斑信息计算过域像素数和光斑信息的灰度质心坐标,计算机采用螺旋扫描方定位高轨卫星的位置,然后计算机根据过阈像素数判断转台是否捕获到高轨卫星,若转台捕获到高轨卫星,则通过计算机控制转台移动,实现准确捕获,完成高轨卫星的闭环实时跟踪,即获得高轨卫星星地激光远场分布情况,实现了获得高轨卫星星地激光远场分布测试。
后续工作人员根据高轨卫星星地激光远场分布测试方法获得的高轨卫星的具体数据,进行后续分析及处理。
相机2由光电雪崩二极管(APD)及配套电路即可将功能由搜索信号源扩展为光通信接收。或者,转台3携带两个望远镜1,其中一个望远镜与相机2机械硬连接,另一个望远镜与光电雪崩二极管(APD)及配套电路机械硬连接,该装置能够允许瞄准、捕获、跟踪以及光通信接收同时实现。
具体实施方式二、本实施方式是对实施方式一所述的高轨卫星星地激光远场分布测试方法的进一步说明,本实施方式中,步骤二中计算机4根据光斑信息计算过阈像素数和光斑信息的灰度质心坐标,根据光斑信息计算过阈像素数的具体过程如下;
设定光斑信息的光斑像素指定阈值d,将光斑信息的光斑像素灰度值与光斑信息的光斑像素指定阈值d进行比较;
若光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d,则计数一次;
若光斑信息的光斑像素灰度值≤光斑信息的光斑像素指定阈值d,则不计数;
对光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d的计数数目进行求和统计,获得的总的计数数目,即为过域像素数。
本实施方式中,计算机根据相机记录的光斑信息,获取光斑信息的光斑像素灰度值,并在计算机中设定光斑信息的光斑像素指定阈值d,若计算机中获得的光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d,则计算机计数一次,所有光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d的总的光斑数目即为过域像素数。
具体实施方式三、本实施方式是对实施方式一所述的高轨卫星星地激光远场分布测试方法的进一步说明,本实施方式中,步骤二中计算机4根据光斑信息计算过阈像素数和光斑信息的灰度质心坐标,根据光斑信息计算光斑信息的灰度质心坐标的具体过程如下;
设定光斑信息的大小为M1×M2维图像,第m1行m2列像元坐标记为(xm1,ym2);M1×M2维图像由m1行m2列像元构成;m1和m2均为整数;
I(xm1,ym2)表示像元的质量密度,即图像上x、y处像素点的灰度值;x、y值均为整数;
则光斑总能量Wm的表达式为:
光斑信息的质心坐标的表达式(xc,yc)为:
具体实施方式四、本实施方式是对实施方式一所述的高轨卫星星地激光远场分布测试方法的进一步说明,本实施方式中,步骤三中的计算机4控制转台3采用螺旋扫描方式扫描高轨卫星,即采用等线速的方式从初始瞄准点开始螺旋向外侧进行扫描,具体过程如下:
步骤三一、以初始瞄准点作为坐标原点(0,0)建立极坐标系;螺旋扫描曲线可表示为:
式中,r——极坐标半径变量;θ——极坐标角变量;Iθ——扫描步长;
步骤三二、从初始瞄准点(0,0)到任意点(r,θ)的扫描时间近似表达为:
式中,△t——每步停留时间,其大小取决于扫描模式,任意点(r,θ)即为高轨卫星的坐标。
本实施方式中,扫描方式有很多,其中最常见的有矩形扫描、螺旋扫描以及矩形螺旋扫描。其中螺旋扫描以有效覆盖和工程易实现等优点而被广泛采用。本发明选用螺旋扫描,采用等线速的方式从瞄准初始位置开始螺旋向外侧进行扫描,螺旋扫描在捕获视场中的扫描轨迹如图4所示。图中虚线代表设定扫描的不确定域的边界,实线代表螺旋扫描轨迹,其中的圆点代表扫描过程中每步停留的地方。图4中5表示激光通信地面终端;6表示星地激光通信链路下行接收光场。
对于单场扫描捕获,每步停留时间△tI取决于发射端控制系统带宽FT,即
具体实施方式五、本实施方式是对实施方式一所述的高轨卫星星地激光远场分布测试方法的进一步说明,本实施方式中,步骤三中计算机4根据过阈像素数判断转台3是否捕获到高轨卫星,若是,则执行步骤四;若否,则返回步骤一;
若过阈像素数﹥4,表示转台3捕获到高轨卫星;
若过阈像素数≤3,表示转台3未捕获到高轨卫星。
本发明提出了高轨卫星星地激光远场分布测试方法及装置。与现有空间激光通信光端机相比,该装置的设计思路简洁、模块单元成熟,价格低廉、性能稳定、功耗低、质量轻、体积小、方便、灵活、隐蔽性强、可机载、舰载和车载。
该高轨卫星星地激光远场分布测试方法中的装置可以实时记录瞄准、捕获、跟踪过程的高轨卫星的相关参数,包括时间、日期、经度、纬度、方位、俯仰、赤经、赤纬、过域像素数、灰度质心坐标等。
具体实施方式六、本实施方式为一个实施例:
(1)望远镜
采用INTANE 80mm APO型号望远镜,宽波段复消色差设计,其主要参数特点为,通光口径为φ80mm;焦距为560mm。
(2)相机
采用AVT Marlin F131B NIR型号相机来接收800nm通信激光(对于1550nm等其他波段通信激光,采用对应光谱范围的其他型号相机即可),其主要参数特点如表1:
表1参数
光谱范围 400nm-1000nm
接口 IEEE 1394a
分辨率 1280×1024
像元尺寸 6.7μm
位深 8位
(3)转台
采用iOptron MiniTower Pro型号转台用于系统提前量补偿瞄准、捕获和跟踪。其主要特点如下,USB接口,内置32通道GPS接收机,俯仰、方位全周手动/电动微调,计算机控制自动寻星、自动跟踪经纬仪。其他硬件包括:800nm窄带滤光片、计算机、1394转接卡及数据线、USB连接线、USB转串口驱动。
(4)扫描参数设置
扫瞄范围3mrad(全角),扫描步长Iθ=1mrad,扫瞄每步停留时间△t=50ms。
(5)空间光场分布测试
利用高轨卫星星地激光远场分布测试装置,结合高轨星地激光链路系统可进行激光远场分布的数据采集试验分析,其原理图如图4所示。
星上光通信终端的发散角为30μrad,到达地面站后的远场光束直径为1200m,因此远场光场分布接收探测测试端机的分布应涵盖上述范围。一般来说,远场光场分布接收探测测试端机(高轨卫星星地激光远场分布测试装置)的数量越多,则分布测试的精度越高。考虑到试验成本问题,一般至少需要6套远场光场分布装置(高轨卫星星地激光远场分布测试装置)进行测试。每次测试涵盖光束远场的某一方向半径,最后通过累加给出光束远场二维分布测试结果。

Claims (5)

1.高轨卫星星地激光远场分布测试方法,其特征在于,该方法是基于高轨卫星星地激光远场分布测试装置实现的,该装置包括望远镜(1)、相机(2)、转台(3)和计算机(4);
望远镜(1)安装在转台(3)上,且随转台(3)移动;相机(2)安装在望远镜(1)的物镜上;转台(3)的转动信号输入输出端连接计算机(4)的转动信号输出输入端;相机(2)的拍照信号输入输出端连接计算机(4)的拍照信号输出输入端;
该方法包括下述步骤:
步骤一、计算机(4)向转台(3)发送控制指令,使控制转台(3)移动到指定方位、指定俯仰角度,对高轨卫星进行提前预瞄准;
步骤二、计算机(4)持续读取相机(2)中的光斑信息,并根据光斑信息计算过阈像素数和光斑信息的灰度质心坐标;
步骤三、计算机(4)读取转台(3)的当前方位和俯仰角度,并且计算机(4)控制转台(3)采用螺旋扫描方式扫描高轨卫星,并计算机(4)根据过阈像素数判断转台(3)是否捕获到高轨卫星,若是,则执行步骤四;若否,则返回步骤一;
步骤四、计算机(4)根据步骤二获得的光斑信息的灰度质心坐标持续反馈控制转台(3)移动,完成高轨卫星闭环跟踪,即获得高轨卫星星地激光远场分布测试。
2.根据权利要求1所述的高轨卫星星地激光远场分布测试方法,其特征在于,步骤二中计算机(4)根据光斑信息计算过阈像素数的具体过程如下;
设定光斑信息的光斑像素指定阈值d,将光斑信息的光斑像素灰度值与光斑信息的光斑像素指定阈值d进行比较;
若光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d,则计数一次;
若光斑信息的光斑像素灰度值≤光斑信息的光斑像素指定阈值d,则不计数;
对光斑信息的光斑像素灰度值>光斑信息的光斑像素指定阈值d的计数数目进行求和统计,获得的总的计数数目,即为过域像素数。
3.根据权利要求1所述的高轨卫星星地激光远场分布测试方法,其特征在于,步骤二中计算机(4)根据光斑信息计算光斑信息的灰度质心坐标的具体过程如下;
设定光斑信息的大小为M1×M2维图像,第m1行m2列像元坐标记为(xm1,ym2);M1×M2维图像由m1行m2列像元构成;m1和m2均为整数;
其中,I(xm1,ym2)表示像元的质量密度,即图像上x、y处像素点的灰度值;x、y值均为整数;
则光斑总能量Wm的表达式为:
则光斑总能量Wm的表达式为:
<mrow> <msub> <mi>W</mi> <mi>m</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>1</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>1</mn> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>2</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>2</mn> </msub> </munderover> <mi>I</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
光斑信息的质心坐标的表达式(xc,yc)为:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>x</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mo>&amp;Sigma;</mo> <mi>x</mi> <mi>I</mi> </mrow> <msub> <mi>W</mi> <mi>m</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>1</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>1</mn> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>2</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>2</mn> </msub> </munderover> <msub> <mi>x</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mi>I</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>1</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>1</mn> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>2</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>2</mn> </msub> </munderover> <mi>I</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>y</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mo>&amp;Sigma;</mo> <mi>y</mi> <mi>I</mi> </mrow> <msub> <mi>W</mi> <mi>m</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>1</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>1</mn> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>2</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>2</mn> </msub> </munderover> <msub> <mi>y</mi> <mrow> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mi>I</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>1</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>1</mn> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mn>2</mn> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mn>2</mn> </msub> </munderover> <mi>I</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
4.根据权利要求1所述的高轨卫星星地激光远场分布测试方法,其特征在于,
步骤三中的计算机(4)控制转台(3)采用螺旋扫描方式扫描高轨卫星,即采用等线速的方式从初始瞄准点开始螺旋向外侧进行扫描,具体过程如下:
步骤三一、以初始瞄准点作为坐标原点(0,0)建立极坐标系;螺旋扫描曲线可表示为:
<mrow> <mi>r</mi> <mo>=</mo> <mfrac> <msub> <mi>I</mi> <mi>&amp;theta;</mi> </msub> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <mi>&amp;theta;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中,r——极坐标半径变量;θ——极坐标角变量;Iθ——扫描步长;
步骤三二、从初始瞄准点(0,0)到任意点(r,θ)的扫描时间近似表达为:
<mrow> <mi>T</mi> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>&amp;ap;</mo> <mfrac> <mrow> <msup> <mi>&amp;pi;r</mi> <mn>2</mn> </msup> </mrow> <msubsup> <mi>I</mi> <mi>&amp;theta;</mi> <mn>2</mn> </msubsup> </mfrac> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中,Δt——每步停留时间,其大小取决于扫描模式,任意点(r,θ)即为高轨卫星的坐标。
5.根据权利要求1所述的高轨卫星星地激光远场分布测试方法,其特征在于,
步骤三中计算机(4)根据过阈像素数判断转台(3)是否捕获到高轨卫星,若是,则执行步骤四;若否,则返回步骤一;
若过阈像素数>4,表示转台(3)捕获到高轨卫星;
若过阈像素数≤3,表示转台(3)未捕获到高轨卫星。
CN201511021877.XA 2015-12-29 2015-12-29 高轨卫星星地激光远场分布测试方法 Active CN105606342B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511021877.XA CN105606342B (zh) 2015-12-29 2015-12-29 高轨卫星星地激光远场分布测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511021877.XA CN105606342B (zh) 2015-12-29 2015-12-29 高轨卫星星地激光远场分布测试方法

Publications (2)

Publication Number Publication Date
CN105606342A CN105606342A (zh) 2016-05-25
CN105606342B true CN105606342B (zh) 2018-03-16

Family

ID=55986455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511021877.XA Active CN105606342B (zh) 2015-12-29 2015-12-29 高轨卫星星地激光远场分布测试方法

Country Status (1)

Country Link
CN (1) CN105606342B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111879299B (zh) * 2020-08-17 2022-04-15 中国科学院上海天文台 一种地基望远镜全自动指星方法
CN115603807B (zh) * 2022-11-28 2023-03-28 长春理工大学 全自动光束捕获扫描跟踪系统、方法、存储介质以及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10044522C2 (de) * 2000-09-04 2002-08-01 Fraunhofer Ges Forschung Optische Anordnung zur Strahlführung
CN201213261Y (zh) * 2008-07-04 2009-03-25 中国科学院上海光学精密机械研究所 卫星激光通信双终端双向远距离传输模拟与地面检测装置
CN102412898A (zh) * 2011-11-16 2012-04-11 中国科学院上海光学精密机械研究所 无波面畸变自由空间远距离激光传输模拟装置
CN104034416A (zh) * 2014-05-21 2014-09-10 中国科学院西安光学精密机械研究所 高动态范围激光远场焦斑测量装置及方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10044522C2 (de) * 2000-09-04 2002-08-01 Fraunhofer Ges Forschung Optische Anordnung zur Strahlführung
CN201213261Y (zh) * 2008-07-04 2009-03-25 中国科学院上海光学精密机械研究所 卫星激光通信双终端双向远距离传输模拟与地面检测装置
CN102412898A (zh) * 2011-11-16 2012-04-11 中国科学院上海光学精密机械研究所 无波面畸变自由空间远距离激光传输模拟装置
CN104034416A (zh) * 2014-05-21 2014-09-10 中国科学院西安光学精密机械研究所 高动态范围激光远场焦斑测量装置及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
卫星激光通信端机远场分布测试分析及测试系统设计;王建民 等;《光电子-激光》;20081130;第19卷(第11期);第1467-1473页 *

Also Published As

Publication number Publication date
CN105606342A (zh) 2016-05-25

Similar Documents

Publication Publication Date Title
CN112525162B (zh) 一种无人机测量输电线路影像距离的系统及方法
US7417717B2 (en) System and method for improving lidar data fidelity using pixel-aligned lidar/electro-optic data
CN111896973B (zh) 基于主被动融合的超远距离目标三维运动轨迹预测方法
KR100965678B1 (ko) 항공 정찰 시스템
US7920182B2 (en) Digital camera with non-uniform image resolution
CN104166142B (zh) 一种多元光子计数激光测距的三维成像系统
US6664529B2 (en) 3D multispectral lidar
EP2097715B1 (en) Three-dimensional optical radar method and device which use three rgb beams modulated by laser diodes, in particular for metrological and fine arts applications
CN106788763A (zh) 机载激光通信设备及其控制方法
KR101214081B1 (ko) 항공촬영이미지와 수치정보를 합성 처리하는 영상도화 시스템
CN103675795A (zh) 激光雷达收发光轴自动匹配装置及方法
CN109073467A (zh) 光谱照相机控制装置、光谱照相机控制程序、光谱照相机控制系统、搭载有该系统的飞行体及光谱图像拍摄方法
CN108139115B (zh) 用于定日镜的校准方法
CN104991255A (zh) 一种基于目视原理的多点激光测距雷达
CN105606342B (zh) 高轨卫星星地激光远场分布测试方法
CN104729482A (zh) 一种基于飞艇的地面微小目标侦测系统及方法
CN112098343A (zh) 基于etl新型自对焦一体化高光谱成像探测系统
CN111121725A (zh) 一种基于无人机的高光谱数据和三维建模数据采集系统
CN112598702B (zh) 一种基于合作信标的水下目标全自动识别与跟踪系统和方法
US20150092179A1 (en) Light ranging with moving sensor array
Arvizu et al. ATP subsystem for Optical Communications on a CubeSat
CN112764052A (zh) 防空导弹飞行监测系统
CN108305290A (zh) 一种动目标的精确测谱方法
CN209070117U (zh) 一种激光自动扫描探测装置
Sergeyev et al. Near the ground laser communication system: anisoplantic studies based on the PSF measurements

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant