如何通过np.meshgrid快速形成网格坐标值（对）？

东方耀

如何通过np.meshgrid快速形成网格坐标值（对）？

1. # -*- coding: utf-8 -*-
   
2. __author__ = u'东方耀 微信：dfy_88888'
   
3. __date__ = '2019/11/18 8:14'
   
4. __product__ = 'PyCharm'
   
5. __filename__ = 'dfy_test_demo'
   
6. 
   
7. import numpy as np
   
8. 
   
9. cx = np.linspace(11, 30, 20)
   
10. cy = np.linspace(1, 10, 10)
    
11. print('所有x轴上的点的值：', cx, '共有多少个：', len(cx))
    
12. print('所有y轴上的点的值：', cy, '共有多少个：', len(cy))
    
13. 
    
14. cx_grid, cy_grid = np.meshgrid(cx, cy)
    
15. print('开始形成网格meshgrid的坐标值（对）:')
    
16. print(cx_grid.shape)
    
17. print(cx_grid[:2, :])
    
18. 
    
19. print(cy_grid.shape)
    
20. print(cy_grid[:2, :])
    
21. print('x轴与y轴点相互可形成的坐标总数为：', cx_grid.size)
    
22. print('x轴与y轴点相互可形成的坐标总数为：', cy_grid.size)
    
23. print('第一种方法：')
    
24. coords = np.dstack((cx_grid.flat, cy_grid.flat))[0]
    
25. print(coords)
    
26. print(coords.shape)
    
27. print('第二种方法：开始打印出每个坐标的具体值（对）：')
    
28. cx_grid_vector = np.reshape(cx_grid, (-1, ))
    
29. cy_grid_vector = np.reshape(cy_grid, (-1, ))
    
30. for x, y in zip(cx_grid_vector, cy_grid_vector):
    
31.     print('横坐标为：%d, 纵坐标为：%d' % (x, y))
    
32. 
    
33. 
    
34. cx_grid = np.expand_dims(cx_grid, -1)  # 为了如下的 np.tile() 做准备
    
35. cy_grid = np.expand_dims(cy_grid, -1)  # 为了如下的 np.tile() 做准备
    
36. print(cx_grid.shape)
    
37. print(cy_grid.shape)
    
38. # np.tile(A,reps)贴瓷砖：通过在对应的维度axis上重复rep给定的次数来构造数组
    
39. 
    
40. cx_res = np.tile(cx_grid, reps=(2, 3, 4))
    
41. cy_res = np.tile(cy_grid, reps=(1, 1, 4))
    
42. print(cx_res.shape)
    
43. print(cy_res.shape)
    
44. 
    
45. 
    
46. 
    

复制代码