买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：武汉大学

摘要：本发明涉及一种基于空间几何比值关系的SAR影像叠掩区域定位方法。首先获取SAR数据及其对应区域的DEM数据；根据SAR影像附带的xml文件，获取影像入射角；读取RPC参数文件，构星载SAR影像的有理函数模型RFM；确定SAR影像四角像点坐标对应的大地坐标，并获取大地坐标组成的最小外接矩形；根据最小外接矩形裁剪DEM数据，并对其重采样；确定叠掩发生时物方像素与像方像素数量的比值关系阈值；在对重采样DEM数据物方分块的基础上，结合多线程并行的方式计算每个物方子块的掩膜，合并掩膜并进行后处理生成最终的叠掩掩膜。本发明提出的方法可用于快速生成SAR叠掩盖掩膜，为SAR影像的进一步应用提供技术保障。

主权项：1.一种基于空间几何比值关系的SAR影像叠掩区域定位方法，其特征在于，包括以下步骤：步骤1，获取SAR原始影像，根据影像中心点范围预估所在位置，获取覆盖影像范围对应区域的DEM数据；步骤2，解析SAR影像附带的xml文件，获取影像入射角；步骤3，读取RPC参数文件，构建星载SAR影像的有理函数模型；步骤4，计算SAR影像四个角点的像点坐标，将四个角点的像点坐标换算成大地坐标，并根据大地坐标计算对应最小外接矩形；步骤5，根据步骤4得到的最小外接矩形，裁剪步骤1中获取的DEM数据，并对裁剪出的DEM数据重采样到指定分辨率；步骤6，根据SAR叠掩成像几何构建地距与斜距的对应关系，并进行空间几何转换从而得到叠掩发生时物方像素与像方像素的比值关系阈值；地距Lground在X方向长度满足下式： 式中，LX表示Lground在X方向的投影长度，θ是距离向与X方向的夹角；考虑到斜距与地距存在如下关系： 式中，Lslant表示斜距；Lground为地距；表示视角，与入射角相等；地距在X方向的投影长度LX由DEM重采样后的分辨率Resg与该长度下对应像素个数CNTX乘积表示：LX＝Resg×CNTX10斜距由斜距分辨率Resg与其对应方向的像素个数CNTs进行表示：Lslant＝Ress×CNTs11式中，CNTs表示SAR斜距向像素个数，即SAR影像宽度W，则式11可以进一步表示为：Lslant＝Ress×W12联立式8、式9式10、式12可得： 当地距与斜距的比值不小于阈值TH时，将地距范围判定为叠掩，即： 实际计算中，地距不单独计算而是通过地距在X方向的投影进行间接表示，将式14按照分辨率与像素个数乘积的方式进行换算得： 式中，ng×Resg为地距在X方向的投影，表示地距，ng表示DEM位于X方向的像素个数，ns表示斜距方向的像素个数；由15得，当叠掩发生时，ng与ns的比值关系满足以下条件： 由式16可知，当在重采样的DEM上有不低于THfinal个像素对应的大地坐标X,Y,Z利用定位模型反算的像点坐标R,C处于同一个像方像素格网时，大地坐标对应的位置判定为叠掩；步骤7，根据步骤6得到的比值关系阈值生成叠掩区域的掩膜，实现叠掩区域的快速定位；步骤7.1，对步骤5重采样后DEM数据按照指定大小进行分块，获取DEM子块；步骤7.2，对DEM子块逐像素获取对应的大地坐标X,Y,Z，并利用步骤3中有理函数模型计算位于SAR影像的像点坐标R,C，根据像点坐标最大值、最小值获取最小外接矩形，在最小外接矩形基础上进行边界扩展，获得像方分块；步骤7.3，在获取DEM子块和其对应的像方子块后，根据步骤6计算得到的阈值THfinal生成子块对应的掩膜；步骤7.4，对步骤7.3得到的子块掩膜进行合并，生成与重采样DEM相同大小的掩膜，得到SAR影像对应的物方初始掩膜；步骤7.5，对初始掩膜利用形态学“闭”运算，消除部分“孔”、“洞”，得到最终的SAR影像叠掩区域掩膜。

全文数据：

权利要求：

百度查询： 武汉大学 一种基于空间几何比值关系的SAR影像叠掩区域定位方法