GK2B-GOCI2数据介绍与MATLAB批量下载

最新推荐文章于 2024-04-19 11:07:25 发布

海研人

最新推荐文章于 2024-04-19 11:07:25 发布

阅读量2.2k

点赞数 8

文章标签： matlab

本文链接：https://blog.csdn.net/m0_51143849/article/details/129498675

版权

GOCI (Geostationary Ocean Color Imager)，是搭载在韩国卫星COMS上的海洋遥感传感器，于2010年成功发射，21年停运。作为世界首颗静止轨道海洋遥感卫星，它具有覆盖可见光到近红外范围的8个波段，500m空间分辨率，1h的时间分辨率，覆盖范围2500km×2500km，涵盖了中国的渤黄东海。优越的性能，特别是一天8景的高时间分辨率使其在海洋环境检测、渔业和航运等领域有着广泛的应用。

官网：http://kosc.kiost.ac.kr/index.nm?lang=en&contentId=82

作为其后继者，GOCI2于2020年2月发射，长江后浪推前浪，GOCI2的性能比他大哥更优越。高达13个波段、翻倍的空间分辨率（250m）、一天内可生成10景（北京时间7-16点）。随着数据量的不断积累，想必GOCI2会逐渐成为海洋遥感领域的又一大杀器。

GOCI2波段参数

在GOCI2的官网可下载从2020年12月至今的L1B级（天顶辐亮度）与二级产品数据。GOCI2数据是分块获取的，其中S009与S010覆盖了中国大部分海岸线。

GOCI2数据区块

GOCI2数据下载界面

下载地址：https://www.nosc.go.kr/eng/program/actionGociDownload.do

为了批量下载所需数据，我们可以用Python、Matlab等来编程实现，下面给大家分享一套我用Matlab编写的使用websave批量下载nc数据与jpg快视图的程序。程序中添加了保姆级注释，下载其他网站的数据其实思路和这个一样，同学们可在这个基础上进行删改。

卫星宫中号“海研人” 回复“GOCI2”领取示例代码

Step1丨“人工智能”筛选少云影像

为了获取研究区域上空无云/少云影像，我们可以通过批量下载快视图，然后通过“人工智能目视解译”（说白了就是自己一张一张看）的方式，来筛选出我们需要的影像日期。以下是一个批量下载指定日期与指定区块的快视图的MATLAB程序quickview.m。

程序功能很简单，更改起始日期和结束日期并运行程序后，就会在m文件所在的目录生成S009qw与S010qw两个文件夹并分别储存两个区块的快视图文件，使用时要注意将dindex.m这个函数添加到路径，且要在良好的网络环境下运行。下载后效果如下。

其中，部分jpg文件打不开，大小也只有2KB，是因为GOCI2在该时间、该区块没有数据，websave下载失败，中国海岸线所在的区块9和10经常会缺失00时的数据，具体原因我也不清楚。

下载好所需的快视图后，手动在excel记录需要的影像的日期编号，如210101，便可运行另一个程序直接继续进行下载。

Step2丨根据统计的日期编号下载对应数据

统计好所需影像的日期编号后便可运行batchdownload.m来批量下载nc数据和对应快视图并储存到相应文件夹中。该程序与获取快视图的程序思想一致，只不过添加了许多用来检查所下载文件是否正常以及预防网络出现波动后，再次运行程序可以在断点继续下载的步骤。

首先，读取储存日期编号的excel文件，设置所需下载的区块slot（别打错了）与时间。

然后，自动生成nc与jpg文件的文件名与下载链接url，并进入第一个检查模块。该模块会判断目标文件夹中是否已经下载了所需文件，若没有就进行下载，有就先判断它的大小正不正常，不正常就删了再下一次。这一步有两个作用，一是当下载进程由于某种原因而中断，需要再次运行程序下载时，不需要额外更改程序便可自动从中断处继续下载；二是重新下载已下载文件中因为网络波动而下载到的大小不正常的破损文件，一景GOCI2影像的大小普遍为422MB，大小低于400MB的nc文件将被视为破损文件并重新下载。

第二个检查模块会在每次下载后检查nc文件的大小，也有两个作用，首先可以删掉空文件，即区块9与10经常缺失的00时的数据，这些数据的url不存在，程序中生成url下载后是一个2KB的无效文件；其次，有时GOCI2的成像分秒不是15分30，而是16分30，后缀也就变成了1630而不是1530，这样的文件在上一步下载时得到的也是无效文件，会在这一步删除并重新下载。