![cover](https://img-blog.csdnimg.cn/img_convert/c871bae5d19d4b64bdd8f6da99df754e.png)
Sentinel-1A-IW-SLC预处理及极化分解全过程(影像处理)
哨兵一号数据极化分解(cloude分解)在SNAP中操作全过程(完整版)
接上回!
影像处理的软件是SNAP软件,当然,也可以在其他软件中进行,但是SNAP是欧空局生产的软件,个人认为与哨兵一号数据更契合一些!
所以第一步就是需要下载SNAP软件。图标如下:
至于下载安装过程可以在CSDN中进行搜索,流程比较全面且成熟!个人建议安装到非C盘,因为在之后软件进行运作的时候回占用C盘的空间,将占内存较大的数据存到其他盘会比较稳妥!不过不必担心,当你进行完一步操作之后,C盘会在一段时间内自动释放内存,如果关掉SNAP的话,释放速度会更快。一般C盘中保存的数据会存到下图位置:
如果后期进行完某一项目后,不需要元数据,可以在此途径进行删除以释放内存。
正式开始预处理操作!
第一步:轨道校正(Apply Orbit File)
将原数据(压缩包)拖至Product Explorer区域。压缩包数据无需解压。
打开Radar操作中的Apply Orbit File
参数设置完成后点击“Run”运行即可,时间大概几分钟吧!运行完成后点击“Close”关掉页面即可。
第二步:S-1 TOPS Split
I/O Parameters面板中需注意数据产品是进行完“轨道校正”后的数据,此过程结束后会添加一个后缀名“split”,之后的操作也会自动在其后面添加相应操作的缩写后缀名。个人建议可以将名称进行修改,因为默认的名称是当时下载的数据的名称,过于长,可以进行适当修改。(我自己也需要重新进行修改,因为忘记添加时间了,如果研究期很长,有很多幅影像,没有时间的话可能会在后期的其他研究中出现时间顺序的错误,所以最好将时间加上。比较保险。)
例如名称可以改为我这样的:包含成像模式(IW)、数据产品类型(SLC)、时间(20240225)以及相应操作之后的后缀名(Orb)。有同学可能会发现我的名称有的没改有的改了,没关系,因为我也是边做边修改这篇文章,但这不影响,因为我在开头说明了这个事情就可以了。不用太过于纠结我的文件名称。
Processing Parameters面板中“Subswath”中三个条带依次选择一次(所以split这一操作需要进行三次)
极化方式默认会全都选中。至于“Bursts”条带的选择,可以根据自己的研究区占哪几个条带选择哪几个就可以。如果不确定的话可以都选择。设置完之后点击“Run”即可。运行时间大概十几秒吧,不是很长时间,不过这得看电脑配置以及网速。
之后还是按照之前的操作重复进行IW2与IW3,如下图:(数据源还是轨道校正完之后的数据源和IW1处理的数据源一样)
注:我在做IW3的相关内容时发现总是无法成功,界面会出现错误提示,通过一番查找我发现我本来的原数据中就是没有IW3内容,所以,当大家在进行操作的时发现无法成功的时候就可以在原数据中进行查看,查看路径如下图;因此,遇到任何问题不需要慌张,多查一查资料就可以啦!
再分享一个我遇到的问题,如果大家也遇到过同样的问题有解决方法的可以留言告诉我哦!谢谢大家啦!
我当时在查看我的原数据的时候发现它与之前下载的数据不一样,但我放置的研究区是同一个。具体情况如下图;研究区一样但是影像时间不一样导致我下载的影像也是不一样的,我目前没有找到解决的办法,想到的可能性就是:可能我下载数据的这个网址产生了问题,可以在哨兵一号的官网上下载,之后对比两副影像;也可能卫星在进行扫描时就是没有产生这一波段。都有可能吧,我会先尝试继续做后续步骤,如果可以完整的将影像裁剪的话,那就没有问题啦!
2024.3.14 终于找到了解决方法,这过程说实话就是非常的曲折。因为我从ASF网站上下载的影像打开后是没有IW3,我下载了好多不同时间,我以为是卫星的问题,我还通过下载数据推断得出从2023年10月份以后他就缺少IW3。然后哨兵一号的官网也关闭了,我又找到了新的官网,下载过程异常的慢,反正也是出了一点问题。最终我发现,只要把压缩包解压了他的图像就是完整的。我真的服了。
第三步:辐射校正(Calibration)
此处的source是我们上一步处理得到的IW1和IW2以及IW3,所以会做三次。
Processing Parameters面板中,“Polarisations”将所有极化方式选中,此处的极化类型取决于你当时数据下载时选择的类型,所以发现不一样不需要担心喔!
默认勾选的是“Output sigma0 band”,但我们需改为"Save as complex output"。点击“Run”。等待大约几分钟吧,不是很长时间,不过这得看电脑配置以及网速。
最终在左侧面板得到的结果如下:
第四步:burst拼接(Deburst)
这一操作同样进行两次或三次(如果有IW3的话)
结束之后我们看一下Deburst操作后的效果:
下面两张是上一步的结果图;
下面两张是进行deburst操作后的影像:
会发现影像拼接的更加完整,没有了条带之间的间隔(就是上一步结果图中的黑道)。
第五步:burst合并(TOPSAR Merge)
点击“Run”即可。
第六步:极化矩阵(Polarimetric-Matrices)
点击“Run”运行即可。
第七步:多视处理(Multilooking)
点击“Run”即可。
第八步:Speckle滤波(polarimetric-speckle-filter)
说明:该方法适用于极化SAR数据。
点击“Run”即可。
第九步:极化分解操作(Polarimetric-Decomposition)
这里极化方式的选择根据个人需求选择即可。目前我知道的是如果选择Cloude decomposition的话会不成功,因为哨兵一号的数据是双极化而Cloude decomposition要求全极化数据。可以自己进行选择和尝试,尝试过后就知道可不可以喽!
点击“Run”即可。
第十步:地理编码(Terrain-Correction)
点击“Run”即可!
至此,在SNAP中进行极化分解的操做全部结束。
![Logo](https://devpress.csdnimg.cn/79de2bf0b7994defa4242ef90d5513fa.jpg)
开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!
更多推荐
所有评论(0)