一种基于离散化处理的行星红外辐射仿真方法技术

本发明专利技术提供了一种基于离散化处理的行星红外辐射仿真方法，包含以下步骤：(1)以行星球心为原点建立平行于地心直角坐标系的行星直角坐标系；(2)将行星表面离散化，并对行星表面的离散面元遍历，判断其是否处于视场范围内；(3)判断各离散面元是否处于日照区域，根据其法线方向和太阳直射方向的夹角，以及行星表面太阳直射点温度和阴影温度，计算该面元温度；(4)根据黑体辐射和小面源辐射模型，计算行星红外辐射量。本发明专利技术在已知行星位置信息的前提下，实现了行星表面温度从太阳直射点到阴影区域的过渡，进而通过黑体辐射模型，计算行星的红外辐射。仿真表明，本发明专利技术能够有效的完成传感器视场内行星红外辐射的仿真。

【技术实现步骤摘要】

本专利技术属于空间目标探测和识别的
，涉及一种基于离散化处理的行星红外辐射仿真方法。
技术介绍
空间背景的红外特性涉及了空间恒星、行星和黄道带的红外辐射，空间背景的红外辐射计算是空间目标探测和识别的重要依据，其中的黄道带影响较弱，恒星由于距离远，并且具有较恒定的辐射光谱，但对于行星的辐射特性，还没有非常有效的仿真计算模型。根据VSOP87C理论可以求得太阳系内行星的初始位置，即行星在地心坐标系、日心坐标系下的位置信息，因此可以借助太阳和行星的几何关系计算行星表面温度，进而计算行星的红外辐射。太阳系的行星在围绕太阳运动的同时，自身也在进行自转，不同的时刻，太阳会照射到行星的不同位置，同时，在地球上观测时，不同的时刻会观测行星不同的部分，由于观测区域和太阳照射区域都是时变的，导致在地球上观测行星的红外辐射时面临较大的困难。本专利技术通过将行星表面离散化为多个面元，由于面元面积较小，可以看作温度均匀的辐射源，通过判断面元法线与太阳直射方向的夹角可以估计面元的温度，通过对全部面元的辐射求和可以计算行星的红外辐射。太阳系行星的红外辐射计算结果将应用于空间目标探测和识别。
技术实现思路
视场范围内的行星温度是行星辐射量仿真的重要参数之一。由于行星表面存在日照区域以及阴影区域，而从太阳直射点到阴影区域的温度是按照一定的规律过渡的，传统的温度计算方法并没有考虑这一点，将日照区域的温度统一认为直射点的温度，从而影响了行星表面温度计算的准确度。本专利技术针对行星表面进行离散化处理，建立平行于地心直角坐标系，以行星球心为原点的直角坐标系，对行星上的离散化面元依次遍历，判断其是否处于传感器视场范围内，对于处于日照区域的面元，计算其法线方向和太阳直射方向的夹角，从而可以计算该面元的温度，对于阴影区域的点，不同的点，由于其温度变化微弱，采用阴影温度统一处理，可以提高温度的计算准确度。由于对行星表面进行了离散化处理，因此每个离散单元均可以看作是小面元,是温度均匀的辐射源，根据黑体辐射模型和小面源辐射模型，通过对全部面元的辐射求和可以计算行星的红外辐射。本专利技术解决其技术问题采用的技术方案为：一种基于离散化处理的行星红外辐射仿真方法，包括以下具体步骤：步骤(1)、建立以行星球心为原点的平行于地心直角坐标系的行星直角坐标系在已知太阳系各行星在地心直角坐标系、日心直角坐标系下的位置信息前提下，计算视场范围内行星表面温度的时候，若采用地心直角坐标系，行星表面的点就很难用坐标直接表示出来，而采用日心直角坐标系存在同样的问题，因此就需要建立新的坐标系，为了计算和表示方便，本专利技术建立了平行于地心直角坐标系，以行星球心为原点的行星直角坐标系。在实际应用中，为了方便计算，首先根据传感器视场指向将地平坐标系转换到沿视线方向坐标系下，得到视线方向坐标系下的行星位置信息，判断该行星是否可以被探测到，若该行星位于视场范围内，进行步骤(2)，否则，改变传感器视轴指向继续遍历。步骤(2)、将行星表面离散化，并对离散面元遍历，判断其是否处于视场观测范围内由于行星直角坐标系，地心直角坐标系，太阳直角坐标系是两两平行的，因此行星坐标系下的地球和太阳位置信息就可以通过简单的坐标转换得到。对行星表面上的点进行离散化处理，依次遍历，求得各离散面元到地球的距离，以及和地球的连线与视轴方向的夹角，根据传感器视场大小以及地球视线和行星相切点到地球的距离，来判断各离散面元是否处于视场内。假设P为行星表面一点，为了方便离散化处理，用球坐标系的形式表示其中r为行星半径，θ为行星上任意一点与Z轴之间的夹角，为行星上任意一点到XOY平面的投影与球心的连线和X轴的夹角。将θ,离散化处理之后，然后将球坐标表示的形式转换到以行星为原点的直角坐标系下的离散面元P(x,y,z)。根据坐标系平行的特性，行星表面离散面元在地心坐标系下的坐标P'(x',y',z')为行星球心在地心坐标系下的坐标加上P相对行星球心的偏移量。根据传感器视轴指向(A,H)将P'(x',y',z')投影到沿视线方向的坐标系下Pr(xr,yr,zr)，A为传感器的方位角，H为俯仰角。然后计算Pr(xr,yr,zr)到地球的连线和视轴之间的夹角α：接下来需要判断当前面元P是否处于地球观测视场范围内，还需要对其到地球观测位置的距离进行判断，由于地球只能观测到面向地球的一侧，当视线与行星相切时，求得可以观测到的点距离地球的最远距离为：其中，为地球到行星球心的距离，r为行星半径。所以仅当P才位于视场范围内，FOV为传感器视场范围。步骤(3)、判断视场内离散面元是否位于日照区域，并进行面元温度计算在对行星上的离散面元位置判断完成后，需要继续判断其是否位于日照区域，同样地，计算各离散面元和太阳之间的距离，根据太阳光线和行星相切点到太阳的距离，判断各离散面元是否处于日照区域，然后计算各离散面元法线方向和太阳直射方向的夹角，求得各离散面元的温度。已知日心坐标系下的行星位置，根据坐标系平行的特性，对日心坐标系下的行星坐标改变符号，即可得到以行星球心为原心的坐标下的太阳坐标，S为太阳位置。求得太阳和行星面向太阳一侧点可以达到的最大距离：为太阳到行星球心的距离，S为太阳的位置，r为行星半径。所以仅当行星表面面元P与太阳的距离满足时，P才位于日照区域。对于日照区域内的面元P，根据余弦定理可以求得行星表面面元P法线方向和太阳直射方向的夹角为β：则行星表面各离散面元的温度tp可以按照下式计算：其中，tsunsub为太阳直射行星点的温度；tshadow为行星阴影温度。步骤(4)、计算行星红外辐射求出行星某个离散小面元温度之后，由于面元面积较小，可以看作温度均匀的辐射源，可利用普朗克函数求解给定波长下行星的光谱辐射出射度，然后求出辐射亮度，根据小面源辐射照度模型，求得视场内的行星面元辐射照度，将各离散面元的辐射求和即可得到给定波长下行星的红外辐射量。最后，若给定传感器的起止波长值，可通过对整个波段上积分求和得到该波段内行星的红外辐射量。本专利技术与现有技术相比的优点在于：(1)本专利技术所建立的直角坐标系以行星球心为原心，并且与地心坐标系平行，在已知行星位置信息前提下，只需要通过简单的向量运算就可以完成行星表面面元的表示和遍历；(2)本专利技术所采用的离散化方法对行星表面进行遍历，日照区域各离散面元的温度根据其法线方向和太阳直射方向夹角求得，从直射点到阴影区域逐渐过渡，而不是采用单一的直射点温度；(3)本专利技术将行星表面离散化为多个面元，由于面元面积较小，可以看作温度均匀的辐射源，通过对全部面元的辐射求和可以完成行星红外辐射的仿真。附图说明图1为本专利技术基于离散化处理计算行星表面温度的流程图；图2为本专利技术建立的平行于地心坐标系的以行星直角坐标系；图3为本专利技术判断行星上的面元是否处于视场范围内的示意图；图4为本专利技术采用的小面源辐射模型。具体实施方式下面结合附图以及具体实施方式进一步说明本专利技术。在视场内行星红外辐射仿真过程中，本专利技术采用行星表面的离散化处理，建立平行于地心坐标系的行星球心直角坐标系，依次遍历行星表面的面元，判断其是否处于视场范围内，对于视场范围内的面元，通过判断面元法线方向与太阳直射方向的夹角可以估计面元的温度，通过对全部面元的辐射求和可以计算行星的红外辐射。本文档来自技高网...

