CN113223136A - 复杂电磁环境下飞机表面场强分布的纹理投影映射方法 - Google Patents
复杂电磁环境下飞机表面场强分布的纹理投影映射方法 Download PDFInfo
- Publication number
- CN113223136A CN113223136A CN202110488390.1A CN202110488390A CN113223136A CN 113223136 A CN113223136 A CN 113223136A CN 202110488390 A CN202110488390 A CN 202110488390A CN 113223136 A CN113223136 A CN 113223136A
- Authority
- CN
- China
- Prior art keywords
- data
- texture
- double
- coordinate
- projection
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000013507 mapping Methods 0.000 title claims abstract description 35
- 238000009826 distribution Methods 0.000 title claims abstract description 15
- 230000004927 fusion Effects 0.000 claims abstract description 5
- 230000000007 visual effect Effects 0.000 claims abstract description 5
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 230000009466 transformation Effects 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000009877 rendering Methods 0.000 claims description 14
- 230000005672 electromagnetic field Effects 0.000 claims description 11
- 102100034598 Angiopoietin-related protein 7 Human genes 0.000 claims description 9
- 101100039010 Caenorhabditis elegans dis-3 gene Proteins 0.000 claims description 9
- 101000924546 Homo sapiens Angiopoietin-related protein 7 Proteins 0.000 claims description 9
- 101100030351 Schizosaccharomyces pombe (strain 972 / ATCC 24843) dis2 gene Proteins 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 101100328884 Caenorhabditis elegans sqt-3 gene Proteins 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000009432 framing Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims description 2
- 238000004590 computer program Methods 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 9
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000012800 visualization Methods 0.000 description 9
- 238000007794 visualization technique Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 238000012952 Resampling Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013079 data visualisation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/04—Texture mapping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
Abstract
本发明涉及一种针对复杂电磁环境下飞机表面合成场强的纹理投影映射方法，属于纹理投影映射技术领域。以采集系统的实验数据或由CST软件计算得到的电磁仿真场强数据作为原始数据，再由经纬度坐标选定的电磁环境数据场的空间范围内，利用数据插值和颜色映射的方法并结合飞机所在的高度将该高度所对应的电磁态势合成场强数据的纹理投影到该飞机表面上，通过调节投影的位置与角度使飞机模型上的各点与场强纹理数据上各点一一对应，从而准确的将纹理投影到飞机表面上，最后进行颜色融合设置透明效果。最终以一种简洁、直观、高效的方式展示飞机在空中飞行时所受周围电磁环境影响下的场强分布情况。
Description
技术领域
本发明属于纹理投影映射领域，具体涉及一种复杂电磁环境下飞机表面场强分布的纹理投影映射方法。
背景技术
随着电子技术的不断发展，越来越多的用频装备被投入到军队作战之中，空间中各种各样的辐射干扰源以及战场所处的自然环境导致战场上的电磁信号密度将十分密集，且信号形式复杂多样。科学、直观地描述电磁态势，有助于指挥员及时了解电子信息装备与主战武器在复杂电磁环境下受到的影响，帮助指挥员进行及时准确的科学决策与武器装备的合理部署。
目前电磁态势可视化表现形式分为二维和三维两种方式。其中周倜、程芳等人对于通过统计图表、频谱图、瀑布图等方法对二维电磁态势可视化方法的研究已趋近完善，但二维电磁态势可视化方法并不能生动、直观的展现战场电磁态势。因此，三维电磁态势可视化方法的研究是目前研究的热点问题。三维电磁态势的构建方法主要是面绘制和体绘制两类，其中面绘制法一般先在三维空间数据场中构造出二维平面、曲面等几何图元，然后利用计算机图形学技术实现绘制展现。主要绘制思路是先提取目标物体的表面信息，再根据光照、明暗模型调用绘制算法进行阴影和渲染得到最后的显示图像。直接体绘制也称为体绘制方法，直接由数据生成三维图像，不用构造中间的几何图元，以体素作为基本单元，能够表示数据的内部信息，反映三维数据的整体信息，绘制图像质量高。但是，相应的计算量大，包括数据的采样、重构、重采样、绘制等操作。
但目前电磁态势可视化的研究均主要集中在用频设备的探测范围受复杂电磁环境影响的可视化展现、战场电磁态势3维符号的设计以及战场大范围空间内电磁云图的实现与可视化表征方面，这些均是从数据场可视化的角度进行描述绘制，而复杂电磁环境下各用频装备的武器效能、装备可靠性以及不规则表面装备平台的动态合成场强分布的展示问题并没有给予太多有关可视化表征方面的研究。而一般的纹理映射可视化方法是通过预先定义的纹理坐标进行纹理映射，由于仿真场景中纹理映射的物体表面一般都是由多个三角面片所组成，而纹理映射的坐标是通过插值运算得出，但是在物体表面为凹凸不规则情况下，纹理插值运算会产生不同的梯度，从而导致物体表面上纹理分布不均，进而出现的纹理无法准确映射到模型的具体位置或出现纹理丢失等不可控的情况，同时使用一般的纹理映射技术需要预先确定4个角点的坐标，而投影纹理技术则只需要知道摄像头的坐标和投影方向2方面的信息，相比一般的纹理映射方法，其工作量减少了一半。所以，能够兼顾快速高效的渲染方式，又能清晰直观地反映用频装备受电磁环境局部与整体影响的表面的场强分布可视化方法，是十分具有意义的。
发明内容
要解决的技术问题
为了避免现有技术的不足之处，本发明提出一种复杂电磁环境下飞机表面场强分布的纹理投影映射方法，能够简洁、直观、高效的方式展示飞机在空中飞行时所受周围电磁环境影响下的场强分布情况。
技术方案
一种复杂电磁环境下飞机表面场强分布的纹理投影映射方法，其特征在于步骤如下：
步骤1：获取并载入原始电磁环境体数据
通过实验采集系统获得原始数据或者通过CST软件计算得到的电磁环境离散数据，其中包含了地理空间标定信息和含有坐标信息的离散电磁场强度数据，然后将这些数据载入地理信息系统中；
步骤2：空间坐标系的变换与统一
输入电磁体数据覆盖范围的经纬度坐标，然后将其变换为地面平面坐标系的坐标，再与对应坐标的高度数据相结合，最后将地面平面坐标系转换为地球物理坐标系；
步骤3：电磁环境数据计算与处理
根据用户自行选定的电磁波传播衰减模型来计算在选定模型下的经过传播衰减的电磁场强数据，并将计算所得的电磁场强度数据结果进行归一化处理；
步骤4：设定投影纹理图片大小及纹理层数
根据分辨率要求设置电磁环境体数据的纹理大小即像素点个数，同时根据体数据的高度设置合理的纹理层数；
步骤5：像素点逐层、逐点处理
遍历每一层的像素点，根据像素点附近的数据点插值计算该像素点的场强数据，设置传递函数，将电磁合成功率映射为颜色属性，将颜色(R,G,B)及坐标位置(x,y)赋予各个像素点；
重复步骤5直至遍历完电磁环境体数据中的每个像素点，并存储进缓存中；
步骤6：投影坐标归一化并映射
利用取景变换矩阵和投影变换矩阵把视点范围内的物体模型点依次转换到视点坐标系下和投影空间里，由于经过投影变换后的图元点所在空间是[-1,1]，而所需的纹理空间坐标是[0,1]，所以再经过变换矩阵将上述转换到投影空间坐标系的图元数据转换成最终所需的纹理坐标，然后再将模型物体表面的每个点根据得到的纹理坐标值映射到飞机高度所对应的相应的纹理空间中，获取对应的纹理像素值，再设置融合方式，进行渲染；最终得到电磁场强度数据投影在飞机上的可视化效果。
本发明进一步的技术方案：步骤2中所述的坐标系转换和统一的具体方法是：
经纬度转化为地面平面坐标系：
X＝(longitude+180)*QuiKuandu/360
Y＝(latitude+90)*QuiGaodu/180
地面平面坐标系转化为地球物理坐标系：
AngX＝(-180+X/QiuKuandu*360)*π/180
AngY＝(-90+Y/QiuGaodu*180)*π/180
x＝cos(AngY)*(QiuRad+Y)*cos(AngX)
y＝cos(AngY)*(QiuRad+Y)*sin(AngX)
z＝(QiuRad+Z)*sin(AngY)
其中，longitude为经度，latitude为纬度，QiuKuandu为地球赤道长度，QiuGaodu为二分之一赤道长度，QiuRad为地球半径，(X,Y)为经过转换后得到的平面坐标系下的坐标，(x,y,z)为经过转换后得到的地球物理坐标系下的坐标。
本发明进一步的技术方案：步骤5中所述的数据插值计算的具体方法是：
设原始体数据Data[level][row][col]共有level层，每层的数据大小是col*row，设定每层的纹理大小为TextureWidth*TextureHeight；
(1)计算像素点到其附近数据点的相对距离
int row＝(int)(y*(row-1)/(TextureHeight-1))
int col＝(int)(x*(col-1)/(TextureWidth-1))
double row＝(double)(y*(row-1)/(TextureHeigh-1))
double col＝(double)(x*(col-1)/(TextureWidth-1))
double dfi＝double row-(double)int row
double dfj＝double col-(double)int col
其中int代表当前数据点对应着整形类型数据点的索引值，double代表当前数据点对应着双精度类型数据点的索引值，double dfi、double dfj表示像素点到其附近数据点的相对距离；
(2)计算像素点与附近数据点的距离权值
double dis1＝[1.0/sqrt(dfi*dfi+dfj*dfj)]2
double dis4＝[1.0/sqrt(dfi*dfi+(1.0-dfj)*(1.0-dfj))]2
double dis3＝[1.0/sqrt((1.0-dfi)*(1.0-dfi)+dfj*dfj)]2
double dis2＝[1.0/sqrt((1.0-dfi)*(1.0-dfi)+(1.0-dfj)*(1.0-dfj))]2
double disall＝dis1+dis2+dis3+dis4
(3)获取该像素点附近的数据值
float v1＝Data[level][int row][int col]
float v2＝Data[level][int row+1][int col]
float v3＝Data[level][int row+1][int col+1]
float v4＝Data[level][int row][int col+1]
(4)加权插值
Data＝v1*dis1/disall+v2*dis2/disall+v3*dis3/disall+v4*dis4/disall
Data即为最终该像素点的电磁数据值。
本发明进一步的技术方案：步骤6中所述的投影纹理具体方法是：
(1)新建一光源，使之一直保持在飞机编队正上方并伴随飞机编队的移动而移动，同时在该光源处再设置一个摄像机作为投影视点；
(2)计算投影纹理坐标
根据视点相机的位置、投影角度，物体坐标，对物体顶点进行坐标空间转换，经过世界变换、取景变换和投影变换后空间坐标x,y,z范围都是(-1,1)，而纹理坐标的原点在屏幕左上角且坐标范围为(0,1)，所以还要用一个平移缩放矩阵将投影变换后的投影顶点坐标归一化到(0,1)空间中，最后将归一化后的坐标转换为对应的纹理坐标，即投影纹理坐标。
一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述的方法。
一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现上述的方法。
一种计算机程序，其特征在于包括计算机可执行指令，所述指令在被执行时用于实现上述的方法。
有益效果
本发明提供了一种针对复杂电磁环境下飞机表面合成场强的纹理投影映射方法。以采集系统的实验数据或由CST软件计算得到的电磁仿真场强数据作为原始数据，再由经纬度坐标选定的电磁环境数据场的空间范围内，利用数据插值和颜色映射的方法并结合飞机所在的高度将该高度所对应的电磁态势合成场强数据的纹理投影到该飞机表面上，通过调节投影的位置与角度使飞机模型上的各点与场强纹理数据上各点一一对应，从而准确的将纹理投影到飞机表面上，最后进行颜色融合设置透明效果。最终以一种简洁、直观、高效的方式展示飞机在空中飞行时所受周围电磁环境影响下的场强分布情况。
本发明方法较一般的纹理映射，不需要预先生成纹理坐标，纹理坐标的计算与纹理渲染是同时进行，是一个流水线过程，它们分别由2个线程控制，在新坐标计算完后便交由渲染线程进行渲染，而新的顶点坐标也会在线程运行时继续计算，提高了计算运行效率，并且由于投影纹理的纹理坐标运算是以图元顶点为单位进行计算的，不会存在一般坐标纹理映射对复杂不规则物体表面映射时产生的梯度问题，所以也就避免二维纹理映射至三维物体上造成的纹理扭曲失真等现象。
附图说明
附图仅用于示出具体实施例的目的，而并不认为是对本发明的限制，在整个附图中，相同的参考符号表示相同的部件。
图1本方法的系统框图；
图2投影纹理映射处理后的可视化效果图1：(a)CST软件图；(b)投影纹理效果图；
图3投影纹理映射处理后的可视化效果图2：(a)原图；(b)投影纹理效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。此外，下面描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
图1是本发明提出的基于投影纹理的电磁环境体数据可视化方法的系统框图，具体包括以下步骤：
步骤一：获取和载入原始电磁环境体数据。
实验采集系统获得的原始数据或者通过DST软件计算得到的离散数据，包括地理空间标定信息，确定电磁环境体数据的空间覆盖范围，含有坐标信息的离散电磁场强度数据(x,y,z,E)按照层数、行数和列数的格式存储Data[level][row][col]。其中level为层数索引，row为行数索引，col为列数索引。
步骤二：坐标系转换和统一。
输入电磁环境空间覆盖范围内的经纬度坐标，将其转换为地球表面铺平的平面坐标，再加以高度数据，转化为地球物理坐标系，最后设置渲染绘制时的基准点，计算各点的相对位置，建立渲染相对坐标系，完成从经纬度坐标到渲染坐标的转换和统一。
经纬度转化为地面平面坐标系：
X＝(longitude+180)*QuiKuandu/360
Y＝(latitude+90)*QuiGaodu/180
地面平面坐标系转化为地球物理坐标系：
AngX＝(-180+X/QiuKuandu*360)*π/180
AngY＝(-90+Y/QiuGaodu*180)*π/180
x＝cos(AngY)*(QiuRad+Y)*cos(AngX)
y＝cos(AngY)*(QiuRad+Y)*sin(AngX)
z＝(QiuRad+Z)*sin(AngY)
其中，longitude为经度，latitude为纬度，QiuKuandu为地球赤道长度，QiuGaodu为二分之一赤道长度，QiuRad为地球半径，(X,Y)为经过转换后得到的平面坐标系下的坐标，(x,y,z)为经过转换后得到的地球物理坐标系下的坐标。
步骤三：设置纹理贴图大小及纹理层数。
电磁环境数据计算并归一化后，根据数据高度范围设置对应纹理层数，采用水平堆叠平面纹理的方式表示原始电磁环境体数据。根据需求设置纹理贴图的大小TextureWidth*TextureHeight，TextureWidth为贴图宽度，TextureHeight为贴图高度，它们决定了每层纹理的像素点个数。
步骤四：逐层处理，生成纹理对象。
为每层纹理创建DIB设备无关位图文件，遍历设定每个像素点，为每个像素点进行步骤五。
设原始体数据Data[level][row][col]共有level层，每层的数据大小是col*row，设定每层的纹理大小为TextureWidth*TextureHeight。
(3)计算像素点到其附近数据点的相对距离
int row＝(int)(y*(row-1)/(TextureHeight-1))
int col＝(int)(x*(col-1)/(TextureWidth-1))
double row＝(double)(y*(row-1)/(TextureHeigh-1))
double col＝(double)(x*(col-1)/(TextureWidth-1))
double dfi＝double row-(double)int row
double dfj＝double col-(double)int col
其中int代表当前数据点对应着整形类型数据点的索引值，double代表当前数据点对应着双精度类型数据点的索引值，double dfi、double dfj表示像素点到其附近数据点的相对距离。
(2)计算像素点与附近数据点的距离权值
double dis1＝[1.0/sqrt(dfi*dfi+dfj*dfj)]2
double dis4＝[1.0/sqrt(dfi*dfi+(1.0-dfj)*(1.0-dfj))]2
double dis3＝[1.0/sqrt((1.0-dfi)*(1.0-dfi)+dfj*dfj)]2
double dis2＝[1.0/sqrt((1.0-dfi)*(1.0-dfi)+(1.0-dfj)*(1.0-dfj))]2
double disall＝dis1+dis2+dis3+dis4
(3)获取该像素点附近的数据值
float v1＝Data[level][int row][int col]
float v2＝Data[level][int row+1][int col]
float v3＝Data[level][int row+1][int col+1]
float v4＝Data[level][int row][int col+1]
(4)加权插值
Data＝v1*dis1/disall+v2*dis2/disall+v3*dis3/disall+v4*dis4/disall
Data即为最终该像素点的电磁数据值。
步骤五：数据插值和颜色映射。
根据像素点位置，载入相邻4个原始数据点，插值计算各像素点处的电磁场强度。设置传递函数，剔除无需可视化点，将物理属性映射为颜色属性，赋予各点颜色值(R,G,B)和数据坐标(x,y)。
重复步骤四、步骤五，生成所有层数的纹理对象并存储在纹理缓存中。
步骤六：投影坐标归一化并映射。
利用取景变换、投影变换以及归一化变换矩阵求得物体模型点所对应的投影纹理坐标，再根据编队飞机所在高度，从纹理缓存中选取对应的电磁数据纹理位图，根据物体模型上每个点的纹理坐标逐坐标的对应到纹理空间中，并将纹理数据映射到飞机表面对应位置，设置融合方式，最终得到电磁场强在飞机表面上的投影可视化效果如图2所示。而如图3所示，对于其它图片该技术也同样可以实现很好的投影映射效果。
根据视点相机的位置、投影角度，物体坐标，对物体顶点进行坐标空间转换，经过世界变换、取景变换和投影变换后空间坐标x,y,z范围都是(-1,1)，而纹理坐标的原点在屏幕左上角且坐标范围为(0,1)，所以还要用一个平移缩放矩阵将投影变换后的投影顶点坐标归一化到(0,1)空间中，如此便得到物体的投影纹理坐标。
取景变换：
D3DXMATRIXA16 lightViewMatrix；
D3DXMatrixLookAtLH(&lightViewMatrix,pLightPosition,&D3DXVECTOR3(0,0,0),&D3DXVECTOR3(0,-1,-))；
投影变换：
D3DXMATRIXA16 projMatrix；
D3DXMatrixPerspectiveFovLH(&projMatrix,50.0f*3.14f/180.0f,1.72f,0.25,20.0)变换矩阵：
D3DXMATRIXA16 RangeMappingMatrix(
0.5,0,0,0,
0,0.5,0,0,
0,0,0.5,0,
0.5,0.5,0.5,1)；
投影坐标矩阵：
texViewProjMatrix＝modelMatrix*lightViewMatrix*projMatrix*
RangeMappingMatrix；
texViewProjCoordinate＝texViewProjMatrix*modelCoordinate；
其中，RangeMappingMatrix为对应平移缩放变换矩阵，texViewProjMatrix为归一化后的投影坐标矩阵，texViewProjCoordinate为最终的投影纹理坐标。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明公开的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。
Claims (7)
1.一种复杂电磁环境下飞机表面场强分布的纹理投影映射方法，其特征在于步骤如下：
步骤1：获取并载入原始电磁环境体数据
通过实验采集系统获得原始数据或者通过CST软件计算得到的电磁环境离散数据，其中包含了地理空间标定信息和含有坐标信息的离散电磁场强度数据，然后将这些数据载入地理信息系统中；
步骤2：空间坐标系的变换与统一
输入电磁体数据覆盖范围的经纬度坐标，然后将其变换为地面平面坐标系的坐标，再与对应坐标的高度数据相结合，最后将地面平面坐标系转换为地球物理坐标系；
步骤3：电磁环境数据计算与处理
根据用户自行选定的电磁波传播衰减模型来计算在选定模型下的经过传播衰减的电磁场强数据，并将计算所得的电磁场强度数据结果进行归一化处理；
步骤4：设定投影纹理图片大小及纹理层数
根据分辨率要求设置电磁环境体数据的纹理大小即像素点个数，同时根据体数据的高度设置合理的纹理层数；
步骤5：像素点逐层、逐点处理
遍历每一层的像素点，根据像素点附近的数据点插值计算该像素点的场强数据，设置传递函数，将电磁合成功率映射为颜色属性，将颜色(R，G，B)及坐标位置(x，y)赋予各个像素点；
重复步骤5直至遍历完电磁环境体数据中的每个像素点，并存储进缓存中；
步骤6：投影坐标归一化并映射
利用取景变换矩阵和投影变换矩阵把视点范围内的物体模型点依次转换到视点坐标系下和投影空间里，由于经过投影变换后的图元点所在空间是[-1，1]，而所需的纹理空间坐标是[0，1]，所以再经过变换矩阵将上述转换到投影空间坐标系的图元数据转换成最终所需的纹理坐标，然后再将模型物体表面的每个点根据得到的纹理坐标值映射到飞机高度所对应的相应的纹理空间中，获取对应的纹理像素值，再设置融合方式，进行渲染；最终得到电磁场强度数据投影在飞机上的可视化效果。
2.根据权利要求1所述的一种复杂电磁环境下飞机表面场强分布的纹理投影映射方法，其特征在于步骤2中所述的坐标系转换和统一的具体方法是：
经纬度转化为地面平面坐标系：
X＝(longitude+180)*QuiKuandu/360
Y＝(latitude+90)*QuiGaodu/180
地面平面坐标系转化为地球物理坐标系：
AngX＝(-180+X/QiuKuandu*360)*π/180
AngY＝(-90+Y/QiuGaodu*180)*π/180
x＝cos(AngY)*(QiuRad+Y)*cos(AngX)
y＝cos(AngY)*(QiuRad+Y)*sin(AngX)
z＝(QiuRad+Z)*sin(AngY)
其中，longitude为经度，latitude为纬度，QiuKuandu为地球赤道长度，QiuGaodu为二分之一赤道长度，QiuRad为地球半径，(X，Y)为经过转换后得到的平面坐标系下的坐标，(x，y，z)为经过转换后得到的地球物理坐标系下的坐标。
3.根据权利要求1所述的一种复杂电磁环境下飞机表面场强分布的纹理投影映射方法，其特征在于步骤5中所述的数据插值计算的具体方法是：
设原始体数据Data[level][row][col]共有level层，每层的数据大小是col*row，设定每层的纹理大小为TextureWidth*TextureHeight；
(1)计算像素点到其附近数据点的相对距离
int row＝(int)(y*(row-1)/(TextureHeight-1))
int col＝(int)(x*(col-1)/(TextureWidth-1))
double row＝(double)(y*(row-1)/(TextureHeigh-1))
double col＝(double)(x*(col-1)/(TextureWidth-1))
double dfi＝double row-(double)int row
double dfj＝double col-(double)int col
其中int代表当前数据点对应着整形类型数据点的索引值，double代表当前数据点对应着双精度类型数据点的索引值，double dfi、double dfj表示像素点到其附近数据点的相对距离；
(2)计算像素点与附近数据点的距离权值
double dis1＝[1.0/sqrt(dfi*dfi+dfi*dfj)]2
double dis4＝[1.0/sqrt(dfi*dfi+(1.0-dfj)*(1.0-dfj))]2
double dis3＝[1.0/sqrt((1.0-dfi)*(1.0-dfi)+dfj*dfj)]2
double dis2＝[1.0/sqrt((1.0-dfi)*(1.0-dfi)+(1.0-dfj)*(1.0-dfj))]2
double disall＝dis1+dis2+dis3+dis4
(3)获取该像素点附近的数据值
float v1＝Data[level][int row][int col]
float v2＝Data[level][int row+1][int col]
float v3＝Data[level][int row+1][int col+1]
float v4＝Data[level][int row][int col+1]
(4)加权插值
Data＝v1*dis1/disall+v2*dis2/disall+v3*dis3/disall+v4*dis4/disall
Data即为最终该像素点的电磁数据值。
4.根据权利要求1所述的一种复杂电磁环境下飞机表面场强分布的纹理投影映射方法，其特征在于步骤6中所述的投影纹理具体方法是：
(1)新建一光源，使之一直保持在飞机编队正上方并伴随飞机编队的移动而移动，同时在该光源处再设置一个摄像机作为投影视点；
(2)计算投影纹理坐标
根据视点相机的位置、投影角度，物体坐标，对物体顶点进行坐标空间转换，经过世界变换、取景变换和投影变换后空间坐标x，y，z范围都是(-1，1)，而纹理坐标的原点在屏幕左上角且坐标范围为(0，1)，所以还要用一个平移缩放矩阵将投影变换后的投影顶点坐标归一化到(0，1)空间中，最后将归一化后的坐标转换为对应的纹理坐标，即投影纹理坐标。
5.一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现权利要求1所述的方法。
6.一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现权利要求1所述的方法。
7.一种计算机程序，其特征在于包括计算机可执行指令，所述指令在被执行时用于实现权利要求1所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110488390.1A CN113223136A (zh) | 2021-05-06 | 2021-05-06 | 复杂电磁环境下飞机表面场强分布的纹理投影映射方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110488390.1A CN113223136A (zh) | 2021-05-06 | 2021-05-06 | 复杂电磁环境下飞机表面场强分布的纹理投影映射方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113223136A true CN113223136A (zh) | 2021-08-06 |
Family
ID=77091125
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110488390.1A Pending CN113223136A (zh) | 2021-05-06 | 2021-05-06 | 复杂电磁环境下飞机表面场强分布的纹理投影映射方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113223136A (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103337095A (zh) * | 2013-06-25 | 2013-10-02 | 桂林理工大学 | 一种真实空间三维地理实体的立体虚拟显示方法 |
| CN105336003A (zh) * | 2015-09-28 | 2016-02-17 | 中国人民解放军空军航空大学 | 结合gpu技术实时流畅绘制出三维地形模型的方法 |
| CN108460232A (zh) * | 2018-03-27 | 2018-08-28 | 西北工业大学 | 一种基于纹理映射的电磁环境体数据可视化方法 |
-
2021
- 2021-05-06 CN CN202110488390.1A patent/CN113223136A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103337095A (zh) * | 2013-06-25 | 2013-10-02 | 桂林理工大学 | 一种真实空间三维地理实体的立体虚拟显示方法 |
| CN105336003A (zh) * | 2015-09-28 | 2016-02-17 | 中国人民解放军空军航空大学 | 结合gpu技术实时流畅绘制出三维地形模型的方法 |
| CN108460232A (zh) * | 2018-03-27 | 2018-08-28 | 西北工业大学 | 一种基于纹理映射的电磁环境体数据可视化方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108648269B (zh) | 三维建筑物模型的单体化方法和系统 | |
| CN110383337A (zh) | 可变速率着色 | |
| CN108460232B (zh) | 一种基于纹理映射的电磁环境体数据可视化方法 | |
| CN106991215B (zh) | 一种基于gpu的高精度图形电磁散射计算方法 | |
| CN113674389B (zh) | 场景渲染方法、装置、电子设备及存储介质 | |
| CN110908510B (zh) | 一种倾斜摄影建模数据在沉浸式显示设备中的应用方法 | |
| CN107330964A (zh) | 一种复杂三维物体的显示方法及系统 | |
| CN113593027B (zh) | 一种三维航电显控界面装置 | |
| CN103544731B (zh) | 一种基于多相机的快速反射绘制方法 | |
| CN110706325A (zh) | 一种三维海底环境实时动态渲染方法和渲染系统 | |
| CN110852952A (zh) | 一种基于gpu的大规模地形实时绘制方法 | |
| US10347034B2 (en) | Out-of-core point rendering with dynamic shapes | |
| CN113223136A (zh) | 复杂电磁环境下飞机表面场强分布的纹理投影映射方法 | |
| CN106408644A (zh) | 三维控制笼构造方法 | |
| Xu et al. | Visualization methods of vector data on a Digital Earth System | |
| CN117274465B (zh) | 匹配真实地理水域环境的水体渲染方法、系统、介质和设备 | |
| US20180005432A1 (en) | Shading Using Multiple Texture Maps | |
| CN117456074B (zh) | 基于数字孪生仿真的海上风电冲刷坑三维渲染方法和设备 | |
| CN114627258B (zh) | 一种重力场双体球同构建模方法与系统 | |
| CN117274527A (zh) | 一种发电机设备三维可视化模型数据集构建方法 | |
| Luo et al. | Spherical projective displacement mesh | |
| Wang et al. | 3D Reconstruction of Plant Based on NeRF | |
| Sourin | From Ancient Greeks to Pixels | |
| CN111563947A (zh) | 全球三维云的互动实时体渲染方法 | |
| Guo et al. | Research on rapid loading and measuring method of panorama based on depth image |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210806 |