Cartopy虽然对地理图形的绘制提供了极大的方便,但是个人感觉还是有很多反人类的地方,尤其是在同一个代码中绘制多副图形的时候,比如说地理坐标轴的设置,每幅子图都要设置一遍,最终搞得代码十分的冗长。为了解决这类问题,自定义的封装函数就起了很大作用,只需在最开始定义好,接下来每次使用直接调用即可,再也不用每次画图都重新设置什么刻度,海岸线之类的了。
其实用到的就是python的def功能,通过def自定义函数,return返回所需值,可以极大的简洁重复代码。
比如说绘制下面这幅图,我分别使用未封装的代码和使用封装后的代码绘制,方便对比。
未使用封装函数的完整代码:
#公共设置(地图投影,地图边界,坐标刻度地理格式)
proj = ccrs.PlateCarree(central_longitude=80)
img_extent = [0,160, 0, 80]
lon_formatter = cticker.LongitudeFormatter()
lat_formatter = cticker.LatitudeFormatter()
fig2 = plt.figure(figsize=(15,15))
#子图1
f2_ax1 = fig2.add_axes([0.1, 0.1, 0.4, 0.3],projection = proj)
#边界,海岸线,湖泊,坐标刻度,坐标格式
f2_ax1.set_extent(img_extent, crs=ccrs.PlateCarree())
f2_ax1.add_feature(cfeature.COASTLINE.with_scale('50m'))
f2_ax1.add_feature(cfeature.LAKES, alpha=0.5)
f2_ax1.set_xticks(np.arange(leftlon,rightlon+20,20), crs=ccrs.PlateCarree())
f2_ax1.set_yticks(np.arange(lowerlat,upperlat+20,20), crs=ccrs.PlateCarree())
f2_ax1.xaxis.set_major_formatter(lon_formatter)
f2_ax1.yaxis.set_major_formatter(lat_formatter)
#图序
f2_ax1.set_title('(a) ',loc='left')
#填色
f2_ax1.contourf(lon,lat, r, levels=np.arange(-0.9,1.0,0.1),extend='both', zorder=0, transform=ccrs.PlateCarree(), cmap=plt.cm.RdBu_r)
#子图2,同上
f2_ax2 = fig2.add_axes([0.6, 0.1, 0.4, 0.3],projection = proj)
f2_ax2.set_extent(img_extent, crs=ccrs.PlateCarree())
f2_ax2.add_feature(cfeature.COASTLINE.with_scale('50m'))
f2_ax2.add_feature(cfeature.LAKES, alpha=0.5)
f2_ax2.set_xticks(np.arange(leftlon,rightlon+20,20), crs=ccrs.PlateCarree())
f2_ax2.set_yticks(np.arange(lowerlat,upperlat+20,20), crs=ccrs.PlateCarree())
f2_ax2.xaxis.set_major_formatter(lon_formatter)
f2_ax2.yaxis.set_major_formatter(lat_formatter)
f2_ax2.set_title('(b) ',loc='left')
f2_ax2.contourf(lon,lat, r, levels=np.arange(-0.9,1.0,0.1),extend='both', zorder=0, transform=ccrs.PlateCarree(), cmap=plt.cm.RdBu_r)
可以看到,对底图设置的越精细所需的代码越长,但是绘制多图时均为重复代码,那么再看看封装后是怎样的:
使用封装函数:
首先将重复部分进行封装,那我们先分析重复部分的内容,包含了边界,海岸线,湖泊,坐标刻度,坐标格式等内容,需要需要传递的主要是地图的边界信息和坐标的间隔信息,以及对应的子图
def contour_map(fig,img_extent,spec):
fig.set_extent(img_extent, crs=ccrs.PlateCarree())
fig.add_feature(cfeature.COASTLINE.with_scale('50m'))
fig.add_feature(cfeature.LAKES, alpha=0.5)
fig.set_xticks(np.arange(leftlon,rightlon+spec,spec), crs=ccrs.PlateCarree())
fig.set_yticks(np.arange(lowerlat,upperlat+spec,spec), crs=ccrs.PlateCarree())
lon_formatter = cticker.LongitudeFormatter()
lat_formatter = cticker.LatitudeFormatter()
fig.xaxis.set_major_formatter(lon_formatter)
fig.yaxis.set_major_formatter(lat_formatter)
封装好后,后面只需直接调用即可。需要注意的是,这里并没有return,是因为函数中直接作用在了传递进去的fig上,因此无需再return。
#公共设置
proj = ccrs.PlateCarree(central_longitude=80)
leftlon, rightlon, lowerlat, upperlat = (0,160,0,80)
img_extent = [leftlon, rightlon, lowerlat, upperlat]
fig2 = plt.figure(figsize=(15,15))
#子图1
f2_ax1 = fig2.add_axes([0.1, 0.1, 0.4, 0.3],projection = proj)
#调用封装函数
contour_map(f2_ax1,img_extent,20)
f2_ax1.set_title('(a) ',loc='left')
f2_ax1.contourf(lon,lat, r, levels=np.arange(-0.9,1.0,0.1),extend='both', zorder=0, transform=ccrs.PlateCarree(), cmap=plt.cm.RdBu_r)
#子图2
f2_ax2 = fig2.add_axes([0.6, 0.1, 0.4, 0.3],projection = proj)
#调用封装函数
contour_map(f2_ax2,img_extent,20)
f2_ax2.set_title('(b) ',loc='left')
f2_ax2.contourf(lon,lat, r, levels=np.arange(-0.9,1.0,0.1),extend='both', zorder=0, transform=ccrs.PlateCarree(), cmap=plt.cm.RdBu_r)
仅仅是两幅子图,并没有缩减太多的代码,但是当子图较多时,这个方法就十分实用了。尤其是在jupyter notebook中,仅需要在最开始的代码块中随着各种库的一起引入,后边无论新建多长的代码块,都能使用,可以说只需要定义一次,在整个科研阶段中都可以随意调用。
实际上这里只是以绘图的方式举例def功能,还有很多地方也可以用到类似的思路,比如说自己写了一个算法,比如说是滑动t检验,你需要对多个序列进行检验,如果每次都写一遍滑动t的代码,那整个脚本的长度就太可怕了,完全可以用def的方法,只定义一次,接下来每个序列只需要一行调用函数就可以得到想要的结果,并且可以通过return的设置,返回所需要的各种中间变量。