ArcGIS制作全球地图并生成纬度统计分布线
-
全球气候变化相关研究，或者发表在高水平期刊的研究，往往需要选择全球大尺度，需要进行世界地图的制作。
-
而世界地图的出图，也有一些技巧。有的研究是相当美观的，如发表在Nature上的一篇文章：
若是只进行普通的出图，结果还是差了一点：
今天这篇文章教大家添加经纬度线。我们仿照的是下图：
地图部分
数据准备
-
通过一系列操作，如裁剪、重分类。
-
再建立布局视图，添加格网线，生成一张世界地图。
修改投影
首先需要改投影，用椭圆形，这是标准的世界地图格式：
选择一个投影坐标系：在World文件夹里找到Robinson(world)或者Natural Earth(world)，这两个效果都可以，这里选择Natural Earth
基本效果如上图
按纬度统计
我们需要按纬度求均值，使用python的gdal包几行代码就可以搞定了
from osgeo import gdal
import numpy as np
import pandas as pd
dataset = gdal.Open("D:/Acdemic/acad211229/data/mean90_mmc.tif")
cols=dataset.RasterXSize#图像长度
rows=(dataset.RasterYSize)#图像宽度
r=dataset.ReadAsArray(0,0,cols,rows)
r[r < -500000] = np.nan # 空值的像元，我这里是一极小值。可以把条件判断改成你空值的判断
x = np.nanmean(r, axis=1) # 按y轴统计
dataframe = pd.DataFrame({
'value': x})
# 将DataFrame存储为csv,index表示是否显示行名，default=True
dataframe.to_csv("test.csv", index=False)
于是我们导出了csv的纬度统计
扫描二维码关注公众号，回复：
14789828 查看本文章
可以去Origin制图了
直接复制粘贴数据：
然后选择Plot——Line
翻转x，y轴、修改样式。
最终效果
本文由 mdnice 多平台发布