CN114003749A - 一种多波段光学二进制遥感数据可视化方法 - Google Patents

一种多波段光学二进制遥感数据可视化方法 Download PDF

Info

Publication number
CN114003749A
CN114003749A CN202111226746.0A CN202111226746A CN114003749A CN 114003749 A CN114003749 A CN 114003749A CN 202111226746 A CN202111226746 A CN 202111226746A CN 114003749 A CN114003749 A CN 114003749A
Authority
CN
China
Prior art keywords
remote sensing
image
data
screen
show
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.)
Granted
Application number
CN202111226746.0A
Other languages
English (en)
Other versions
CN114003749B (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.)
CETC 54 Research Institute
Original Assignee
CETC 54 Research Institute
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 CETC 54 Research Institute filed Critical CETC 54 Research Institute
Priority to CN202111226746.0A priority Critical patent/CN114003749B/zh
Publication of CN114003749A publication Critical patent/CN114003749A/zh
Application granted granted Critical
Publication of CN114003749B publication Critical patent/CN114003749B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
    • G06F16/54—Browsing; Visualisation therefor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/048—Interaction techniques based on graphical user interfaces [GUI]
    • G06F3/0484—Interaction techniques based on graphical user interfaces [GUI] for the control of specific functions or operations, e.g. selecting or manipulating an object, an image or a displayed text element, setting a parameter value or selecting a range
    • G06F3/04845—Interaction techniques based on graphical user interfaces [GUI] for the control of specific functions or operations, e.g. selecting or manipulating an object, an image or a displayed text element, setting a parameter value or selecting a range for image manipulation, e.g. dragging, rotation, expansion or change of colour
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/048—Interaction techniques based on graphical user interfaces [GUI]
    • G06F3/0484—Interaction techniques based on graphical user interfaces [GUI] for the control of specific functions or operations, e.g. selecting or manipulating an object, an image or a displayed text element, setting a parameter value or selecting a range
    • G06F3/0485—Scrolling or panning
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/048—Interaction techniques based on graphical user interfaces [GUI]
    • G06F3/0487—Interaction techniques based on graphical user interfaces [GUI] using specific features provided by the input device, e.g. functions controlled by the rotation of a mouse with dual sensing arrangements, or of the nature of the input device, e.g. tap gestures based on pressure sensed by a digitiser

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Computer Interaction (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Controls And Circuits For Display Device (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种多波段光学二进制遥感数据可视化方法，涉及遥感卫星影像数据可视化领域。本发明包括：获取数据描述信息；计算显示缩放比例；计算有效显示区显示时的横向和纵向像素个数；计算待显示的实体数据；创建二进制遥感数据缓存区，依据显示需要进行可视化；监听屏幕鼠标操作事件；计算遥感影像缩放后显示区域对应的多波段光学二进制遥感数据区域；计算遥感影像平移后显示区域对应的多波段光学二进制遥感数据区域。本发明在不需要很高的计算机硬件配置条件下，可实现多波段光学遥感影像的流畅实时交互显示。

Description

一种多波段光学二进制遥感数据可视化方法
技术领域
本发明涉及遥感卫星影像数据可视化领域，尤其涉及一种多波段光学二进制遥感数据可视化方法。
背景技术
随着遥感成像技术日益进步，遥感影像的分辨率逐渐提高，数据量也逐渐增大。例如国内的高分系列卫星、吉林一号系列卫星、高景一号卫星的空间分辨率可达或优于亚米级，其中高景一号图像的空间分辨率可达0.5米，单景图像的数据量达数个G字节级别，且独立通道数据的位深通常大于8比特，无法像普通图片一样在操作系统中直接显示。若直接将整个影像文件读入内存中，利用CPU进行位深压缩计算，然后将计算的结果绘制到显示设备中，则会对计算机内存容量的构成挑战，导致大量占用内存，影响其他程序运行，导致图像的显示实时性和流畅性差，不能满足遥感应用处理业务快速人机交互的要求。
目前，用于解决多波段光学遥感卫星影像可视化的问题，采用的主要是图像金字塔和GPU加速技术。
图像金字塔是一系列以金字塔形状排列的分辨率和尺寸自下而上逐渐降低的图像集合。金字塔结构可分为静态金字塔结构和动态金字塔两种。静态金字塔具有算法简单，时间复杂度低的优点，但是需要占用较大的硬盘存储空间；动态金字塔由于是在内存中完成创建过程，所以不占用硬盘空间，但是算法复杂，时间复杂度高。
GPU具有强大的浮点计算能力和并行处理能力，在遥感影像的可视化过程中，可以使用GPU来提高图像渲染效率以及优化显示流畅性。但图像显示效率也往往依赖于GPU自身性能，购置性能强大的GPU设备，为用户带来较高的使用成本。
发明内容
有鉴于此，本发明提出一种多波段光学二进制遥感数据可视化方法，该方法具备计算效率高、图像显示流畅度高、不依赖GPU设备、占用CPU、内存等资源少等优点，既能支持C/S(客户机/服务器)模式的遥感影像可视化业务需求，也能支持B/S(浏览器/服务器)模式的Web网页端的遥感影像可视化业务需求。
本发明采用的技术方案为：
一种多波段光学二进制遥感数据可视化方法，包括以下步骤：
(1)输入光学二进制遥感数据地址，获取数据描述信息；
(2)根据光学二进制遥感数据D的横向像素个数dataw、纵向像素个数datah和屏幕上图像的有效显示区Rscreen的横向像素个数sw、纵向像素个数sh，计算整图显示缩放比例α；
(3)计算整图在屏幕上图像的有效显示区Rscreen显示时的横向像素个数imgw、纵向实际像素个数imgh；
(4)多波段光学二进制遥感数据D中，待显示的实体数据为Dread，在全图显示的情况下，Dread的每个波段的起点行列号(xstart,ystart)为(0,0)，每个波段的列数readw为dataw、每个波段的行数readh为datah，波段数为3；
(5)创建二进制遥感数据缓存区buffer，依据显示需要，从多波段光学二进制遥感数据D中，读取对应区域的实体数据Dread，并进行可视化；
(6)监听屏幕鼠标操作事件，如果为滚动鼠标滚轮，则执行步骤(7)，如果为按下左键并移动，则执行步骤(8)；
(7)计算遥感影像缩放后显示区域Ishow对应的多波段光学二进制遥感数据区域Dread，返回步骤(5)；
(8)计算遥感影像平移后显示区域Ishow对应的多波段光学二进制遥感数据区域Dread，返回步骤(5)。
进一步的，步骤(5)具体包括以下步骤：
(501)为实体数据Dread创建数据缓存区buffer，buffer的大小为datalength个字节，datalength为遥感影像显示区域Ishow的宽度showw和高度showh以及像素位深depth的乘积，即：
datalength＝showw·showh·depth (1)
(502)依据实际显示的需要，指定红色通道、绿色通道、蓝色通道对应的二进制数据波段区域
Figure BDA0003314591170000031
k＝1,2,3；
(503)令k＝1；
(504)对区域
Figure BDA0003314591170000032
内的数据使用最近邻插值法进行空间重采样处理，生成空间重采样影像数据
Figure BDA0003314591170000033
并存入数据缓冲区buffer；
(505)对
Figure BDA0003314591170000034
进行逐像素线性拉伸，将其像素数值取值范围映射至0-255之间，生成待显示图像数据
Figure BDA0003314591170000035
(506)令k＝k+1，当k≤3时执行步骤(504)至(505)，否则执行步骤(507)；
(507)将
Figure BDA0003314591170000041
进行波段合成，生成待显示图像数据Dshow；
(508)将待显示图像数据Dshow，在屏幕上图像的有效显示区Rscreen内进行绘制与显示。
进一步的，步骤(7)具体包括以下步骤：
(701)读取鼠标在Rscreen中的位置P(xview,yview)，其中xview≤sw，yview≤sh；
(702)获取鼠标滚轮滚动前的当前整图逻辑影像I的显示缩放比例η，计算P(xview,yview)与I的左边界、右边界、上边界、下边界的距离li，其中i＝1,2,3,4；
(703)计算鼠标滚轮滚动后的图像显示缩放比例η′；
(704)计算P(xview,yview)与缩放η′倍后的整图逻辑影像I′的四个边界的距离li′，以及I′的左上角点A′在Rscreen中的位置(x′item,y′item)：
li′＝li·η′/η,i＝1,2,3,4 (2)
Figure BDA0003314591170000042
(705)计算I′与Rscreen的交集，获得遥感影像缩放显示区域，并将其记为Ishow；
(706)计算P与Ishow左边界、右边界、上边界、下边界的距离si,i＝1,2,3,4；
(707)计算Ishow的宽度showw和高度showh，以及Ishow的左上角点B在Rscreen中的位置(xshow,yshow)；
(708)计算Ishow对应在多波段光学二进制遥感数据D中的数据区域Dread，每个波段中的起点行列号(xstart,ystart)为：
Figure BDA0003314591170000051
(709)计算多波段光学二进制遥感数据区域Dread每个波段的列数readw和行数readh：
Figure BDA0003314591170000052
进一步的，步骤(702)中计算点P与I四边的距离l1、l2、l3和l4：
Figure BDA0003314591170000053
其中，η为当前整图逻辑影像I的显示缩放比例，其初始值为多波段光学二进制遥感数据D在整图显示的情况下的缩放比例α；(xitem,yitem)为I的左上角点A在Rscreen中的位置，其初始值为D整图显示整图显示的情况下，整图逻辑影像左上点在在Rscreen中的位置(xITEM,yITEM)；dataw和datah分别为光学二进制遥感数据D的横向像素个数、纵向像素个数。
进一步的，步骤(8)具体包括以下步骤：
(801)获取鼠标左键按下前的当前整图逻辑影像I的显示缩放比例η，I左上角点A在屏幕上图像的有效显示区Rscreen位置(xitem,yitem)；
(802)在屏幕上图像的有效显示区Rscreen中，计算鼠标左键按下时的位置P(xview,yview)与鼠标移动至松开左键时对应的位置P′(x′view,y′view)之间的位移(Δx,Δy)；
(803)计算平移后的整图逻辑影像I′左上角点A′在屏幕上图像的有效显示区Rscreen中的位置(x′item,y′item)；
(804)计算P′(x′view,y′view)与I′的左边界、右边界、上边界、下边界的距离l′i，其中i＝1,2,3,4；
(805)计算I′与Rscreen的交集，获得遥感影像平移显示区域，并将其记为Ishow；
(806)计算P′与Ishow左边界、右边界、上边界、下边界的距离si,i＝1,2,3,4；
(807)计算Ishow的宽度showw，高度showh，左上角点B在Rscreen中的位置(xshow,yshow)；
(808)计算Ishow对应的多波段光学二进制遥感数据D中的数据区域Dread，每个波段中的起点行列号(xstart,ystart)为：
Figure BDA0003314591170000061
(809)计算多波段光学二进制遥感数据区域Dread每个波段的列数readw和行数readh：
Figure BDA0003314591170000062
进一步的，步骤(804)中计算点P′与I′四边的距离l1′、l2′、l3′和l4′的方式为：
Figure BDA0003314591170000063
其中，(xitem,yitem)为当前整图逻辑影像I左上角点A在屏幕上图像的有效显示区Rscreen位置，dataw和datah分别为光学二进制遥感数据D的横向像素个数、纵向像素个数。
本发明相比背景技术具有如下优点：
1、本发明采用基于图像显示区域的影像实体数据提取方法，实时提取数据并刷新显示，大幅降低计算机后台处理的数据量，在不依赖GPU和无需生成图像金字塔的条件下，有效解决数据量在1GB以上的遥感影像显示卡顿问题。
2、本发明采用基于人机交互的实时坐标转换方法，通过实时获取鼠标所在位置与屏幕上图像的有效显示区四角点以及当前整图逻辑影像四角点的几何关系，实现物理屏幕坐标与实体数据行列号之间的快速精准转换，使缩放和平移操作更加流畅迅速。
3、本发明使用基于图像显示区域的内存映射使用策略，根据实际显示区域的大小动态分配内存，并采取分段读取的方法，大大降低对计算机内存的需求，有效防止出现内存崩溃的情况。
附图说明
图1是本发明实施例方法的整体流程图。
图2是遥感影像缩放功能实现流程图。
图3是缩放后显示区域尺寸计算示意图。其中图(a)为缩放前点P与当前整图逻辑影像位置关系示意图，图(b)为计算缩放后显示区域示意图。
图4是遥感影像平移功能实现流程图。
图5是平移后显示区域尺寸计算示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，下面结合附图和具体实施步骤对本发明作进一步详细说明。
如图1所示，一种多波段光学二进制遥感数据可视化方法，具体包括以下步骤：
步骤1：读取多波段光学二进制遥感数据D的描述信息，获取数据尺寸、数据类型、波段数、各波段的最大值和最小值；
步骤2：根据屏幕上图像的有效显示区Rscreen的横向像素个数sw、纵向像素个数sh，光学二进制遥感数据D的横向像素个数dataw、纵向像素个数datah，计算整图显示缩放比例α：
α＝min{dataw/sw,datah/sh} (1)
步骤3：计算整图在屏幕上图像的有效显示区Rscreen显示时的横向像素个数imgw、纵向实际像素个数imgh，左上角点在Rscreen的位置(xITEM,yITEM)，其中imgh、imgw和(xITEM,yITEM)的计算公式如下：
Figure BDA0003314591170000081
Figure BDA0003314591170000082
步骤4：计算整图显示时的遥感影像显示区域Ishow的宽度showw和高度showh。此时有showw＝imgw，showh＝imgh；
步骤5：计算在全图显示时，多波段光学二进制遥感数据D中，待显示的实体数据为Dread的波段数、每个波段的起点行列号(xstart,ystart)、每个波段的列数readw、每个波段的行数readh。此时，波段数为3，(xstart,ystart)为(0,0),readw＝dataw，readh＝datah；
步骤6：创建大小为datalength个字节的二进制遥感数据缓存区buffer，datalength的计算公式如下：
datalength＝showw·showh·depth (4)
其中，showw和showh分别为遥感影像显示区域Ishow的宽度和高度，depth为像素位深；
步骤7：依据实际显示的需要，指定红色通道、绿色通道、蓝色通道对应的二进制数据波段区域
Figure BDA0003314591170000091
k＝1,2,3；
步骤8：从Dread第一个波段开始进行可视化前的数据处理，令k＝1；
步骤9：使用最近邻插值法，对区域
Figure BDA0003314591170000092
内的数据进行空间重采样处理，生成空间重采样影像数据
Figure BDA0003314591170000093
并存入数据缓冲区buffer；
步骤10：对
Figure BDA0003314591170000094
进行逐像素线性拉伸，将其像素数值取值范围映射至0-255之间，生成待显示图像数据
Figure BDA0003314591170000095
线性拉伸的计算公式如下：
Figure BDA0003314591170000096
其中，maxk为
Figure BDA0003314591170000097
数据中的最大值，mink为
Figure BDA0003314591170000098
数据中最小值，f(x)为拉伸前的数值，F(x)为拉伸后的数值；
步骤11：令k＝k+1，当k≤3时执行步骤8至步骤10，否则执行步骤12；
步骤12：将
Figure BDA0003314591170000099
进行波段合成，生成待显示图像数据Dshow；
步骤13：将待显示图像数据Dshow，在屏幕上图像的有效显示区Rscreen内进行绘制与显示，并释放缓冲区buffer；
步骤14：监听鼠标交互事件，当鼠标滚轮滚动时，向后台发送缩放显示信号，执行步骤15；当鼠标左键按下并移动时，向后台发送平移信号，执行步骤17。
步骤15：对前端的信号做出响应，调用缩放函数，计算遥感影像缩放后显示区域Ishow对应的多波段光学二进制遥感数据区域Dread，实现原理如图3所示，其具体步骤如图2所示，包括：
步骤1501：读取鼠标在Rscreen中的位置P(xview,yview)，其中xview≤sw，yview≤sh；
步骤1502：获取鼠标滚轮滚动前的当前整图逻辑影像I的显示缩放比例η，如图3(a)所示，计算P(xview,yview)与I的左边界、右边界、上边界、下边界的距离li，其中i＝1,2,3,4：
Figure BDA0003314591170000101
其中，η的初始值为多波段光学二进制遥感数据D在整图显示的情况下的缩放比例α。(xitem,yitem)为I的左上角点A在Rscreen中的位置，其初始值为D整图显示整图显示的情况下，整图逻辑影像左上点在Rscreen中的位置(xITEM,yITEM)。dataw和datah分别为光学二进制遥感数据D的横向像素个数、纵向像素个数。
步骤1503：计算鼠标滚轮滚动后的图像显示缩放比例η′：
η′＝η·ωi k (7)
其中，ωi,i＝1,2为每次滚动时的缩放因子，滚轮向前滚动时为ω1∈(0,1)，向后滚动时为ω2∈(1,∞)，k为滚动次数；
步骤1504：计算P(xview,yview)与缩放η′倍后的整图逻辑影像I′的左边界、右边界、上边界、下边界的距离li′，其中i＝1,2,3,4，以及I′的左上角点A′在Rscreen中的位置(x′item,y′item)；
li′＝li·η′/η,i＝1,2,3,4 (8)
Figure BDA0003314591170000111
步骤1505：计算I′与Rscreen的交集，获得遥感影像缩放显示区域，并将其记为Ishow；
步骤1506：计算计算P与Ishow左边界、右边界、上边界、下边界的距离si,i＝1,2,3,4；
Figure BDA0003314591170000112
步骤1507：计算Ishow的宽度showw和高度showh，以及Ishow的左上角点B在Rscreen中的位置(xshow,yshow)；
Figure BDA0003314591170000113
Figure BDA0003314591170000121
步骤1508：计算Ishow对应在多波段光学二进制遥感数据D中的数据区域Dread，每个波段中的起点行列号(xstart,ystart)；
Figure BDA0003314591170000122
步骤1509：计算多波段光学二进制遥感数据区域Dread每个波段的列数readw和行数readh。
Figure BDA0003314591170000123
步骤16：返回步骤6；
步骤17：对前端的信号做出响应，调用平移函数，计算遥感影像平移后显示区域Ishow对应的多波段光学二进制遥感数据区域Dread，实现流程如图4所示，具体过程如下：
步骤1701：获取鼠标左键按下前的当前整图逻辑影像I的显示缩放比例η，I左上角点A在屏幕上图像的有效显示区Rscreen位置(xitem,yitem)；
步骤1702：在屏幕上图像的有效显示区Rscreen中，计算鼠标左键按下时的位置P(xview,yview)与鼠标移动至松开左键时对应的位置P′(x′view,y′view)之间的位移(Δx,Δy)；
Figure BDA0003314591170000124
步骤1703：计算平移后的整图逻辑影像I′左上角点A′在屏幕上图像的有效显示区Rscreen中的位置(x′item,y′item)；
Figure BDA0003314591170000131
步骤1704：如图5所示，计算P′(x′view,y′view)与I′的左边界、右边界、上边界、下边界的距离l′i，其中i＝1,2,3,4；
Figure BDA0003314591170000132
其中，dataw和datah分别为光学二进制遥感数据D的横向像素个数、纵向像素个数；
步骤1705：计算计算I′与Rscreen的交集，获得遥感影像平移显示区域，并将其记为Ishow；
步骤1706：计算P′与Ishow左边界、右边界、上边界、下边界的距离si,i＝1,2,3,4，计算公式如公式(10)所示；
步骤1707：如图5所示，计算Ishow的宽度showw，高度showh，左上角点B在Rscreen中的位置(xshow,yshow)，计算公式如公式(11)、(12)所示；
步骤1708：计算Ishow对应的多波段光学二进制遥感数据D中的数据区域Dread，每个波段中的起点行列号(xstart,ystart)；
Figure BDA0003314591170000133
步骤1709：计算多波段光学二进制遥感数据区域Dread每个波段的列数readw和行数readh。
Figure BDA0003314591170000141
步骤18：返回步骤6。
总之，本发明在不需要很高的计算机硬件配置条件下，可实现多波段光学遥感影像的流畅实时交互显示。
以上内容是结合具体的实施方式，对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (6)

1.一种多波段光学二进制遥感数据可视化方法，其特征在于，包括以下步骤：
(1)输入光学二进制遥感数据地址，获取数据描述信息；
(2)根据光学二进制遥感数据D的横向像素个数dataw、纵向像素个数datah和屏幕上图像的有效显示区Rscreen的横向像素个数sw、纵向像素个数sh，计算整图显示缩放比例α；
(3)计算整图在屏幕上图像的有效显示区Rscreen显示时的横向像素个数imgw、纵向实际像素个数imgh；
(4)多波段光学二进制遥感数据D中，待显示的实体数据为Dread，在全图显示的情况下，Dread的每个波段的起点行列号(xstart,ystart)为(0,0)，每个波段的列数readw为dataw、每个波段的行数readh为datah，波段数为3；
(5)创建二进制遥感数据缓存区buffer，依据显示需要，从多波段光学二进制遥感数据D中，读取对应区域的实体数据Dread，并进行可视化；
(6)监听屏幕鼠标操作事件，如果为滚动鼠标滚轮，则执行步骤(7)，如果为按下左键并移动，则执行步骤(8)；
(7)计算遥感影像缩放后显示区域Ishow对应的多波段光学二进制遥感数据区域Dread，返回步骤(5)；
(8)计算遥感影像平移后显示区域Ishow对应的多波段光学二进制遥感数据区域Dread，返回步骤(5)。
2.根据权利要求1所述的一种多波段光学二进制遥感数据可视化方法，其特征在于，步骤(5)具体包括以下步骤：
(501)为实体数据Dread创建数据缓存区buffer，buffer的大小为datalength个字节，datalength为遥感影像显示区域Ishow的宽度showw和高度showh以及像素位深depth的乘积，即：
datalength＝showw·showh·depth (1)
(502)依据实际显示的需要，指定红色通道、绿色通道、蓝色通道对应的二进制数据波段区域
Figure FDA0003314591160000021
(503)令k＝1；
(504)对区域
Figure FDA0003314591160000022
内的数据使用最近邻插值法进行空间重采样处理，生成空间重采样影像数据
Figure FDA0003314591160000023
并存入数据缓冲区buffer；
(505)对
Figure FDA0003314591160000024
进行逐像素线性拉伸，将其像素数值取值范围映射至0-255之间，生成待显示图像数据
Figure FDA0003314591160000025
(506)令k＝k+1，当k≤3时执行步骤(504)至(505)，否则执行步骤(507)；
(507)将
Figure FDA0003314591160000026
进行波段合成，生成待显示图像数据Dshow；
(508)将待显示图像数据Dshow，在屏幕上图像的有效显示区Rscreen内进行绘制与显示。
3.根据权利要求1所述的一种多波段光学二进制遥感数据可视化方法，其特征在于，步骤(7)具体包括以下步骤：
(701)读取鼠标在Rscreen中的位置P(xview,yview)，其中xview≤sw，yview≤sh；
(702)获取鼠标滚轮滚动前的当前整图逻辑影像I的显示缩放比例η，计算P(xview,yview)与I的左边界、右边界、上边界、下边界的距离li，其中i＝1,2,3,4；
(703)计算鼠标滚轮滚动后的图像显示缩放比例η′；
(704)计算P(xview,yview)与缩放η′倍后的整图逻辑影像I′的四个边界的距离li′，以及I′的左上角点A′在Rscreen中的位置(x′item,y′item)：
li′＝li·η′/η,i＝1,2,3,4 (2)
Figure FDA0003314591160000031
(705)计算I′与Rscreen的交集，获得遥感影像缩放显示区域，并将其记为Ishow；
(706)计算P与Ishow左边界、右边界、上边界、下边界的距离si,i＝1,2,3,4；
(707)计算Ishow的宽度showw和高度showh，以及Ishow的左上角点B在Rscreen中的位置(xshow,yshow)；
(708)计算Ishow对应在多波段光学二进制遥感数据D中的数据区域Dread，每个波段中的起点行列号(xstart,ystart)为：
Figure FDA0003314591160000032
(709)计算多波段光学二进制遥感数据区域Dread每个波段的列数readw和行数readh：
Figure FDA0003314591160000033
4.根据权利要求3所述的一种多波段光学二进制遥感数据可视化方法，其特征在于，步骤(702)中计算点P与I四边的距离l1、l2、l3和l4：
Figure FDA0003314591160000041
其中，η为当前整图逻辑影像I的显示缩放比例，其初始值为多波段光学二进制遥感数据D在整图显示的情况下的缩放比例α；(xitem,yitem)为I的左上角点A在Rscreen中的位置，其初始值为D整图显示整图显示的情况下，整图逻辑影像左上点在在Rscreen中的位置(xITEM,yITEM)；dataw和datah分别为光学二进制遥感数据D的横向像素个数、纵向像素个数。
5.根据权利要求1所述的一种多波段光学二进制遥感数据可视化方法，其特征在于，步骤(8)具体包括以下步骤：
(801)获取鼠标左键按下前的当前整图逻辑影像I的显示缩放比例η，I左上角点A在屏幕上图像的有效显示区Rscreen位置(xitem,yitem)；
(802)在屏幕上图像的有效显示区Rscreen中，计算鼠标左键按下时的位置P(xview,yview)与鼠标移动至松开左键时对应的位置P′(x′view,y′view)之间的位移(Δx,Δy)；
(803)计算平移后的整图逻辑影像I′左上角点A′在屏幕上图像的有效显示区Rscreen中的位置(x′item,y′item)；
(804)计算P′(x′view,y′view)与I′的左边界、右边界、上边界、下边界的距离li′，其中i＝1,2,3,4；
(805)计算I′与Rscreen的交集，获得遥感影像平移显示区域，并将其记为Ishow；
(806)计算P′与Ishow左边界、右边界、上边界、下边界的距离si,i＝1,2,3,4；
(807)计算Ishow的宽度showw，高度showh，左上角点B在Rscreen中的位置(xshow,yshow)；
(808)计算Ishow对应的多波段光学二进制遥感数据D中的数据区域Dread，每个波段中的起点行列号(xstart,ystart)为：
Figure FDA0003314591160000051
(809)计算多波段光学二进制遥感数据区域Dread每个波段的列数readw和行数readh：
Figure FDA0003314591160000052
6.根据权利要求5所述的一种多波段光学二进制遥感数据可视化技术，其特征在于，步骤(804)中计算点P′与I′四边的距离l1′、l2′、l3′和l4′的方式为：
Figure FDA0003314591160000053
其中，(xitem,yitem)为当前整图逻辑影像I左上角点A在屏幕上图像的有效显示区Rscreen位置，dataw和datah分别为光学二进制遥感数据D的横向像素个数、纵向像素个数。
CN202111226746.0A 2021-10-21 2021-10-21 一种多波段光学二进制遥感数据可视化方法 Active CN114003749B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111226746.0A CN114003749B (zh) 2021-10-21 2021-10-21 一种多波段光学二进制遥感数据可视化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111226746.0A CN114003749B (zh) 2021-10-21 2021-10-21 一种多波段光学二进制遥感数据可视化方法

Publications (2)

Publication Number Publication Date
CN114003749A true CN114003749A (zh) 2022-02-01
CN114003749B CN114003749B (zh) 2022-12-09

Family

ID=79923394

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111226746.0A Active CN114003749B (zh) 2021-10-21 2021-10-21 一种多波段光学二进制遥感数据可视化方法

Country Status (1)

Country Link
CN (1) CN114003749B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105469434A (zh) * 2015-12-31 2016-04-06 中南林业科技大学 一种在Android系统上浏览大型遥感影像的方法
CN113190195A (zh) * 2021-04-25 2021-07-30 中国电子科技集团公司第五十四研究所 一种大幅遥感影像可视化方法及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105469434A (zh) * 2015-12-31 2016-04-06 中南林业科技大学 一种在Android系统上浏览大型遥感影像的方法
CN113190195A (zh) * 2021-04-25 2021-07-30 中国电子科技集团公司第五十四研究所 一种大幅遥感影像可视化方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
扎西央宗等: "积雪遥感监测实时计算显示关键技术研究", 《河北省科学院学报》 *
查东平等: "基于GDAL的遥感影像数据快速读取与显示方法的研究", 《中南林业科技大学学报》 *

Also Published As

Publication number Publication date
CN114003749B (zh) 2022-12-09

Similar Documents

Publication Publication Date Title
JP4696635B2 (ja) 画像領域の高凝縮要約画像を生成する方法、装置およびプログラム
CN102737243B (zh) 获取多幅图像的描述信息的方法及装置与图像匹配方法
CN110533594B (zh) 模型训练方法、图像重建方法、存储介质及相关设备
EP4207051A1 (en) Image super-resolution method and electronic device
WO2020098721A1 (zh) 虚拟现实设备的图像呈现方法、装置、设备和计算机存储介质
CN104660951A (zh) 一种高清转超高清视频图像的超分辨率放大方法
CN108537718B (zh) Gpu中一种实现点图元反走样的装置和方法
CN104537697A (zh) 一种渐进模糊动画的实现方法
CN108694733B (zh) 一种基于mcu的背景图上叠加动态图显示系统及其显示方法
CN111462136B (zh) 一种基于边界分割算法的图像矢量化系统和方法
JP2882465B2 (ja) 画像生成方法およびその装置
CN114003749B (zh) 一种多波段光学二进制遥感数据可视化方法
CN114219716B (zh) 多层图像显示方法和显示引擎
CN105427235A (zh) 一种图像浏览方法及系统
Jia et al. Super-resolution with deep adaptive image resampling
CN104732491A (zh) 一种边缘先验引导的单帧遥感图像超分辨率处理方法
Ding et al. Single image super-resolution via dynamic lightweight database with local-feature based interpolation
CN103489268A (zh) 一种用于pos平台的阿拉伯语显示方法
CN105741335A (zh) 一种基于分块组织的多层次兵棋地图快速绘制方法
CN108629734B (zh) 图像几何变换方法、装置及终端
CN115205115A (zh) 一种基于双线性插值算法的图像放大方法
CN111161340B (zh) 基于深度特征提取的图像重定向方法及系统、终端
US10387997B2 (en) Information processing device, information processing method, and storage medium
Safinaz et al. VLSI realization of Lanczos interpolation for a generic video scaling algorithm
CN110580274A (zh) 一种gis数据渲染方法

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