Topic: callibrate a 3dPolyline

[cadplayer]

Hello!

Iīm wondering if somebody have thought about how I can  modify a 3dPolyline like picture "before" to two 3dPolyline like picture "after".
It coming from surfeyor as he gones sections of a road with same pointCodes. My function draws a 3dPoly after direction of pointNumber. But in this point I need propose improvements.
Anybody a idea how to begin ?

[Lee Mac]

I don't see any images Dirk...

[cadplayer]

oh sorry I forgot it ...
white line shows the direction as Polyline was measured.

[CAB]

Still on my second cup of coffee so maybe I'm not awake yet, but are we looking at a before profile (bottom purple line)
and an after profile (top purple line)?

What pline are you trying to create?

[cadplayer]

No itīs not a profile itīs in plan view I mean the man out measured one 3dPolyline in this form as white line shows. Now I want to convert this Polyline (white) in two Polylines (magenta) so it will be two straight Polylines.
I have a try with select Polyline and compare angle between first and next segment. If the angle is bigger than 10grads third Polylinevertex erased. Can you follow me?

PS.(thanks Gile for your great stringSplitter)

Code: [Select]

;;; (sublist '(1 2 3 4 5 6) 2 2) -> (3 4)
(defun sublist (lst start leng / n r)
  (if (or (not leng) (< (- (length lst) start) leng))
    (setq leng (- (length lst) start))
  )
  (setq n (+ start leng))
  (repeat leng
    (setq r (cons (nth (setq n (1- n)) lst) r))
  )
)
;; (split-list '(1 2 3 4 5 6 7 8) 2) -> ((1 2) (3 4) (5 6) (7 8))
(defun split-list (lst n)
  (if lst
    (cons (sublist lst 0 n)
  (split-list (sublist lst n nil) n)
    )
  )
)
(defun c:foo (/ z ent)
  (while (= z nil)
    (if (setq ent (entsel "\nSelect Polyline" ))
      (if (or (= "POLYLINE" (cdr (assoc 0 (entget (car ent)))))
      (= "LWPOLYLINE" (cdr (assoc 0 (entget (car ent)))))
   )
(progn
  (setq z 0)
  (setq coordList (split-list (vlax-get (vlax-ename->vla-object (car ent)) 'Coordinates) 3))
)
(princ "\n No Polyline selected")

      )
      (setq z nil)
    )
   )
  (setq z nil)
  (princ)
  )

(defun checkSegments ( p1 p2 p3)
  (if (equal (angle p1 p2) (angle p2 p3) (* pi (/ 10 180)))
    (setq lista  (append (list lista) (list p1) (list p2) (list p3)))
    (setq lista  (append (list lista) (list p1) (list p2)))
    )
  )
(defun c:test ()
  (while (/= coordList nil)
    (checkSegments
      (car coordList)
      (cadr coordList)
      (caddr coordList)
      )
    (setq coordList (cdr coordList))(setq coordList (cdr coordList))(setq coordList (cdr coordList))
    )
  )


[cadplayer]

Sorry was a error

Code: [Select]

(defun checkSegments ( p1 p2 p3)
  (if (equal (angle p1 p2) (angle p2 p3) (* pi (/ 10 180)))
    (setq lista  (append lista (list p1) (list p2) (list p3)))
    (setq lista  (append lista (list p1) (list p2)))
    )
  (setq mspace (vla-get-ModelSpace
(vla-get-ActiveDocument
   (vlax-get-acad-object)
   )
)
  )
  (vlax-invoke mspace 'Add3dPoly (apply 'append (reverse lista)))
  (princ)
  )


[cadplayer]

Take a look on this video, maby better than my explanation

[cadplayer]

my idea is to check segment (a) and follow segment (b), if the follow segment(b) goes ca. same angle as segment(a) it should draw a new 3dPolyline otherwise it should stop drawing and check next follow segment(c). Endresult will be a 3dPolyline that goes straight along so I did it in video.

[ribarm]

This worked on my example - see attachment...

Code - Auto/Visual Lisp: [Select]

1. (defun c:reconst3dpoly (/ 3dpl fuzzang fuzzangr initang pt ptlst pt1lst pt2lst ptlst1 ptlst11 ptlst2 ss tst tt vert)
2.   (prompt "\nPick zig-zag reference 3d polyline")
3.   (while (not ss)
4.     (setq ss (ssget "_+.:E:S:L" (list '(0 . "POLYLINE") '(70 . 8))))
5.     (if (not ss) (prompt "\nMissed, try again"))
6.   )
7.   (setq 3dpl (ssname ss 0))
8.   (setq vert 3dpl)
9.   (while (and (setq vert (entnext vert)) (eq (cdr (assoc 0 (entget vert))) "VERTEX"))
10.     (setq pt (cdr (assoc 10 (entget vert))))
11.     (setq ptlst (cons pt ptlst))
12.   )
13.   (setq ptlst (reverse ptlst))
14.   (initget 7)
15.   (setq fuzzang (getreal "\nInput fuzz tolerance angle in decimal degrees for recreating new 3d polyines of reference 3d polyline vertices: "))
16.   (setq fuzzangr (* (/ fuzzang 180.0) pi))
17.   (setq initang (angle (car ptlst) (cadr ptlst)))
18.   (setq pt2lst ptlst)
19.   (foreach pt1 ptlst
20.     (setq pt1lst (cons pt1 pt1lst))
21.     (foreach pt pt1lst (setq pt2lst (vl-remove pt pt2lst)))
22.     (if (and ptlst1 (equal pt1 (car ptlst1) 1e-6)) (setq tt nil))
23.     (foreach pt2 pt2lst
24.       (if (and (not tt) (equal (angle pt1 pt2) initang fuzzangr)) (progn (setq initang (angle pt1 pt2)) (setq ptlst1 (cons pt1 ptlst1) ptlst1 (cons pt2 ptlst1) tt t)))
25.     )
26.   )
27.   (setq ptlst1 (reverse ptlst1))
28.   (setq ptlst1 (acet-list-remove-duplicates ptlst1 1e-6))
29.   (setq initang (angle (car ptlst1) (cadr ptlst1)))
30.   (setq tst (mapcar '(lambda (a b) (if (equal (angle a b) initang fuzzangr) (progn (setq initang (angle a b)) T) nil)) ptlst1 (cdr ptlst1)))
31.   (while (eq (car tst) T)
32.     (setq ptlst11 (cons (car ptlst1) ptlst11))
33.     (setq ptlst1 (cdr ptlst1) tst (cdr tst))
34.   )
35.   (setq ptlst11 (cons (car ptlst1) ptlst11))
36.   (setq ptlst1 ptlst11)
37.   (setq ptlst2 ptlst)
38.   (foreach pt ptlst1
39.     (setq ptlst2 (vl-remove pt ptlst2))
40.   )
41.   (entmake
42.     (list
43.       '(0 . "POLYLINE")
44.       '(100 . "AcDbEntity")
45.       '(100 . "AcDb3dPolyline")
46.       '(66 . 1)
47.       '(62 . 1)
48.       '(10 0.0 0.0 0.0)
49.       '(70 . 8)
50.       '(210 0.0 0.0 1.0)
51.     )
52.   )
53.   (foreach pt ptlst1
54.     (entmake
55.       (list
56.         '(0 . "VERTEX")
57.         '(100 . "AcDbEntity")
58.         '(100 . "AcDbVertex")
59.         '(100 . "AcDb3dPolylineVertex")
60.         (cons 10 pt)
61.         '(70 . 32)
62.       )
63.     )
64.   )
65.   (entmake
66.     (list
67.       '(0 . "SEQEND")
68.       '(100 . "AcDbEntity")
69.     )
70.   )
71.   (entmake
72.     (list
73.       '(0 . "POLYLINE")
74.       '(100 . "AcDbEntity")
75.       '(100 . "AcDb3dPolyline")
76.       '(66 . 1)
77.       '(62 . 1)
78.       '(10 0.0 0.0 0.0)
79.       '(70 . 8)
80.       '(210 0.0 0.0 1.0)
81.     )
82.   )
83.   (foreach pt ptlst2
84.     (entmake
85.       (list
86.         '(0 . "VERTEX")
87.         '(100 . "AcDbEntity")
88.         '(100 . "AcDbVertex")
89.         '(100 . "AcDb3dPolylineVertex")
90.         (cons 10 pt)
91.         '(70 . 32)
92.       )
93.     )
94.   )
95.   (entmake
96.     (list
97.       '(0 . "SEQEND")
98.       '(100 . "AcDbEntity")
99.     )
100.   )
101.   (princ)
102. )
103.  

M.R.

[cadplayer]

Thanks for big effort Marco, I test on your example, but I donīt get it. I donīt realy understand which fuzzy I should typing.

[ribarm]

For my example it's 15-20 degree, but if top projection of imaginary 3d poly connections are close to line the fuzz angle should be smaller...

[cadplayer]

How do I know, which angle I must typing - "try and error" or can I get help. Do you mean angle of first Polylinesegment ?

My try to get angle of segment like this... itīs not helping

Code: [Select]

(setq fuzzang (getreal (strcat "\nInput fuzz tolerance angle in decimal degrees for recreating new 3d polyines of reference 3d polyline vertices: " (rtos (* (/ 180 pi) (angle (car ptlst) (cadr ptlst))) 2 (getvar "AUPREC")))))

[ribarm]

Look cadplayer, this angle is tolerance angle, not angle of starting segment... As I said if top projections of imaginary connections of resulting polylines are close to straight line, that angle should be near 0.0 degree - 0.0 isn't allowed so you must enter 1.0 or some small angle, and if projections are curving, like in my posted example then you should enter angle close to difference angles of adjacent segments of imaginary projection polyline...

[CAB]

Marko,
Running the routine I get a copy of the original less the north end segments.
Not an outside trace as requested. I used one & 5 & then 10 as the angle tolerance.

[cadplayer]

Marko, can you explain, what did it do your code. You collect all vertex and than your working with a toleranz angle from what. Is it the angle between seg1 to seg2 and seg2 to seg3 and seg3 to seg 4 ... ?
My example in last drawing didnīt work if I typed 1, 5, 10 or 15 degrees