ArcGIS制作全球地图并生成纬度统计分布线
-
全球气候变化相关研究,或者发表在高水平期刊的研究,往往需要选择全球大尺度,需要进行世界地图的制作。
-
而世界地图的出图,也有一些技巧。有的研究是相当美观的,如发表在Nature上的一篇文章:
![image-20220314155140246](https://i-blog.csdnimg.cn/blog_migrate/7f7a6fc9df81977a05a5a8e148548b35.png)
若是只进行普通的出图,结果还是差了一点:
![image-20220313105145058](https://i-blog.csdnimg.cn/blog_migrate/95cf685d4e2bb0339128e5a83632a3eb.png)
今天这篇文章教大家添加经纬度线。我们仿照的是下图:
![image-20220313105256488](https://i-blog.csdnimg.cn/blog_migrate/bdbd9edb47dc792db2df7f8f763a4ae3.png)
地图部分
数据准备
-
通过一系列操作,如裁剪、重分类。
-
再建立布局视图,添加格网线,生成一张世界地图。
![image-20220313115417989](https://i-blog.csdnimg.cn/blog_migrate/ec243c13b2cca2b766b5783815d53039.png)
修改投影
首先需要改投影,用椭圆形,这是标准的世界地图格式:
![image-20220313115551906](https://i-blog.csdnimg.cn/blog_migrate/84db4fa1742a63acde878388bece682c.png)
![alt](https://i-blog.csdnimg.cn/blog_migrate/3bf9b9a6f2dbbce02ff3a5dc9c6f2ad0.png)
选择一个投影坐标系:在World文件夹里找到Robinson(world)或者Natural Earth(world),这两个效果都可以,这里选择Natural Earth
![image-20220313120004142](https://i-blog.csdnimg.cn/blog_migrate/fb2d0eff278fd98e618c8ea279b1b490.png)
基本效果如上图
按纬度统计
我们需要按纬度求均值,使用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的纬度统计
![image-20220313195513866](https://i-blog.csdnimg.cn/blog_migrate/15696b0ecc79dc9bf3ae2180e360b227.png)
可以去Origin制图了
直接复制粘贴数据:
![image-20220313195952648](https://i-blog.csdnimg.cn/blog_migrate/a1d27b7a3b45a17c3733f73c2c7d43e6.png)
然后选择Plot——Line
![image-20220313200730690](https://i-blog.csdnimg.cn/blog_migrate/062cf2f39fd43baddc5705348864e3cc.png)
翻转x,y轴、修改样式。
![image-20220313195920713](https://i-blog.csdnimg.cn/blog_migrate/7670a9948b28ac2764bff87661e65e8f.png)
最终效果
![image-20220313202109489](https://i-blog.csdnimg.cn/blog_migrate/3c92e0725ef8487da72af8e623052be5.png)
本文由 mdnice 多平台发布