【技术保护点】
一种基于离散化处理的行星红外辐射仿真方法，其特征在于：该方法的实现步骤如下：步骤(1)、建立以行星球心为原点的平行于地心直角坐标系的行星直角坐标系在已知太阳系各行星在地心直角坐标系、日心直角坐标系下的位置信息前提下，为了计算和表示方便，建立了平行于地心直角坐标系，以行星球心为原点的行星直角坐标系，首先根据传感器视轴指向将地平坐标系转换到沿视线方向坐标系下，得到视线方向坐标系下的行星位置信息，判断该行星是否可以被探测到，若该行星位于视场范围内，则进行步骤(2)，否则，改变视轴指向继续遍历；步骤(2)、将行星表面离散化，并对离散面元遍历，判断其是否处于视场观测范围内由于行星直角坐标系，地心直角坐标系，太阳直角坐标系是两两平行的，因此行星坐标系下的地球和太阳位置信息就可以通过简单的坐标转换得到，对行星表面进行离散化处理，依次遍历各离散面元，求得各离散面元到地球的距离，以及和地球的连线与视轴方向的夹角，根据传感器视场大小以及地球视线和行星相切点到地球的距离，来判断各离散面元是否处于视场内，若处于视场内，则进行步骤(3)，否则，继续遍历各离散面元；步骤(3)、判断视场内离散面元是否位于日照区域，并进行面元温度计算在对行星上的离散面元位置判断完成后，需要继续判断各离散面元是否位于日照区域，同样地，计算各离散面元和太阳之间的距离，根据太阳光线和行星相切点到太阳的距离，判断各离散面元是否处于日照区域，通过判断面元法线方向与太阳直射方向的夹角可以估计面元的温度；步骤(4)、计算行星红外辐射求出行星某个离散小面元温度之后，由于面元面积较小，可以看作温度均匀的辐射源，通过对全部面元的辐射求和可以计算行星的红外辐射，可以根据普朗克函数求解给定波长下行星的光谱辐射出射度，然后求出辐射亮度，根据小面源辐射照度模型，求得视场内的行星面元辐射照度，将各离散面元的辐射求和即可得到给定波长下行星的红外辐射量；最后，若给定传感器的起止波长值，可通过对整个波段上积分求和得到该波段内行星的红外辐射量。...

