[swkx]

Контур - замкнутая 2D-полилиния неправильной формы.
Как программно определить, лежит ли точка с координатами X и Y внутри контура ?
Нужна не готовая функция, а просто подходы к решению.
А может уже есть стандартная из серии VLA-....?

[Дима_]

Домой пора - короче нефига не понял как по одному вектору можно у НЕВЫПУКЛОЙ фигуры направление понять. В понедельник отпишусь.

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

Сообщение от Дима_ Домой пора - короче нефига не понял как по одному вектору можно у НЕВЫПУКЛОЙ фигуры направление понять. В понедельник отпишусь.

Направление обхода точек можно узнать программно - заходи на мой сайт и бери или ищи на этом форуме.
Направление сегмента можно получить из вектора направленности - результат работы функции vlax-curve-getFirstDeriv, т.е. угол между двумя точками '(0. 0.) и вектором
Направление на ближайшую точку контура, вообще просто - угол между двумя точками - тестовой и ближайшей
Угол между двумя углами, вообще просто - вычел из одного другой и получил нужный угол, только проверить чтоб был в диапазоне -180/+180 градусов...

[VVA]

>Евгений, все таки наверное алгоритм для выпуклых контуров
Ну и алгоритм поста #20 в коде (для выпуклых контуров)

Код:

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

