pygac-fdr

使用pygac创建AVHRR GAC数据基本数据记录（FDR）的Python包

安装

要安装最新版本

pip install pygac-fdr

要安装最新开发版本

pip install git+https://github.com/pytroll/pygac-fdr

用法

要读取和校准AVHRR GAC第1b级数据，请修改etc/pygac-fdr.yaml中的配置模板，然后运行

pygac-fdr-run --cfg=my_config.yaml /data/avhrr_gac/NSS.GHRR.M1.D20021.S0*

结果以netCDF格式写入指定的输出目录。之后，收集并补充生成的netCDF文件的元数据

pygac-fdr-mda-collect --dbfile=test.sqlite3 /data/avhrr_gac/output/*

这可能需要一些时间，因此结果将保存到数据库中。您可以指定来自多个平台上的文件；将分别对每个平台的元数据进行分析。在大量文件中，您可能会遇到命令行参数大小限制（“参数列表太长”）。在这种情况下，请使用以下命令从文件中读取文件名列表（每行一个）

pygac-fdr-mda-collect --dbfile=test.sqlite3 @myfiles.txt

最后，就地更新netCDF元数据

pygac-fdr-mda-update --dbfile=test.sqlite3

AVHRR GAC FDR用户提示

检查全局质量标志

全局质量标志可以从以下命令行进行检查

ncks -CH -v global_quality_flag -s "%d" myfile.nc

裁剪重叠

由于数据接收机制，连续的AVHRR GAC文件通常部分包含相同的信息。这就是我们所说的重叠。例如，文件A末尾的一些扫描行也出现在文件B的开头。在pygac-fdr输出文件中的overlap_free_start和overlap_free_end属性指示重叠。有两种方法可以去除它

• 与后续文件裁剪重叠：选择扫描行0:overlap_free_end
• 与前一文件裁剪重叠：选择扫描行overlap_free_start:-1

如果用户还想创建每日合成图，则必须使用包含两天观测数据的文件两次：一次只使用UTC 00:00之前的部分，另一次只使用UTC 00:00之后的部分。同时裁剪重叠部分和整个一天会稍微复杂一些，因为重叠部分可能会覆盖UTC 00:00。这就是为什么提供了pygac-fdr-crop工具的原因。

$ pygac-fdr-crop AVHRR-GAC_FDR_1C_N06_19810330T225108Z_19810331T003506Z_...nc --date 19810330
0 8260
$ pygac-fdr-crop AVHRR-GAC_FDR_1C_N06_19810330T225108Z_19810331T003506Z_...nc --date 19810331
8261 12472

返回的数字是起始和结束扫描线（基于0）。