【技术特征摘要】
1.一种基于离散化处理的行星红外辐射仿真方法，其特征在于：该方法的实现步骤如下：步骤(1)、建立以行星球心为原点的平行于地心直角坐标系的行星直角坐标系在已知太阳系各行星在地心直角坐标系、日心直角坐标系下的位置信息前提下，为了计算和表示方便，建立了平行于地心直角坐标系，以行星球心为原点的行星直角坐标系，首先根据传感器视轴指向将地平坐标系转换到沿视线方向坐标系下，得到视线方向坐标系下的行星位置信息，判断该行星是否可以被探测到，若该行星位于视场范围内，则进行步骤(2)，否则，改变视轴指向继续遍历；步骤(2)、将行星表面离散化，并对离散面元遍历，判断其是否处于视场观测范围内由于行星直角坐标系，地心直角坐标系，太阳直角坐标系是两两平行的，因此行星坐标系下的地球和太阳位置信息就可以通过简单的坐标转换得到，对行星表面进行离散化处理，依次遍历各离散面元，求得各离散面元到地球的距离，以及和地球的连线与视轴方向的夹角，根据传感器视场大小以及地球视线和行星相切点到地球的距离，来判断各离散面元是否处于视场内，若处于视场内，则进行步骤(3)，否则，继续遍历各离散面元；步骤(3)、判断视场内离散面元是否位于日照区域，并进行面元温度计算在对行星上的离散面元位置判断完成后，需要继续判断各离散面元是否位于日照区域，同样地，计算各离散面元和太阳之间的距离，根据太阳光线和行星相切点到太阳的距离，判断各离散面元是否处于日照区域，通过判断面元法线方向与太阳直射方向的夹角可以估计面元的温度；步骤(4)、计算行星红外辐射求出行星某个离散小面元温度之后，由于面元面积较小，可以看作温度均匀的辐射源，通过对全部面元的辐射求和可以计算行星的红外辐射，可以根据普朗克函数求解给定波长下行星的光谱辐射出射度，然后求出辐射亮度，根据小面源辐射照度模型，求得视场内的行星面元辐射照度，将各离散面元的辐射求和即可得到给定波长下行星的红外辐射量；最后，若给定传感器的起止波长值，可通过对整个波段上积分求和得到该波段内行星的红外辐射量。2.根据权利要求1中所述的一种基于离散化处理的行星红外辐射仿真方法，其特征在于：所述步骤(2)中的根据行星表面各离散面元和地球连线与视轴方向的夹角，以及地球视线和行星相切点到地球的距离来判断各离散面元是否处于视场范围内的方法：假设P表示行星表面一点，用球坐标系的形式表示该点：其中r为行星半径，θ为行星上任意一点与Z轴之间的夹角，为行星上任意一点到XOY平面的投影与球...

【专利技术属性】
技术研发人员：王文光，屈高龙，李晨鸣，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11