通过随机数,生成48角度的隐式三维空间的对称平面w=ws()
# -*- coding: utf-8 -*-
"""
Created on Mon Jan 11 14:19:19 2021
@author: LZQ
"""
import numpy as np
import math
v1=np.array([-1,0,0])
v2=np.array([-1,-1,1])
v=np.array([0,0,1])
v,v1,v2=v[:,None],v1[:,None],v2[:,None]
print("------------indices--------------")
indices = np.linspace(1, 48, 48)
print(indices)
print("------------phi--------------")
alpha=0.5*np.pi
phi = np.arccos(1 + (math.cos(alpha) - 1) * indices / 48)
print(phi)
print("------------theta--------------")
theta = np.pi * (1 + 5 ** 0.5) * indices
print(theta)
print("------------r--------------")
r = np.sin(phi)
print(r)
print("------------w--------------")
w=(v*np.cos(phi)+r*(v1*np.cos(theta)+v2*np.sin(theta))).T
print(w)
print(w.shape)
print("-----------ws--------------")
ws=w/w[:,2:]
print(ws)
print(ws.shape)
输出结果展示:
runfile('C:/Users/mayn/Desktop/test0.py', wdir='C:/Users/mayn/Desktop')
------------indices--------------
[ 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18.
19. 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36.
37. 38. 39. 40. 41. 42. 43. 44. 45. 46. 47. 48.]
------------phi--------------
[0.2044802 0.28968699 0.3554212 0.41113786 0.46049343 0.50536051
0.54685017 0.58568554 0.62236849 0.65726415 0.69064807 0.72273425
0.75369289 0.78366205 0.81275556 0.84106867 0.86868205 0.89566479
0.92207665 0.94796974 0.97338991 0.99837775 1.02296948 1.04719755
1.07109129 1.09467727 1.11797973 1.1410209 1.1638212 1.18639955
1.2087735 1.23095942 1.25297262 1.27482752 1.29653771 1.31811607
1.33957485 1.36092573 1.38217994 1.40334825 1.42444107 1.4454685
1.46644035 1.48736624 1.50825556 1.52911759 1.54996149 1.57079633]
------------theta--------------
[ 10.16640738 20.33281477 30.49922215 40.66562954 50.83203692
60.99844431 71.16485169 81.33125908 91.49766646 101.66407385
111.83048123 121.99688862 132.163296 142.32970338 152.49611077
162.66251815 172.82892554 182.99533292 193.16174031 203.32814769
213.49455508 223.66096246 233.82736985 243.99377723 254.16018462
264.326592 274.49299939 284.65940677 294.82581415 304.99222154
315.15862892 325.32503631 335.49144369 345.65785108 355.82425846
365.99066585 376.15707323 386.32348062 396.489888 406.65629539
416.82270277 426.98911015 437.15551754 447.32192492 457.48833231
467.65473969 477.82114708 487.98755446]
------------r--------------
[0.20305822 0.28565228 0.34798527 0.39965263 0.44439019 0.48412292
0.51999933 0.5527708 0.58296119 0.61095326 0.63703687 0.66143783
0.68433615 0.7058781 0.72618438 0.74535599 0.76347843 0.78062475
0.79685798 0.81223286 0.82679728 0.84059338 0.85365839 0.8660254
0.87772393 0.88878038 0.89921841 0.90905934 0.91832237 0.92702481
0.93518232 0.94280904 0.94991776 0.95652002 0.96262625 0.96824584
0.9733872 0.97805788 0.9822646 0.9860133 0.98930917 0.99215674
0.99455986 0.99652173 0.99804496 0.99913157 0.99978296 1. ]
------------w--------------
[[ 0.28689266 0.13716385 0.84200281]
[-0.30953188 -0.28455852 1.24289186]
[ 0.06443361 0.27616137 0.66133863]
[ 0.32393106 -0.06961227 0.98627894]
[-0.61347325 -0.23851668 1.13435001]
[ 0.59320519 0.4675248 0.4074752 ]
[-0.22180125 -0.46147259 1.31563926]
[-0.32960765 0.18962174 0.6437116 ]
[ 0.76129016 0.22243258 0.59006742]
[-0.8123114 -0.55336131 1.34502798]
[ 0.41718287 0.60783772 0.16299561]
[ 0.24063348 -0.33164994 1.08164994]
[-0.81531481 -0.14694027 0.87610694]
[ 0.98342362 0.57745235 0.13088098]
[-0.62684049 -0.72016295 1.40766295]
[-0.08959329 0.48034242 0.18632425]
[ 0.79437172 0.03154527 0.61428806]
[-1.10396682 -0.55063704 1.17563704]
[ 0.83281544 0.79600742 -0.19184075]
[-0.10321412 -0.62361914 1.20695248]
[-0.70916257 0.11025093 0.45224907]
[ 1.17009773 0.48008965 0.06157702]
[-1.02020548 -0.83284337 1.3536767 ]
[ 0.32083301 0.75140427 -0.25140427]
[ 0.56943025 -0.26677101 0.74593767]
[-1.17960146 -0.37277451 0.83110785]
[ 1.17665814 0.8295049 -0.3920049 ]
[-0.54773669 -0.85540989 1.27207656]
[-0.38565745 0.42723181 -0.03139848]
[ 1.13269472 0.23628953 0.13871047]
[-1.29238894 -0.78660705 1.14077372]
[ 0.76949746 0.92918684 -0.59585351]
[ 0.16938948 -0.58163752 0.89413752]
[-1.03222931 -0.07897516 0.37064182]
[ 1.36032973 0.70663012 -0.43579679]
[-0.97295769 -0.96823432 1.21823432]
[ 0.06694 0.72094442 -0.49177775]
[ 0.88362485 -0.09025937 0.2985927 ]
[-1.37627024 -0.5938381 0.7813381 ]
[ 1.14657141 0.97009121 -0.80342455]
[-0.31040364 -0.83731524 0.98314857]
[-0.69463366 0.26223917 -0.13723917]
[ 1.33893478 0.4540885 -0.34992183]
[-1.28078304 -0.93437533 1.01770866]
[ 0.54807963 0.92438515 -0.86188515]
[ 0.47498363 -0.42788746 0.46955413]
[-1.24996703 -0.29456168 0.31539501]
[ 1.36837726 0.86275491 -0.86275491]]
(48, 3)
-----------ws--------------
[[ 0.34072649 0.16290189 1. ]
[ -0.24904168 -0.22894874 1. ]
[ 0.09742908 0.41757938 1. ]
[ 0.32843757 -0.07058072 1. ]
[ -0.54081478 -0.21026727 1. ]
[ 1.45580687 1.14736995 1. ]
[ -0.16858819 -0.35075921 1. ]
[ -0.51204242 0.29457561 1. ]
[ 1.29017488 0.37696129 1. ]
[ -0.60393644 -0.41141249 1. ]
[ 2.55947303 3.72916623 1. ]
[ 0.22246891 -0.30661486 1. ]
[ -0.93061107 -0.16771956 1. ]
[ 7.51387719 4.41204174 1. ]
[ -0.44530581 -0.51160183 1. ]
[ -0.48084612 2.57799202 1. ]
[ 1.29315833 0.05135257 1. ]
[ -0.93903712 -0.46837333 1. ]
[ -4.34118109 -4.1493135 1. ]
[ -0.08551631 -0.51668906 1. ]
[ -1.56807966 0.24378365 1. ]
[ 19.00218279 7.79657207 1. ]
[ -0.7536552 -0.6152454 1. ]
[ -1.27616369 -2.98882856 1. ]
[ 0.7633751 -0.35763177 1. ]
[ -1.41931214 -0.44852725 1. ]
[ -3.00164138 -2.11605747 1. ]
[ -0.43058469 -0.67245158 1. ]
[ 12.28267893 -13.60676742 1. ]
[ 8.16589201 1.70347293 1. ]
[ -1.13290561 -0.6895382 1. ]
[ -1.29142053 -1.55942162 1. ]
[ 0.18944455 -0.65050119 1. ]
[ -2.78497795 -0.21307676 1. ]
[ -3.12147718 -1.62146702 1. ]
[ -0.79866219 -0.79478496 1. ]
[ -0.1361184 -1.46599641 1. ]
[ 2.95929819 -0.30228257 1. ]
[ -1.76142728 -0.76002706 1. ]
[ -1.42710527 -1.20744533 1. ]
[ -0.31572404 -0.85166704 1. ]
[ 5.06148242 -1.91081867 1. ]
[ -3.82638249 -1.29768553 1. ]
[ -1.25849675 -0.91811671 1. ]
[ -0.63590796 -1.07251546 1. ]
[ 1.01156311 -0.91126334 1. ]
[ -3.96317943 -0.93394527 1. ]
[ -1.58605561 -1. 1. ]]
(48, 3)
In [107]:
三维重建及计算机视觉交流群:
代码技术相关交流群,不定期分享pdf技术相关书籍(全部免费) !