Cartopy画地图第七天（python画浮雕地图和比例尺）

本文链接：https://blog.youkuaiyun.com/weixin_42372313/article/details/119885922

本文介绍了如何使用Python的Cartopy库绘制浮雕地图，并手动创建比例尺。作者分享了下载浮雕地图文件的来源，以及详细步骤，包括设置地图边界、导入浮雕图像、绘制经纬度网格、画比例尺和各种散点图例的方法。此外，还提供了自定义函数以实现比例尺和不同形状的散点图。

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

Cartopy画地图第七天（python画浮雕地图和比例尺）

本文利用了python、cartopy进行了浮雕地图的绘制，同时还画了比例尺。
先上图为敬，一些图例符号不对请不要介意，随便表示的
在这里插入图片描述

第一、下载浮雕地图

想画浮雕地图的朋友们，第一步首先是要下载到浮雕地图文件，网站是Natural Earth，网址是：https://www.naturalearthdata.com/downloads/
在这里插入图片描述
只有1：10和1：50的有浮雕地图，就是第三个选项“栅格”（别问为啥我是中文的，问我就告诉你我是浏览器翻译的。）
点进去之后找自己需要的下载就好了，反正我的电脑1：10里面最大的那个地图画不了，报错显示地图太大，其他都没问题，大家自己按实际情况下载就好了。
在这里插入图片描述

下载完成以后是一个压缩包，解压得到里面的xxxxx.tif文件就好了，放到代码目录下就能使用。

第二、画浮雕地图

其实没啥难度，就是告诉大家cartopy是可以画浮雕地图的

    #地图四周经纬度
    west = 110
    east = 125
    south = 20
    north = 30
    #解决中文乱码问题
    plt.rcParams['font.sans-serif'] = ['SimHei']
    #画布设置，投影设置
    fig = plt.figure(figsize=(16,9.6))
    ax = fig.add_subplot(111,projection = ccrs.PlateCarree())
    #边界设置
    img_extent = [west, east, south, north]
    ax.set_extent(img_extent,crs = ccrs.PlateCarree())

    #浮雕地图导入
    fname = 'HYP_LR_SR_OB_DR.tif'
    ax.imshow(plt.imread(fname), origin='upper', transform=ccrs.PlateCarree(),extent=[-180, 180, -90, 90])

在这里插入图片描述
有些朋友到这里就OK了。

第三、地图美画

经纬度网格设置

	#设置经纬度网格和标签，透明度改成1，就会出现网格 alpha=1
    gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, linewidth=1, color='k', alpha=1,linestyle='--')#坐标轴设置
    gl.xformatter = LONGITUDE_FORMATTER ##坐标刻度转换为经纬度样式
    gl.yformatter = LATITUDE_FORMATTER
    #设置经纬度的显示，例如经度从70到135，间隔5显示
    gl.xlocator = mticker.FixedLocator(np.arange(70,135,5))
    gl.ylocator = mticker.FixedLocator(np.arange(15,55,5))

在这里插入图片描述

画比例尺和黑白边框

def drow_the_scale(y,x,text,length = 1.5,lw = 5):
    #画比例尺函数
    # y代表比例尺所在纬度
    # x代表比例尺开始的经度
    # text代表比例尺最后刻度值
    # length代表比例尺的长度，单位为多少个经度
    # lw代表比例尺的宽度
    step = length/5#计算步长，画五格
    #画黑白线五条
    plt.hlines(y=y,xmin=x,xmax=x + step,colors="black", ls="-", lw=lw)
    plt.hlines(y=y,xmin=x + step,xmax=x + step*2,colors="white", ls="-", lw=lw)
    plt.hlines(y=y,xmin=x + step*2,xmax=x + step*3,colors="black", ls="-", lw=lw)
    plt.hlines(y=y,xmin=x + step*3,xmax=x + step*4,colors="white", ls="-", lw=lw)
    plt.hlines(y=y,xmin=x + step*4,xmax=x + step*5,colors="black", ls="-", lw=lw)
    #画长刻度两个
    plt.vlines(x = x, ymin = y - (lw/100) *3, ymax = y + lw/100, colors="black", ls="-", lw=1)
    plt.vlines(x = x + length, ymin = y - (lw/100) *3, ymax = y + lw/100, colors="black", ls="-", lw=1)
    #画段刻度四个
    plt.vlines(x = x + step, ymin = y - (lw/100) *2, ymax = y + lw/100, colors="black", ls="-", lw=1)
    plt.vlines(x = x + step*2, ymin = y - (lw/100) *2, ymax = y + lw/100, colors="black", ls="-", lw=1)
    plt.vlines(x = x + step*3, ymin = y - (lw/100) *2, ymax = y + lw/100, colors="black", ls="-", lw=1)
    plt.vlines(x = x + step*4, ymin = y - (lw/100) *2, ymax = y + lw/100, colors="black", ls="-", lw=1)
    #写字，0，500，km
    plt.text