CN116680871B - 一种全球背景辐射数据及丰度的获取方法、装置及设备 - Google Patents
一种全球背景辐射数据及丰度的获取方法、装置及设备 Download PDFInfo
- Publication number
- CN116680871B CN116680871B CN202310529370.3A CN202310529370A CN116680871B CN 116680871 B CN116680871 B CN 116680871B CN 202310529370 A CN202310529370 A CN 202310529370A CN 116680871 B CN116680871 B CN 116680871B
- Authority
- CN
- China
- Prior art keywords
- polygon
- pixel
- point
- background
- value
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 106
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000004088 simulation Methods 0.000 claims abstract description 92
- 238000004364 calculation method Methods 0.000 claims abstract description 57
- 230000005484 gravity Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000002310 reflectometry Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 abstract description 6
- 230000003287 optical effect Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 14
- 240000004282 Grewia occidentalis Species 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012163 sequencing technique Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 238000013475 authorization Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
- G06F16/9017—Indexing; Data structures therefor; Storage structures using directory or table look-up
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/20—Configuration CAD, e.g. designing by assembling or positioning modules selected from libraries of predesigned modules
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Image Generation (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种全球背景辐射数据及丰度的获取方法、装置及设备,涉及光学卫星成像仿真及计算机图形学领域。本发明通过计算仿真像元与背景影像格网的交点,按照顺时针方向遍历交点,生成边缘裁剪像元对应的多边形及内部完整像元对应的多边形,通过计算各多边形的几何重心坐标及面积获取全部覆盖像元的行列号及丰度。由于基于空间解析几何方法较栅格裁剪方法的计算量小,同时像元坐标系下几何图形的面积即为对应完整像元的占比,而内部完整像元无需计算面积,此外几何计算与栅格数据分离,通过行列号索引从缓存中获取仿真像元背景辐射数据,可减少磁盘I/O操作,从而提高了获取全球背景辐射数据及丰度的速度。
Description
技术领域
本发明涉及光学卫星成像仿真领域,特别涉及一种全球背景辐射数据及丰度的获取方法、装置及设备。
背景技术
全球背景辐射数据包括地表反射率和发射率,通常采用栅格数据格式存储,在进行光学卫星成像仿真时,需要逐像元获取仿真图像各像元覆盖的背景辐射数据像元,从而获取对应的背景辐射数据及所占面积比例。
相关技术中,通常使用卫星影像各像元的覆盖范围作为裁剪范围,对背景辐射数据影像进行栅格裁剪,得到各像元覆盖的背景像元对应的辐射数据,同时基于各背景辐射数据像元的四角坐标构建背景像元的多边形,对背景像元多边形进行矢量裁剪,计算各像元覆盖的背景像元对应的面积比例。
在采用上述方法时,逐像元生成仿真图像需要用每一个仿真像元对应的矢量范围对背景辐射数据做栅格裁剪及矢量裁剪,计算需要耗费较多的时间,并且在做栅格数据裁剪计算时需要先读取栅格数据才能进行裁剪计算,对于百万或千万甚至数亿像元级别的栅格裁剪计算及矢量裁剪计算效率极低,导致成像仿真过程中获取背景辐射数据耗费的时间较多。
发明内容
为解决上述技术问题,本发明提供一种全球背景辐射数据及丰度的获取方法、装置及设备,能够提高获取全球背景辐射数据及丰度的速度。
为达到上述目的,本发明采用的技术方案如下。
第一方面,提供了一种全球背景辐射数据的获取方法,方法包括:
将全球背景辐射数据存储为栅格数据格式,使用像元行列号构建像元唯一标识P_ID(i,j),其中i为行号,j为列号;
将仿真像元覆盖范围记录为多边形,该多边形包含仿真像元的四个顶点,用P1,P2,P3,P4分别表示各顶点,将各点的坐标转换为像元坐标,按照顺时针方向逐点记录多边形各点;
遍历多边形各点,计算多边形与背景影像格网的交点,生成包含背景格网交点的仿真像元多边形点集列表L1;
遍历包含背景格网交点的仿真像元多边形点集列表L1,生成边缘背景像元多边形,计算每个多边形的面积,并将这些多边形加入背景像元多边形列表L2;
遍历内部完整背景像元,生成内部背景像元多边形,其中每个多边形的面积为1,将这些多边形加入背景像元多边形列表L2,得到完整的背景像元多边形列表。
对上述背景像元多边形列表L2,计算每个多边形的几何重心坐标并取整,得到各背景像元多边形对应的背景辐射数据的行列号,根据行列号从全球背景辐射数据中获取对应行列号位置的背景辐射数据值Vi,其中i表示背景像元多边形在背景像元多边形列表L2中的序号;使用各多边形在像元坐标系下的几何面积做为丰度值Ki,计算丰度值Ki与背景辐射数据值Vi的乘积并求和,得到混合后的仿真像元值V。
在一些实施方式中,遍历多边形各点,计算多边形与背景影像格网的交点,方法包括:
创建一个列表作为仿真像元多边形点集列表L1,按照顺时针方向,将多边形第一个点P1作为当前点,加入仿真像元多边形点集列表L1;
取多边形第二个点P2作为下一个点,确定线段P1P2在背景影像格网上的走向及线段P1P2在背景影像格网上的X方向与Y方向的增量dx及dy的符号;
根据线段P1P2的x坐标,确定X方向列坐标的整数起始值和终止值,即线段P1P2在背景影像格网上跨越的列号起始值Sx和列号终止值Ex;
以Cx为X方向游标,取大于P1点x坐标值的最小整数作为当前X方向的游标值,计算线段P1P2在x等于Cx时的y值Cy,得到线段P1P2与背景格网竖线的第一个交点P(Cx,Cy);
根据点P1和P的y坐标,确定Y方向行坐标的整数起始值和终止值,即线段P1P在背景影像格网上跨越的行号的起始值Fy和行号终止值Ty;
以Cy为Y方向游标,根据线段P1P2在Y方向的走向,取第一个整数y坐标作为当前Y方向的游标值,移动Y方向游标Cy,使Cy增加dy,从而遍历Fy到Ty之间的每一个整数值Cy,计算线段P1P2在每一个y等于Cy值时对应的x值,得到线段P1P与背景格网横线的全部交点,并依次按顺序将每一个交点加入仿真像元多边形点集列表L1;
移动X方向游标Cx,使Cx增加dx,重复以上所述步骤,从而遍历Sx到Ex之间的每一个x值,实现线段P1P2与背景格网竖线交点及横线交点的计算。
将P2作为当前点,P3作为下一点,重复以上所述步骤,直到处理完全部仿真像元顶点,得到包含背景格网交点的仿真像元多边形点集列表L1。
在一些实施方式中,遍历包含背景格网交点的仿真像元多边形点集列表L1,生成边缘背景像元多边形,方法包括:
生成背景像元多边形列表L2,按照顺时针方向,从仿真像元多边形点集列表L1中第一个仿真像元多边形与背景格网的交点C1开始,依次构建边缘裁剪像元多边形,并将构建出的多边形加入背景像元多边形列表L2,具体方法为:
取仿真像元多边形点集列表L1中的第2个点和第3个点作为Ca和Cb,根据Cb是否是仿真像元多边形与背景格网的交点,按照不同情形生成裁剪像元多边形。
所述不同情形包括情形1与情形2:
情形1为Cb是仿真像元多边形与背景格网横线或竖线的交点,此时分两种情况:如果Cb的x坐标是整数,则Cb是仿真像元多边形与背景格网竖线的交点,如果Cb的y坐标是整数,则Cb是仿真像元多边形与背景格网横线的交点。对于这两种情况,按顺时针方向,依次添加背景像元坐标格点构成边缘被裁剪的背景像元多边形,并将多边形加入背景像元多边形列表L2中。
情形2为Cb不是仿真像元多边形与背景格网横线或竖线的交点,此时添加仿真像元多边形点集列表L1中的下一个点并按顺时针方向添加背景像元坐标格点,构成多边形顶点处被裁减的背景像元多边形,并将多边形加入背景像元多边形列表L2中。
在一些实施方式中,遍历内部完整背景像元,生成内部背景像元多边形,方法包括:
遍历已经生成的包含边缘被裁剪的背景像元多边形列表L2,计算每个多边形的几何重心,对几何重心坐标取整后得到该多边形所在像元的行列号坐标(R,C),其中R表示行号,C表示列号;
基于行列号坐标,以行号为键,以列号列表为值,构建行号-列号列表键值对Key-Map=(R,[C1,C2]);
根据Key-Map的键值,确定背景像元多边形列表L2所有多边形跨越的行号的起止范围;
对于每一行,对列号列表Key-Map[R]排序,得到每一行列号的起止范围;
遍历每一行的列号,查找列号列表Key-Map[R]中缺失的列值,即可得到对应内部完整像元的行列号坐标(r,c),其中r表示行号,c表示列号;
用行列号坐标(r,c)对应的像元的四角坐标{(r,c),(r,c+1),(r+1,c+1),(r+1,c)}构成一个多边形,将此多边形加入背景像元多边形列表L2。
重复上述步骤,遍历全部行号,即可完成内部完整像元多边形的构建及添加。
在一些实施方式中,使用背景像元多边形列表L2中各多边形的行列号获取对应行列号位置的背景辐射数据Vi,所述多边形的行列号通过计算多边形几何重心坐标并对坐标值取整得到。
在一些实施方式中,使用背景像元多边形列表L2中各多边形在像元坐标系下的几何面积做为丰度值Ki,在像元坐标系下,完整像元面积为1,因此内部完整像元无需计算面积,在具体实施中,只需要计算边缘被裁剪的像元的多边形面积。
第二方面,提供了全球背景辐射数据的获取装置,装置包括:
几何计算模块,用于完成所述第一方面的方法包含的计算过程,包括根据仿真像元范围,计算仿真像元多边形与全球背景影像格网的交点,根据交点构造边缘裁剪多边形与内部完整多边形;
读取及缓存模块,用于读取全球背景辐射数据并缓存在计算机内存中,计算机内存中的图像数据构成一个栅格图像矩阵,像元值可基于全球背景辐射数据像元行列号进行唯一值标识,当进行图像仿真计算时,可根据像元唯一值标识P_ID(i,j)直接索引获取对应背景像元点的数据。
像元混合计算模块,用于对仿真像元裁剪的背景辐射像元进行像元值与丰度值的加权求和,权值为像元占完整像元的面积比例,亦即丰度,背景辐射像元为对应像元的反射率或发射率。
第三方面,提供一种计算机设备,所述设备包括处理器以及存储器:
所述存储器用于存储程序代码,并将所述程序代码传输给所述处理器;
所述处理器用于根据所述程序代码中的指令执行上述方面所述的方法。
本发明上述方法的优点在于:
采用解析几何方法计算仿真像元覆盖的背景辐射数据像元,可实现几何计算与实际数据相分离,几何计算采用解析几何计算方法,其计算量较栅格裁剪计算要少;计算得到裁剪及完整的背景像元多边形后,可通过计算多边形几何重心坐标并取整,直接得到该像元在全球背景辐射数据中对应的行列号,从而便于直接获取背景像元辐射数据;此外多边形的几何图形在像元坐标系下的面积,就是该多边形所占完整像元的面积比例,对于内部完整多边形,面积为1不用计算,因此只需要计算边缘裁剪多边形的面积,从而以较少的计算量得到相应背景辐射像元的丰度值;在全图逐像元仿真计算过程中,可通过几何计算得到的全球背景辐射数据像元的行列号列表,从预先缓存的全球背景辐射数据中获取背景辐射数据值,从而减少重复的磁盘I/O操作。基于以上优点,可实现仿真像元背景辐射数据及丰度的快速获取及仿真像元混合计算,以2000x2000的图像仿真计算为例,可将图像仿真计算耗费的时间从数天降低到数秒。
附图说明
图1是本发明实施例提供的一种全球背景辐射数据及丰度的获取方法的流程图;
图2是本发明实施例提供的一种仿真像元范围裁剪背景格网像元的解析几何计算方法的示意图;
图3是本发明实施例提供的一种计算仿真像元范围与背景格网交点的方法示意图;
图4是本发明实施例提供的一种生成边缘裁剪像元对应的多边形的示意图;
图5是本发明实施例提供的一种生成边缘裁剪像元对应的多边形时对于格网交点如何生成多边形的规则的示意图;
图6是本发明实施例提供的一种生成边缘裁剪像元时对于仿真像元顶点如何生成多边形的规则的示意图;
图7是本发明实施例提供的一种获取全球背景辐射数据及丰度进行像元混合计算的流程图;
图8是本发明实施例提供的一种全球背景辐射数据及丰度的获取装置800的结构示意图;
图9是本发明实施例提供的一种计算机设备900的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
下面对本发明实施例涉及的一些术语概念做解释说明。
内存:计算机主存储器的简称,包括只读存储器(read-only memory,ROM)和随机存取存储器(random access memory,RAM),本实施例中例如是随机存取存储器。
栅格数据:以二维矩阵的形式来表示空间地物或现象分布的数据组织方式,每个矩阵的单元称为一个像元,每个像元表示一个信息值。
波段:在遥感技术中,通常按遥感器的响应范围将电磁波谱划分为多个波段,每个波段对应一个二维矩阵,多个波段构成一个栅格数据文件,在本发明中泛指栅格数据文件中的波段。
在基于全球背景辐射数据进行卫星视场范围内图像逐像元仿真计算的过程中,需要根据仿真像元在背景辐射数据上的投影范围,获取落入仿真像元范围内的背景辐射数据像元及每个像元的丰都值,这个过程实际上是一个逐块对背景辐射数据影像进行裁剪的过程,由于栅格数据裁剪计算本身需要耗费较多的时间,而对于一景诸如2000x2000的图像则需要进行400万次裁剪计算,这种规模的计算量耗费时间在数天之久;此外每次进行整景全图像仿真的过程中,都需要读取背景辐射数据,产生的磁盘I/O操作也需要耗费较多的时间。因此,如果采用栅格裁剪计算方法确定落入仿真像元范围内的背景辐射数据及丰度,同时不对磁盘I/O操作进行优化,对于像元数达到数亿级的大幅宽影像的仿真计算,单景影像的计算时间可达数月之久,耗费时间较长,严重影像了图像仿真效率。
基于此,本发明的一些实施例中,提供了一种基于解析几何计算获取背景辐射数据的方法,通过计算仿真像元范围与背景辐射数据格网的交点,逐点构建边缘裁剪像元多边形,之后通过计算过程中记录的行列范围,添加内部完整像元多边形,从而快速生成裁剪影像的几何数据,通过计算各像元多边形的几何重心坐标值并且取整得到对应的行列号索引,通过计算各像元多边形的面积,得到对应像元的丰度值,从而将计算量降低到最小;同时采用像元行列号做索引将几何计算与实际数据分离开,通过像元行列号从缓存的影像数据中获取像元数据值,从而避免了栅格裁剪计算中栅格数据与几何信息相耦合的问题及多次仿真计算中的磁盘重复I/O问题。基于以上优势,从而极大的节省了背景辐射数据及丰度的获取时间,以2000x2000的图像仿真为例,可将图像仿真计算过程中获取背景辐射数据及丰度的时间耗费量从数天降低到数秒级。
下面结合图1与图2,对本发明实施例提供的一种全球背景辐射数据及丰度的获取方法进行介绍。参见图1,该图是本发明实施例提供的一种全球背景辐射数据及丰度的获取方法的流程图,图1所示方法包括步骤S101至步骤S104,图2是本发明实施例提供的一种仿真像元范围裁剪背景像元格网的解析几何计算方法的示意图,是图1所示步骤S101至S104的具体案例的示意图。
步骤S101,仿真像元范围点集经纬度坐标转像元坐标。
仿真像元范围为根据卫星轨道位置及相机成像参数计算得到的影像像元对应地面的范围,包含四个顶点坐标,用经纬度表示。仿真计算需要获取的,正是仿真像元范围对应地面范围内的背景辐射数据及丰度。地面背景辐射数据为包含全球范围的栅格数据集。全球背景辐射数据为可见光各波段的反射率、红外波段的发射率等,包含至少一个波段。仿真图像也是栅格数据,其中的某一个像元为图像上的一个格点,其范围根据成像几何关系,投影到地面后,受相机观测角度及视场角度影响,表现为发生一定畸变的多边形,该多边形采用经纬度坐标系表示。全球背景辐射数据也采用经纬度坐标系,可以根据像元行列号及图像分辨率计算每一个像元中心点的经纬度坐标。对于栅格数据,除了表示其真实世界位置的经纬度坐标外,还包含栅格行列坐标,也称为像元坐标,通常在图像处理及计算过程中,使用像元坐标较为方便。像元坐标采用像元坐标系表示,取栅格矩阵左上角为(0,0)点,向右为列坐标,向下为行坐标。经纬度坐标与像元坐标的转换,通过仿射变换实现,转换方法为:
X=X0+C*Rx+R*Ty
Y=Y0+C*Tx+R*Ry
其中,X,Y表示点的经纬度坐标值,R表示行坐标,C表示列坐标,Rx表示X方向分辨率,Ry表示Y方向分辨率,Tx表示X方向旋转系数,Ty表示Y方向旋转系数。
根据以上坐标转换关系,用仿真像元四角坐标构成仿真像元多边形,将仿真像元多边形各点的坐标从经纬度坐标转为像元坐标系下的坐标,按顺时针方向,依次记录为点P1,P2,P3,P4。
在一些实施方式中,将全球背景辐射数据存储为栅格数据格式,使用像元行列号构建像元唯一标识P_ID(i,j),其中i为行号,j为列号;由空间某一点的经纬度坐标,计算得到行列号,从而直接取得对应矩阵位置的背景辐射数据。
步骤S102,遍历多边形各点,计算多边形与背景影像格网的交点,生成包含交点的多边形点集L1。如图3所示,图3是本发明实施例提供的一种计算仿真像元范围与背景格网交点的方法示意图,具体方法为:
步骤S102.1,创建一个列表作为仿真像元多边形点集列表L1,按照顺时针方向,将多边形第一个点P1作为当前点,加入仿真像元多边形点集列表L1中;
步骤S102.2,取多边形第二个点P2,计算P2与P1在X方向与Y方向的差,确定X方向增量dx与Y方向增量dy的符号:
dx=sign(P2.x-P1.x)
dy=sign(P2.y-P1.y)
其中,sign(t)为符号函数,当变量t大于0时值为1,当t等于0时值为0,当t小于0时值为-1;P1.x及P2.x分别表示点P1与P2的x坐标,P1.y及P2.y分别表示点P1与P2的y坐标。
步骤S102.3,确定当前点P1在背景格网上的列号Cx及P1到P2在X方向的起始列号Sx和终止列号Ex,其中:
其中,Cx为P1在背景格网上的列号,作为X方向游标的当前值;Int(t)为取整函数,返回不大于t的最大整数值。
步骤S102.4,确定当前点P1在背景格网上的行号Cy及P1到P2在Y方向的起始行号Sy和终止行号Ey,其中:
其中,Cy为P1在背景格网上的行号,作为Y方向游标的当前值。
步骤S102.5,以X方向为主方向,以Cx作为X方向游标,Cy为Y方向游标,以Sx为X方向起始值,Ex为X方向终止值,dx为X方向增量,计算线段P1P2与背景格网竖线的交点,具体做法如下:
步骤S102.5A,取Cx作为x值,计算x值对应的y值,得到线段P1P2与竖线的交点C2,计算方法如下:
步骤S102.5B,计算P1与C2之间线段P1P2在Y方向穿过时与横线产生的交点C1,计算方法为:
步骤S102.5B1,取Fy为P1点的y坐标值,Ty为C2点的y值,当Cy值介于Fy与Ty之间时,取y’为Cy,计算对应的x值x’,计算方法为:
步骤S102.5B2,记点(x’,y’)为C1,将C1加入仿真像元多边形点集列表L1;
步骤S102.5B3,移动Y方向游标,使Cy增加dy;
重复步骤S102.5B,直到游标Cy超越Ty值;
步骤S102.5C,将交点C2(x,y)加入仿真像元多边形点集列表L1;
取多边形第二个点为当前点,第三个点为下一点,重复步骤S102.1-S102.5,直到处理完全部多边形点。
图3所示方法是上述步骤S102中计算线段P1P2与第一个格网竖线交点C2及线段P1C2之间线段P1P2与格网横线交点的具体示例。图3中P1,P2,P3,P4为仿真像元的四个角点,C1-C20依次为仿真像元范围与背景格网的交点。
步骤S103,遍历仿真像元多边形点集列表L1,生成边缘背景像元多边形,计算每个多边形的面积,并将这些多边形加入背景像元多边形列表L2,图4是此步骤的处理过程示意图,其中S401-S420为按照顺时针方向依次构建边缘背景像元多边形的过程,具体方法为:
生成背景像元多边形列表L2,按顺时针方向,从多边形点集列表L1中第一个多边形与背景格网的交点C1开始,依次构建边缘背景像元多边形,并将构建出的多边形加入背景像元多边形列表L2,具体方法如下:
步骤S103.1,构建多边形点集列表B1,该列表用于记录当前多边形;
步骤S103.2,取仿真像元多边形点集列表L1中的第2个点C1作为当前点Ca,取仿真像元多边形点集列表L1中第三个点C2作为下一点Cb,将Ca和Cb加入多边形点集列表B1;
步骤S103.3,判断Cb是否是仿真像元多边形与背景影像格网的交点,判断方法为:
如果Cb的x坐标为整数,则Cb是仿真像元多边形与背景影像格网竖线Int(Cb.x)的交点;
如果Cb的y坐标为整数,则Cb是仿真像元多边形与背景影像格网横线Int(Cb.y)的交点;
如果Cb的x坐标与y坐标均不为整数,则Cb不是仿真像元多边形与背景影像格网的交点,而是仿真像元多边形的顶点;
根据Cb是否是仿真像元多边形与背景影像格网的交点,产生两种情形:
情形1,Cb是仿真像元多边形与背景影像格网的交点。此时,如图5所示,图5提供了一种生成边缘裁剪像元对应的多边形时对于格网交点如何生成多边形的规则,对于情形1,又分情形C501和情形C502两种情况:
情形C501,Cb是仿真像元多边形与背景影像格网竖线的交点,此时通过点Cb的横坐标Cb.x是否大于点Ca的横坐标Ca.x判断线段CaCb在X轴方向的走向,根据线段CaCb在X轴方向的走向,分C501A和C501B两种情形:
情形C501A,如果线段CaCb的走向为X轴正方向,此时通过点Cb的纵坐标Cb.y是否大于点Ca的纵坐标Ca.y判断线段CaCb在Y轴方向的走向,根据线段CaCb在Y轴方向的走向,分C501A1和C501A2两种情形;
情形C501A1,如果线段CaCb的走向为Y轴反方向,则取y=Int(Cb.y)+1,这种情况下,分两种情形:
情形C501A1-1,如果y等于Ca的y坐标,则添加点Cc=(Int(Cb.x),Int(Cb.y)+1)到多边形点集列表B1;
情形C501A1-2,如果y不等于Ca的y坐标,则添加点Cc=Int(Cb.x,Int(Cb.y)+1)及Cd=(Cb.x-1,Int(Cb.y)+1)到多边形点集列表B1;
情形C501A2,如果线段CaCb的走向为Y轴正方向,则取x=Int(Cb.x)-1,这种情况下,分两种情形:
情形C501A2-1,如果x不等于Ca的x坐标,则添加点Cc=(Cb.x,Int(Cb.y)+1),Cd=(Cb.x-1,Int(Cb.y)+1),Ce=(Cb.x-1,Int(Cb.y))到多边形点集列表B1;
情形C501A2-2,如果x等于Ca的x坐标,则添加点Cc=(Cb.x,Int(Cb.y)+1),Cd=(Cb.x-1,Int(Cb.y)+1)到多边形点集列表B1;
情形C501B,如果线段CaCb的走向为X轴反方向,此时通过点Cb的纵坐标Cb.y是否大于点Ca的纵坐标Ca.y判断线段CaCb在Y轴方向的走向,此时有C501B1和C501B2两种情形;
情形C501B1,如果线段CaCb的走向为Y轴正方向,则取y=Int(Cb.y),这种情况下,分两种情形:
情形C501B1-1,如果y等于Ca的y坐标,则添加点Cc=(Int(Cb.y),Cb.x)到多边形点集列表B1;
情形C501B1-2,如果y不等于Ca的y坐标,则添加点Cc=(Cb.x,Int(Cb.y))及Cd=(Cb.x+1,Int(Cb.y))到多边形点集列表B1;
情形C501B2,如果线段CaCb的走向为Y轴反方向,则取x=Int(Cb.x)+1,这种情况下,分两种情形:
情形C501B2-1,如果x等于Ca的x坐标,则添加点Cc=(Cb.x,Int(Cb.y)),Cd=(Cb.x+1,Int(Cb.y))到多边形点集列表B1;
情形C501B2-2,如果x不等于Ca的x坐标,则添加点Cc=(Cb.x,Int(Cb.y)),Cd=(Cb.x+1,Int(Cb.y)),Ce=(Cb.x+1,Int(Cb.y)+1)到多边形点集列表B1;
情形C502,如果Cb是仿真像元多边形与背景影像格网横线的交点,此时通过点Cb的横坐标Cb.x是否大于点Ca的横坐标Ca.x判断线段CaCb在X轴方向的走向,此时分C502A和C502B两种情形:
情形C502A,如果线段CaCb的走向为X轴正方向,此时通过点Cb的纵坐标Cb.y是否大于点Ca的纵坐标Ca.y判断线段CaCb在Y轴方向的走向,此时分C502A1和C502A2两种情形;
情形C502A1,如果线段CaCb的走向为Y轴反方向,则取x=Int(Cb.x)-1,这种情况下,分两种情形:
情形C502A1-1,如果x等于Ca的x坐标,则添加点Cc=(Int(Cb.x)+1,Cb.y),Cd=(Int(Cb.x)+1,Cb.y+1),Ce=(Int(Cb.x,Cb.y+1)到多边形点集列表B1;
情形C502A1-2,如果x不等于Ca的x坐标,则添加点Cc=(Int(Cb.x)+1,Cb.y),Cd=(Int(Cb.x)+1,Cb.y+1)到多边形点集列表B1;
情形C502A2,如果线段CaCb的走向为Y轴正方向,取y=Cb.y-1,这种情况下,分两种情形:
情形C502A2-1,如果y等于Ca的y坐标,则添加点Cc=(Int(Cb.x),Cb.y),Cd=(Int(Cb.x),Cb.y-1)到多边形点集列表B1;
情形C502A2-2,如果y不等于Ca的y坐标,则添加点Cc=(Int(Cb.x),Cb.y)到多边形点集列表B1;
情形C502B,如果线段CaCb的走向为X轴反方向,此时通过点Cb的纵坐标Cb.y是否大于点Ca的纵坐标Ca.y判断线段CaCb在Y轴方向的走向,此时有C502B1和C502B2两种情形;
情形C502B1,如果线段CaCb的走向为Y轴正方向,则取y=Cb.y-1,这种情况下,分两种情形:
情形C502B1-1,如果y等于Ca的y坐标,则添加点Cc=(Int(Cb.x),Cb.y),Cd=(Int(Cb.x),Cb.y-1)到多边形点集列表B1;
情形C502B1-2,如果y不等于Ca的y坐标,则添加点Cc=(Int(Cb.x),Cb.y),Cd=(Int(Cb.x),Cb.y-1),Ce=(Int(Cb.x)+1,Cb.y-1)到多边形点集列表B1;
情形C502B2,如果线段CaCb的走向为Y轴反方向,则取x=Int(Cb.x)+1,这种情况下,分两种情形:
情形C502B2-1,如果x等于Ca的x坐标,则添加点Cc=(Int(Cb.x)+1,Cb.y)到多边形点集列表B1;
情形C502B2-2,如果x不等于Ca的x坐标,则添加点Cc=(Int(Cb.x)+1,Cb.y),Cd=(Int(Cb.x)+1,Cb.y+1)到多边形点集列表B1;
情形2,Cb不是仿真像元多边形与背景影像格网交点,而是仿真像元多边形的顶点,此时取仿真像元多边形点集列表L1中的下一个点Cc加入多边形点集列表B1。
对于情形2,图6提供了一种生成边缘裁剪像元时对于仿真像元顶点如何生成多边形的规则的示意图,如图6所示,根据点Cc所在格网的位置,又分情形C601和情形C602两种情况:
情形C601,Cc是仿真像元多边形与背景影像格网竖线的交点,此时通过判断点Cc的纵坐标Cc.y是否大于点Cb的纵坐标Cb.y确定线段CbCc在Y轴方向的走向,此时又可分C601A和C601B两种情况:
情形C601A,如果线段CbCc在Y轴方向的走向为Y轴正方向,则取y=Int(Cc.y)+1,这种情况下分两种情形:
情形C601A-1,如果y等于Ca的y坐标,则添加点Cd=(Cc.x,Int(Cc.y)+1到多边形点集列表B1;情形C601A-2,如果y不等于Ca的y坐标,则添加点Cd=(Cc.x,Int(Cc.y)+1,Ce=(Cc.x-1,Int(Cc.y)+1)到多边形点集列表B1;
情形C601B,如果线段CbCc在Y轴方向的走向为Y轴反方向,则取y=Int(Cc.y),这种情况下分两种情形:
情形C601B-1,如果y等于Ca的y坐标,则添加点Cd=(Cc.x,Int(Cc.y))到多边形点集列表B1中;
情形C601B-2,如果y不等于Ca的y坐标,则添加点Cd=(Cc.x,Int(Cc.y)),Ce=(Cc.x+1,Int(Cc.y))到多边形点集列表B1中;
情形C602,如果Cc是仿真像元多边形与背景影像格网横线的交点,此时通过点Cc的横坐标Cc.x是否大于点Cb的横坐标Cb.x判断CbCc在X轴方向的走向,此时分C602A和C602B两种情况:
情形C602A,如果线段CbCb在X轴方向的走向为X轴正方向,则取x=Int(Cc.x),这种情况下分两种情形:
情形C602A-1,如果x等于Ca的x坐标,则添加点Cd=(Int(Cc.x),Cc.y)到多边形点集列表B1中;
情形C602A-2,如果x不等于Ca的x坐标,则添加点Cd=(Int(Cc.x)+1,Cc.y),Ce=(Int(Cc.x)+1,Cc.y+1)到多边形点集列表B1中;
情形C602B,如果线段CbCb在X轴方向的走向为X轴反方向,则取x=Int(Cc.x)+1,这种情况下分两种情形:
情形C602B-1,如果x等于Ca的x坐标,则添加点Cd=(Int(Cc.x)+1,Cc.y)到多边形点集列表B1中;
情形C602B-2,如果x不等于Ca的x坐标,则添加点Cd=(Int(Cc.x),Cc.y),Ce=(Int(Cc.x),Cc.y-1)到多边形点集列表B1中;
此时多边形点集列表B1构成一个完整的边缘多边形,计算多边形B1的面积,并将多边形点集列表B1加入背景像元多边形列表L2;
重复步骤S103.2-S103.3,直到处理完仿真像元多边形点集列表L1中的全部点;
步骤S104,遍历内部完整背景像元,生成内部背景像元多边形,其中每个多边形的面积为1,将这些多边形加入背景像元多边形列表L2,得到完整的背景像元多边形列表。构造内部完整像元对应的多边形,其方法为:
步骤S104.1,遍历步骤S103产生的边缘多边形,计算每一个多边形的几何重心,计算方法为:
其中,Cx,Cy分别表示当前边缘多边形的几何重心的x坐标和y坐标,xi和yi表示当前边缘多边形的第i个点的x坐标和y坐标,n表示构成当前边缘多边形的点的个数。
步骤S104.2,对上述边缘多边形几何重心坐标值Cx,Cy分别取整,得到该多边形对应像元的行列号R,C,即:
R=Int(Cy)
C=Int(Cx)
其中Int(t)表示取整函数,该函数取得不大于变量t的最大整数。
在上述处理过程中,使用行号做键,使用列号列表做值,构建行号-列号列表键值对Key-Map=(R,[C1,C2]),其中R表示行号,[C1,C2]表示列号列表,C1,C2为列号;
步骤S104.3,根据上述Key-Map的键值范围,确定仿真像元覆盖的行号范围,对每一行的列号列表Key-Map[R]进行排序后确定列号的最大值与最小值,查找最大值与最小值之间缺失的列号值,即可得到内部完整像元的行列号坐标(r,c),其中r为行号,c为列号;
步骤S104.4,对上述每一个内部完整像元,构造其对应的多边形,其多边形由像元的四个格点构成,分别是(r,c),(r,c+1),(r+1,c+1),(r+1,c),将此多边形加入背景像元多边形列表L2,从而得到仿真像元覆盖的全部背景像元。
在本实施例提供的方法中,考虑到使用全球背景辐射数据对应的栅格坐标作为栅格数据的唯一标识P_ID(i,j)(其中i表示行号,j表示列号),在计算过程中,无需读取背景辐射数据,只需要根据仿真像元范围计算其覆盖及裁剪的背景格网,通过计算背景像元多边形的几何重心坐标并取整,即可得到背景像元的行列号,从而确定该像元的唯一标识P_ID(i,j),在完成仿真图像全部像元的计算后,得到整个仿真图像的P_ID(i,j)唯一值列表,从而可根据该列表从缓存的背景辐射栅格数据中得到相应的背景辐射数据值;此外在做逐像元仿真时,每一个仿真像元裁剪多个背景辐射像元,仿真像元值最终由多个像元混合而成,在进行混合计算时,使用背景辐射数据像元落入仿真像元的面积占其自身面积的比例做丰度,采用本实施例提供的方法,在像元坐标系下,一个完整像元其面积为1,对于边缘被裁剪的像元,其面积就是需要计算的面积比例,对于内部完整像元,其面积为1无需计算,因此本实施例提供的方法可省去计算面积比例的步骤,同时省去计算内部像元面积的操作,可快速获取每个像元的丰度值。
图7是本发明实施例提供的一种获取全球背景辐射数据及丰度进行像元混合计算的流程图。在进行图像仿真计算时,对于第i行第j列的仿真像元(i,j),根据上述S101-S104所述方法获取该像元覆盖的全部背景辐射数据像元的辐射数据值及丰度值,采用加权求和方法进行像元混合计算,得到仿真像元值,最终生成完整的仿真图像。
像元混合计算时,对于某一个仿真像元,取其对应的背景像元多边形列表L2中第i个多边形对应像元的行列号,从计算机内存中缓存的全球背景辐射数据中读取对应的背景辐射数据Vi,同时像元的面积值为丰度Ki,通过加权求和,即可计算得到仿真像元值V:
V=∑Ki*Vi
其中,i表示背景像元多边形列表L2中的第i个多边形,Ki为该多边形的丰度值,Vi为该多边形的背景辐射值。
图8是本发明实施例提供的一种全球背景辐射数据及丰度获取装置800的结构示意图。
全球背景辐射数据及丰度获取装置800包括:
几何计算模块801,用于根据仿真像元范围进行左边转换,计算仿真像元多边形与全球背景数据格网的交点,根据交点构造边缘裁剪多边形与内部完整多边形;
读取及缓存模块802,用于读取全球背景辐射数据并缓存在计算机内存中,计算机内存中的图像数据构成一个栅格图像矩阵,像元值可基于全球背景辐射数据像元行列号进行唯一值标识,当进行图像仿真计算时,可根据像元唯一值P_ID(i,j)标识作为索引取得对应的数据值。
像元混合计算模块803,用于对仿真像元裁剪的背景辐射像元进行加权求和,权值为像元面积比例,亦即丰度值,背景辐射像元值为对应像元的反射率或发射率。
图9是本发明实施例提供的一种计算机设备900的结构示意图,如图9所示,所述计算机设备900包括处理器901以及存储器902:
所述处理器901用于根据所述程序代码中的指令执行上述实施例提供的一种全球背景辐射数据及丰度的获取方法。
所述存储器902用于存储程序代码,并将所述程序代码传输给所述处理器;
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分可互相参考,对于实施例公开的系统或装置而言,由于其与实施例公开的方法相对应,相关之处参见方法部分说明即可。
本发明实施例所涉及的数据(包括但不限于用于分析的数据、存储的数据、展示的数据等),均为经用户授权或者经过各方充分授权的,且相关数据的收集、使用和处理需要遵守相关国家和地区的相关法律法规和标准。例如,本发明中涉及到的全球背景辐射数据都是在充分授权的情况下获取的。
上述实施例可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行计算机程序指令时,全部或部分地产生按照本发明实施例描述的流程或功能。计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。可用介质可以是磁性介质,(例如,软盘、硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘Solid State Disk(SSD))等。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (7)
1.一种全球背景辐射数据及丰度的获取方法,其特征在于,所述方法包括:
将仿真像元覆盖范围记录为多边形,将多边形坐标转换为像元坐标,按照顺时针方向逐点记录多边形各点;
遍历多边形各点,计算多边形与背景影像格网的交点,生成包含背景影像格网交点的多边形点集列表,所述背景影像格网为全球背景辐射数据的栅格行列号格网,所述栅格行列号格网基于栅格行列号构建,所述栅格行列号为像元唯一标识、作为获取数据的索引;
遍历包含背景影像格网交点的多边形点集列表,生成边缘多边形,计算每个边缘多边形的面积,并将这些边缘多边形加入背景像元多边形列表;
遍历内部完整像元,生成内部完整像元多边形,其中每个内部完整像元多边形的面积为1,将这些内部完整像元多边形加入背景像元多边形列表,得到完整的背景像元多边形列表;
遍历上述背景像元多边形列表,计算每一个多边形的几何重心坐标并取整,得到各多边形对应的背景像元的行列号;
读取全球背景辐射数据并缓存在内存中,基于上述背景像元多边形列表,使用其中各多边形对应的行列号从缓存的全球背景辐射数据中获取背景辐射数据值,使用其中各多边形的几何面积作为丰度值,计算丰度值与背景辐射数据值的乘积并求和,得到混合后的仿真像元值;
所述遍历包含背景影像格网交点的多边形点集列表,生成边缘多边形,包括:
从包含背景影像格网交点的多边形点集列表的第1个点开始,按顺时针方向依次构建边缘裁剪多边形,具体方法为:
按照顺时针方向,取所述包含背景影像格网交点的多边形点集列表的第2个点作为当前点,取第3个点作为下一点,根据下一点是否是背景影像格网交点,分不同情形构造完整边缘裁剪多边形;
如果所述下一点是多边形与背景格网的交点,则根据上一点与下一点的走向可以确定有限个相交情形,根据具体相交情形,按照顺时针方向依次添加像元的顶点坐标构造边缘裁剪多边形;
如果所述下一点不是多边形与背景格网的交点,则该点是原始多边形的顶点,此时取多边形点集列表的再下一个点作为下一点,同样可以根据前后两点的走向确定有限个相交情形,根据具体的相交情形,按顺时针方向依次添加像元的顶点坐标构造边缘裁剪多变形。
2.根据权利要求1所述的一种全球背景辐射数据及丰度的获取方法,其特征在于,所述遍历多边形各点,计算多边形与背景影像格网的交点,包括:
创建多边形点集列表,将多边形第1个点作为当前点加入列表;
按顺时针方向,将多边形第二个点作为下一个点,计算有向线段穿越的整数横坐标值;
对于每一个整数横坐标值,计算对应的纵坐标得到相应的格网竖线交点;
对于每一个格网竖线交点,计算前一点与该点之间穿越的整数纵坐标值;
对于每一个整数纵坐标值,计算对应的横坐标得到相应的格网横线交点;
重复以上操作,直到处理完毕全部多边形点,得到包含背景影像格网交点的多边形点集列表。
3.根据权利要求2所述的一种全球背景辐射数据及丰度的获取方法,其特征在于,所述遍历内部完整像元,生成内部完整像元多边形,方法包括:
对所述的边缘裁剪多边形,计算几何重心坐标,对几何重心坐标取整后得到对应像元的行列号;
使用行号为键,使用列号列表为值,构造键值对列表;
根据键值范围确定行号起止范围,对于每一行,根据列号列表中记录的列号最大值与最小值确定列号起止范围;
根据行号及列号起止范围,遍历内部完整像元,依次添加完整像元各顶点坐标,生成内部完整像元多边形。
4.根据权利要求3所述的一种全球背景辐射数据及丰度的获取方法,其特征在于,使用背景像元多边形列表中各多边形对应的行列号从缓存的全球背景辐射数据中获取背景辐射数据值,包括:
读取全球背景辐射数据,并在内存中进行缓存,每个像元的数据使用行列号做唯一标识;
计算每个背景像元多边形的几何重心坐标并取整,得到相应像元的行列号坐标;
使用所述行列号,从缓存的背景辐射数据中读取相应的背景辐射数据值。
5.根据权利要求4所述的一种全球背景辐射数据及丰度的获取方法,其特征在于,使用背景像元多边形列表中各多边形的几何面积作为丰度值,计算丰度值与背景辐射数据值的乘积并求和,得到混合后的仿真像元值,包括:
计算边缘裁剪像元对应的多边形及内部完整像元对应的多边形的面积,得到每个背景像元多边形的面积,将各背景像元多边形的面积作为丰度;
基于像元的行列号获取的背景辐射数据值与所述丰度值相乘并求和,得到混合后的仿真像元值。
6.根据权利要求1-5之一的所述一种全球背景辐射数据及丰度的获取方法的获取装置,其特征在于,包括:
几何计算模块,用于根据仿真像元范围,计算仿真像元多边形与全球背景数据格网的交点,根据交点构造边缘裁剪多边形与内部完整多边形;
读取及缓存模块,用于读取全球背景辐射数据并缓存在计算机内存中,计算机内存中的图像数据构成一个栅格图像矩阵,像元值基于全球背景辐射数据像元行列号设置唯一值标识,当进行图像仿真计算时,根据像元唯一值标识直接索引获取对应背景像元点的数据;
像元混合计算模块,用于仿真像元裁剪的全部背景辐射像元对应的像元值与丰度值进行加权求和,权值为像元面积比例,即丰度,背景辐射像元为对应像元的反射率或发射率。
7.一种计算机设备,其特征在于,所述计算机设备包括:处理器,所述处理器与存储器耦合,所述存储器中存储有至少一条计算机程序指令,所述至少一条计算机程序指令由所述处理器加载并执行,以使所述计算机设备实现权利要求1-5中任一项所述的一种全球背景辐射数据及丰度的获取方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310529370.3A CN116680871B (zh) | 2023-05-11 | 2023-05-11 | 一种全球背景辐射数据及丰度的获取方法、装置及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310529370.3A CN116680871B (zh) | 2023-05-11 | 2023-05-11 | 一种全球背景辐射数据及丰度的获取方法、装置及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116680871A CN116680871A (zh) | 2023-09-01 |
CN116680871B true CN116680871B (zh) | 2024-03-12 |
Family
ID=87788127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310529370.3A Active CN116680871B (zh) | 2023-05-11 | 2023-05-11 | 一种全球背景辐射数据及丰度的获取方法、装置及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116680871B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542035A (zh) * | 2011-12-20 | 2012-07-04 | 南京大学 | 基于扫描线法的多边形栅格化并行转换方法 |
CN104821013A (zh) * | 2015-05-11 | 2015-08-05 | 武汉大学 | 基于大地坐标系数字高程模型的地表面积提取方法及系统 |
CN106204606A (zh) * | 2016-07-15 | 2016-12-07 | 中国石油天然气集团公司 | 一种地理背景图像的裁剪方法及系统 |
CN110322556A (zh) * | 2019-04-29 | 2019-10-11 | 武汉大学 | 一种基于边界裁剪的高速高精度矢栅叠置分析方法 |
CN111563962A (zh) * | 2020-04-09 | 2020-08-21 | 中国科学院空天信息创新研究院 | 一种基于几何辐射一体化采样的遥感图像仿真方法 |
CN113568994A (zh) * | 2021-06-16 | 2021-10-29 | 数简科技(苏州)有限公司 | 一种基于栅格瓦片地图的真实数据获取方法 |
CN115841512A (zh) * | 2022-10-31 | 2023-03-24 | 太原理工大学 | 一种复杂地表真实面积的计算方法 |
-
2023
- 2023-05-11 CN CN202310529370.3A patent/CN116680871B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542035A (zh) * | 2011-12-20 | 2012-07-04 | 南京大学 | 基于扫描线法的多边形栅格化并行转换方法 |
CN104821013A (zh) * | 2015-05-11 | 2015-08-05 | 武汉大学 | 基于大地坐标系数字高程模型的地表面积提取方法及系统 |
CN106204606A (zh) * | 2016-07-15 | 2016-12-07 | 中国石油天然气集团公司 | 一种地理背景图像的裁剪方法及系统 |
CN110322556A (zh) * | 2019-04-29 | 2019-10-11 | 武汉大学 | 一种基于边界裁剪的高速高精度矢栅叠置分析方法 |
CN111563962A (zh) * | 2020-04-09 | 2020-08-21 | 中国科学院空天信息创新研究院 | 一种基于几何辐射一体化采样的遥感图像仿真方法 |
CN113568994A (zh) * | 2021-06-16 | 2021-10-29 | 数简科技(苏州)有限公司 | 一种基于栅格瓦片地图的真实数据获取方法 |
CN115841512A (zh) * | 2022-10-31 | 2023-03-24 | 太原理工大学 | 一种复杂地表真实面积的计算方法 |
Non-Patent Citations (1)
Title |
---|
DDA算法在地表面积计算中的应用研究;麻雪丽 等;测绘地理信息;第44卷(第6期);第79-82页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116680871A (zh) | 2023-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115359195B (zh) | 一种正射影像生成方法、装置、存储介质和电子设备 | |
CN112419460B (zh) | 烘焙模型贴图的方法、装置、计算机设备和存储介质 | |
CN111709876B (zh) | 一种图像拼接的方法及装置、设备、存储介质 | |
CN111340882A (zh) | 基于图像的室内定位方法及装置 | |
CN113724391A (zh) | 三维模型构建方法、装置、电子设备和计算机可读介质 | |
CN113963072B (zh) | 双目摄像头标定方法、装置、计算机设备和存储介质 | |
JP7432793B1 (ja) | 三次元点群に基づくマッピング方法、装置、チップ及びモジュール機器 | |
CN116680871B (zh) | 一种全球背景辐射数据及丰度的获取方法、装置及设备 | |
CN116385622B (zh) | 云图处理方法、装置、计算机及可读存储介质 | |
CN115082322B (zh) | 图像处理方法和装置、图像重建模型的训练方法和装置 | |
CN114118367B (zh) | 增量式神经辐射场构建的方法及设备 | |
CN115760641A (zh) | 基于多尺度特征注意力网络的遥感影像去云雾方法和设备 | |
CN116503474A (zh) | 位姿获取方法、装置、电子设备、存储介质及程序产品 | |
CN115393423A (zh) | 目标检测方法和装置 | |
CN114519679A (zh) | 一种sar目标图像数据智能增强方法 | |
CN114549650A (zh) | 相机标定方法、装置、电子设备及可读存储介质 | |
CN116503464B (zh) | 基于遥感图像的农田建筑物高度预测方法 | |
CN116295031B (zh) | 弧垂的测量方法、装置、计算机设备和存储介质 | |
RU2538319C1 (ru) | Устройство поиска дубликатов изображений | |
CN117557711B (zh) | 可视域的确定方法、装置、计算机设备、存储介质 | |
WO2024087137A1 (zh) | 一种渲染方法及相关装置 | |
CN115953544B (zh) | 一种三维重建方法、装置、电子设备及可读存储介质 | |
CN118097566B (zh) | 基于深度学习的场景变动检测方法、装置、介质及设备 | |
CN116468611B (zh) | 图像拼接方法、装置、设备和存储介质 | |
CN118247149A (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 |