科学计算和可视化

1、numpy ,matplotlib的笔记：

 

 2、python123成绩

import numpy as np
import matplotlib.pyplot as plt
 
# 中文和负号的正常显示
plt.rcParams['font.sans-serif'] = 'Microsoft YaHei'
plt.rcParams['axes.unicode_minus'] = False
 
# 使用ggplot的绘图风格
plt.style.use('ggplot')
 
# 构造数据
values = [5,5,5,5,5,5,4]
feature = ['第一周','第二周','第三周','第四周','第五周','第六周','第七周']
 
N = len(values)
# 设置雷达图的角度，用于平分切开一个圆面
angles=np.linspace(0, 2*np.pi, N, endpoint=False)
# 为了使雷达图一圈封闭起来，需要下面的步骤
values=np.concatenate((values,[values[0]]))
angles=np.concatenate((angles,[angles[0]]))
 
# 绘图
fig=plt.figure()
ax = fig.add_subplot(111, polar=True)
# 绘制折线图
ax.plot(angles, values, 'bo-', linewidth=2, label = '学号2019310143001')
# 填充颜色
ax.fill(angles, values, color='m',alpha=0.35)
 
# 添加每个特征的标签
ax.set_thetagrids(angles * 180/np.pi, feature)
# 设置雷达图的范围
ax.set_ylim(0,6)
# 添加标题
plt.title('ghh的成绩单')
 
# 添加网格线
ax.grid(True)
# 设置图例
plt.legend(loc = 'best')

# 显示图形
plt.show()

结果：

3、手绘：

from PIL import Image
import numpy as np
vec_el=np.pi/2.2
vec_az=np.pi/4
depth=10
im=Image.open('ooj.jpg').convert('L')
a=np.asarray(im).astype('float')
grad=np.gradient(a)
grad_x,grad_y=grad
grad_x=grad_x*depth/100
grad_y=grad_y*depth/100
dx=np.cos(vec_el)*np.cos(vec_az)
dy=np.cos(vec_el)*np.sin(vec_az)
dz=np.sin(vec_el)
A=np.sqrt(grad_x**2+grad_y**2+1.)
uni_x=grad_x/A
uni_y=grad_y/A
uni_z=1./A
a2=255*(dx*uni_x+dy*uni_y+dz*uni_z)
a2=a2.clip(0,255)
im2=Image.fromarray(a2.astype('uint8'))
im2.save('ddd.jpg')

原图：

 

 结果：