当前位置： 首页 > news >正文

网站快速收录软件大连企业网站排名

news 2025/9/17 19:53:20
网站快速收录软件,大连企业网站排名,成都有几个区,个人求职网页制作python绘制年平均海表温度、盐度、ph图 文章目录 python绘制年平均海表温度、盐度、ph分布图前言一、数据准备二、代码编写2.1. python绘制年平均海表温度（主要）2.2. python绘制年平均海表盐度（选看）2.3. python绘制年平均海表ph&…

python绘制年平均海表温度、盐度、ph图

文章目录

  • python绘制年平均海表温度、盐度、ph分布图
    • 前言
    • 一、数据准备
    • 二、代码编写
      • 2.1. python绘制年平均海表温度（主要）
      • 2.2. python绘制年平均海表盐度（选看）
      • 2.3. python绘制年平均海表ph（选看）
    • 总结

python绘制年平均海表温度、盐度、ph分布图

所属目录：紫菜创建时间：2025/2/18更新时间：2025/2/19URL：https://blog.csdn.net/2301_78630677/article/details/145716784

前言

本文主要使用python绘制年平均海表温度、盐度、ph分布图，所用数据来源于Bio-ORACLE
参考文章：
Python绘制海表温度
【python海洋专题十二】年平均的南海海表面温度图

所用到的中国地图shp文件：
链接：https://pan.baidu.com/s/1q9hitI11CCYDWvBTWbAevg
提取码：9ju8

一、数据准备

在Bio-ORACLE官网 下载所需的 2010-2020平均海表Ocean temperature、Salinity、pH数据
点击前往下载地址
在这里插入图片描述

下载下来的环境数据为.nc文件，也就是NetCDF格式。

‌NetCDF（Network Common Data Form）格式是一种用于存储和共享科学数据的标准格式，广泛应用于气象学、海洋学、地球科学等领域‌。.nc文件是NetCDF文件的扩展名，主要用于存储大型科学和工程数据集。

.nc文件的基本结构和特点
‌自描述性‌：.nc文件包含关于数据集的元数据，这些元数据描述了数据集的结构和内容，使得用户无需其他文档即可理解数据。
‌可移植性‌：.nc文件是二进制格式，能够在不同的平台上无缝迁移和使用。
‌多维数组结构‌：.nc文件通常包含多个维度，如时间、经度和纬度，适用于存储复杂的多维数据集。

补充：
所用的Bio-ORACLE环境数据合集（有需要就下载吧）
Bio-ORACLE数据分享[decade 2010-2020] [Surface layers]

二、代码编写

接下来主要讲述 python绘制年平均海表温度的代码，另外两个类似，只需要稍加修改

2.1. python绘制年平均海表温度（主要）

该代码用于绘制中国周边海域的海表温度（SST）分布图，并添加了省份边界、等温线和网格线等细节，最后保存pdf文件

