System.ArgumentOutOfRangeException: 索引和长度必须引用该字符串内的位置。 参数名: length 在 System.String.Substring(Int32 startIndex, Int32 length) 在 zhuanliShow.Bind()
【技术实现步骤摘要】
本专利技术属于微波遥感,尤其涉及一种ka-sar卫星标准产品生成方法及计算机可读介质。
技术介绍
1、随着卫星规模化应用热潮兴起,微小型sar卫星迅速发展,为在短时间内获得重复观测提供了新的途径。ka-sar卫星可以实现应用于海洋环境监测和洪水探测等方面,具备广泛应用前景。但当前的ka-sar卫星主要是以微小卫星为主受限于卫星尺寸,与标准卫星相比,微小卫星的几何性能相对较差,限制了其在摄影测量领域的应用。传统方法可以通过增加分布良好的地面控制点来提高卫星图像的几何定位精度,但需要大量财力及人力资源。为此,提出一种ka-sar卫星标准产品生成方法。结合标准sar卫星(如高分三号)高几何定位精度优点,通过与标准sar卫星生成的dom匹配高精度控制点而实现ka-sar影像定位精度提升,进而生成高精度ka-sar卫星标准产品,为其在水利等行业精细化应用提供支撑。本文方法可为多源sar卫星影像数据协同处理提供一种新思路。
技术实现思路
1、针对多源sar卫星初始定位精度差的问题,本专利技术提出了一种ka-sar卫星标准产品生成方法及计算机可读介质。
2、本专利技术方法的技术方案为一种ka-sar卫星标准产品生成方法,具体如下:
3、进行强度图转换,生成对应的强度图影像;
4、对强度图影像进行子块划分;
5、根据强度图影像信息,计算几何纠正的输出分辨率;
6、对强度图影像进行局部几何纠正得到强度图影像的多个几何纠正子块,根据每个子块范
7、采用异源匹配算法获取多组同名点对位于dom重采样子块地理坐标与强度图影像几何纠正子块地理坐标;依次映射得到dom上每个dom重采样子块对应的控制点地理三维坐标、强度图影像上每个几何纠正子块地理坐标对应的像素坐标;
8、采用区域网平差进行粗差剔除,得到精确的dom上控制点地理三维坐标、强度图影像控制点的像素坐标;
9、采用精确的dom上控制点地理坐标、强度图影像控制点的像素坐标对rfm模型进行更新,得到更新的rfm模型;
10、迭代更新得到最终的rfm模型;采用最终的rfm模型,结合辅助dem数据对强度图影像进行间接校正,得到ka-sar卫星影像的标准产品。
11、本专利技术的具体步骤如下:
12、步骤1:获取原始ka-sar单视复数影像、辅助的有理函数模型参数文件、c波段的高分三号dom、数字高程模型;
13、步骤2:对原始ka-sar单视复数影像进行强度图转换,生成对应的强度图影像;
14、步骤3:获取强度图影像的像素宽、像素高,根据像素宽、像素高对强度图影像进行子块划分;
15、步骤4:根据强度图影像信息,计算几何纠正的输出分辨率;
16、步骤5:根据每个地理子块,结合dem辅助数据,对强度图影像进行局部几何纠正得到强度图影像的多个几何纠正子块,根据每个子块范围,在dom上重采样对应区域的像素值,得到dom的多个重采样子块;
17、步骤6:将每个dom重采样子块作为参考影像、强度图影像的几何纠正子块作为配准影像,采用异源匹配算法获取多组同名点对位于dom重采样子块地理坐标与强度图影像几何纠正子块地理坐标;将dom重采样子块地理坐标结合辅助dem映射成dom上每个dom重采样子块对应的控制点地理三维坐标;将强度图影像的每个几何纠正子块地理坐标映射成强度图影像上每个几何纠正子块地理坐标对应的像素坐标;
18、步骤7:将dom上每个dom重采样子块对应的控制点地理三维坐标与强度图影像上每个几何纠正子块地理坐标对应的像素坐标采用区域网平差进行粗差剔除,得到精确的dom上控制点地理三维坐标、强度图影像控制点的像素坐标;
19、步骤8:采用精确的dom上控制点地理坐标、强度图影像控制点的像素坐标对rfm模型进行更新,得到更新的rfm模型;
20、步骤9:将强度图影像的rfm模型采用步骤8中更新后的精确的rfm模型进行替代,并循环步骤3—步骤8,得到最终的rfm模型。
21、步骤10:采用最终的rfm模型,结合辅助dem数据对强度图影像进行间接校正,得到ka-sar卫星影像的标准产品。
22、作为优选,步骤2所述生成对应的强度图影像,具体如下:
23、所述ka-sar单视复数影像定义为
24、所述强度图影像定义为iintensity;
25、第p行、第q列强度图对应的像素值计算如下:
26、
27、其中,
28、表示ka-sar单视复数影像第1波段第p行、第q列对应的像素值;
29、表示ka-sar单视复数影像第2波段第p行、第q列对应的像素值;
30、iintensity(p,q)表示第p行、第q列强度图影像对应的像素值;
31、作为优选,步骤3所述对强度图影像进行子块划分,具体过程如下:
32、设定列方向分子块划分个数为ngx,设定行方向分子块划分个数为ngy;
33、分别计算每个子块在列方向的间隔、每个子块在行方向的间隔,如下:
34、
35、其中,cenx表示每个子块在列方向的间隔,ceny表示每个子块在行方向的间隔,width表示强度图影像的像素宽,height表示强度图影像的像素高,ngx为设定列方向子块划分个数,ngy为设定行方向子块划分个数;
36、第j行第i列的子块,定义为:
37、celi,j
38、i∈[1,ngx],j∈[1,ngy]
39、其中,i表示列方向的子块序号,j表示行方向的子块序号,ngx为设定列方向子块划分个数,ngy为设定行方向子块划分个数;
40、设定子块尺寸为kernel;
41、计算第j行第i列的子块的左上角像素坐标,具体如下:
42、
43、
44、
45、其中,cenx表示每个子块在列方向的间隔,ceny表示每个子块在行方向的间隔,表示第j行第i列的子块的左上角像素坐标,表示第j行第i列的子块的左上角像素坐标的列坐标,表示第j行第i列的子块的左上角像素坐标的行坐标,kernel表示设定子块尺寸大小;
46、第j行第i列的子块的左上角像素坐标对应的地理坐标,定义为:
47、
48、其中,表示第j行第i列的子块的左上角像素坐标对应的地理坐标的经度;表示第j行第i列的子块的左上角像素坐标对应的地理坐标的纬度;表示第j行第i列的子块的左上角像素坐标对应的地理坐标的高程;
49、计算第j行第i列的子块的右上角像素坐标,具体如下:
50、
51、
52、
53、其中,cenx表示每个子块在列方向的间本文档来自技高网...
【技术保护点】
1.一种Ka-SAR卫星标准产品生成方法,其特征在于,包括:
2.根据权利要求1所述的卫星标准产品生成方法,其特征在于,包括以下步骤:
3.根据权利要求2所述的卫星标准产品生成方法,其特征在于,步骤2所述生成对应的强度图影像,具体如下:
4.根据权利要求3所述的卫星标准产品生成方法,其特征在于:
5.根据权利要求4所述的卫星标准产品生成方法,其特征在于:
6.根据权利要求5所述的卫星标准产品生成方法,其特征在于:
7.根据权利要求6所述的卫星标准产品生成方法,其特征在于:
8.根据权利要求7所述的卫星标准产品生成方法,其特征在于:
9.根据权利要求8所述的卫星标准产品生成方法,其特征在于:
10.一种计算机可读介质,其特征在于,其存储电子设备执行的计算机程序,当所述计算机程序在电子设备上运行时,使得所述电子设备执行如权利要求1-9任一项所述方法的步骤。
【技术特征摘要】
1.一种ka-sar卫星标准产品生成方法,其特征在于,包括:
2.根据权利要求1所述的卫星标准产品生成方法,其特征在于,包括以下步骤:
3.根据权利要求2所述的卫星标准产品生成方法,其特征在于,步骤2所述生成对应的强度图影像,具体如下:
4.根据权利要求3所述的卫星标准产品生成方法,其特征在于:
5.根据权利要求4所述的卫星标准产品生成方法,其特征在于:
6.根据权利要求5...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。