ArcGIS制作全球地图并生成纬度统计分布线

最新推荐文章于 2024-01-31 18:00:00 发布

地学万事屋

最新推荐文章于 2024-01-31 18:00:00 发布

阅读量2.9k

点赞数 4

文章标签：程序人生

本文链接：https://blog.csdn.net/wlh2067/article/details/129901804

版权

本文介绍了如何使用ArcGIS进行地图制作，包括数据准备、修改投影到NaturalEarth，并利用Python的gdal包进行纬度统计。之后，数据被转换为CSV并在Origin中创建线条图表进行可视化展示。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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的纬度统计

可以去Origin制图了

直接复制粘贴数据：

然后选择Plot——Line

翻转x，y轴、修改样式。

最终效果

本文由 mdnice 多平台发布