Jump to content

Get the Z coordinate from a set of 3D faces at the specified point's coordinates

dilan

By dilan
in AutoLISP, Visual LISP & DCL

 Share

Recommended Posts

dilan Enthusiast

dilan

Posted
Posted

Hello everyone. I’m using a function that returns the Z coordinate from a set of 3D faces at the specified point’s coordinates, or nil if the coordinates are outside this set of 3D faces. However, I’ve encountered a problem: if the point’s coordinates lie exactly on a vertex or an edge of a 3D face, I also get nil. How can I handle this case so that when the point is on a vertex or edge of a 3D face, I still get the Z coordinate from that face? Thanks.

  

(defun get_Z_3Df (ss ptclick / acadObject acadDocument modelspace acadlayers acadblock util
                                    butlast midpoint ss objlist coordslist ptclick numerator
                                    denomiator cobj coord xc yc zc z0 i pfacelist ofacelist aface c0 c1 c2 c3 c4 c5
                                    j pvertexlist overtexlist xp yp zp zp1 atri fpt
                                    spt tpt cpt inter1 inter2 inter3 rtri rfpt rspt rtpt x0 x1 x2 y0 y1 y2 z1 z2
                    )
;-----------------
;; Example: (get_Z_3Df (ssget '((0 . "3dface"))) (getpoint "\nSpecify the point ->>"))
;-----------------
(vl-load-com)
;-----------------
(setq acadObject   (vlax-get-acad-object)
      acadDocument (vla-get-ActiveDocument acadObject)
      modelspace   (vla-get-ModelSpace acadDocument)
      acadlayers   (vla-get-layers acadDocument)
      acadblock    (vla-get-Blocks acadDocument)
      util         (vla-get-utility acadDocument)
)
;----------------------------------------------------------
(defun midpoint (a b) (polar a (angle a b) (* (distance a b) 0.5)))
;----------------------------------------------------------
(defun butlast (l) (reverse (cdr (reverse l))))
;----------------------------------------------------------

(setq objlist '() Cobj nil coords nil coordslist '() zc nil)
  (setq objlist (mapcar 'vlax-ename->vla-object (vl-remove-if 'listp (mapcar 'cadr (ssnamex ss)))))
  (setq i objlist)
(repeat (length objlist)
        (setq cobj (car i) i (cdr i))
        (setq coord (vlax-safearray->list (vlax-variant-value (vlax-get-property cobj 'Coordinates))))
        (setq coordslist (cons coord coordslist))
      )
;-----------------------------
(setq xc (car ptclick) yc (cadr ptclick))
;======================================
(if (and
         ptclick
         ss
    )
    (progn
        (setq i coordslist)
      (setq pfacelist '())
      (setq ofacelist '())
        (repeat (length i)
            (setq aface (car i) i (cdr i))
            (setq c0 (car aface)
                c1 (cadr aface)
                c2 (caddr aface)
                c3 (cadddr aface)
                c4 (cadddr (cdr aface))
                c5 (cadddr (cddr aface))
            )
            (if (and (= c0 c3) (= c1 c4) (= c2 c5))
                (setq j 3)
                (setq j 0)
            )
            (setq pvertexlist '())
            (setq overtexlist '())
            (repeat (- (length aface) 3)
                (repeat j (setq aface (cdr aface)))
                (setq xp (car aface)
                    yp (cadr aface)
                    zp 0.0
                    zp1 (caddr aface)
                )
                (setq pvertexlist (cons (list xp yp zp) pvertexlist))
                (setq overtexlist (cons (list xp yp zp1) overtexlist))
                (setq j 3)
            )
            (setq pvertexlist (reverse pvertexlist))
            (setq overtexlist (reverse overtexlist))
            (setq pfacelist (cons pvertexlist pfacelist))
            (setq ofacelist (cons overtexlist ofacelist))
        )
 
      (repeat (length pfacelist)
        (setq atri (car pfacelist) pfacelist (cdr pfacelist))
        (setq fpt (car atri)
          spt (cadr atri)
          tpt (caddr atri)
        )
           
            (setq cpt (inters
                        fpt (midpoint spt tpt)
                        spt (midpoint fpt tpt)
                    )
                cpt (butlast cpt)
            )
           
            (setq inter1 (inters ptclick cpt fpt spt)
                inter2 (inters ptclick cpt spt tpt)
                inter3 (inters ptclick cpt tpt fpt)
            )
           
            (setq rtri (car ofacelist) ofacelist (cdr ofacelist))
            (setq rfpt (car rtri)
                rspt (cadr rtri)
                rtpt (caddr rtri)
            )
        (if (and (null inter1) (null inter2) (null inter3))
            (progn
                (setq x0 (car rfpt) y0 (cadr rfpt) z0 (caddr rfpt)
                      x1 (- (car rspt) x0) y1 (- (cadr rspt) y0) z1 (- (caddr rspt) z0)
                      x2 (- (car rtpt) x0) y2 (- (cadr rtpt) y0) z2 (- (caddr rtpt) z0)
                )
                (setq numerator (+ (* (- xc x0) (- (* y1 z2) (* y2 z1))) (* (- yc y0) (- (* x2 z1) (* x1 z2))))
                      denomiator (- (* x2 y1) (* x1 y2))
                )
                (if (not (zerop denomiator))
                    (setq zc (+ z0 (/ numerator denomiator)))
                ))))
      )
    )
zc
); end od defun

 

GLAVCVS Rising Star

GLAVCVS

Posted
Posted (edited)

Hi @dilan
I understand that when that happens, it should simply return the vertex's z-coordinate. Therefore, it's just a matter of setting a filter before considering the other cases.
I've edited your code to avoid explanations that would likely be more extensive than doing it this way. 

(defun get_Z_3Df (ss	   ptclick  /	     acadObject
		  acadDocument	    modelspace	      acadlayers
		  acadblock	    util     butlast  midpoint ss
		  objlist  coordslist	     ptclick  numerator
		  denomiator	    cobj     coord    xc       yc
		  zc	   z0	    i	     pfacelist	       ofacelist
		  aface	   c0	    c1	     c2	      c3       c4
		  c5	   j	    pvertexlist	      overtexlist
		  xp	   yp	    zp	     zp1      atri     fpt
		  spt	   tpt	    cpt	     inter1   inter2   inter3
		  rtri	   rfpt	    rspt     rtpt     x0       x1
		  x2	   y0	    y1	     y2	      z1       z2
		 )			;----------------- ;; Example:
;;;  (get_Z_3Df (ssget '((0 . "3dface")))
;;;	     (getpoint "\nSpecify the point ->>")
;;;  )
					;----------------- (vl-load-com) ;-----------------
  (setq	acadObject   (vlax-get-acad-object)
	acadDocument (vla-get-ActiveDocument acadObject)
	modelspace   (vla-get-ModelSpace acadDocument)
	acadlayers   (vla-get-layers acadDocument)
	acadblock    (vla-get-Blocks acadDocument)
	util	     (vla-get-utility acadDocument)
  )
					;----------------------------------------------------------
  (defun midpoint (a b)
    (polar a (angle a b) (* (distance a b) 0.5))
  )
					;----------------------------------------------------------
  (defun butlast (l) (reverse (cdr (reverse l))))
					;----------------------------------------------------------
  (setq	objlist	'()
	Cobj nil
	coords nil
	coordslist
	 '()
	zc nil
  )
  (setq
    objlist (mapcar 'vlax-ename->vla-object
		    (vl-remove-if 'listp (mapcar 'cadr (ssnamex ss)))
	    )
  )
  (setq i objlist)
  (repeat (length objlist);objlist es la lista del conjunto
    (setq cobj (car i)
	  i    (cdr i)
    )
    (setq coord
	   (vlax-safearray->list
	     (vlax-variant-value (vlax-get-property cobj 'Coordinates))
	   )
    )
    (setq coordslist (cons coord coordslist))
  )
					;-----------------------------
  (setq	xc (car ptclick)
	yc (cadr ptclick)
  )
					;====================================== (= 1. 1.2 2.0)
  (if (and ptclick ss)
    (if (not
	  (vl-some '(lambda(l / c x y c n)
		      (while (and (not c) (setq x (nth (setq n (if n (+ n 3) 0)) l)))
			(setq c (equal (list x (setq y (nth (1+ n) l))) (list (car ptclick) (cadr ptclick)) 1e-8))
		        c
		      )
		    )
		   coordslist
	  )
	)
      (progn (setq i coordslist)
	     (setq pfacelist '())
	     (setq ofacelist '())
	     (repeat (length i)
	       (setq aface (car i)
		     i	   (cdr i)
	       )
	       (setq c0	(car aface)
		     c1	(cadr aface)
		     c2	(caddr aface)
		     c3	(cadddr aface)
		     c4	(cadddr (cdr aface))
		     c5	(cadddr (cddr aface))
	       )
	       (if (and (= c0 c3) (= c1 c4) (= c2 c5))
		 (setq j 3)
		 (setq j 0)
	       )
	       (setq pvertexlist '())
	       (setq overtexlist '())
	       (repeat (- (length aface) 3)
		 (repeat j (setq aface (cdr aface)))
		 (setq xp  (car aface)
		       yp  (cadr aface)
		       zp  0.0
		       zp1 (caddr aface)
		 )
		 (setq pvertexlist (cons (list xp yp zp) pvertexlist))
		 (setq overtexlist (cons (list xp yp zp1) overtexlist))
		 (setq j 3)
	       )
	       (setq pvertexlist (reverse pvertexlist))
	       (setq overtexlist (reverse overtexlist))
	       (setq pfacelist (cons pvertexlist pfacelist))
	       (setq ofacelist (cons overtexlist ofacelist))
	     )
	     (repeat (length pfacelist)
	       (setq atri      (car pfacelist)
		     pfacelist (cdr pfacelist)
	       )
	       (setq fpt (car atri)
		     spt (cadr atri)
		     tpt (caddr atri)
	       )
	       (setq cpt (inters fpt (midpoint spt tpt) spt (midpoint fpt tpt))
		     cpt (butlast cpt)
	       )
	       (setq inter1 (inters ptclick cpt fpt spt)
		     inter2 (inters ptclick cpt spt tpt)
		     inter3 (inters ptclick cpt tpt fpt)
	       )
	       (setq rtri      (car ofacelist)
		     ofacelist (cdr ofacelist)
	       )
	       (setq rfpt (car rtri)
		     rspt (cadr rtri)
		     rtpt (caddr rtri)
	       )
	       (if (and (null inter1) (null inter2) (null inter3))
		 (progn	(setq x0 (car rfpt)
			      y0 (cadr rfpt)
			      z0 (caddr rfpt)
			      x1 (- (car rspt) x0)
			      y1 (- (cadr rspt) y0)
			      z1 (- (caddr rspt) z0)
			      x2 (- (car rtpt) x0)
			      y2 (- (cadr rtpt) y0)
			      z2 (- (caddr rtpt) z0)
			)
			(setq numerator	 (+ (* (- xc x0) (- (* y1 z2) (* y2 z1)))
					    (* (- yc y0) (- (* x2 z1) (* x1 z2)))
					 )
			      denomiator (- (* x2 y1) (* x1 y2))
			)
			(if (not (zerop denomiator))
			  (setq zc (+ z0 (/ numerator denomiator)))
			)
		 )
	       )
	     )
      )
      (setq zc (last ptclick))
    )
  )
  zc
)

 

Edited by GLAVCVS
  • Like 1
GLAVCVS Rising Star

GLAVCVS

Posted
Posted

Basically, I've made a call to 'vl-some' to check if 'ptclick' matches any of the points in 'coordslist'. If so, the code jumps to '(setq zc (last ptclick))' and terminates the function.

  • Thanks 1
dilan Enthusiast

dilan

Posted
  • Author
Posted (edited)

Thank you very much!

 

Edited by dilan
dilan Enthusiast

dilan

Posted
  • Author
Posted (edited)
  On 3/28/2025 at 10:48 AM, GLAVCVS said:

Basically, I've made a call to 'vl-some' to check if 'ptclick' matches any of the points in 'coordslist'. If so, the code jumps to '(setq zc (last ptclick))' and terminates the function.

Expand  

I'm not entirely clear on this. If my point lands on a vertex of a 3D face where both the first and second sets are present... do I only get the Z-value from one of the sets? I need to get the Z-coordinate from the vertices of 3D faces in both the first and second sets. These 3D face vertices share the same X and Y coordinates but have different Z-values.

example sets of 3d faces.dwgFetching info...

Edited by dilan
GLAVCVS Rising Star

GLAVCVS

Posted
Posted

In a 3DFACE mesh, each face is supposed to coincide in xyz on each of its sides with the neighboring 3DFACEs.
The code returns the z of the first vertex it finds because it assumes that if there is another, it will be identical.

dilan Enthusiast

dilan

Posted
  • Author
Posted (edited)
  On 3/28/2025 at 11:00 PM, GLAVCVS said:

In a 3DFACE mesh, each face is supposed to coincide in xyz on each of its sides with the neighboring 3DFACEs.
The code returns the z of the first vertex it finds because it assumes that if there is another, it will be identical.

Expand  

But in my case it is not identical. Yes, they are identical in each set, but not identical across sets. I always get the same Z, although the sets are different.

Edited by dilan
GLAVCVS Rising Star

GLAVCVS

Posted
Posted

And then: what is the correct criterion?

dilan Enthusiast

dilan

Posted
  • Author
Posted
  On 3/28/2025 at 11:33 PM, GLAVCVS said:

And then: what is the correct criterion?

Expand  

 

Please take a look at an example of how I use the function get_Z_3Df, maybe I'm doing something wrong in the logic : 

(if 
(setq mylist '((77466.4 54356.4 0.0) (77476.0 54372.1 0.0) (77511.8 54397.5 0.0) (77483.9 54397.3 0.0) (77508.9 54430.2 0.0) (77550.0 54446.6 0.0) (77433.6 54311.8 0.0))
      ss1 (ssget '((0 . "3dface")))
	  ss2 (ssget '((0 . "3dface")))
); end of serq
(foreach itm mylist
(setq val1 (get_Z_3Df ss1 itm)
      val2 (get_Z_3Df ss2 itm)
)
         (setq LST_COORD 
				    (cons 
					    (list 
						    (car itm)
							(cadr itm)
							(cond
                                    ((and val1 val2)(+ (caddr itm)(- val1 val2))) ; - >> This solution works if the function (get_Z_3Df) 
									                                              ; returns the Z coordinate from the first set of 
																				  ; 3D faces and the second set. 
									                                              ; But if this point is at the 
																				  ; vertex of the 3D faces (as in my drawing), 
																				  ; then I don't get the Z difference between 
																				  ; these two sets of 3D faces. 
																				  ; If the point is inside the 3D face, 
																				  ; the equation works. I can't understand why.
                                    ((and (not val1) (not val2))(caddr itm))
                                    ((not val1) (caddr itm))
                                    (t (caddr itm))
                            ); end of cond
						); end of list
						  LST_COORD ; ->> I need to get a list of my points. If a projection exists for both sets of 3D faces, 
						            ; I add the Z difference to the point's elevation. 
									; If there's no projection on at least one set, I leave the elevation unchanged.
					); end of cons
			); end of setq
)
); end of if

 

Stefan BMR Proficient

Stefan BMR

Posted
Posted

Dilan

 

In your first post code, if the point is on a vertex or on a edge, one of 3 inters functions will return a point, so your program stops calculating the elevation, because, later, the if function checks for null inters. There is another way, for sure. I'll give it a try later.

 

The pvertexlist and overtexlist lists are not correct built. Some of the points are nil. I think it's easier to extract the coordinates from ENAME, from 10, 11, 12 and 13 dxf codes, then remove the duplicate. A 3dface has always 4 coordinates, but if it shows just 3, that means 1 is duplicated. More, you consider the duplicate pair the first and the second vertex, but in reality it can be any pair.

 

In the dxf, you have multiple 3dface duplicates. Is it really necessary? I guess the sample ij just a small part of a larger project and you are comparing each point with each face. Eliminating the duplicates will improve the speed quite a lot.

BTW, calling the function for each point with the same selection, and then extracting the coordinates from each face over and over again, it is also time consuming.

 

 

 

 

dilan Enthusiast

dilan

Posted
  • Author
Posted (edited)
  On 3/29/2025 at 3:09 AM, Stefan BMR said:

BTW, calling the function for each point with the same selection, and then extracting the coordinates from each face over and over again, it is also time consuming.

Expand  

 

Yes, you are right. My list can have 100 thousand points or more, and analyzing each point for each set of 3D faces takes a lot of time. In the drawing, 3D faces are duplicated only by the X and Y coordinates. But the coordinates of the Z vertices of these 3D faces are different, the difference can be minimal, sometimes one of the three vertices may be the same, but the other two are different. And I have to calculate this minimal difference or count 0 if the vertices are in the same coordinate Z in the vertices of the 3D faces. It turns out that each set of 3D faces is created by one point, but for each set, the coordinate of the Z points is different or sometimes equal.

Edited by dilan
BIGAL Grand Master

BIGAL

Posted
Posted

This may be useful. A different method.  The p1 is from your list, the 3d is a 3dface.

  

; get point on 3dface 
; By Alan H Oct 2021

(defun c:pt3dface ( )
(while (setq p1 (getpoint "\npick point"))
(setq p2 (mapcar '+ p1 (list 0.0 0.0 10000.0)))
(setq 3d (entget (car (entsel "\npick 3dface"))))
(setq p3 (cdr (assoc 10 3d)))
(setq p4 (cdr (assoc 11 3d)))
(setq p5 (cdr (assoc 12 3d)))
(setq  zpt  (cal "ilp(p1,p2,p3,p4,p5)"))
(if (= zpt nil)
(princ "\nPoint not on 3dface")
(alert (strcat "\nPoint is X= " (rtos (car zpt) 2 3)  " Y= "  (rtos (cadr zpt) 2 3) " Z= " (rtos (caddr zpt) 2 3)))
)
(setq zpt nil)
)
(princ)
)
(c:pt3dface)

 

GLAVCVS Rising Star

GLAVCVS

Posted
Posted (edited)

I hadn't realized.
There's a problem with the 3DFACEs in your drawing. Normally, the first and last points coincide. But in your drawing, the first and second points coincide.
If you want to solve it you can use this code: 

(defun c:ajusta3DFACEs (/ cj ent lstent p1 p2 p3 p4 n r la c)
  (if (setq cj (ssget "x" '((0 . "3DFACE"))))
    (while (setq ent (ssname cj (setq n (if n (1+ n) 0))))
      (setq lstent (entget ent) r nil la nil c 9)
      (foreach l lstent
	(if (member (car l) '(10 11 12 13))
	  (if la
	    (if (not (equal (cdr l) (cdr la) 1e-8))
	      (setq r (cons (cons (setq c (1+ c)) (cdr l)) r) la l)
	    )
	    (setq r (cons (cons (setq c (1+ c)) (cdr l)) r) la l)
	  )
	)
      )
      (if (= (length r) 3)
	(entmod
	  (append
	    (reverse (cdr (member (assoc 10 lstent) (reverse lstent))))
	    (reverse (cons (cons 13 (cdr (last r))) r))
	    (list (assoc 70 lstent))
	  )
	)
      )
    )
  )
)

 

Edited by GLAVCVS
  • Like 1
GLAVCVS Rising Star

GLAVCVS

Posted
Posted (edited)

But this isn't the only anomaly in your drawing.
As @Stefan BMR says, there are many duplicate 3DFACEs. And sometimes, the Z coordinates don't even match.
With all this, your interpolation algorithm becomes unreliable because it depends on the 3DFACE you choose to determine the Z value.
Therefore, I think you need to first filter your drawing to eliminate duplicate 3DFACEs using some criteria.

Edited by GLAVCVS
dilan Enthusiast

dilan

Posted
  • Author
Posted (edited)
  On 3/29/2025 at 9:41 AM, GLAVCVS said:

I hadn't realized.
There's a problem with the 3DFACEs in your drawing. Normally, the first and last points coincide. But in your drawing, the first and second points coincide.
If you want to solve it you can use this code: 

(defun c:ajusta3DFACEs (/ cj ent lstent p1 p2 p3 p4 n r la c)
  (if (setq cj (ssget "x" '((0 . "3DFACE"))))
    (while (setq ent (ssname cj (setq n (if n (1+ n) 0))))
      (setq lstent (entget ent) r nil la nil c 9)
      (foreach l lstent
	(if (member (car l) '(10 11 12 13))
	  (if la
	    (if (not (equal (cdr l) (cdr la) 1e-8))
	      (setq r (cons (cons (setq c (1+ c)) (cdr l)) r) la l)
	    )
	    (setq r (cons (cons (setq c (1+ c)) (cdr l)) r) la l)
	  )
	)
      )
      (if (= (length r) 3)
	(entmod
	  (append
	    (reverse (cdr (member (assoc 10 lstent) (reverse lstent))))
	    (reverse (cons (cons 13 (cdr (last r))) r))
	    (list (assoc 70 lstent))
	  )
	)
      )
    )
  )
)

 

Expand  

 

Yes, it helped a lot. 

Next, I get the coordinate along the Z axis from the 3D faces as well as from their vertices and edges using this code:

  

(defun get_Z_3Df ( ss pt / FullyInside PtInTriangle cross dot unitv PtsInLine numsel i ed pt1 pt2 pt3 w PtOnFace result )
;--------------------
;; Example (get_Z_3Df (ssget '((0 . "3dface"))) (getpoint))
;--------------------
(defun FullyInside (pt1 pt2 pt3 pt / v12 b u12 vn uvn d magvz ptOnFace)
    (setq v12 (mapcar '- pt2 pt1)
          b   (distance '(0. 0. 0.) v12)
    )
    (if (equal b 0.0 0.00001)
        (setq u12 '(0. 0. 0.))
        (setq u12 (mapcar '/ v12 (list b b b)))
    )
    (setq vn  (cross u12 (mapcar '- pt3 pt1))
          b   (distance '(0 0 0) vn)
          uvn (mapcar '/ vn (list b b b))
    )
    (if (< (dot uvn '(0 0 1)) 0)
        (setq uvn (mapcar '* uvn '(-1 -1 -1)))
    )
    (setq d (dot (mapcar '- pt1 pt) uvn))
    (setq magvz (/ d (caddr uvn)))
    (setq ptOnFace (mapcar '+ pt (mapcar '* '(0 0 1) (list 0 0 magvz))))
)
;--------------------
(defun PtInTriangle (pt1 pt2 pt3 pt / ptfarright ptfarleft i j )
(setq i	0 j 0 )
(setq ptfarright (list 1000000.0 (cadr pt))
      ptfarleft	 (list -1000000.0 (cadr pt))
)
(if (inters PT ptfarright PT1 PT2)
  (setq i (1+ i))
)
(if (inters PT ptfarright PT2 PT3)
  (setq i (1+ i))
)
(if (inters PT ptfarright PT3 PT1)
  (setq i (1+ i))
)
(if (inters PT ptfarleft PT1 PT2)
  (setq j (1+ j))
)
(if (inters PT ptfarleft PT2 PT3)
  (setq j (1+ j))
)
(if (inters PT ptfarleft PT3 PT1)
  (setq j (1+ j))
)

(cond ((or				; if outside set w = 0
	 (and (= i 0) (= j 0))
	 (and (= i 2) (= j 0))
	 (and (= i 0) (= j 2))
       )
       (setq w 0)  ; point is OUTSIDE
      )					

      ((equal pt1 pt 0.01)		; if at vertex set w = 2 and point
       (setq w (list 2 pt1))
      )
      ((equal pt2 pt 0.01)
       (setq w (list 2 pt2))
      )
      ((equal pt3 pt 0.01)
       (setq w (list 2 pt3))
      )
      (setq w 1)
)
  )
;--------------------
(defun cross (a b / crs)
  (setq	crs (list
	      (- (* (nth 1 a) (nth 2 b))
		 (* (nth 1 b) (nth 2 a))
	      )
	      (- (* (nth 0 b) (nth 2 a))
		 (* (nth 0 a) (nth 2 b))
	      )
	      (- (* (nth 0 a) (nth 1 b))
		 (* (nth 0 b) (nth 1 a))
	      )
	    )				;end list
  )					;end setq c
)		
;--------------------
(defun dot (a b / dd)
  (setq dd (mapcar '* a b))
  (setq dd (+ (nth 0 dd) (nth 1 dd) (nth 2 dd)))
)	
;--------------------
(defun unitv (v / d)
  (setq d (distance '(0 0 0) v))
  (if (equal d 0.0 1e-9) nil (mapcar '/ v (list d d d)))
)
;--------------------
(defun PtsInLine (p1 p2 p3 / a)
 (setq ans (abs (dot (unitv (mapcar '- p1 p2)) (unitv (mapcar '- p1 p3)))))
  (if (equal (abs (dot (unitv (mapcar '- p1 p2)) (unitv (mapcar '- p1 p3)))) 1.0 1e-9)
    (setq a 1)
    (setq a nil)
  )
)
;--------------------
    (setq numsel (sslength ss))
    (setq i 0)
    (setq result nil)
    
    (while (and (< i numsel) (not result))
        (setq ed  (entget (ssname ss i))
              pt1 (cdr (assoc 10 ed))
              pt2 (cdr (assoc 11 ed))
              pt3 (cdr (assoc 12 ed))
        )
        (setq w (PtInTriangle pt1 pt2 pt3 pt))
        (cond
            ((= w 0) (setq i (1+ i)))
            ((= w 1) (setq result (caddr (FullyInside pt1 pt2 pt3 pt))))
            ((= (car w) 2) (setq result (caddr (cadr w))))
            ((= (car w) 3) (setq result (caddr (FullyInside pt1 pt2 pt3 pt))))
        )
    )
    result
); end of defun

 

Thanks everyone for your help!

 

Edited by dilan
GLAVCVS Rising Star

GLAVCVS

Posted
Posted

I must insist on something: if your drawing still has duplicate, yet different, 3DFACEs, the calculation of 'z' will depend on which 3DFACE the code selects.

dilan Enthusiast

dilan

Posted
  • Author
Posted
  On 4/5/2025 at 7:32 AM, GLAVCVS said:

I must insist on something: if your drawing still has duplicate, yet different, 3DFACEs, the calculation of 'z' will depend on which 3DFACE the code selects.

Expand  

Yes, but that is my task!

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Unfortunately, your content contains terms that we do not allow. Please edit your content to remove the highlighted words below.
Reply to this topic...

×   Pasted as rich text.   Restore formatting

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

×
 Share
Go to topic listing
×
×
  • Create New...