import xarray as xr  
import matplotlib.pyplot as plt  
import cartopy.crs as ccrs  
import cartopy.feature as cfeature  
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter  # 导入经纬度格式器 
import numpy as np
from cartopy.io.shapereader import Reader
from cartopy.feature import ShapelyFeature
# 加载NetCDF格式的SST数据（替换为你的SST数据文件路径）  
ds = xr.open_dataset(r"D:\oceandata\Bio-ORACLE\Temperature [mean].nc")  # 假设你的SST数据在该文件中  print(ds.variables) #打印SST数据的所有变量名及其相关信息，通过查看这些信息，你可以确定要使用哪个变量进行绘图和分析# 选择SST变量（替换为你的SST变量名）  
sst = ds['thetao_mean']# 计算时间轴上的平均值（如果时间是一个维度）  
sst_mean = sst.mean(dim='time')  # 假设'time'是时间维度  # 创建一个地图并设置投影  
fig = plt.figure(figsize=(10, 5))  
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())  # 添加陆地和海洋特征  
ax.add_feature(cfeature.LAND, color='lightgray')  
ax.add_feature(cfeature.OCEAN, color='w', edgecolor='lightgray')  
ax.coastlines(color='black')  
# 添加省份边界
shapefile = r"C:\Users\www\Desktop\china_map\china_SHP\省界_Project.shp"  # 替换为你的Shapefile文件路径
china_provinces = ShapelyFeature(Reader(shapefile).geometries(), ccrs.PlateCarree(), edgecolor='black', facecolor='none')
ax.add_feature(china_provinces)# 绘制SST平均值数据  
sst_plot = sst_mean.plot.contourf(ax=ax, transform=ccrs.PlateCarree(), cmap='coolwarm', levels=25, extend='both', add_colorbar=False,vmin=5, vmax=30)  
# levels参数可以调整等值线的数量  （具体来说，levels=25 表示将数据范围分成25个间隔，并绘制出相应的等值线。这些等值线将数据集的值范围（在此例中是5到30°C）平均分成25个部分，每个部分的上限和下限定义了一条等值线。）# 添加颜色条  
cbar = fig.colorbar(sst_plot, drawedges=True, ax=ax, location='right', shrink=0.95, pad=0.08, spacing='uniform', label='Average Sea Surface Temperature (°C)')  
cbar.ax.tick_params(labelsize=10)  # 设置色标尺标签大小 # 设置颜色条的刻度标签
cbar.set_ticks(np.arange(5, 31, 5))# 添加等温线  
sst_contour = sst_mean.plot.contour(ax=ax, transform=ccrs.PlateCarree(), colors='gray', levels=130,linewidths=0.5)   # levels参数可以调整等温线的数量  # 为等值线添加标签  
plt.clabel(sst_contour, inline=True, fontsize=10, fmt='%1.1f') # 设置地图的经纬度范围（可选）  
ax.set_extent([110, 135, 20, 40], crs=ccrs.PlateCarree())  # 添加网格线  
gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, xlocs=np.arange(110, 135, 5), ylocs=np.arange(20, 40, 5),linewidth=0.5, linestyle='--', color='k', alpha=0.8)  # 添加网格线
gl.xlabels_top = False  
gl.ylabels_right = False  
gl.xformatter = LongitudeFormatter()  # 使用默认的经度格式器  
gl.yformatter = LatitudeFormatter()   # 使用默认的纬度格式器  
gl.xlabel_style = {'size': 10, 'color': 'black'}  
gl.ylabel_style = {'size': 10, 'color': 'black'}  
print("Map created successfully!")
# 保存地图为PDF文件
plt.savefig('scs_sst_1.pdf', dpi=600, bbox_inches='tight', pad_inches=0.1)
# 显示地图  
plt.show()

代码大概包括以下流程：

  1. 加载数据：读取 NetCDF 文件并提取 SST 变量。
  2. 数据处理：对时间维度取平均。
  3. 创建地图：设置投影、添加陆地和海洋特征、绘制省份边界。
  4. 绘制 SST 数据：绘制填充等温线和等温线，添加颜色条和标签。
  5. 设置地图范围与网格线：调整地图范围，添加网格线并格式化标签。
  6. 保存与显示：保存地图为 PDF 文件并显示。

print(ds.variables) #打印SST数据的所有变量名及其相关信息，通过查看这些信息，你可以确定要使用哪个变量进行绘图和分析
例如该SST数据包括四个变量：time、latitude、longitude、thetao_mean， 其中thetao_mean就是SST变量名。

在这里插入图片描述

结果显示
在这里插入图片描述

2.2. python绘制年平均海表盐度（选看）

以下代码与年平均海表温度的代码类似（一些注释信息就没改过来了，知道意思即可）


