用矩阵旋转图元的方法，这样么写对么

brainstorm

刚接触矩阵，试着写了一个动态拖动旋转的联系，总是感觉写的太繁琐求指导一下，有什么其他好方法。我知道有封装的函数xdrx_drag_rotate
只是针对这个练习

1. (defun c:tt ()
   
2.   (defun callback (p)
   
3.     (if        (not p1)
   
4.       (setq p1 p)
   
5.     )
   
6.     (if        (not an0)
   
7.       (setq an0 (angle p p0))
   
8.     )
   
9.     (if        (not (equal p p1 1e-3))
   
10.       (progn
    
11.         (setq an1 (angle p1 p0))
    
12.         (setq an0 (- (angle p p0) an1))
    
13.         (setq p1 p)
    
14. 
    
15.       )                                        ; point moved
    
16.     )
    
17. 
    
18.     (setq mat (xdrx_matrix_setrotation
    
19.                 mat0
    
20.                 an0
    
21.                 '(0 0 1)
    
22.                 p0
    
23.               )
    
24.     )
    
25.     (xdrx_entity_transform e mat)
    
26.     (princ an0)
    
27.   )
    
28.   (setq mat0 (xdrx_matrix_identity 3))
    
29.   (setq        e  (car (entsel))
    
30.         p0 (getpoint)
    
31.   )
    
32.   (xdrx_pointmonitor "callback")
    
33.   (getpoint p0)
    
34.   (xdrx_pointmonitor)
    
35. )
    

