之前写过一次关于提取全国空气质量站点的,共计是全国2000个左右的站点数据,2013年到最近的,小时的。不过现在需要绘制出这些数据来看看长什么样子,于是决定写下来
首先数据是分散的,还没有按各省各站点分布,关于怎么按各省分站点单独提取出来的内容之前已经写过了,可以查看一下阅读


最后提取得到的就是分省按站点存储的数据了

因此在之前提取的基础上,我想做一个空间分析的简单操作,看看要素在全国某一个时间点上的分布情况,于是简单提取一下,预期的效果就是第一列表示全国的站点,然后是经纬度和字段名字
由于字段比较多,这里仅仅使用了aqi,pm2.5和pm10三个字段就行了

然后遍历二级文件夹,读取csv文件,提取第一个时间点,将列表转成numpy数组再转成DataFrame,最后将数据导出,准备用ArcGIS来看看全国空间尺度上的空气质量分布情况。当然了也可以不用提取第一个时间的,可以弄成月的,年的都可以

然后通过ArcGIS展示一下全国站点的分布,会发现有一个站点跑出去了

然后插值出来看看提取的这个时间点的aqi空间分布情况

最后导出图形去就好了,这个还是比较简单的
