[[N.S]]

Если нарисовать большое количество 3DSold-ов и часть из них заключить в сферу, то выполняя команду "Interfere" (ACAD2006, англ.) за какие-то доли секунды можно получить количество попадающих в эту сферу объектов. Однако, набор элементов, попадающих в сферу, таким образом получить не удается.
Господа, подскажите пожалуйста:
1) Есть ли возможность вытащить такой набор при выполнении команды?
2) Какой есть сопоставимо быстрый способ получить набор объектов, попадающих в сферу из заданной точки и с заданным радиусом?
Спасибо

[Дима_]

Незнаю насчет быстрый, но первое что приходт в голову, проверить набор внутри окружности по любой из осей на предмет vla-boolean c созданной сферой.

[Евгений А.]

Есть для тел функция vla-CheckInterference.

[[N.S]]

Дима_: Можно пояснить поподробней?
Евгений А.: vla-CheckInterference не подойдет. Она работает только с двумя объектами и может возвращать объект-пересечение двух тел. Можно последовательно проверять каждое тело чертежа (из нескольких тысяч) на пересечение со сферой , но это, думаю, займет много времени.

[Puroshev]

Интересная задача, но если можно опиши всю проблему целиком.

[Елпанов Евгений]

Сообщение от [N.S] Дима_: Можно пояснить поподробней? Евгений А.: vla-CheckInterference не подойдет. Она работает только с двумя объектами и может возвращать объект-пересечение двух тел. Можно последовательно проверять каждое тело чертежа (из нескольких тысяч) на пересечение со сферой , но это, думаю, займет много времени.

В любом случае, автокад не может получить набор для сферы. Здесь придется проходить по объектам и проверять расстояние. Примерный алгоритм:
1. Получаем набор всех объектов, попадающих в куб со стороной, равной диаметру сферы.
2. проходим по всем объектам и проверяем возможность их пересечений.

Как вариант, можно пытаться ускорить этот алгоритм.
Если объекты не слишком сложные или объекты предсказуемой конфигурации, я бы предложил написать свою функцию получения расстояний от объекта до точки.

[[N.S]]

Елпанов Евгений: Про куб то я не подумал! Спасибо. Ведь можно наверно отсечь ненужные элементы кубом, а потом со сферой разбираться..)
Вообщето мне нужно отобрать попадающие в сферу элементы типа линий (0 . "*LINE"), и блоков (0 . "INSERT"). Через интерференцию тел - это извращенный способ, который мне первый в голову пришел.)
Так. Какой-бы теперь такой хитрый фильтр для ssget придумать, чтобы сразу выбрать примитивы попадающие в куб?

Собственно говоря задача про попадание в сферу полилиний всего чертежа уже мной решена:

Вот кривенькая ф-ция координат полилайна

Код:

[Выделить все]

;|
Получить координаты полилайна 
Параметры:
	ent	| ENAME	| имя примитива, координаты которого нужно получить
Возвращает:
	res	| LIST	| список координат в WCS вида ((x1 y1 z1) (x2 y2 z2) ...)
