如何同时使用MEIC和MIX清单?
作者:王浩帆
MEIC清单仅为中国境内的排放清单,但是在模拟全国污染场的案例中,中国周边国家的排放是不容忽视的,因此需要通过MIX清单来对MEIC进行一个补充。
不论是模拟网格分辨率大于等于清单网格分辨率,还是模拟网格分辨率小于清单网格分辨率的情况,同时使用MEIC和MIX清单的关键步骤都是如何将MEIC清单镶嵌到MIX中, 作为一系列新的GeoTIFF文件来作为coarse_emission_2_fine_emission.py[1]和 fine_emission_2_coarse_emission.py[2]的输入。
因此本部分将重点讲解如何使用工具来完成两个系列GeoTIFF的镶嵌工作。
1.将MIX清单和MEIC清单都转换为GeoTiff格式。
•使用mix_2_GeoTiff.py[3]将MIX清单转换为GeoTiff格式。•使用meic_2_GeoTiff.py[4]将MEIC清单转换为GeoTiff格式。由于MIX清单中没有PMC,因此需要通过calculate-pmc.py[5]将其计算出来。
1.配置combine.py[6]中的输入参数。
•upper_raster_dir:上层GeoTiff所在目录路径。
•bottom_raster_dir:下层GeoTiff所在目录路径。
•output_dir:输出GeoTiff目录路径。
•upper_raster_pollutants:上层GeoTiff需要合并的污染物名称。
•bottom_raster_pollutants:下层GeoTiff所对应的污染物名称。
•output_pollutants:一一对应到的输出污染物的名称。
•upper_label:上层GeoTiff标签。
•bottom_label:下层GeoTiff标签。
•output_label:输出GeoTiff标签。
•upper_raster_template:任意一个上层GeoTiff文件路径。
•bottom_raster_template:任意一个下层GeoTiff文件路径。
•upper_resolution:上层GeoTiff的分辨率。
•bottom_resolution:下层GeoTiff的分辨率。
•sectors:需要合并的部门。
•upper_year:上层GeoTiff的年份。
•bottom_year:下层GeoTiff的年份。
•output_year:定义输出GeoTiff的年份。
1.运行combine.py[7]
在终端中输入:
代码语言:javascript复制python ./combine.py
便可以开始运行程序,程序结束后将在output_dir
中产生合并后的系列GeoTiff。
1.进行空间分配、物种分配和时间分配。
此步骤和第一个教程[8]或第二个教程中的步骤完全相同,不再赘述。
References
[1]
coarse_emission_2_fine_emission.py: ../coarse_emission_2_fine_emission.py
[2]
fine_emission_2_coarse_emission.py: ../fine_emission_2_coarse_emission.py
[3]
mix_2_GeoTiff.py: ../PREP/mix_2_GeoTiff.py
[4]
meic_2_GeoTiff.py: ../PREP/meic_2_GeoTiff.py
[5]
calculate-pmc.py: ../calculate-pmc.py
[6]
combine.py: ../UTIL/combine/combine.py
[7]
combine.py: ../UTIL/combine/combine.py
[8]
第一个教程: 1-adopt_meic_for_prd_emission_file.md