CN109872268B - 一种高码速率遥感卫星原始数据实时快视方法 - Google Patents
一种高码速率遥感卫星原始数据实时快视方法 Download PDFInfo
- Publication number
- CN109872268B CN109872268B CN201811561450.2A CN201811561450A CN109872268B CN 109872268 B CN109872268 B CN 109872268B CN 201811561450 A CN201811561450 A CN 201811561450A CN 109872268 B CN109872268 B CN 109872268B
- Authority
- CN
- China
- Prior art keywords
- data
- remote sensing
- bit
- cuda
- pixel
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种高码速率遥感卫星原始数据实时快视方法,首先,在GPU设备端建立CUDA和OpenGL互操作的初始化环境,然后,CPU主机端将遥感数据传入GPU设备端显存,然后,采用CUDA并行提取遥感数据中的有效图像数据,并通过CUDA和OpenGL互操作完成遥感影像快视,最后,判断遥感数据是否传输完毕,未传输完毕则返回第二步,传输完毕则结束;本发明能够提高GPU设备端的利用率,提高数据并行处理速度和实现快视性能,能够更加高效和可靠地满足目前高分辨率遥感卫星地面系统对高码速率遥感卫星数传数据实时快视的需求。
Description
技术领域
本发明涉及卫星遥感技术领域,尤其涉及一种高码速率遥感卫星原始数据实时快视方法。
背景技术
目前,随着电子元器件和通信技术水平的提高,遥感卫星有效载荷分辨率越来越高,CCD像元数量越来越多,数传速率越来越高,如吉林一号光学A星空间分辨率可达亚米级,全色为0.72米,多光谱为2.88米,全色CCD像元数量为16384个,双通道编码后数传速率可达688Mbps。而快视系统可在接收遥感数据的同时或接收完成之后即刻生成未经任何校正处理的可视图像,通过它用户可直观地了解有关遥感图像质量、云层覆盖和评估、目标图像快速判别等方面的基本情况,以及进行相机成像质量评价和星地间数传链路测试等。
但是,传统的快视方法目前主要存在以下几点缺陷:
(1)传统的快视方法主要采用CPU来完成像素数据由高阶向低阶的截取处理(如10位截取为8位),然后由GDI+或DirectDraw技术来显示影像。
在数据截取处理方面,遥感数据快视需要从大量的连续数据码流中提取出有效图像数据,并将提取出的数据转化为像素值,最后在相应的窗口上绘制显示,以吉林一号光学A星数据为例,星上将图像数据量化为10位,处理时需要对原始数据码流每10位分割,然后从每一个10位数据中截取出8位作为一个像素点显示。常规的提取图像数据方法是采用CPU遍历全部图像数据,利用移位和位与运算提取出图像数据,这种方法处理耗时较长且长时间占用CPU;
在影像显示方面,基于GDI+的影像快视系统效率很低,不能充分发挥现代计算机的硬件能力,目前已逐渐淘汰,而基于DirectDraw技术的影响快视系统,虽然显示效率高于基于GDI+技术显示的快视系统,但其所实现的硬件并行加速能力不可控,缺乏真正的线程级并行处理能力,且该接口已经被微软废弃。
(2)传统快视方法在硬件设备方面需要配置大内存(通常为128GB)服务器来完成数据处理和显示,造成内存资源的极大浪费,不能充分发挥现代计算机的硬件能力,尤其是在GPU在并行处理后需将并行处理的结果数据移动回或拷贝回计算机主内存,需要进行多次数据移动和复制,容易出现数据丢失,且效率低,未能充分利用GPU在图形图像绘制方面的能力。
发明内容
本发明的目的是提供一种高码速率遥感卫星原始数据实时快视方法,能够提高GPU设备端的利用率,提高数据并行处理速度和实现快视性能,能够更加高效和可靠地满足目前高分辨率遥感卫星地面系统对高码速率遥感卫星数传数据实时快视的需求。
本发明采用的技术方案为:
包括以下步骤:
A、 在GPU设备端建立CUDA和OpenGL互操作的初始化环境;具体过程如下:
a1:初始化OpenGL,选择CUDA设备并允许CUDA设备与OpenGL互操作;
a2:创建OpenGL像素缓冲对象PBO并注册为CUDA资源;
a3:创建OpenGL纹理Texture;
a4:创建矩形绘制窗口顶点和纹理坐标;
B、 CPU主机端将遥感数据传入GPU设备端显存;遥感数据来源于CPU主机端实时接收的原始遥感数据或回放的CPU主机端硬盘中存储的历史遥感数据;
C、 采用CUDA并行提取遥感数据中的有效图像数据,并通过CUDA和OpenGL互操作完成遥感影像快视;具体过程如下:
c1:映射CUDA资源并返回显存地址:将已注册的OpenGL像素缓冲对象PBO映射到CUDA存储空间,获取已映射像素缓冲对象PBO对应的设备指针,然后将获取的设备指针作为参数传入CUDA内核函数中;
c2:执行 CUDA内核函数完成数据并行提取;
c3:解除映射CUDA资源;
c4:互操作影像快视显示;
D、 判断遥感数据是否传输完毕,未传输完毕则返回步骤B,传输完毕则结束。
进一步地,所述步骤c2所述的数据并行提取过程如下:
c2-1.1:从原始数据中连续提取M位原始像素数据;M位像素数据中包括t个8位字节,且M等于8与遥感卫星对应成像载荷的每一位像素的量化位数Q的最小公倍数,
c2-1.2:从M位字节中提取s个Q位像素数据,M=s*Q;
c2-1.3:将提取的每个Q位像素数据分别截取为8位量化数据,获取s个8位量化数据。
进一步地,所述遥感卫星采用吉林一号A星全色载荷遥感卫星,所述步骤c2的数据并行提取过程如下:
c2-2.1:将原始数据每40位对齐;对齐的40位原始像素数据包括五个8位字节;
c2-2.2:从对齐的40位原始像素数据中提取四个10位像素;具体方法为:第一个10位像素取第一个字节+第二个字节的高两位;第二个10位像素取第二个字节的低六位+第三个字节的高四位;第三个10位像素取第三个字节的低四位+第四个字节的高六位;第四个10位像素取第四个字节的低两位+第五个字节;
c2-2.3:将提取的每一个10位像素通过右移两位的方式截取成为8位像素字节。
进一步地,所述步骤b8中CUDA和OpenGL互操作方式显示快视影像采用的方法为:将像素缓冲对象PBO与纹理Texture绑定,将像素缓冲对象PBO作为纹理数据源,以绘制纹理的方式来绘制像素缓冲对象PBO。
进一步地,所述像素缓冲对象PBO和纹理Texture之间通过DMA执行像素数据异步传送。
本发明具有以下有益效果:
(1)通过采用CUDA和OpenGL互操作的方式来完成遥感影像快视,既能够充分发挥GPU在通用计算(GPGPU)方面的能力,也能充分利用GPU在图形图像渲染绘制方面的能力,还可以减少计算机主存和显存之间的数据拷贝次数,在GPU完成数据并行处理后不必再将像素数据传送回主存,而是基于互操作接口将显存中数据作为纹理数据源直接递送给OpenGL纹理绘制接口完成影像快视显示,整个过程不需要任何数据移动或复制;
(2)通过采用GPU完成遥感数据的并行提取,相比于传统的CPU数据提取方案,并行处理时间仅为传统CPU并行提取方案时间的1/4,大大提高数据提取效率,同时,由于采用了基于CUDA和OpenGL互操作的遥感影像快视绘制方案,避免了再将并行处理后的结果数据传回计算机主存的移动或拷贝,可进一步提高遥感卫星原始数据并行处理和实时快视性能。
附图说明
图1为本发明的流程图;
图2为具体实施方式中吉林一号A星全色载荷遥感卫星的CUDA并行提取示意图。
具体实施方式
为了更好地理解本发明,下面结合附图对本发明的技术方案做进一步说明。
本发明所公开的一种高码速率遥感卫星原始数据实时快视方法包括以下步骤:
A、在GPU设备端建立CUDA和OpenGL互操作的初始化环境;具体过程如下:
a1:初始化OpenGL,选择CUDA设备并允许CUDA设备与OpenGL互操作;
a2:创建OpenGL像素缓冲对象PBO并注册为CUDA资源;
a3:创建OpenGL纹理Texture;
a4:创建矩形绘制窗口顶点和纹理坐标。
B、CPU主机端将遥感数据传入GPU设备端显存;遥感数据来源于CPU主机端实时接收的原始遥感数据或回放的CPU主机端硬盘中存储的历史遥感数据。
C、采用CUDA并行提取遥感数据中的有效图像数据,并通过CUDA和OpenGL互操作完成遥感影像快视。
CUDA和OpenGL互操作的基本方式就是先将OpenGL资源映射到CUDA地址空间,再使用CUDA向该地址空间写入所需数据,最后利用OpenGL在屏幕上绘制出图形。可以注册到CUDA的OpenGL资源有缓冲对象、纹理和渲染对象等。
具体过程如下:
c1:映射CUDA资源并返回显存地址:将已注册的OpenGL像素缓冲对象PBO映射到CUDA存储空间,获取已映射像素缓冲对象PBO对应的设备指针,然后将获取的设备指针作为参数传入CUDA内核函数中;
在快视实现过程中,CUDA内核函数不断地从原始码流数据(即遥感数据)中并行地提取出有效图像数据,并将提取的有效图像数据存入经过互操作注册和映射后的像素缓冲对象PBO中。存入像素缓冲对象PBO中后,传统的处理方式是将像素数据拷贝回计算机主存,然后再用OpenGL图形接口传回显存完成二维图像绘制,该方法虽然简单但效率不高。本发明采用将像素缓冲对象PBO与二维纹理绑定,将像素缓冲对象PBO作为二维纹理数据源,以绘制纹理的方式来绘制PBO,具体方法是将PBO与纹理texture绑定,打开二维纹理开关,以四边形的方式将纹理绘制出来。
使用像素缓冲对象PBO的主要优势在于OpenGL中可以在像素缓冲对象PBO和纹理对象之间通过DMA(Direct Memory Access)执行像素数据异步传送,显著提高像素缓冲对象PBO数据上载到纹理对象的性能,而不需要CPU的参与。
c2:执行 CUDA内核函数完成数据并行提取;
由于不同型号卫星对应成像载荷中每一像素的量化位数Q不一样,因此,并行提取不同类型遥感卫星的成像载荷像素数据,并确定原始码流数据中连续提取的像素数据的位数极为重要。本发明采用的并行提取方式具体如下:
c2-1.1:从原始数据中连续提取M位原始像素数据;M位像素数据中包括t个8位字节,且M等于8与遥感卫星对应成像载荷的每一位像素的量化位数Q的最小公倍数;
c2-1.2:从M位字节中提取s个Q位像素数据,M=s*Q;
c2-1.3:将提取的每个Q位像素数据分别截取为8位量化数据,获取s个8位量化数据。
下面以以吉林一号A星全色载荷遥感卫星为例进行并行提取数据过程的进一步描述,具体过程如下:
c2-2.1:将原始数据每40位对齐;对齐的40位原始像素数据包括五个8位字节;
c2-2.2:从对齐的40位原始像素数据中提取四个10位像素;具体方法为:第一个10位像素取第一个字节+第二个字节的高两位;第二个10位像素取第二个字节的低六位+第三个字节的高四位;第三个10位像素取第三个字节的低四位+第四个字节的高六位;第四个10位像素取第四个字节的低两位+第五个字节;
c2-2.3:将提取的每一个10位像素通过右移两位的方式截取成为8位像素字节。
本发明所采取的并行提取数据方式,能够适用于不同类型卫星的数据提取,在提高提取效率的同时,应用更为灵活和广泛;
c3:解除映射CUDA资源;
c4:互操作影像快视显示;
D、判断遥感数据是否传输完毕,未传输完毕则返回步骤B,传输完毕则结束。
本发明能够加速遥感影像并行提取,提高快视效率,同时降低了CPU资源利用率,满足目前高分辨率遥感卫星地面系统对高码速率遥感卫星数传数据实时快视的需求,在遥感卫星地面应用系统中具有极大应用价值和发展前景。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解,其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换,而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的范围。
Claims (5)
1.一种高码速率遥感卫星原始数据实时快视方法,其特征在于:包括以下步骤:
A、在GPU设备端建立CUDA和OpenGL互操作的初始化环境;具体过程如下:
a1:初始化OpenGL,选择CUDA设备并允许CUDA设备与OpenGL互操作;
a2:创建OpenGL像素缓冲对象PBO并注册为CUDA资源;
a3:创建OpenGL纹理Texture;
a4:创建矩形绘制窗口顶点和纹理坐标;
B、CPU主机端将遥感数据传入GPU设备端显存;遥感数据来源于CPU主机端实时接收的原始遥感数据或回放的CPU主机端硬盘中存储的历史遥感数据;
C、采用CUDA并行提取遥感数据中的有效图像数据,并通过CUDA和OpenGL互操作完成遥感影像快视;具体过程如下:
c1:映射CUDA资源并返回显存地址:将已注册的OpenGL像素缓冲对象PBO映射到CUDA存储空间,获取已映射像素缓冲对象PBO对应的设备指针,然后将获取的设备指针作为参数传入CUDA内核函数中;
c2:执行 CUDA内核函数完成数据并行提取;
c3:解除映射CUDA资源;
c4:互操作影像快视显示;
D、判断遥感数据是否传输完毕,未传输完毕则返回步骤B,传输完毕则结束。
2.根据权利要求1所述的高码速率遥感卫星原始数据实时快视方法,其特征在于:所述步骤c2所述的数据并行提取过程如下:
c2-1.1:从原始数据中连续提取M位原始像素数据;M位像素数据中包括t个8位字节,且M等于8与遥感卫星对应成像载荷的每一位像素的量化位数Q的最小公倍数,
c2-1.2:从M位字节中提取s个Q位像素数据,M=s*Q;
c2-1.3:将提取的每个Q位像素数据分别截取为8位量化数据,获取s个8位量化数据。
3.根据权利要求2所述的高码速率遥感卫星原始数据实时快视方法,其特征在于:所述遥感卫星采用吉林一号A星全色载荷遥感卫星,所述步骤c2的数据并行提取过程如下:
c2-2.1:将原始数据每40位对齐;对齐的40位原始像素数据包括五个8位字节;
c2-2.2:从对齐的40位原始像素数据中提取四个10位像素;具体方法为:第一个10位像素取第一个字节+第二个字节的高两位;第二个10位像素取第二个字节的低六位+第三个字节的高四位;第三个10位像素取第三个字节的低四位+第四个字节的高六位;第四个10位像素取第四个字节的低两位+第五个字节;
c2-2.3:将提取的每一个10位像素通过右移两位的方式截取成为8位像素字节。
4.根据权利要求1所述的高码速率遥感卫星原始数据实时快视方法,其特征在于:所述步骤b8中CUDA和OpenGL互操作方式显示快视影像采用的方法为:将像素缓冲对象PBO与纹理Texture绑定,将像素缓冲对象PBO作为纹理数据源,以绘制纹理的方式来绘制像素缓冲对象PBO。
5.根据权利要求4所述的高码速率遥感卫星原始数据实时快视方法,其特征在于:所述像素缓冲对象PBO和纹理Texture之间通过DMA执行像素数据异步传送。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811561450.2A CN109872268B (zh) | 2018-12-20 | 2018-12-20 | 一种高码速率遥感卫星原始数据实时快视方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811561450.2A CN109872268B (zh) | 2018-12-20 | 2018-12-20 | 一种高码速率遥感卫星原始数据实时快视方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109872268A CN109872268A (zh) | 2019-06-11 |
CN109872268B true CN109872268B (zh) | 2022-09-16 |
Family
ID=66917183
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811561450.2A Active CN109872268B (zh) | 2018-12-20 | 2018-12-20 | 一种高码速率遥感卫星原始数据实时快视方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109872268B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114494136A (zh) * | 2021-12-27 | 2022-05-13 | 北京吉威数源信息技术有限公司 | 面向卫星影像覆盖率检测方法、装置、设备及存储介质 |
CN116506621B (zh) * | 2023-06-26 | 2023-09-01 | 中国科学院空天信息创新研究院 | 遥感卫星全分辨率快视图像实时传输方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106023101B (zh) * | 2016-05-16 | 2018-12-18 | 中国资源卫星应用中心 | 一种基于视觉保真的卫星遥感图像处理方法 |
CN106251291A (zh) * | 2016-07-28 | 2016-12-21 | 武汉斗鱼网络科技有限公司 | 利用OpenGL与OpenCL协作实现图像缩放的方法及系统 |
-
2018
- 2018-12-20 CN CN201811561450.2A patent/CN109872268B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109872268A (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3129974B1 (en) | Gradient adjustment for texture mapping to non-orthonormal grid | |
US7612783B2 (en) | Advanced anti-aliasing with multiple graphics processing units | |
GB2467049A (en) | Processing Graphics Using Primitive Distribution Lists | |
US20160203635A1 (en) | Frustum tests for sub-pixel shadows | |
CN109872268B (zh) | 一种高码速率遥感卫星原始数据实时快视方法 | |
US8988434B1 (en) | Text rendering for browsers and mobile based applications | |
EP3964949B1 (en) | Graphics processing method and apparatus | |
US20210343066A1 (en) | Residency map descriptors | |
US11501467B2 (en) | Streaming a light field compressed utilizing lossless or lossy compression | |
CN104360319A (zh) | 一种基于cuda平台的雷达b显处理方法 | |
WO2017048357A1 (en) | Facilitating quantization and compression of three-dimensional graphics data using screen space metrics at computing devices | |
CN114002701A (zh) | 一种实时渲染点云的方法、装置、电子设备及系统 | |
EP3251081B1 (en) | Graphics processing unit with bayer mapping | |
CN114004972A (zh) | 一种图像语义分割方法、装置、设备和存储介质 | |
US10269168B2 (en) | Graphics processing systems | |
CN113393555B (zh) | 一种基于着色器的屏幕填充率统计方法及其系统 | |
CN104268169A (zh) | 一种基于ps软件的遥感影像数据快速处理方法 | |
EP4379647A1 (en) | Render format selection method and device related thereto | |
US9251557B2 (en) | System, method, and computer program product for recovering from a memory underflow condition associated with generating video signals | |
US8427500B1 (en) | Spatially aware sub-pixel rendering | |
EP4350605A1 (en) | Image processing method and apparatus, device, and medium | |
CN112598566A (zh) | 基于gpu的cuda的图像处理方法、装置和存储介质 | |
CN117788609A (zh) | 界面图元的拾取方法、装置、设备及存储介质 | |
TW202322043A (zh) | 小網格著色圖譜 | |
CN115904371A (zh) | 一种基于图像截取的仪表显示方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |