다각형의 무작위 샘플링


9

다각형에서 균일하게 임의의 점을 샘플링하고 싶습니다 ...

많은 수의 표본을 추출하는 경우 면적이 같은 경우 두 지역에 속할 가능성이 높습니다.

좌표로 [0,1]에 두 개의 난수를 사용할 것이기 때문에 정사각형이면 매우 간단합니다.

내가 가진 모양은 일반 다각형이지만 모든 다각형에서 작동하고 싶습니다.

/programming/3058150/how-to-find-a-random-point-in-a-quadrangle

답변:


9
  1. 다각형 삼각 분할
  2. 점이 어느 삼각형에 놓여 야하는지 결정 (삼각형 영역에 가중치 부여)
  3. 이 게시물에 설명 대로 삼각형의 점을 샘플링하십시오.

이 질문이 이전 질문과 중복되지 않습니까?
Raphael

@Raphael : 관련이 있지만 더 일반적입니다.
A.Schulz

4

쉬운 방법 중 하나는 다각형의 경계 상자를 찾고 거부 샘플링을 사용하는 것입니다. 경계 상자에서 샘플을 가져 와서 다각형 내에 들어간 경우 수락합니다 . 적어도 확률로 발생합니다 (생각합니다).1/2

다른 가능성은 다각형을 삼각 측량하는 것입니다. 먼저 비례 적으로 삼각형을 샘플링 한 다음 삼각형에서 임의의 점을 샘플링하십시오. 후자는 간단합니다. 변환을 세밀하게하기 위해 모든 삼각형은 입니다. 분포에서 점을 균일하게 샘플링하려면 먼저 밀도 에 따라 을 샘플링합니다 (즉, 균일 한 샘플링하고 계산합니다) ) 다음 균일하게 샘플링합니다 (즉, 균일 한 샘플링하고 계산합니다 ). 더 간단한 방법은 을 샘플링 하는 것입니다.{(x,y):x,y0,x+y1}x[0,1]2(1x)r[0,1]x=11ry[0,1x]s[0,1]y=(1x)sx,y[0,1]x+y>1 대체 가진 .(x,y)(1x,1y)


거부 샘플링은 2 차원에서 최대 1/2의 확률로 거부하지만, 높은 차원에서는 거부 확률이 훨씬 나빠질 수 있습니다.
DW

거부 샘플링의 거부율은 1/2보다 클 수 있습니다. 약간 돌출 된 나선을 생각해보십시오.
A.Schulz

다각형이 볼록한 것이 보장되면 어떻게합니까?
Yuval Filmus

경계 상자가 축 정렬되어 있으면 볼록성이 도움이되지 않습니다. 앞의 질문에 대한 답변에서 알 수 있듯이 매우 큰 x에 대해 (0, 1), (1, 0) 및 (x, x)에 정점이있는 삼각형을 고려하십시오. 이는 경계 상자의 사라지는 작은 비율을 차지합니다 x는 무한대로갑니다. 가장 작은 경계 상자에 대해 이야기하고 있다면 아마도 볼록한 모양이 차지하는 볼륨에서 경계를 도출 할 수 는 있지만 상자를 찾아야합니다 ...
Steven Stadnicki

4

이것은 약간 미친 짓이지만 다각형이 매우 이상하더라도 잘 작동합니다.

Reimann 매핑 정리 를 사용하여 단위 디스크에서 다각형까지의 등각 맵을 찾아 의 하위 세트로 봅니다. 예를 들어 다음의 참조를 참조하십시오.C

http://siam.org/pdf/news/1297.pdf

그런 다음 Metropolis-Hastings MCMC 샘플링 에서 제안 밀도로 디스크에 균일 한 밀도를 적용하십시오 .


컨 포멀 맵은 반드시 지역을 보존하는 것은 아닙니다. 그것들은 각도 보존이지만 , 다각형을 균일하게 샘플링 하지 않도록 거의 보장 됩니다.
Steven Stadnicki

따라서 실제 샘플러가 아닌 MCMC의 제안으로 사용할 필요가 있습니다. Poincare 불평등을 사용하면 등고선 맵의 변형이 균일 함을 경계로 표시 할 수 있습니다.
Nick Alger

나는 아직도 그것을 놓치고있다; 위키 백과 뾰족한-에 페이지의 논의는 '재판 분포는'아직 할 필요가 있다고 직접 원하는 분포에 비례; 즉, 일부 상수 및 경우 가 아니라 입니다. 매핑 된 밀도의 로컬 분산은 여전히 ​​샘플링에서 로컬 분산으로 이어집니다. aP(x)<f(x)<bP(x)abf(x)=cP(x)x
Steven Stadnicki

Metropolis Hastings MCMC의 요점은 제안이 실제 배포가 아니라는 것입니다. MCMC 체인의 수렴 속도는 제안이 실제 분포에 얼마나 근접한 지에 달려 있습니다. 가장 일반적인 제안은 ... 당신이 샘플하려고하는 분포에 관계없이 현재 시점에서 가우스를두고있다
닉 알제
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.