相同卫星产品的较高值范围(NO2_OMI_TEMIS)

Question

我在处理对流层的NO2柱密度，我的数据源是特米．

NO2 level3数据可以从OMI仪器的原始信息中得到。

我下载了两类数据，.grd而且.kml同月。

我下载的数据在这里上传1,2

KML数据

我在谷歌Earth打开它，图是这样的:

我们可以看到NO2热带的数据范围。列为0~20，与网站上的模板图相同。

接地的数据

我没有找到关于这个数据的详细资料。在它的内容中，-999年视为no_data位置。

我用python来读取和绘制空间分布。

filename = './ cho - no2 /no2_201306. txt . txt . filename = './ cho - no2 /no2_201306. txt接地的' def read_grd(filename): ncols = np.array(linecache.getline(filename, 1)[6:10]).astype(float) nrows = np.array(linecache.getline(filename, 2)[6:10]).astype(float) xllcorner = np.array(linecache.getline(filename, 3)[10:14]).astype(float) yllcorner = np.array(linecache.getline(filename, 4)[10:14]).astype(float) cellsize = np.array(linecache.getline(filename, 5)[9:14]).astype(float) nan_value = np.array(linecache.getline(filename, 6)[13:17]).astype(float) longitude = xllcorner + cellsize * np.arange(ncols) latitude = yllcorner + cellsize * np.arange(nrows) value = np.loadtxt(filename, skiprows=7) value = value[::-1] return value, longitude, latitude, nan_value no2,lon_no2, lat_no2, nan_value = read_grd(filename) no2[no2 == nan_value] = np.nan def isnt_NaN(num): return num == num no2[isnt_NaN(no2)].max() > 9999.0

看来高价值在接地的格式数据不符合常规条件，这可以从以往的研究中了解到(NO2柱的热点约为15~20 10^15 molec/cm2)。

有人熟悉OMI-NO2数据吗?我不知道如何处理不规则值，这比实际值高得多。

Answer 1

似乎您担心数据中的值超过kml文件中显示的颜色条上的最大值。我不担心这个。您使用的是月平均“3级”产品，根据定义，这是高度质量控制的产品。对于全球大多数地区，典型的数值将在1 x 10^14和5 x 10^15摩尔/平方厘米之间变化。超过20 x 10^15 molec/cm2的值只出现在高度污染的地区。OMI NO2产品中非常大的值是现实的，特别是在东南亚，那里获得了最大的NO2柱。我有时甚至看到值超过1 x 10^17。当信号变大时，检索的可信度就会提高。实际上，应该谨慎对待的是较低的值(例如小于8 x 10^14 mol /cm2)，因为任何低于这个值的值都可能是“噪音”。

为了显示全球范围内值的方差(可以跨越4个数量级)，选择的颜色刻度效果很好。请注意，颜色条是对数刻度，这有助于定义具有典型值的区域中的图。不过，如果您对东南亚的更好的绘图感兴趣，您可以重写kml文件，并选择不同的最大值。

如果你有兴趣了解更多关于产品的信息，我强烈建议你看看Boersma (2011)也可以参考用户手册．