(defun c:tt ()<br />   (defun callback (p)<br />     (if        (not p1)<br />       (setq p1 p)<br />     )<br />     (if        (not an0)<br />       (setq an0 (angle p p0))<br />     )<br />     (if        (not (equal p p1 1e-3))<br />       (progn<br />         (setq an1 (angle p1 p0))<br />         (setq an0 (- (angle p p0) an1))<br />         (setq p1 p)<br /> <br />       )                                        ; point moved<br />     )<br /> <br />     (setq mat (xdrx_matrix_setrotation<br />                 mat0<br />                 an0<br />                 '(0 0 1)<br />                 p0<br />               )<br />     )<br />     (xdrx_entity_transform e mat)<br />     (princ an0)<br />   )<br />   (setq mat0 (xdrx_matrix_identity 3))<br />   (setq        e  (car (entsel))<br />         p0 (getpoint)<br />   )<br />   (xdrx_pointmonitor "callback")<br />   (getpoint p0)<br />   (xdrx_pointmonitor)<br /> )<br />

newer

你这个 (setq an1 (angle p1 p0)) 有点复杂了，点顺序反了 用 (angle p0 p)) 就简洁多了，因为p0是基点。

可以简化下，另外，你的代码在UCS下还不对，试试下面代码，UCS也也行

1. (defun c:tt ()
   
2.   (defun callback (p)
   
3.     (if        (not (equal p p0 1e-3))
   
4.       (progn
   
5.         (setq an0 (angle p0 p))
   
6.         (setq an1 (- an0 oldan))
   
7.         (setq mat
   
8.                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)
   
9.         )
   
10.         (xdrx_entity_transform e mat)
    
11.         (setq oldan an0)
    
12.       )
    
13.     )
    
14.   )
    
15.   (setq mat0 (xdrx_matrix_identity 3))
    
16.   (setq        e     (car (entsel))
    
17.         xdir  (getvar "ucsxdir")
    
18.         ucsan (xdrx_vector_angle xdir)
    
19.         oldan ucsan
    
20.         p0    (getpoint "\n基点：")
    
21.         p0 (trans p0 1 0)
    
22.   )
    
23.   (xdrx_pointmonitor "callback")
    
24.   (getpoint (trans p0 0 1))
    
25.   (xdrx_pointmonitor)
    
26. )
    
27. 
    

(defun c:tt ()<br />   (defun callback (p)<br />     (if        (not (equal p p0 1e-3))<br />       (progn<br />         (setq an0 (angle p0 p))<br />         (setq an1 (- an0 oldan))<br />         (setq mat<br />                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)<br />         )<br />         (xdrx_entity_transform e mat)<br />         (setq oldan an0)<br />       )<br />     )<br />   )<br />   (setq mat0 (xdrx_matrix_identity 3))<br />   (setq        e     (car (entsel))<br />         xdir  (getvar "ucsxdir")<br />         ucsan (xdrx_vector_angle xdir)<br />         oldan ucsan<br />         p0    (getpoint "\n基点：")<br />         p0 (trans p0 1 0)<br />   )<br />   (xdrx_pointmonitor "callback")<br />   (getpoint (trans p0 0 1))<br />   (xdrx_pointmonitor)<br /> )<br /> <br />

brainstorm

谢谢 里面有些东西没明白 说了就清楚多了

newer

brainstorm 发表于 2017-9-3 12:40

1、

(defun callback(dynpt)
   ...
)

点监视器回调函数的参数，dynpt是WCS点




2、getpoint 是UCS点


3、angle 求的是和当前UCS X轴的夹角




都统一到WCS就解决了UCS下的问题了。

newer

brainstorm 发表于 2017-9-3 12:40

你试试实现下下面这个动态，地球绕太阳公转，自身自转。

1、多段线跟着鼠标走，位置是多段线的重心和鼠标点重合
2、走的同时，随着鼠标和基点连线的方向，多段线自身绕重心点旋转，每次10度
3、旋转方向同鼠标移动的时针方向，鼠标连线顺时针，自转也顺时针，逆时针，自转也逆时针。

看看矩阵怎么实现。

brainstorm

newer 发表于 2017-9-3 13:51

好 有空练习练习

brainstorm

1. (defun c:tt ()
   
2.   (defun callback (p)
   
3.     (if        (not (equal p p0 1e-3))
   
4.       (progn
   
5.         (setq an0 (angle p0 p))
   
6.         (setq an1 (- an0 oldan))
   
7.         (setq mat
   
8.                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)
   
9.         )
   
10.         (xdrx_entity_transform e mat)
    
11.         (setq pcent (xd::polyline:centroid e)
    
12.               pcent (trans pcent 1 0)
    
13.         )
    
14.         (setq mat
    
15.                (xdrx_matrix_setrotation mat0 (* 10.0 an1) '(0 0 1) pcent)
    
16.         )
    
17.         (xdrx_entity_transform e mat)
    
18.         (setq oldan an0)
    
19.       )
    
20.     )
    
21.   )
    
22.   (setq mat0 (xdrx_matrix_identity 3))
    
23.   (setq        e     (car (entsel))
    
24.         p0    (getpoint "\n基点：")
    
25.         p0    (trans p0 1 0)
    
26.         pcent (xd::polyline:centroid e)
    
27.         pcent (trans pcent 1 0)
    
28.         oldan (angle p0 pcent)
    
29.   )
    
30.   (xdrx_pointmonitor "callback")
    
31.   (getpoint (trans p0 0 1))
    
32.   (xdrx_pointmonitor)
    
33. )
    

(defun c:tt ()<br />   (defun callback (p)<br />     (if        (not (equal p p0 1e-3))<br />       (progn<br />         (setq an0 (angle p0 p))<br />         (setq an1 (- an0 oldan))<br />         (setq mat<br />                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)<br />         )<br />         (xdrx_entity_transform e mat)<br />         (setq pcent (xd::polyline:centroid e)<br />               pcent (trans pcent 1 0)<br />         )<br />         (setq mat<br />                (xdrx_matrix_setrotation mat0 (* 10.0 an1) '(0 0 1) pcent)<br />         )<br />         (xdrx_entity_transform e mat)<br />         (setq oldan an0)<br />       )<br />     )<br />   )<br />   (setq mat0 (xdrx_matrix_identity 3))<br />   (setq        e     (car (entsel))<br />         p0    (getpoint "\n基点：")<br />         p0    (trans p0 1 0)<br />         pcent (xd::polyline:centroid e)<br />         pcent (trans pcent 1 0)<br />         oldan (angle p0 pcent)<br />   )<br />   (xdrx_pointmonitor "callback")<br />   (getpoint (trans p0 0 1))<br />   (xdrx_pointmonitor)<br /> )<br />

newer

brainstorm 发表于 2017-9-4 09:06

能转了，但是中心没移动到鼠标点



在试试改改， 平移和旋转矩阵的组合应用

brainstorm

1. (defun c:tt ()
   
2.   (defun callback (p)
   
3.     (if        (not (equal p p0 1e-3))
   
4.       (progn
   
5.         (setq an0 (angle p0 p))
   
6.         (setq an1 (- an0 oldan))
   
7.         (setq mat
   
8.                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)
   
9.         )
   
10.         (xdrx_entity_transform e mat)
    
11.         (setq pcent (xd::polyline:centroid e)
    
12.               pcent (trans pcent 1 0)
    
13.         )
    
14.         (setq mat
    
15.                (xdrx_matrix_setrotation mat0 (* 10.0 an1) '(0 0 1) pcent)
    
16.         )
    
17.        
    
18.         (xdrx_entity_transform e mat)
    
19.         (setq mat
    
20.              (xdrx_matrix_settranslation mat0 (mapcar '- p pcent ))
    
21.         )
    
22.         (xdrx_entity_transform e mat)
    
23.         (setq oldan an0)
    
24.       )
    
25.     )
    
26.   )
    
27.   (setq mat0 (xdrx_matrix_identity 3))
    
28.   (setq        e     (car (entsel))
    
29.         p0    (getpoint "\n基点：")
    
30.         p0    (trans p0 1 0)
    
31.         pcent (xd::polyline:centroid e)
    
32.         pcent (trans pcent 1 0)
    
33.         oldan (angle p0 pcent)
    
34.   )
    
35.   (xdrx_pointmonitor "callback")
    
36.   (getpoint (trans p0 0 1))
    
37.   (xdrx_pointmonitor)
    
38. )
    

(defun c:tt ()<br />   (defun callback (p)<br />     (if        (not (equal p p0 1e-3))<br />       (progn<br />         (setq an0 (angle p0 p))<br />         (setq an1 (- an0 oldan))<br />         (setq mat<br />                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)<br />         )<br />         (xdrx_entity_transform e mat)<br />         (setq pcent (xd::polyline:centroid e)<br />               pcent (trans pcent 1 0)<br />         )<br />         (setq mat<br />                (xdrx_matrix_setrotation mat0 (* 10.0 an1) '(0 0 1) pcent)<br />         )<br />         <br />         (xdrx_entity_transform e mat)<br />         (setq mat<br />              (xdrx_matrix_settranslation mat0 (mapcar '- p pcent ))<br />         )<br />         (xdrx_entity_transform e mat)<br />         (setq oldan an0)<br />       )<br />     )<br />   )<br />   (setq mat0 (xdrx_matrix_identity 3))<br />   (setq        e     (car (entsel))<br />         p0    (getpoint "\n基点：")<br />         p0    (trans p0 1 0)<br />         pcent (xd::polyline:centroid e)<br />         pcent (trans pcent 1 0)<br />         oldan (angle p0 pcent)<br />   )<br />   (xdrx_pointmonitor "callback")<br />   (getpoint (trans p0 0 1))<br />   (xdrx_pointmonitor)<br /> )<br />

newer

brainstorm 发表于 2017-9-4 10:07

很完美了，就一个地方修改下，现在的回调函数里面，用了三次的矩阵变换实体，其实用一次就可以，这也是矩阵神奇的地方，也是数学神奇的地方。不管中间有多少次变换矩阵，这些变换矩阵可以矩阵左乘，组合成一个最终的变换矩阵，实体作用这个最终的变换矩阵，变换一次，实体就变换到最终的位置，这节省了很多的时间和消耗。看过不少专业的软件，中间变换好几次，每次都重新生成实体，本来几秒完成的事，用了几十倍上百倍的时间。

你用了三个矩阵，按照最终变换顺序的反向，使用xdrx_matrix_product 矩阵左乘组合成最终的变换矩阵

(setq mat (xdrx_matrix_product mat3 mat2 mat1))

试试下面的。

1. (defun c:tt ()
   
2.   (defun callback (p)
   
3.     (if        (not (equal p p0 1e-3))
   
4.       (progn
   
5.         (setq an0 (angle p0 p))
   
6.         (setq an1 (- an0 oldan))
   
7.         (setq mat1
   
8.                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)
   
9.         )
   
10.         (setq pcent (xd::polyline:centroid e)
    
11.               pcent (trans pcent 1 0)
    
12.         )
    
13.         (setq mat2
    
14.                (xdrx_matrix_setrotation mat0 (* 10.0 an1) '(0 0 1) pcent)
    
15.         )
    
16.         (setq mat3
    
17.              (xdrx_matrix_settranslation mat0 (mapcar '- p pcent ))
    
18.         )
    
19.         (setq mat (xdrx_matrix_product mat3 mat2 mat1))
    
20.         (xdrx_entity_transform e mat)
    
21.         (setq oldan an0)
    
22.       )
    
23.     )
    
24.   )
    
25.   (setq mat0 (xdrx_matrix_identity 3))
    
26.   (setq        e     (car (entsel))
    
27.         p0    (getpoint "\n基点：")
    
28.         p0    (trans p0 1 0)
    
29.         pcent (xd::polyline:centroid e)
    
30.         pcent (trans pcent 1 0)
    
31.         oldan (angle p0 pcent)
    
32.   )
    
33.   (xdrx_pointmonitor "callback")
    
34.   (getpoint (trans p0 0 1))
    
35.   (xdrx_pointmonitor)
    
36. )
    
37. 
    

(defun c:tt ()<br />   (defun callback (p)<br />     (if        (not (equal p p0 1e-3))<br />       (progn<br />         (setq an0 (angle p0 p))<br />         (setq an1 (- an0 oldan))<br />         (setq mat1<br />                (xdrx_matrix_setrotation mat0 an1 '(0 0 1) p0)<br />         )<br />         (setq pcent (xd::polyline:centroid e)<br />               pcent (trans pcent 1 0)<br />         )<br />         (setq mat2<br />                (xdrx_matrix_setrotation mat0 (* 10.0 an1) '(0 0 1) pcent)<br />         )<br />         (setq mat3<br />              (xdrx_matrix_settranslation mat0 (mapcar '- p pcent ))<br />         )<br />         (setq mat (xdrx_matrix_product mat3 mat2 mat1))<br />         (xdrx_entity_transform e mat)<br />         (setq oldan an0)<br />       )<br />     )<br />   )<br />   (setq mat0 (xdrx_matrix_identity 3))<br />   (setq        e     (car (entsel))<br />         p0    (getpoint "\n基点：")<br />         p0    (trans p0 1 0)<br />         pcent (xd::polyline:centroid e)<br />         pcent (trans pcent 1 0)<br />         oldan (angle p0 pcent)<br />   )<br />   (xdrx_pointmonitor "callback")<br />   (getpoint (trans p0 0 1))<br />   (xdrx_pointmonitor)<br /> )<br /> <br />

brainstorm

唉 我刚才正在试怎么先把矩阵变换好 然后依次tranform，我没有数学基础，基本都是靠看别人代码一点点实验 你这么一指点 我就明白了 太感谢你了

newer

brainstorm 发表于 2017-9-4 10:49

这也不太需要数学基础，知道怎么用就行了，数学的计算都有现成的函数。知道怎么算出矩阵，就那么几种，平移，缩放，旋转，镜像，对齐。然后左乘就够了。精力放到算法的实现上。

brainstorm

newer 发表于 2017-9-4 10:59

是的。
无论怎么说，这段时间你给我的帮助很大，真心感谢！

xinxirong

晓东api可以非等比吗？

brainstorm

求助，如果当前在用户坐标系下，进行矩阵变换，是要把矩阵、点都转换到用户坐标系下面么