TROPOMI卫星甲醛数据和模式模拟对比评估及数据处理

​
目的：使用卫星数据评估模式模拟结果

数据来源：下载方法

TROPOMI数据下载方法
这里下载了TROPOMI源数据，因为是nc格式的，又没办法在arcgis中打开，因此最后不得不选择GEE下载l3水平的产品，并将每一天的数据和月均的数据下载和计算出来，并保存为tiff格式

处理：

1.tif——>nc

因为GEE数据的分辨率 是1113.2m，而我模式模拟的是3km的，为了绘图方便，因此需要统一分辨率。但是这个tiff文件我没有试过直接转为nc，也许在arcgis的“多维”——“栅格转NetCDF”可以转，我的是师兄直接通过冯师兄的账户帮我转的，非常方便。

2.nc——>tif

因为我不会用沈师兄的插值代码，所以只能用笨办法，就是将计算好的月均的甲醛柱密度导出到nc文件里面，然后在arcgis中处理这个nc文件。

var_0 =col_avg_1st/1.0E15  ;这个是最终计算得出的变量
;-----
diml=dimsizes(lat2d)
data = new((/diml(0),diml(1)/),float)
copy_VarMeta(lat2d, data)
copy_VarAtts(lat2d, data)
copy_VarCoords(lat2d, data)
do i=0,173
  do j=0,167
   x=var_0(i,j) ;需要一个中间变量来传递，而不能直接赋值
   data(i,j)=x  ;这里data[south_north | 174] x [west_east | 168]
  end do
end do
hcho= data  ;这里分别把要用的甲醛、经纬度赋值到变量里面，他们的维度都是一样的
lat1=lat2d
lon1=lon2d
outfile  = addfile("/data1/loucx/COPY/CNAQ/data1/cctm/outtest.nc","c") ;导出变量
outfile->HCHO=hcho
outfile->XLAT=lat1 ;别忘了经纬度
outfile->XLONG=lon1

存在的问题：因为直接导出之后到arcgis里面它的单位是位置单位，而且是按照网格来划分的，完全就不是想要的按照经纬度划分的，无法绘图。
​​

1.将nc文件导入arcgis

“多维工具”——“创建NetCDF栅格图层”
但是会报如下的错误：

接下来需要按照这3个文档作为参考：
1.ArcGIS10.2读取NetCDF格式的气象数据含风向
学习它如何设置参数，这个设置完之后会出现一个点图层：

2.将矢量的图层变成栅格图片

参考这个链接提到的IDW插值，所以“Spatial Analyst工具”——“插值分析”——“反距离权重法”
这里设置像元大小（重采样后的栅格）、环境-处理范围
点转矢量图

3.生成符合目标要求的像元大小的tif 文件，实现了将模式的结果提到arcgis里面

​​​