|;
(defun polyline-ends
       		(ent / 	obj acad_doc utility e n lst ptlst trpt res
		 ;|local func|; _group-by-num _pline-get-verts)

  
;|
Сгруппировать элементы lst по num элементов
Пример:
	(_group-by-num '(1 2 3 4) 2)
	((1 2) (3 4))
|;
  	(defun _group-by-num (lst num / ls ret)
		(if
	  		(= (rem (length lst) num ) 0)
 			(repeat
				(/ (length lst) num)
  				(repeat
					num
					(setq 	ls	(cons (car lst) ls)
						lst 	(cdr lst)
					)
				) ; _end repeat
  				(setq 	ret
					(append	ret
						(list (reverse ls))
					)
					ls 	nil
				) ; _end setq
			) ; _end repeat
		) ; _end if
  		ret
	) ; _end defun _group-by-num
;|
Получить список 2D или 3D точек в зависимости от типа полилинии
|;
	(defun _pline-get-verts (pline_obj / verts)
		(setq 	verts 	(vlax-get pline_obj 'Coordinates)
	     		verts
      			(cond
		  	  (	(wcmatch	
			 	 	(vlax-get pline_obj 'Objectname)
					"AcDb2dPolyline,AcDb3dPolyline"
				)
         			(_group-by-num verts 3)
		  	  )
        	 	  (	(eq
			  		(vlax-get pline_obj 'Objectname)
					"AcDbPolyline"
				)
         			(_group-by-num verts 2)
		  	  )
		  	  (1 nil)
			) ; _end cond
		) ; _end setq
	) ; _end defun _pline-get-verts

  	(setq
	  	elist	(entget ent)
		obj	(vlax-ename->vla-object ent)
		acad_doc 	(vla-get-ActiveDocument (vlax-get-acad-object))
		utility		(vla-get-Utility	acad_doc)
		
		e	(vlax-get obj 'Elevation ) ;|расстояние от цента по нормали|;
		n	(vlax-get obj 'Normal)	;|вектор нормали|;
		lst	(_pline-get-verts obj)	;|первоначальный список точек|;
		ptlst	nil	
	) ; _end setq
  	(if	;|Если 2-х мерн. координаты - добавляем третью = Elevation|;
	  	(= 2 (length (car lst)))
	  	(progn
		  (foreach pt lst
		    	(setq 	ptlst
				(append	ptlst
					(list (append	pt (list e)))
				)
			) ; _end setq
		  ) ; _end foreach
		)
	) ; _end if
  	(foreach pt	;|Переводим каждую координату из OCS в WCS|;
		 ptlst
 		(setq	trpt 
			(vlax-variant-value
				(vla-TranslateCoordinates
			      		utility
					(vlax-3d-point pt)
					acOCS
					acWorld
					0
					(vlax-3d-point n)
				) ; _end vla-TranslateCoordinates
			) ; _end vlax-variant-value
			res
			(append	res
				(list (list 
					(vlax-safearray-get-element trpt 0)
					(vlax-safearray-get-element trpt 1)
					(vlax-safearray-get-element trpt 2)
				)) ; _end list list
			)
		) ; _end setq
	) ; _end foreach
  	res
) ; _end defun 
;|

(polyline-ends (car (entsel)))

|;

Вот преобразование набора в ((имя-примитива список точек) ...):

Код:

[Выделить все]

;|
Преобразовать набор в список
((<ename> <1/nil маркер блока> <pt1> [<pt2> ...]) ... )
|;
(defun ss-to-coordlist (ss / i lst ent)
 (setq i 0)
 (repeat (sslength ss)
  (setq lst
        (append lst
                (list
                 (append
                  (list (setq ent (ssname ss i)))
                  ;|Если блок - берем точку вставки|;
                  (if (= "INSERT" (cdr (assoc 0 (entget ent))))
                   (list
                    'fake	;|маркер блока|;
                    (vlax-get
                     (vlax-ename->vla-object ent)
                     'InsertionPoint
                     )
                    ) ; _end list
                   ;|Если полилайн - берем координаты концов|;
                   (append
                    (list nil)
                    (polyline-ends ent)
                    ) ; _end append
                   ) ; _end if
                  ) ; _end append
                 ) ; _end list
                ) ; _end append lst
        i	(1+ i)
        ) ; _end setq
  ) ; end repeat
 lst
 ) ; _end defun
;|
(ss-to-coordlist (ssget))
|;

Сам алгоритм проверки попадания в сферу и тест времени выполнения:

Код:

[Выделить все]

(defun in-sphere
       (r cpt B / in onborder out in-p _radius-vector)
  (defun _radius-vector (pt / x y z)
    (setq x (- (car pt) (car cpt))
	  y (- (cadr pt) (cadr cpt))
	  z (- (caddr pt) (caddr cpt))
	  )
    (sqrt (+ (* x x) (* y y) (* z z)))
    ) ; _end defun
  (foreach e B
    (setq in-p nil)
    (foreach pt (cddr e)
      (setq in-p
	     (append in-p
		     (list
		       (if (<= (_radius-vector pt) r) 1 nil)
		       ) ; _end list
		     ) ; _end append
	    ) ; _end setq
      ) ; _end foreach
    (if (apply 'and in-p)
      (setq in (append in (list e)))
      (if (apply 'or in-p)
	(setq onborder (append onborder (list e)))
	(setq out (append out (list e)))
	)
      )
    ) ; _end foreach
  (list in onborder out)
  ) ; _end defun

;; Тест времени
(defun ts (l pt / ss time)
 (setq time (getvar "CDATE"))
 (setq ss (ssget "_X" '((0 . "LWPOLYLINE"))))
 (setq lst (ss-to-coordlist ss))
 (in-sphere l pt lst)
 (* (- (getvar "CDATE") time) 1e6)
 ) ; _end defun

;|
(setq pt (getpoint))
(ts 10000 pt)
|;

У меня работает на вложенном файле 1.5 сек. - что конечно много. Сама in-sphere работает 0.3 сек., а тормозят кривые преобразования координат.. Надо ускорять!

[kp+]

"Что тут хитрить?.." (М. Ю. Лермонтов)

Найти с помощью GETBOUNDINGBOX ограничивающие "кубы" сферы и всех остальных объектов в модели, сравнить результаты. По-моему, этот метод работает достаточно быстро.
Теми, чей "куб" хоть частично попадает в "куб" сферы, заняться вплотную (с помощью довольно тормозных методов CheckInterference и Intersectwith)

[Евгений А.]

Сообщение от [N.S] Евгений А.: vla-CheckInterference не подойдет. Она работает только с двумя объектами и может возвращать объект-пересечение двух тел. ...

она ещё может и не возвращать объект пересечения, а только флаг пересечения, так на порядок быстрее (читай справку по команде).

[Supermax]

Метод IntersectWith. Все графические элементы кроме Pviewport и PoligonMesh.
В книге Полещука Н.Н. "AutoLisp и Visual Lisp в среде AutoCAD" Есть функция ent-int (стр. 878). С ее помощью все и делается.
Строите хоть сферу (3D солид), хоть куб, хоть звезду и находите пересечение с этим телом этой функцией.
А PoligonMesh можно и по частям проверять.

[Евгений А.]

С ее помощью все и делается.

Ан нет!
Метод IntersectWith для тел хоть и приписан, но на самом деле не реализован (может в 2010 версии только, но вряд ли), так что не так всё просто с телами.