百度360必应搜狗淘宝本站头条
当前位置:网站首页 > 技术资源 > 正文

Python利用imshow制作自定义渐变填充柱状图/colorbar

off999 2024-10-27 11:48 41 浏览 0 评论

1. 目的

在各种各样的理论计算中,常常需要绘制各种填充图,绘制完后需要加渐变填充的colorbar。可是有些软件如VMD,colorbar渲染后颜色分布有些失真,不能较准确的表达各颜色对应的数值。用ps中的渐变填充可以解决该问题,但很多电脑配置较低,不能很好的运行ps。Python也可以直接绘制colorbar,填充颜色就好。如cmap中的bwr渐变本人就比较常用。然而,有时候颜色范围是负数范围多于正数范围(如:colorbar需要表示 [-60,40]这段,蓝色表示负数,红色表示正数,白色应该在colorbar由下往上60%处),bwr渐变将white置于50%处显得不够合理,因此需要自定义填充。本文以imshow 函数来进行填充柱状图达到自定义colorbar的目的。interpolation=‘bicubic’ 可以很好的做出渐变效果。

2. 代码

# -*- coding: utf-8 -*-"""Created on Wed Dec 9 10:36:54 2020
@author: fya"""
import matplotlib.pyplot as pltimport numpy as npfrom matplotlib.colors import ListedColormap,LinearSegmentedColormapimport matplotlib as mpl
fig, ax = plt.subplots(dpi=96)ax.set(xlim=(1,10), ylim=(-0.1,101), autoscale_on=False) #创建图像范围
a = np.array([[1, 1], [2, 2], [3, 3], [4, 4], [5, 5]]) #每种渐变色分成五段(array五行),数字表示在colormap对应的深浅print(a.shape)
clist=['white','blue'] #线性变化颜色由上面array值 小到大,越小,越白,达到上白下蓝的渐变效果clist2=['red','white'] #渐变色2,用于白色到红色填充,array越小,越红,达到上红下白的效果newcmp = LinearSegmentedColormap.from_list('chaos',clist)newcmp2 = LinearSegmentedColormap.from_list('chaos',clist2)

plt.imshow(a,cmap=newcmp,interpolation='bicubic',extent=(1,10,0,60))#60%都是蓝色到白色渐变plt.imshow(a,cmap=newcmp2,interpolation='bicubic',extent=(1,10,60,100)) #白色设置在60%处
frame = plt.gca #读取当前图层ax.yaxis.tick_right #纵坐标移到右边ax.set_yticklabels(('-80','-60','-40','-20','0','20','40')) #自定义yticks显示的值,第一个label不显示frame.spines['top'].set_visible(False) #上框线不显示frame.spines['bottom'].set_visible(False)frame.spines['right'].set_visible(False)frame.spines['left'].set_visible(False)plt.xticks([]) #x坐标不要

plt.showfig.savefig('colorbar.tif',dpi=600,format='tif')print('Done!')
#N = 10#x = np.arange(N) + 0.15#y = np.random.rand(N)
#width = 0.4#for x, y in zip(x, y): #ax.imshow(a, interpolation='bicubic', extent=(x, x+width, 0, y), cmap=plt.cm.Blues_r)
#ax.set_aspect('auto')#plt.show

代码2,渐变色分100段

# -*- coding: utf-8 -*-"""Created on Wed Dec 9 10:36:54 2020
@author: fanyiang"""
import matplotlib.pyplot as pltimport numpy as npfrom matplotlib.colors import ListedColormap,LinearSegmentedColormapimport matplotlib as mplimport pandas as pdimport os
fig, ax = plt.subplots(dpi=96)ax.set(xlim=(1,10), ylim=(-0.1,101), autoscale_on=False)
#a = np.array([[1, 1], #[2, 2], #[3, 3], #[4, 4], #[5, 5]]) #每种渐变色分成五段(array五行),数字表示在colormap对应的深浅avalue=locals dfvalue=locals for i in range(1,101): avalue['a'+str(i)]=np.array([[i,i]]) #渐变色分为100段,分的更细 dfvalue['df'+str(i)]=pd.DataFrame(avalue['a'+str(i)]) #转dataframe df=dfvalue['df'+str(i)] df.to_csv("temp.csv", mode='a',header=None) #暂存csv文件,第一列会把每一次循环的index放进去df3=pd.read_csv('temp.csv',header=None)#读取csvdf3.columns=['序号','x','y']#column命名,第一列废弃df3=df3.drop('序号',axis=1)#删除第一列a=np.array(df3) #转arrayprint(df3.head)
#a=np.vstack((a1,a2,a3,a4,a5,a6,a7,a8,a9,a10))
print(a)
clist=['white','blue'] #线性变化颜色由上面array值 小到大clist2=['red','white']newcmp = LinearSegmentedColormap.from_list('chaos',clist)newcmp2 = LinearSegmentedColormap.from_list('chaos',clist2)

plt.imshow(a,cmap=newcmp,interpolation='bicubic',extent=(1,10,0,60))plt.imshow(a,cmap=newcmp2,interpolation='bicubic',extent=(1,10,60,100)) #白色设置在60%处
frame = plt.gca #读取当前图层ax.yaxis.tick_right #纵坐标移到右边ax.set_yticklabels(('-80','-60','-40','-20','0','20','40')) #自定义yticks显示的值,第一个label不显示frame.spines['top'].set_visible(False) #上框线不显示frame.spines['bottom'].set_visible(False)frame.spines['right'].set_visible(False)frame.spines['left'].set_visible(False)plt.xticks([]) #x坐标不要

plt.showfig.savefig('colorbar.tif',dpi=600,format='tif')os.remove("temp.csv") #删除临时的csv文件print('Done!')
#N = 10#x = np.arange(N) + 0.15#y = np.random.rand(N)
#width = 0.4#for x, y in zip(x, y): #ax.imshow(a, interpolation='bicubic', extent=(x, x+width, 0, y), cmap=plt.cm.Blues_r)
#ax.set_aspect('auto')#plt.show

代码3,更改方法2中要暂时存到dataframe的问题

# -*- coding: utf-8 -*-"""Created on Fri Dec 11 10:40:53 2020
@author: fanyiang"""
import matplotlib.pyplot as pltimport numpy as npfrom matplotlib.colors import ListedColormap,LinearSegmentedColormapimport matplotlib as mplimport pandas as pdimport os
fig, ax = plt.subplots(dpi=96)ax.set(xlim=(1,10), ylim=(-0.1,101), autoscale_on=False)
#a = np.array([[1, 1], #[2, 2], #[3, 3], #[4, 4], #[5, 5]]) #每种渐变色分成五段(array五行),数字表示在colormap对应的深浅avalue=locals a=[[1,1]] for i in range(2,1001): avalue['a'+str(i)]=np.array([[i,i]]) #渐变色分为100段,分的更细 a=np.vstack((a,avalue['a'+str(i)])) #直接用vstack来vertical叠加数组print(a)clist=['white','blue'] #线性变化颜色由上面array值 小到大clist2=['red','white']newcmp = LinearSegmentedColormap.from_list('chaos',clist)newcmp2 = LinearSegmentedColormap.from_list('chaos',clist2)

plt.imshow(a,cmap=newcmp,interpolation='bicubic',extent=(1,10,0,60))plt.imshow(a,cmap=newcmp2,interpolation='bicubic',extent=(1,10,60,100)) #白色设置在60%处
frame = plt.gca #读取当前图层ax.yaxis.tick_right #纵坐标移到右边ax.set_yticklabels(('-80','-60','-40','-20','0','20','40')) #自定义yticks显示的值,第一个label不显示frame.spines['top'].set_visible(False) #上框线不显示frame.spines['bottom'].set_visible(False)frame.spines['right'].set_visible(False)frame.spines['left'].set_visible(False)plt.xticks([]) #x坐标不要

plt.showfig.savefig('colorbar.tif',dpi=600,format='tif')#os.remove("temp.csv") #删除临时的csvprint('Done!')

注:该方法主要改变在于

a=[[1,1]] for i in range(2,1001): avalue['a'+str(i)]=np.array([[i,i]]) #渐变色分为100段,分的更细 a=np.vstack((a,avalue['a'+str(i)])) #直接用vstack来vertical叠加数组

3. 代码

效果1

效果2&3

相关推荐

无线网登录(无线网登录入口)

1.打开手机浏览器,访问192.168.0.1,输入用户名和密码就可以进入路由器管理界面。2.手机wifi连接路由器信号之后,会弹出对话框,需要输入用户名和密码(有的路由器只需要密码)。3.正确的输入...

恢复精灵免费版(恢复精灵app软件下载)

一般都是收费的。建议选择正规的软件。推荐你用《安易》数据恢复软件。聊天恢复精灵不需要电脑。因为聊天恢复精灵是属于线上的应用程序,只需要在电子设备里面下载安装聊天恢复精灵软件就可以直接使用不靠谱手机数据...

苹果系统官网下载地址(iphone官网ios下载)

PP助手、同步推等手机助手都可以下载已经下架的应用,这类助手有很多,一搜一大把,而且就我知道的PP助手还能选择下载历史版本,当然也有部分应用是没在商店上架的,他们是通过企业证书公布自己的应用,需要到官...

win10桌面切换(w10系统桌面切换)
win10桌面切换(w10系统桌面切换)

     win10桌面1和桌面2快速切换方法:1.按住win按键按键盘的【win】按键。2.再按Tab按键同时按下【tab】按键。3.点击桌面选项点击界面的桌面选项完成切换。具体方...

2025-12-19 23:51 off999

下载qq号安装(下载qq安装包)
下载qq号安装(下载qq安装包)

QQ不能更新和下载,可尝试一下办法:1,使用数据线,将手机和电脑连接,然后从电脑端下载,保证你不闪退,不掉线。2,将你手机桌面的QQ图标删除,彻底删除.从AppStore里下载.注:在你不联机的状态下,你可能用别的助手会闪退问题,推荐从...

2025-12-19 23:03 off999

系统之家在线一键重装系统(系统之家官方网站)
  • 系统之家在线一键重装系统(系统之家官方网站)
  • 系统之家在线一键重装系统(系统之家官方网站)
  • 系统之家在线一键重装系统(系统之家官方网站)
  • 系统之家在线一键重装系统(系统之家官方网站)
笔记本电脑上网卡多少钱一个月

收费有好几种!1.按流量收费,这适合不怎么上网看电影,下载的人.(想用就用,而且价格便宜)2.按小时收费,用电脑的无线连移动的无线,那有说明。(10元40小时/每月,20元100小时/每月等等,还是手...

电脑无法进入系统如何重装(电脑无法进入系统如何重装系统)
  • 电脑无法进入系统如何重装(电脑无法进入系统如何重装系统)
  • 电脑无法进入系统如何重装(电脑无法进入系统如何重装系统)
  • 电脑无法进入系统如何重装(电脑无法进入系统如何重装系统)
  • 电脑无法进入系统如何重装(电脑无法进入系统如何重装系统)
一键重装系统后一直黑屏(电脑一键重装系统后黑屏)

黑屏故障一般有四种情况:一,全黑。主机、显示器(包括指示灯)均不亮二,显示器的指示灯亮,主机的指示灯不亮三,显示器与主机的指示灯均亮且开关电源冷却风扇也正常旋转,但显示器不显示。四,动态黑屏,开机时显...

激活windows10企业版ltsc(激活Windows10企业版转到设置以激活Windows)

windows10企业版ltsc永久激活方法/步骤1、下载密钥采集器,打开采集器,选择密钥的版本,选好后点开始采集。2、采集完毕后,点击采集到的密钥进行复制,粘贴到密钥擢爻充种的输入窗口里,点击下一步...

win10设置每天定时关机命令(win10设置每天自动关机时间)

首先按【Win和R】键打开运行框,输入【shutdown-t-s600】;-s-t及600前面均有一个空格,其中的数字代表的是时间,单位为秒;如600即代表10分钟,点击【确定】完成设置,此时...

win7旗舰版永久激活码怎么获取

一、在线获取激活密钥1、访问官方网站:打开浏览器,访问微软官方网站。2、注册账号:如果没有微软账号,需要先注册一个账号。3、登录账号:使用注册的账号登录微软官方网站。4、获取密钥:在官方网站上找到wi...

路由器恢复出厂设置怎么办

家里路由器重置以后,需要重新设置宽带网络连接和建立WiFi网络,安装设置方法∶然后打开浏览器,输入192.168.1.1,一般是这个网站,不是的话就看路由器说明。输入用户名admin,密码admin登...

电脑怎么恢复到上一次系统(电脑怎么恢复到之前的系统)

前提你得有备份,用备份还原就可以了,  电脑还原初始状态的步骤如下:  1、将电脑关机然后开机或者直接点击重启,然后按住"DELETE"键,这时,电脑会自动进入到BIOS  2、电脑屏幕上会显示两...

路由器品牌型号(路由器品牌型号在哪查)

其实关于路由器的排名,随便百度一下大把都是,在此我就不再赘述了。但是关于路由器的选择上,我个人的观点是如果家里对不怎么打游戏,房子户型也不太复杂,那么200快钱的小米,华为,TP等等市面上所有这个价位...

取消回复欢迎 发表评论: