一种基于gpu的线阵遥感影像dsm快速提取方法
【专利摘要】本发明公开了一种基于GPU的线阵遥感影像DSM快速提取方法取方法和装置，该方法包括：获取待处理的遥感影像以及遥感影像对应的RPC参数信息；对遥感影像进行Wallis滤波处理；利用GPU以及RPC参数信息对滤波处理后的遥感影像进行核线校正，并生成遥感影像的左核线影像和右核线影像；根据左核线影像和右核线影像的第一同名像点计算出左右核线影像的视差值范围；根据第一同名像点从视差值范围中得到左核线影像和右核线影像的视差值，并根据得到的视差值计算左右核线影像的第二同名像点，并根据该第二同名像点得到原始影像的同名像点，并根据原始影像的同名像点和RPC参数信息，进行基于RPC模型的前方交会生成遥感影像DSM，以解决从线阵遥感影像中获取DSM效率低和精度低的问题。
【专利说明】
-种基于GPU的线阵遥感影像DSM快速提取方法
技术领域
[0001] 本发明设及遥感影像信息提取技术领域，尤其设及一种基于GPU的线阵遥感影像 DSM快速提取方法取方法。
【背景技术】
[0002] 数字表面模型(Digital Surface Model，DSM)是地图产品生产过程中，生成等高 线、真正射影像制作、建筑物提取与重建、数字高程模型(Digital Elevation Model,DEM) 提取、地理信息更新等工作的重要信息源。DSM在城市变化监测、±地管理、交通管理、矿产 资源调查、文物保护等民用领域W及作战指挥自动化控制系统和精确武器制导等军事领域 中都具有重要作用。
[0003] 利用摄影测量的方法获取DSM其本质在于寻找立体影像中的同名像点，核屯、技术 是影像的全自动匹配。然而，现有利用现有的影像匹配技术，从影像中获取的DSM密度与精 度较低，处理效率较低，DSM粗差剔除等后续生产过程需要大量的人工编辑工作，造成DSM获 取的自动化水平较低。
【发明内容】
[0004] 本发明实施例提供一种基于GPU的线阵遥感影像DSM快速提取方法取方法和装置， W至少解决现有从线阵遥感影像中获取的DSM效率低、精度低的技术问题。
[0005] 为达到W上目的，本发明的技术方案为：
[0006] -方面，提供一种基于GPU的线阵遥感影像DSM快速提取方法取方法，包括：
[0007] 获取待处理的遥感影像W及所述遥感影像所对应的有理多项式系数RPC参数信 息;对所述遥感影像进行Wa 11 i S滤波处理;利用图形处理器GPU W及所述RPC参数信息，对滤 波处理后的遥感影像利用投影轨迹法进行核线校正，并通过GPU并行加速方法生成所述遥 感影像的左核线影像和右核线影像;采用基于GPU的SIFT匹配算法快速提取所述左核线影 像和所述右核线影像的第一同名像点，进而计算出左右核线影像的视差值范围；利用SIMD+ 化enMP算法对所述第一同名像点进行并行快速逐像素半全局密集匹配分别从所述视差值 范围中得到所述左核线影像和所述右核线影像的视差值;根据得到的视差值计算所述左核 线影像和所述右核线影像的第二同名像点，并由所述左核线影像和所述右核线影像与原始 影像的关系反算得到原始影像的同名像点；根据所述原始影像的同名像点和所述RPC参数 信息，利用GPU进行基于RP对莫型的前方交会生成遥感影像数字表面模型DSM。
[000引另一方面，提供一种基于GPU的线阵遥感影像DSM快速提取方法取装置，包括:获取 单元，用于获取待处理的遥感影像W及所述遥感影像所对应的有理多项式系数RPC参数信 息；第一处理单元，用于对所述遥感影像进行Wal 1 i S滤波处理;第二处理单元，用于利用图 形处理器GPUW及所述RPC参数信息，对滤波处理后的遥感影像利用投影轨迹法进行核线校 正，并通过GPU并行加速方法生成所述遥感影像的左核线影像和右核线影像；第一计算单 元，用于采用基于GPU的SIFT匹配算法快速提取所述左核线影像和所述右核线影像的第一 同名像点，进而计算出左右核线影像的视差值范围；匹配单元，用于利用SIMD+OpenMP算法 对所述第一同名像点进行并行快速逐像素半全局密集匹配分别从所述视差值范围中得到 所述左核线影像和所述右核线影像的视差值;第二计算单元，用于根据得到的视差值计算 所述左核线影像和所述右核线影像的第二同名像点，并由所述左核线影像和所述右核线影 像与原始影像的关系反算得到原始影像的同名像点；生成单元，用于根据所述原始影像的 同名像点和所述RPC参数信息，利用GPU进行基于RPC模型的前方交会生成遥感影像数字表 面模型DSM。
[0009] 通过上述方法，采用逐像素匹配的策略，利用遥感影像的半全局密集匹配方法，实 现了遥感影像高精度密集点匹配，获取了足够密集的匹配点，无需内插即可直接前方交会 获取密集的DSM信息，从而快速获取高精度的DSM。
【附图说明】
[0010] 图1是本发明实施例提供的一种可选地基于GPU的线阵遥感影像DSM快速提取方法 的流程示意图；
[0011] 图2是本发明实施例提供的一种可选地利用GPU快速实现基于RP对莫型核线影像生 成的流程示意图；
[0012] 图3是本发明实施例提供的一种可选地利用基于GPU的SIFT匹配算法获取同名点 对的流程示意图；
[OOK]图4是本发明实施例提供的一种可选地基于GRJ的线阵遥感影像DSM快速提取装置 的结构示意图；
[0014] 图5是本发明实施例提供的另一种可选地基于GPU的线阵遥感影像DSM快速提取装 置的结构示意图；
[0015] 图6是本发明实施例提供的第S种可选地基于GPU的线阵遥感影像DSM快速提取装 置的结构示意图；
[0016] 图7是本发明实施例提供的第四种可选地基于GPU的线阵遥感影像DSM快速提取 装置的结构示意图。
【具体实施方式】
[0017] 下面结合附图1和具体实施例对本发明进行详细说明。本实施例W本发明技术方 案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限 于下属的实施例。
[001引如图1所示，本发明实施例提供了一种基于GPlKGra地ics Processing化it,图形 处理器)的线阵遥感影像DSM快速提取方法，包括：
[0019] SlOl、获取待处理的遥感影像W及该遥感影像所对应的RPC(rational polynomial coefficient,有理多项式系数)参数信息。
[0020] S102、对该遥感影像进行Wal 1 is滤波处理。
[0021] 其中，通过GPU中Wallis滤波对遥感影像进行预处理，增强遥感影像对比度并且降 低噪声;本步骤可W通过编写Wallis滤波算子核函数实现。
[0022] S103、利用图形处理器GPUW及该RPC参数信息，对滤波处理后的遥感影像利用投 影轨迹法进行核线校正，并通过GPU并行加速方法生成该遥感影像的左核线影像和右核线 影像。
[0023] 在本实施例一种可能的实现方式中，如图2所示，将预处理后的遥感影像从主机端 内存拷贝至GPU的全局存储器中，将RPC参数拷贝至GPU的共享存储器中，根据计算机性能参 数，动态调整片上缓存;设置GPU线程块大小，确定GPU线程数目；在GPU上利用核线参数答解 算子求解核线参数;将求解得到的核线参数从GPU全局存储器中拷贝回主机端内存。
[0024] 其中W生成线阵影像的左像（即左核线影像)和右像(即右核线影像)对为例，基于 RP对莫型进行核线校正的实施过程如下：
[0025] Sl031、从左像的第1列开始，选择列方向的中央像点m(0，H/2 )，利用投影轨迹法 将m和左像投影中屯、S确定的光线投影到右像上，考虑运算量，投影该光线上高程为hmax和 hmin的两个物方点化和化到右像上，确定两个像点qi和Q2 .
[00%] S1032、判断qi和Q2是否在右像的像幅范围内：如果不在则返回S1031，继续从下一 列开始；如果在，则根据该两点qi和Q2确定一条直线1'，记录该直线方程参数a ' x+b ' y+c ' = 0。
[0027] S1033、确定右像上qi和Q2的中点Q3,利用投影轨迹法将Q3和右像投影中屯、S'和其 确定的光线投影到右像上，考虑运算量，投影该光线上高程为hmax和hmin的两个物方点化和化 至Ij左像，确定两个像点Q3和Q4 ;
[0028] S1034、判断Q3和Q4是否在左像的像幅范围内：如果不在则返回S1031，继续从下一 列开始;如果在，则根据该两点Q3和Q4确定一条直线1，记录直线方程ax+by+c = 0的参数。
[0029] S1035、利用左像上的直线1，在左像内沿直线进行重采样，具体是对给定的从0到H 图像的高度变化的系列yi坐标，根据直线方程计算相应的系列Xi坐标，根据系列（xi，yi)坐 标从原始影像通过双线性内插确定W系列的像素灰度值，所得灰度值即为左像核线影像上 的第一行。
[0030] S1036、根据直线1'确定右像核线影像的第一行，返回步骤S1031重复进行，直到确 定核线影像上的最后一行。
[0031] S1037、对左像，记录每一条直线方程的参数，核线有多少行，就有多少组方程参 数，将运些直线参数写入一个文本参数文件，利用该参数文件完成核线影像坐标到原始影 像坐标的转化;对右像执行同样的处理，也形成另一个参数文件，核线影像与原始影像的对 应关系为：
[0032]
[0033] 其中，X和y分别表示原始影像的列方向、行方向的像素坐标值;X'和/表示核线影 像的像素坐标值，分别为列数和行数;hei曲t表示原始影像的高度;邱iwi化h表示核线影像 的宽度，n表示核线方程的个数，i表示直线的编号，ai、bi、ci表示直线方程的参数。
[0034] 对于核线影像的一个点，根据核线参数可求得其在原始影像上对应的坐标，由于 坐标不是整数，采用双线性内插方法进行处理;该过程通过编写基于RPC模型的核线参数答 解算子核函数实现GPU并行，设置计算每个block中thread的维度为（16，16) ,block的维度 为(wi化h/16+l，height/16+l)，每个线程处理I个像素。
[OO%] S104、采用基于GPU的SIFT(Scale-inva;riant fea1:ure transform,尺度不变特征 变换）匹配算法快速提取该左核线影像和该右核线影像的第一同名像点，进而计算出左右 核线影像的视差值范围。
[0036] 其中，如图3所示，为了加快后续步骤图像匹配速度，缩小视差捜索范围，同时减少 对计算机内存的使用大小，可W获取视差值范围，从而在后续计算中从该视差值范围中确 定视差值，具体地，可W采用基于GPU的SIFT匹配算法快速提取左右核线影像的第一同名像 点，剔除粗差后保留若干可靠的同名点。
[0037] S105、利用SIMD(Single Instruction Multiple Data,单指令多数据流）+0penMP 算法对该左核线影像和该右核线影像的第一同名像点进行并行快速逐像素半全局密集匹 配分别从该视差值范围中得到该左核线影像和该右核线影像的视差值。
[0038] 其中，本步骤可W包括W下步骤：
[0039] S1、对该第一同名像点进行Census变换得到该同名像点的匹配代价，并根据该匹 配代价建立视差空间影像DSI。
[0040] 具体地，Census变换W待计算的第一同名像点为中屯、开辟一个矩形窗口，将窗口 中除中屯、像素 W外的其他像素变换为一个比特串；若窗口中一个像素的灰度值比中屯、像素 的灰度值大，则相应位置为1，反之，置为0;相应的Census变换C(P)用比特串的形式表示为：
[0041]
[0042] 参数Ck的值为：
[0043]
[0044] 其中，I(P)为图像P点的像素灰度值，I(q)为图像q点的像素灰度值，P和q分别表示 图像像素位置。函数T[]为化ttes模型，当P点的像素灰度值大于q点的像素灰度值时，返回 1，否则返回0。计算Census变换后参考影像与捜索影像对应窗口之间比特串的Hamming距离 作为任意一个像素点P的匹配代价，即：
[0045] 0(口，(1)=紅3?山8(口山其中，(1为视差；
[0046] 计算完匹配代价后，建立DSI。
[0047] S2、在该DSI内采用半全局密集匹配对该第一同名像点的匹配代价进行匹配代价 聚合，求解最优值处理得到该第一同名像点的视差值。
[004引具体地，在DSI内用多个一维方向的代价聚合代替二维方式捜索，任意一个像素点 P在视差d匹配代价聚合可W通过像素点P点所有一维方向代价的累加求和来计算，即：
[0049]
[0050] 其中，Lr(p，d)是待匹配像素点P在视差d处r方向上的匹配代价，像素点P的视差d =mindS(p,d) 〇
[0051 ]运样，在本步骤中，在分析内存空间和路径数量优化的基础上，针对遥感影像数据 量较大的特点，对影像进行分块生成视差图处理;分析Census变换运算特点，采用化enMP并 行实现Census变换;采用SIMD并行实现匹配聚合W及视差值计算。
[0052] 可选地，在得到视差值后，根据分别得到的该左核线影像和该右核线影像的视差 值进行视差一致性检查，W提高匹配的可靠性，抑制由于遮挡造成的误匹配。
[0053] 其中，视差一致性检查又名双向匹配，将"左右影像匹配得到的视差值作为右左影 像匹配的视差初值，从而提高左右视差一致性检查的运行速度。
[0054] 进一步地，在该视差一致性检查通过后，通过中值滤波器对该视差值进行中值滤 波。
[0055] 其中，该中值滤波中滤波器的尺寸可W为7X7的窗口，利用非线性滤波抑制因噪 声造成的孤立视差值点，在原始影像分辨率匹配时，采用抛物线拟合的方法得到亚像素视 差，从而得到精化后的视差信息。
[0056] S106、根据得到的视差值计算该左核线影像和该右核线影像的第二同名像点，并 由该左核线影像和该右核线影像与原始影像的关系反算得到原始影像的同名像点。
[0057] 其中，W左核线影像为基准，左核线影像点P对应的右核线影像点q的计算公式为 q = P+d，逐像素计算即可得到左核线影像对应的右核线影像的第二同名像点；核线影像与 原始影像的关系参考核线影像生成过程，具体参见步骤S103,利用GPU加速方法计算出左右 原始影像的同名像点。
[005引S107、根据原始影像的同名像点和该RPC参数信息，利用GPU进行基于RPC模型的前 方交会生成遥感影像数字表面模型DSM。
[0059] 其中，用GPU并行运算基于RPC模型的前方交会的核函数，将数据从内存传输到显 卡存储器，调用基于RPC模型的前方交会算子，再将结果从显卡存储器写回内存。核函数将 GPU存储器上的数据分割成块，每块被分到一个block线程中，每块的大小和block中的线程 数对应;一个block中的线程按照接合的特点从GPU存储器中读取数据，然后每个线程执行 一个前方交会操作；同一个block中的不同线程执行的模板操作指令时一样的，让一个 block中的不同线程共享指令前段;在调用前方交会算子计算结束后，计算结果会写回GPU 存储器中，然后再从GPU存储器写回内存的缓冲区中。
[0060] 可选地，在用GPU并行运算基于RPC模型的前方交会的核函数，生成遥感影像DSM之 后，基于GPU对生成的该遥感影像DSM进行移动曲面滤波，获取滤波后的遥感影像DSM。
[0061] 其中，为求取某个=维点处的拟合高程时，该点所在趋势面模型参数由距该点一 定范围内的=维点按最小二乘法求解，即每个=维点处的趋势面都需要根据附近点重新求 解，W保证趋势面的参数与实际地形较好地符合。
[0062] 运样，通过采用逐像素匹配的策略，利用遥感影像的半全局密集匹配方法，实现了 遥感影像高精度密集点匹配，获取了足够密集的匹配点，无需内插即可直接前方交会获取 密集的DSM信息，从而获取高精度的DSM。另外，通过利用GPU大规模并行计算的特点，W及 GPU上任务分解、大规模计算核屯、的分解方法，结合使用了共享存储器、全局存储器对算法 进行加速，实现了快速Wallis图像滤波算法、快速基于RP对莫型的核线影像生成W及快速基 于RPC模型的前方交会。解决了现有从线阵遥感影像中获取的DSM效率低、精度低的技术问 题。
[0063] 图4为本发明实施例提供的一种基于GHJ的线阵遥感影像DSM快速提取方法取装 置，包括：
[0064] 获取单元401，用于获取待处理的遥感影像W及该遥感影像所对应的有理多项式 系数RPC参数信息；
[0(?日]第一处理单元402，用于对该遥感影像进行Wal 1 is滤波处理；
[0066] 第二处理单元403,用于利用图形处理器GPUW及该RPC参数信息，对滤波处理后的 遥感影像利用投影轨迹法进行核线校正，并通过GPU并行加速方法生成该遥感影像的左核 线影像和右核线影像；
[0067] 第一计算单元404,用于采用基于GPU的SIFT匹配算法快速提取该左核线影像和该 右核线影像的第一同名像点，进而计算出左右核线影像的视差值范围；
[0068] 匹配单元405,用于利用SIMD+OpenMP算法对该第一同名像点进行并行快速逐像素 半全局密集匹配分别从该视差值范围中得到该左核线影像和该右核线影像的视差值；
[0069] 第二计算单元406,用于根据得到的视差值计算该左核线影像和该右核线影像的 第二同名像点，并由该左核线影像和该右核线影像与原始影像的关系反算得到原始影像的 同名像点；
[0070] 生成单元407,用于根据得到的原始影像的同名像点利用GPU进行基于RPC模型的 前方交会生成遥感影像DSM。
[0071] 可选地，如图5所示，该装置还包括：
[0072] 检查单元408,用于根据分别得到的该左核线影像和该右核线影像的视差值进行 视差一致性检查。
[0073] 可选地，如图6所示，该装置还包括：
[0074] 第一滤波单元409,用于在该视差一致性检查通过后，通过中值滤波器对该视差值 进行中值滤波。
[0075] 可选地，该匹配单元405,具体用于对该同名像点进行Census变换得到该第一同名 像点的匹配代价，并根据该匹配代价建立视差空间影像DSI;在该DSI内采用半全局密集匹 配对该第一同名像点的匹配代价进行匹配代价聚合，求解最优值处理得到该第一同名像点 的视差值。
[0076] 可选地，如图7所示，该装置还包括：
[0077] 第二滤波单元4010,用于基于GPU对生成的该遥感影像DSM进行移动曲面滤波，获 取滤波后的遥感影像DSM。
[0078] 采用上述装置，通过采用逐像素匹配的策略，利用遥感影像的半全局密集匹配方 法，实现了遥感影像高精度密集点匹配，获取了足够密集的匹配点，无需内插即可直接前方 交会获取密集的DSM信息，从而获取高精度的DSM。另外，通过利用GPU大规模并行计算的特 点，W及GPU上任务分解、大规模计算核屯、的分解方法，结合使用了共享存储器、全局存储器 对算法进行加速，实现了快速Wallis图像滤波算法、快速基于RP对莫型的核线影像生成W及 快速基于RPC模型的前方交会。解决了现有从线阵遥感影像中获取的DSM效率低、精度低的 技术问题。
[0079] 上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
[0080] 在本发明的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有 详述的部分，可W参见其他实施例的相关描述。
[0081] 在本申请所提供的几个实施例中，应该理解到，所掲露的技术内容，可通过其它的 方式实现。其中，W上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，可W为 一种逻辑功能划分，实际实现时可W有另外的划分方式，例如多个单元或组件可W结合或 者可W集成到另一个系统，或一些特征可W忽略，或不执行。另一点，所显示或讨论的相互 之间的禪合或直接禪合或通信连接可W是通过一些接口，单元或模块的间接禪合或通信连 接，可W是电性或其它的形式。
[0082] 所述作为分离部件说明的单元可W是或者也可W不是物理上分开的，作为单元显 示的部件可W是或者也可W不是物理单元，即可W位于一个地方，或者也可W分布到多个 单元上。可W根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0083] 另外，在本发明各个实施例中的各功能单元可W集成在一个处理单元中，也可W 是各个单元单独物理存在，也可W两个或两个W上单元集成在一个单元中。上述集成的单 元既可W采用硬件的形式实现，也可W采用软件功能单元的形式实现。
[0084] 所述集成的单元如果W软件功能单元的形式实现并作为独立的产品销售或使用 时，可W存储在一个计算机可读取存储介质中。基于运样的理解，本发明的技术方案本质上 或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可WW软件产品的形式 体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用W使得一台计算机 设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或 部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memo巧）、随机存取存 储器(RAM,Random Access Memo巧）、移动硬盘、磁碟或者光盘等各种可W存储程序代码的 介质。
[0085] W上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可W做出若干改进和润饰，运些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 一种基于GPU的线阵遥感影像DSM快速提取方法取方法，其特征在于，包括： 获取待处理的遥感影像以及所述遥感影像所对应的有理多项式系数RPC参数信息； 对所述遥感影像进行Wal 1 i s滤波处理； 利用图形处理器GPU以及所述RPC参数信息，对滤波处理后的遥感影像利用投影轨迹法 进行核线校正，并通过GHJ并行加速方法生成所述遥感影像的左核线影像和右核线影像； 采用基于GPU的SIFT匹配算法快速提取所述左核线影像和所述右核线影像的第一同名 像点，进而计算出左右核线影像的视差值范围； 利用SMD+OpenMP算法对所述第一同名像点进行并行快速逐像素半全局密集匹配分别 从所述视差值范围中得到所述左核线影像和所述右核线影像的视差值； 根据得到的视差值计算所述左核线影像和所述右核线影像的第二同名像点，并由所述 左核线影像和所述右核线影像与原始影像的关系反算得到原始影像的同名像点； 根据所述原始影像的同名像点和所述RPC参数信息，利用GPU进行基于RPC模型的前方 交会生成遥感影像数字表面模型DSM。2. 根据权利要求1所述的方法，其特征在于，在所述利用SMD+OpenMP算法对所述第一 同名像点进行并行快速逐像素半全局密集匹配分别从所述视差值范围中得到所述左核线 影像和所述右核线影像的视差值之后，所述方法还包括： 根据分别得到的所述左核线影像和所述右核线影像的视差值进行视差一致性检查。3. 根据权利要求2所述的方法，其特征在于，在所述根据所述原始影像的同名像点和所 述RPC参数信息，利用GHJ进行基于RPC模型的前方交会生成遥感影像数字表面模型DSM之 前，所述方法还包括： 在所述视差一致性检查通过后，通过中值滤波器对所述视差值进行中值滤波。4. 根据权利要求1所述的方法，其特征在于，所述利用SMD+OpenMP算法对所述第一同 名像点进行并行快速逐像素半全局密集匹配分别从所述视差值范围中得到所述左核线影 像和所述右核线影像的视差值包括： 对所述第一同名像点进行Census变换得到所述同名像点的匹配代价，并根据所述匹配 代价建立视差空间影像DSI; 在所述DSI内采用半全局密集匹配对所述第一同名像点的匹配代价进行匹配代价聚 合，求解最优值处理得到所述第一同名像点的视差值。5. 根据权利要求1所述的方法，其特征在于，在用GPU并行运算基于RPC模型的前方交会 的核函数，生成遥感影像DSM之后，所述方法还包括： 基于GPU对生成的所述遥感影像DSM进行移动曲面滤波，获取滤波后的遥感影像DSM。6. -种基于GPU的线阵遥感影像DSM快速提取方法取装置，其特征在于，包括： 获取单元，用于获取待处理的遥感影像以及所述遥感影像所对应的有理多项式系数 RPC参数信息； 第一处理单元，用于对所述遥感影像进行Wal 1 i s滤波处理； 第二处理单元，用于利用图形处理器GPU以及所述RPC参数信息，对滤波处理后的遥感 影像利用投影轨迹法进行核线校正，并通过GPU并行加速方法生成所述遥感影像的左核线 影像和右核线影像； 第一计算单元，用于采用基于GPU的SIFT匹配算法快速提取所述左核线影像和所述右 核线影像的第一同名像点，进而计算出左右核线影像的视差值范围； 匹配单元，用于利用SMD+OpenMP算法对所述第一同名像点进行并行快速逐像素半全 局密集匹配分别从所述视差值范围中得到所述左核线影像和所述右核线影像的视差值； 第二计算单元，用于根据得到的视差值计算所述左核线影像和所述右核线影像的第二 同名像点，并由所述左核线影像和所述右核线影像与原始影像的关系反算得到原始影像的 同名像点； 生成单元，用于根据所述原始影像的同名像点和所述RPC参数信息，利用GPU进行基于 RPC模型的前方交会生成遥感影像数字表面模型DSM。7. 根据权利要求6所述的装置，其特征在于，所述装置还包括： 检查单元，用于根据分别得到的所述左核线影像和所述右核线影像的视差值进行视差 一致性检查。8. 根据权利要求7所述的装置，其特征在于，所述装置还包括： 第一滤波单元，用于在所述视差一致性检查通过后，通过中值滤波器对所述视差值进 行中值滤波。9. 根据权利要求6所述的装置，其特征在于，所述匹配单元，具体用于对所述第一同名 像点进行Census变换得到所述同名像点的匹配代价，并根据所述匹配代价建立视差空间影 像DSI;在所述DSI内采用半全局密集匹配对所述第一同名像点的匹配代价进行匹配代价聚 合，求解最优值处理得到所述第一同名像点的视差值。10. 根据权利要求6所述的装置，其特征在于，所述装置还包括： 第二滤波单元，用于基于GPU对生成的所述遥感影像DSM进行移动曲面滤波，获取滤波 后的遥感影像DSM。
【文档编号】G06K9/00GK106023147SQ201610298084
【公开日】2016年10月12日
【申请日】2016年5月6日
【发明人】折晓宇, 范晓敏, 巩丹超, 陈巍然, 靳笑琳
【申请人】西安航天天绘数据技术有限公司