Cartopy画地图第七天(python画浮雕地图和比例尺) |
您所在的位置:网站首页 › 浮雕图怎么画出来的 › Cartopy画地图第七天(python画浮雕地图和比例尺) |
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(x,y - (lw/100) *7,'0',horizontalalignment = 'center') plt.text(x + length,y - (lw/100) *7,text,horizontalalignment = 'center') plt.text(x + length/2,y + (lw/100)*2,'km',horizontalalignment = 'center') def drowscale(extent,scale_y,scale_x,scale_text,step = 5,lw = 10,scale_length = 1.5,scale_lw = 5): # 画地图黑白边框和比例尺 # extent:表示四周经纬度[west, east, south, north] # scale_y,scale_x,scale_text:代表比例尺的位置,纬度,经度,刻度值 # step:表示步长,一格代表几个经纬度 # lw:代表边框宽度 # scale_length:代表比例尺长度(单位为经度例如1.5个经度) # scale_lw:代表比例尺宽度 for y in [extent[2],extent[3]] :#画上下两边框 xmin = extent[0] while (xmin |
今日新闻 |
点击排行 |
|
推荐新闻 |
图片新闻 |
|
专题文章 |
CopyRight 2018-2019 实验室设备网 版权所有 win10的实时保护怎么永久关闭 |