(defun C:TEST ( )
  (vl-load-COM)
  (and
    (setq e1 (car(entsel "\nВыбери полилинию-контур: ")))
    (setq pt1 (getpoint "\nУкажи точку (не на контуре): "))
    (setq pt2 (vlax-curve-getclosestpointto e1 pt1)) ;_ см п.2 пост #20
    (setq pt3 (vlax-curve-getFirstDeriv e1 (vlax-curve-getParamAtPoint e1 pt2)));_ см п.3 пост #20
    (setq pt1 (mapcar '- pt1 pt2))
    (setq ang (3d_Wnorm pt1 pt3)) ;_ см п.4 пост #20
    (if (lib:pline_clockwise e1) ;;; если контур по часовой стрелке
      (progn
        ;;;Если бы точка была снаружи, то угол был бы отрицательным
        (if (minusp (last ang))
          (alert "Точка снаружи")
          (alert "Точка внутри")
          )
        )
      (progn  ;;; если контур против часовой стрелки
        (if (minusp (last ang))
          (alert "Точка внутри")
          (alert "Точка снаружи")
          )
        )
      )
    )
  )

;********************************
; Векторное произведение векторов
;********************************
; W1, W2 - вектора
; Возвращает: вектор нормали к плоскости заданной векторами  в правой системе координат.
;W1 и W2 не должны лежать на одной прямой).
(defun 3d_Wnorm (W1 W2)
  (if (< (length W1) 3)(setq W1 (list (car W1)(cadr W1) 0)))
  (if (< (length W2) 3)(setq W2 (list (car W2)(cadr W2) 0)))
  (list (- (* (cadr W1)(caddr W2))(* (caddr W1)(cadr W2)))
        (- (* (caddr W1)(car W2)) (* (car W1)(caddr W2)))(- (* (car W1)(cadr W2)) (* (cadr W1)(car W2)))))


(defun lib:pline_clockwise ( lw  / LST MAXP MINP)
;;;Eugeni Elpanov
;;; Arguments
;;;         lw - ename or vla object of polyline
;;; Return
;;;     t - clockwise
;;;    nil - counter-clockwise
  
(if (= (type lw) 'ENAME)
    (setq lw (vlax-ename->vla-object lw)))  
(vla-GetBoundingBox lw 'MinP 'MaxP)
(setq
minp(vlax-safearray->list minp)
MaxP(vlax-safearray->list MaxP)
lst(mapcar(function(lambda(x)
(vlax-curve-getParamAtPoint lw
(vlax-curve-getClosestPointTo lw x))))
(list minp(list(car minp)(cadr MaxP))
MaxP(list(car MaxP)(cadr minp)))))
(if(or
(<=(car lst)(cadr lst)(caddr lst)(cadddr lst))
(<=(cadr lst)(caddr lst)(cadddr lst)(car lst))
(<=(caddr lst)(cadddr lst)(car lst)(cadr lst))
(<=(cadddr lst)(car lst)(cadr lst)(caddr lst))) t nil))

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

VVA, я же описал, если ближайшая точка в вершине, то нужно рассматривать вектор в конце предыдущего сегмента и в начале следующего после ближайшей вершины...
Те же проблемы будут и с выпуклым контуром, если ты укажешь точку снаружи, напротив вершины, т.е. ближайшей будет именно вершина контура.

[zamtmn]

Евгений этот алгорим будет работать только для выпухлых контуров, но он быстр как и все твои алгоритмы =)

Этот алгоритм совсем не быстрый, т.к. поиск ближайшей точки на контуре подразумевает поиск пересечений со всеми линиями образующими контур. Видимо его реализация быстрее за счет

находишь ближайшую точку на полилинии к твоей точке vlax-curve-*

[Дима_]

То Евгений - для тупых - объясни как работает твоя функция на невыпуклой (правой) фигуре?

[zamtmn]

понятие обхода вершин применимо к невыпуклым контурам. в первом случае направление обхода (касательной к контуру) в ближай точке будет в одну сторону, во втором в противоположную.

[Li6-D]

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

Код:

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

(defun conclusion ( / n ent PN P P0 a)
;|Функция запрашивает контур (2D-полининию из отрезков) и точку P.
  Если контур замкнут и точка не находится на его границе,
  а внутри контура, фунция возвращает T. Иначе - nil.
  Функция также печатает в текстовое окно соответствующее сообщение.
  Может некорректно работать, если полилиния содержит дуги|;
  (setq n 0) ;Начальное значение числа пересечений контуром оси PX.
  (princ "\nУкажите контур (2D-полилинию).")
  (and
    (setq ent (ssget "_:S" '((0 . "LWPOLYLINE"))))
    (setq PN (vl-remove-if '(lambda (x) (/= (car x) 10))
                (setq ent (entget (ssname ent 0)))
    )         )
    (or (= (cdr (assoc 70 ent)) 1) (prompt "\nКонтур не замкнут."))
    (setq P (getpoint "\nУкажите точку: "))
    (or
      (not (setq P0 (ssget P)))
      (not (ssmemb (cdar ent) P0))
      (prompt "\nТочка на границе контура.")
    )
    (setq PN (mapcar '(lambda (x) (angle P (cdr x))) PN)
          P0 (last PN))
    (foreach P PN
      (setq a (- P P0) P0 P)
      (or (<= (abs a) Pi) (setq n (1+ n)))
    )
    (or (/= (rem n 2) 0) (prompt "\nТочка снаружи контура."))
    (princ "\nТочка внутри контура.")
    t
) )

Долго не проверял

[Disney]

А у меня другая дилемма, как для команды-boundary найти точку внутри двух контуров красного и синего, чтоб получить зелёный?

[trir]

http://alglib.sources.ru/articles/convpol.php

[Disney]

Сообщение от trir http://alglib.sources.ru/articles/convpol.php

Прочитал, но не понял ни чего

[VVA]

Disney,
Как вариант
1. Найти центроид одного контура
2. Проверить, находится ли он внутри другого
3. Если нет - поменять контура местами
4. Если опять нет - ?

----- добавлено через ~17 мин. -----
Еще один алгоритм
1. 1-й контур преобразовать в region
2. 2-й контур преобразовать в регион
3. Найти пересечение областей
4. У полученного контура найти центроид
Что-то похожее было здесь
Там надо нарисовать замкнутую полилинию, задать размеры отверствия и полученное отверствие "посадить" на один из сегментов

[Disney]

Сообщение от VVA Disney, Как вариант ...Найти центроид ...

С центроидами не вариант, в моём примере он будет лежать вне контура примерно по центру картинки
Я думал о подобном, взять центр из GetBoundingBox.

[Кулик Алексей aka kpblc]

Не прокатит. Возьми фигуру типа "полумесяц"....

[Disney]

Обычно ситуация такая

и там вполне "центр" прокатит, но не исключен вариант с полумесяцами, или например такой

[VVA]

Сообщение от Disney С центроидами не вариант

Сообщение от Кулик Алексей aka kpblc Не прокатит. Возьми фигуру типа "полумесяц"....

Я в #32 добавил еще один вариант, преобразовать в области, найти пересечение контуров и ... мне кажется этого достаточно ?

[Disney]

Сообщение от VVA преобразовать в области, найти пересечение контуров и ... мне кажется этого достаточно ?

Да точно, просто я понять не мог, чем пересечение(IntersectWith) областей отличается от пересечения полилиний.
А оказалось ты о команде Пересечение , я привык её к 3D телам применять, а то что она с областями сработает и не знал. Вообще ни когда области не использовал и не понимал для чего они... я и сейчас не понимаю, но в моём случаи пригодиться.
Спасибо!

----- добавлено через ~18 ч. -----
Всё получилось.
Но теперь как мне достать координаты вершин получившейся области?

[fasadel]

Была у меня такая задачка.

[Disney]

Сообщение от Disney Но теперь как мне достать координаты вершин получившейся области?

Разбил область на линии, и уже из них выудил координаты.

----- добавлено через 20 сек. -----

Сообщение от fasadel Была у меня такая задачка.

Ну и как ты её решил?

[VVA]

Сообщение от Disney Разбил область на линии, и уже из них выудил координаты.

Все правильно. Ссылку на пример я приводил в #32

Сообщение от VVA Что-то похожее было здесь

Фрагмент кода для иллюстрации

Код:

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

  (progn ;_ пересекаются
    (vl-cmdf "_.Region" pline "")(setq pline (entlast))
    (vl-cmdf "_.Region" otv "")(setq otv (entlast))
    (vl-cmdf "_subtract" pline "" otv "")
    (setq otv (entlast) nab (ssadd))
    (foreach item (vlax-invoke (vlax-ename->vla-object otv) 'Explode)
      (ssadd (vlax-vla-object->ename item) nab))
    (if (and (getvar "PEDITACCEPT") (= (getvar "PEDITACCEPT") 1))
    (vl-cmdf "_pedit" "_Multiple" nab "" "_Join" 0 "")
    (vl-cmdf "_pedit" "_Multiple" nab "" "_Y" "_Join" 0 ""))
    (setq pline (entlast) nab nil)
    (entdel otv)
    )