import xarray as xr  
import matplotlib.pyplot as plt  
import cartopy.crs as ccrs  
import cartopy.feature as cfeature  
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter  # 导入经纬度格式器 
import numpy as np
from cartopy.io.shapereader import Reader
from cartopy.feature import ShapelyFeature
# 加载NetCDF格式的SST数据（替换为你的数据文件路径）  ds = xr.open_dataset(r"D:\oceandata\so_baseline_2000_2019_depthsurf_49ed_5fc4_602c_U1739344920620_yandu.nc")  # 假设你的SST数据在该文件中  print(ds.variables)
# 选择SST变量（替换为你的SST变量名）  
sst = ds['so_mean']# 计算时间轴上的平均值（如果时间是一个维度）  
sst_mean = sst.mean(dim='time')  # 假设'time'是时间维度  # 创建一个地图并设置投影  
fig = plt.figure(figsize=(10, 5))  
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())  # 添加陆地和海洋特征  
ax.add_feature(cfeature.LAND, color='lightgray')  
ax.add_feature(cfeature.OCEAN, color='w', edgecolor='lightgray')  
ax.coastlines(color='black')  
# 添加省份边界
shapefile = r"C:\Users\www\Desktop\a5bc0-main\china_SHP\省界_Project.shp"  # 替换为你的Shapefile文件路径
china_provinces = ShapelyFeature(Reader(shapefile).geometries(), ccrs.PlateCarree(), edgecolor='black', facecolor='none')
ax.add_feature(china_provinces)# 绘制SST平均值数据  
sst_plot = sst_mean.plot.contourf(ax=ax, transform=ccrs.PlateCarree(), cmap='coolwarm', levels=15, extend='both', add_colorbar=False,vmin=20, vmax=35)  # levels参数可以调整等温线的数量  # 添加颜色条  
cbar = fig.colorbar(sst_plot, drawedges=True, ax=ax, location='right', shrink=0.95, pad=0.08, spacing='uniform', label='Average Sea Surface Salinity (psu)')  
cbar.ax.tick_params(labelsize=10)  # 设置色标尺标签大小 # 设置颜色条的刻度标签
cbar.set_ticks(np.arange(20, 36, 3))# 添加等温线  
sst_contour = sst_mean.plot.contour(ax=ax, transform=ccrs.PlateCarree(), colors='gray', levels=50,linewidths=0.5)  # 为等值线添加标签  
plt.clabel(sst_contour, inline=True, fontsize=10, fmt='%1.1f') # 设置地图的经纬度范围（可选）  
ax.set_extent([110, 135, 20, 40], crs=ccrs.PlateCarree())  # 添加网格线  
gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, xlocs=np.arange(110, 135, 5), ylocs=np.arange(20, 40, 5),linewidth=0.5, linestyle='--', color='k', alpha=0.8)  # 添加网格线
gl.xlabels_top = False  
gl.ylabels_right = False  
gl.xformatter = LongitudeFormatter()  # 使用默认的经度格式器  
gl.yformatter = LatitudeFormatter()   # 使用默认的纬度格式器  
gl.xlabel_style = {'size': 10, 'color': 'black'}  
gl.ylabel_style = {'size': 10, 'color': 'black'}  
print("Map created successfully!")# 保存地图为PDF文件
plt.savefig('scs_yandu_1.pdf', dpi=600, bbox_inches='tight', pad_inches=0.1)
# 显示地图  
plt.show()

在这里插入图片描述

2.3. python绘制年平均海表ph（选看）

以下代码与年平均海表温度的代码类似，只不过因为ph的值较小，我们可以先计算一下数据中的最大值与最小值，方便确定ph大小


