Python海洋专题五之水深地形图海岸填充

Python海洋专题五之水深地形图海岸填充
海洋与大气科学
上期读取nc水深文件，并出图

但是存在一些不完美，本期修饰

本期内容

障眼法：把大于零的数据填充为陆地的灰色；

把等于零的数据画等深线为陆地和海洋的分界线！为海岸线。

效果图

参考文献及其在本文中的作用

1：matplotlib的colorbar设置显示的刻度个数和指定的刻度值_plt.bar显示个数-CSDN博客

其作用设置colorbar的刻度。

全文代码

# -*- coding: utf-8 -*-
# %%
# Importing related function packages
import pandas as pd
# import os
# from PIL import Image

import matplotlib.pyplot as plt  # https://matplotlib.org/
import numpy as np
from netCDF4 import Dataset
import matplotlib as mpl
# import palettes
from palettable.cartocolors.sequential import DarkMint_7
from palettable.colorbrewer.sequential import GnBu_9
from palettable.colorbrewer.sequential import Blues_9
from pylab import *
import matplotlib.ticker as ticker

def reverse_colourmap(cmap, name='my_cmap_r'):
    reverse = []
    k = []

    for key in cmap._segmentdata:
        k.append(key)
        channel = cmap._segmentdata[key]
        data = []

        for t in channel:
            data.append((1 - t[0], t[2], t[1]))
        reverse.append(sorted(data))

    LinearL = dict(zip(k, reverse))
    my_cmap_r = mpl.colors.LinearSegmentedColormap(name, LinearL)
    return my_cmap_r


# read data
a = Dataset('D:\pycharm_work\data\scs_etopo.nc')
print(a)
lon = a.variables['lon'][:]
lat = a.variables['lat'][:]
ele = a.variables['elevation'][:]
# plot figure
cmap = GnBu_9.mpl_colormap
cmap_r = reverse_colourmap(cmap)
print(Blues_9.colors)
ele[ele > 0] = 1000
#%% 图五
fig = plt.figure(dpi=300, figsize=(2, 2), facecolor='w', edgecolor='blue')
mngr = plt.get_current_fig_manager()  # 获取当前figure manager
mngr.window.wm_geometry("+58+51")  # 调整窗口在屏幕上弹出的位置
ax = fig.add_axes([0.08, 0.15, 0.85, 0.84])  # [距离左边，下边，坐标轴宽度，坐标轴高度] 范围(0, 1)
# 改变文字大小参数-fontsize
plt.xticks(np.arange(105, 125, 2), fontsize= 4)
plt.yticks(np.arange(0, 25, 2), fontsize= 4)
 # 绘制地图
#画0的等值线为海岸线。
cmarker=plt.contour(lon, lat, ele[:, :], levels=0, colors='k',linewidths=0.25)
levels = [-8000,  -6000,  -4000,  -2000, -1000, -500,  -300, -100,-50, 0, 1000]  # 创建分级
#color = [[0,51/255,102/255],[0,76/255,153/255],[0,102/255,204/255],[0,128/255,255/255],[51/255,153/255,255/255],[102/255,178/255,255/255],[153/255,204/255,255/255],[204/255,229/255,255/255],[160/255,160/255,160/255]] # 设置色带
color = [[0.0313725490196078,   0.211764705882353, 0.466666666666667],[0.0313725490196078,    0.274509803921569, 0.556862745098039],[0.0627450980392157,    0.376470588235294, 0.662745098039216],[0.203921568627451, 0.521568627450980, 0.749019607843137],[0.349019607843137, 0.639215686274510, 0.815686274509804],[0.454901960784314, 0.701960784313725, 0.847058823529412],[0.611764705882353, 0.788235294117647, 0.882352941176471],[0.784313725490196, 0.862745098039216, 0.941176470588235],[0.819607843137255, 0.886274509803922, 0.949019607843137],[0.7,0.7,0.7]]
plt.contourf(lon, lat, ele[:, :], levels=levels, colors=color)
# ------tick向内
tick_params(direction='in')
# -------colorbar调整
position = fig.add_axes([0.075, 0.05, 0.85, 0.03])#位置[左,下,右,上]
cb = plt.colorbar(label="depth", cax=position, orientation='horizontal')#方向
cb.ax.tick_params(labelsize=4) #设置colorbar刻度字体大小。
cb.set_label('depth',fontsize= 4,color='k' )#设置colorbar的标签字体及其大小
# 设置colorbar的刻度
cb.set_ticks([-8000,  -6000,  -4000,  -2000, -1000, -500,  -300, -100,-50, 0, 1000] )
cb.update_ticks()
# ------tick向内
tick_params(direction='in')
plt.savefig('scs_elevation.jpg', dpi=600, bbox_inches='tight', pad_inches=0.1)  # 输出地图，并设置边框空白紧密
plt.show()
图片
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77