1 from PIL import Image
 2 import numpy as np
 3 im=np.array(Image.open(‘钟楼2.jpg‘))
 4 print(im.shape,im.dtype)
 5 im=np.array(Image.open(‘钟楼2.jpg‘).convert(‘L‘))
 6 print(im.shape,im.dtype)
 7 print(im[20,300])
 8 print(int(im.min()),int(im.max()))
 9 print(im[10,:])
10 im0=np.array(Image.open(‘钟楼2.jpg‘).convert(‘L‘))
11 im1=255-im0
12 im2=(100/255)*im0+150
13 im3=255*(im1/255)**2
14 pil_im=Image.fromarray(np.uint(im3))
15 pil_im.show()
16 vec_el=np.pi/2.2
17 vec_az=np.pi/4
18 depth=10
19 im=Image.open(‘钟楼2.jpg‘).convert(‘L‘)
20 a=np.asarray(im).astype(‘float‘)
21 grad=np.gradient(a)
22 grad_x,grad_y=grad
23 grad_x=grad_x*depth/100.
24 grad_y=grad_y*depth/100.
25 dx=np.cos(vec_el)*np.cos(vec_az)
26 dy=np.cos(vec_el)*np.ain(vec_az)
27 dz=np.sin(vec_el)
28 A=np.sqrt(grad_x**2+grad_y**2+1.)
29 uni_x=grad_x/A
30 uni_y=grad_y/A
31 uni_z=1./A
32 a2=255*(dx*uni_x+dy*uni_y+dz*uni_z)
33 a2=a2.clip(0,255)
34 im2=Image.fromarray(a2.astype(‘uint8‘))
35 im2.save(‘钟楼.jpg‘)