import xarray as xr  
import matplotlib.pyplot as plt  
import cartopy.crs as ccrs  
import cartopy.feature as cfeature  
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter  # 导入经纬度格式器 
import numpy as np
from cartopy.io.shapereader import Reader
from cartopy.feature import ShapelyFeature
# 加载NetCDF格式的SST数据（替换为你的SST数据文件路径）  ds = xr.open_dataset(r"D:\oceandata\ph_baseline_2000_2018_depthsurf_f606_6dc8_6180_U1739344995788_ph.nc")  # 假设你的SST数据在'sst.nc'文件中  # 选择pH变量（替换为你的pH变量名）  
ph = ds['ph_mean']# 计算pH数据的最高和最低值
ph_max = ph.max(dim=['time', 'latitude', 'longitude'])
ph_min = ph.min(dim=['time', 'latitude', 'longitude'])print("Maximum pH value:", ph_max)
print("Minimum pH value:", ph_min)print(ds.variables)
# 选择SST变量（替换为你的SST变量名）  
sst = ds['ph_mean']# 计算时间轴上的平均值（如果时间是一个维度）  
sst_mean = sst.mean(dim='time')  # 假设'time'是时间维度  # 创建一个地图并设置投影  
fig = plt.figure(figsize=(10, 5))  
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())  # 添加陆地和海洋特征  
ax.add_feature(cfeature.LAND, color='lightgray')  
ax.add_feature(cfeature.OCEAN, color='w', edgecolor='lightgray')  
ax.coastlines(color='black')  
# 添加省份边界
shapefile = r"C:\Users\王浩天\Desktop\a5bc0-main\china_SHP\省界_Project.shp"  # 替换为你的Shapefile文件路径
china_provinces = ShapelyFeature(Reader(shapefile).geometries(), ccrs.PlateCarree(), edgecolor='black', facecolor='none')
ax.add_feature(china_provinces)# 绘制SST平均值数据  
sst_plot = sst_mean.plot.contourf(ax=ax, transform=ccrs.PlateCarree(), cmap='coolwarm', levels=15, extend='both', add_colorbar=False,vmin=8, vmax=8.2)  # levels参数可以调整等温线的数量  # 添加颜色条  
cbar = fig.colorbar(sst_plot, drawedges=True, ax=ax, location='right', shrink=0.95, pad=0.08, spacing='uniform', label='Average Sea Surface pH')  
cbar.ax.tick_params(labelsize=10)  # 设置色标尺标签大小 # 设置颜色条的刻度标签
cbar.set_ticks(np.arange(8, 8.25, 0.05))# 添加等温线  
sst_contour = sst_mean.plot.contour(ax=ax, transform=ccrs.PlateCarree(), colors='gray', levels=200,linewidths=0.5)  # 为等值线添加标签  
plt.clabel(sst_contour, inline=True, fontsize=10, fmt='%1.3f')  #保留小数点后三位# 设置地图的经纬度范围（可选）  
ax.set_extent([110, 135, 20, 40], crs=ccrs.PlateCarree())  # 添加网格线  
gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, xlocs=np.arange(110, 135, 5), ylocs=np.arange(20, 40, 5),linewidth=0.5, linestyle='--', color='k', alpha=0.8)  # 添加网格线
gl.xlabels_top = False  
gl.ylabels_right = False  
gl.xformatter = LongitudeFormatter()  # 使用默认的经度格式器  
gl.yformatter = LatitudeFormatter()   # 使用默认的纬度格式器  
gl.xlabel_style = {'size': 10, 'color': 'black'}  
gl.ylabel_style = {'size': 10, 'color': 'black'}  
print("Map created successfully!")
# 保存地图为PDF文件
plt.savefig('scs_ph_1.pdf', dpi=600, bbox_inches='tight', pad_inches=0.1)
# 显示地图  
plt.show()

在这里插入图片描述

总结

主要是以绘制年平均海表温度分布图为例，其余环境数据也类似。

2025/2/19

http://www.yayakq.cn/news/705489/

相关文章：

  • 餐厅网站建设文案书两学一做 专题网站
  • 如何做旅游计划的网站少林寺网站谁做的
  • 番禺区移动端网站制作美图秀秀在线制作
  • 云县网站建设找那家网站开发维护公司
  • 集客营销软件官方网站网站建设过程中什么最重要
  • wordpress站点美化免费图片制作网站模板
  • 微信支付 公司网站菏泽网站建设推广
  • 六安开发区网站建站技术博客
  • 文山网站建设公司如何做网站主页
  • 2019年做网站还有机会吗石碣镇网站建设
  • 个人怎么开通微信小程序成都外贸seo
  • 个人网站号备案吗杭州网站推广方式
  • 桂林做网站电话号码软广告经典案例
  • 平顶山河南网站建设农家乐网站模板
  • 淘宝网站制作公司哪家好wordpress用户邀请注册
  • 网站设计软件开发百度推广关键词规划师
  • 网站建设栏目流程全民体育世界杯
  • 开发东莞网站制作公司莱芜网吧恢复营业
  • 90设计网站创建时间在线报名网站建设
  • 备案要关闭网站吗上海市工商局企业查询
  • 教师做爰网站互联网站备案登记表
  • 建设铝合金窗网站公司网站建设制作全包
  • 好点子网站建设网站建设与维护实训总结
  • 网站模板怎么弄django做的网站有哪些
  • yellow网站推广联盟网店推广软件有哪些
  • 建网站要多少钱 优帮云网页开发培训网
  • 外贸网站推广教程ens域名注册
  • 做网站一个月30ip烟台当地网站
  • 医疗网站建设教程项目开发流程8个步骤流程图
  • 旅游网站后台模板商务网站内容维护范围