구면 다각형 중심 값 계산


11

구의 다각형에 대한 중심을 계산하는 일반적인 방법을 원합니다.

지금까지 가장 좋은 온라인 참조는 다음과 같습니다.

Jeff Jenness의 그래픽 및 도형 도구 .

여기에 기술 된 방법은 다각형을 다수의 구형 삼각형으로 분해하고, 구형 삼각형 면적에 의해 가중 된 구형 삼각형 중심의 평균을 계산하는 것을 제안한다.

구형 다각형 중심을 정의하는 방법에는 여러 가지가 있지만 점과 폴리선에 대한 다음 정의와 유사한 것을 찾고 있습니다.

  • : 점을 나타내는 데카르트 벡터의 산술 평균.
  • 폴리 라인 : 각 선분의 중간 점을 나타내는 데카르트 벡터의 가중 평균으로, 각 선분의 (구형) 길이로 가중됩니다.

다각형 중심이 면적에 따라 가중 된 삼각형 분해의 가중 평균으로 정의되는 것이 합리적으로 계속되는 것 같습니다 .

내 질문은 사용 된 삼각형 분해에 관계없이 위 참조의 방법이 작동하는지 여부입니다. 특히, 다각형 외부에서도 임의의 점을 기준으로 삼각형으로 분해하여 일부 삼각형은 음의 가중치를 갖는 음의 영역을 갖습니다.

관련 : 어떻게 물체의 형상의 중심을 찾는 방법은?

답변:


9

고정 된 단일 점을 기준으로 모든 삼각 분할을 수행하더라도 일관되게 작동하지 않습니다 . 문제는 그것이 의미하는 바를 고려하지 않고 구형 및 유클리드 계산이 혼합되고 있다는 것입니다.

이를 명백하게하는 한 가지 방법은 반구의 거의 절반과 같은 다소 극단적 인 삼각형을 고려하는 것입니다. 예를 들어, (lon, lat) = (-179, 0)에서 시작하여 적도를 따라 (0, 0)으로 이동 한 다음 (0, 90)에서 북극까지, (- 179, 0). 이것은 서반구의 대부분의 북반부를 포함하는 90-179-90 삼각형입니다. 문제는 끝점 (그림에서 흰색 점으로 표시)이 실제로 평면에 있다는 것입니다. 하나는 기둥에 있고 다른 하나는 거의 반대편에 있습니다. 따라서 구형 (빨간색 점)으로 다시 투영 된 평균은 거의 극점에 도달하지만 가능한 한 합리적인 중심에서 멀어 집니다.

큰 구형 삼각형

또 다른 예로서, 중심, 북극에 대한 상반 구를 나타내는 다각형을 삼각 분할합니다. 우리는 항상 서반구를 두 개의 동일한 반쪽으로 나눌 것입니다. 각 반구는 90-90-90 삼각형입니다. 그러나 동반구는 n 개의 동일한 반음 으로 나뉩니다 . lune k ( k = 1, 2, ..., n ) 의 꼭짓점 에는 (lon, lat) 좌표가 있습니다

((k-1) * 180/n, 0),  (k * 180/n, 0),  (k * 180/n, 90).

k = 8의 Lunes

이 그림은 k = 8 설정을 보여줍니다. 빨간 점은 "그래픽 및 도형 도구"문서, pp 65-67에 따라 계산 된 개별 삼각형 ​​"중심"입니다.

계산을 수행하면 k = 2 인 경우 면적 가중 중심이 실제로 북극에 있지만 (대칭 고려 사항으로 표시됨) n이 증가하면 결과가 서반구로 빠르게 이동하고 한계는 경도 -90도를 따라 위도 89.556도에 도달합니다. 이것은 북극 자체에서 남쪽으로 약 50km 떨어져 있습니다.

물론 20,000km에 이르는 다각형에 대한 +/- 50km 오차는 작습니다. 이 경우 다른 삼각 분할로 인한 임의 변동 의 총량 은 0.5 %에 불과합니다. 분명히 음수 삼각형을 포함하여 상대 오차를 임의로 크게 만들 수 있습니다 (작은 삼각형에 비해 실제로 큰 삼각형을 더하고 빼기 만하면됩니다). 그럼에도 불구하고 구면 계산을 시도하는 사람은 프로젝션 오류를 피하려고 노력하고 있으므로 높은 정확도를 찾고 있습니다. 이 삼각 측량 방법은 권장 할 수 없습니다.


큰 n에 대해 오류가 누적 될 수 있음을 증명했지만 접근 방식에 결함이 있는지 확실하지 않습니다. 한계 값을 달성하기 위해 어떤 n 값을 사용 했습니까?
Jason Davies

또한 계산을 하고이 심층을 살펴 주셔서 대단히 감사합니다. 문제를 해결하기 전에 좀 더 설명하고 싶습니다. :)
Jason Davies

Jason, 나는 당신에게 직관을 제공하기 위해 예비 예를 추가했습니다. 한계 자체가 빠르게 다가옵니다. 수십 개의 루네는 몇 개의 유효 숫자를 얻을 것입니다. 그러나 새로운 가중은이 가중 삼각 분할이 작은 삼각형을 제외하고는 합리적으로 어떤 일을하는지에 대한 의문을 제기해야합니다. 구면 계산을 수행하는 유일한 이유는 분석 영역이 실제로 전역일 때 모든 투영에 많은 왜곡이 발생하기 때문입니다.
whuber

1
정말 고마워요 따라서 올바르게 이해하면 직교 벡터를 평균화해도 구형 삼각형 (특히 첫 번째 예와 같은 큰 삼각형)에 대해 합리적인 중심이 생성되지 않습니다. 예를 들어 큰 원 중앙값의 교차점을 찾는 더 나은 방법을 조사하겠습니다.
Jason Davies

BTW, 나는 위와 비슷한 구면 가중 중심이 작동하기를 희망합니다. 구의 원점에 꼭짓점을 추가하여 각 다각형에 3D 볼륨이 주어진다고 상상해보십시오. 그런 다음 원점에 연결된 보이지 않는 줄로 구를 중단하고 안정된 평형을 찾으십시오. 중심은 최하점입니다 (구면에 질량 중심 투영). 예를 들어 적도 주위를 돌고있는 스트립과 같은 모호한 경우를 제외하고는 합리적인 지점을 선택할 수 있습니다. 가치가 있다고 생각되면 새로운 질문으로 기꺼이 토론하십시오.
Jason Davies

3

다각형의 중심에 있어야하는 속성을 열거하는 것이 좋습니다. 내 기준은 다음과 같습니다.

(a) 정점이나 가장자리 대신 다각형 내부의 속성입니다. 따라서 추가 정점을 삽입하여 모서리를 두 개로 나누면 중심 위치가 바뀌지 않아야합니다. 중심 의 위치는 다각형이 삼각형으로 분할되는 방식에 따라 달라지기 때문에 Jenness의 중심에 대한 정의는 이 기준에서 실패합니다 .

(b) 다각형의 모양을 조금씩 교란 시키면 중심이 조금 움직여야합니다. 여기서 다각형의 전체 범위 (예 : 단일 반구)에 대한 제한을 적용해야합니다. 이 제한이 없으면 중심이 약간의 정점 이동으로 갑자기 지구 반대편으로 흔들리는 경우를 쉽게 구성 할 수 있습니다. 이 조건은 중심이 다각형 안에 있어야하는 방법을 제외합니다.

(c) 작은 다각형에 대한 중심의 평면 정의로 축소되어야합니다.

다음은 이러한 기준을 충족하는 두 가지 접근 방식입니다.

(1) 타원 다각형의 중심을 3 차원으로 계산하고 타원체 표면으로 다시 투영합니다 (타원체의 법선을 따라). 큰 장점 : 다각형을 간단한 모양으로 분할하여 중심을 계산할 수 있습니다.

(2) 중심은 다각형 내부의 모든 점까지 최소 RMS 측지 거리를 가진 점입니다. 버스 및 필 모어, "구면 스플라인 및 보간에 대한 구면 평균 및 응용 프로그램", 그래픽 20 , 95–126 (2001)의 ACM 트랜잭션을 참조하십시오. 큰 장점 : 결과 포인트는 표면이 R 3에 포함되는 방식에 의존하지 않습니다 .

불행히도, 이러한 정의 중 어느 것도 실제로 적용하기에 간단하지 않습니다. 그러나 첫 번째 방법은 단순히 구를 위해 수행 될 수 있습니다. 사용하기에 가장 좋은 "기본"영역은 다각형의 가장자리로 둘러싸인 사변형, 가장자리의 끝점을 통한 두 자오선 및 적도입니다. 전체 다각형의 결과는 가장자리에 대한 기여도를 합산합니다. (다각형이 극을 둘러싸면 추가 단계를 수행해야합니다.)

가장자리의 끝 점이 (φ 1 , λ 1 ) 및 (φ 2 , λ 2 )라고 가정하십시오. 가장자리와 끝점의 방위각을 α 1 과 α 2로 하자 . 구의 반지름이 1이라고 가정하면 사변형의 면적은

  A = α 21
      = 2 tan -1 [tan ½ (λ 21 ) sin ½ (φ 2 + φ 1 ) / cos ½ (φ 2 + φ 1 )]

베셀로 인해이 영역에 대한 공식은 삼각형 영역에 대해 일반적으로 사용되는 L' Huilier의 공식보다 수치 적으로 더 잘 작동합니다.

이 사변형에 대한 중심의 구성 요소는

  2 ⟨ X ⟩ = φ 2 죄 (λ 2 - λ 0 ) - φ 1 죄 (λ 1 - λ 0 )   2 ⟨ Y ⟩ = COS α 02 - σ 1 ) - (φ 2 COS (λ 2 - λ 0 ) - φ 1 COS (λ 1 - λ 0 ))   2 ⟨ Z ⟩ = (λ 2 - λ 1 ) - α 죄 02 - σ

1 )

여기서 σ 21 은 모서리의 길이이고 λ 0 및 α 0 은 적도를 가로 지르는 모서리의 경도와 방위각이며 xy 축은 적도 교차가 x =에 있도록 방향이 지정됩니다. 1, y = 0. ( z 는 물론 극을 통과하는 축입니다.)


Jenness 중심의 위치가 다각형을 삼각형으로 나누는 방법에 의존하는 이유를 설명 할 수 있습니까? @whuber의 예에서 Jenness의 중심 계산이 구형 삼각형에 대해 잘못되었다는 것을 알고 있지만 구형 삼각형 중앙값을 기반으로 한 중심이 대신 사용되면 어떻게됩니까? 여전히 실패합니까?
Jason Davies

Jenness는 구형 다각형을 평면 삼각형 세트로 효과적으로 대체하고 중심을 계산합니다. 분명히 (?) 결과는 분할에 따라 다릅니다. 구형 삼각형의 중심을 사용하여 설명한 계산을 수행하는 것이 좋습니다. JE Brock, 구면 삼각형의 관성 텐서, J. Applied Mechanics 42, 239 (1975)에서 중심에 대한 공식을 찾을 수 있습니다. dx.doi.org/10.1115/1.3423535
cffk

나는 브록의 논문을 다시 보았습니다. 구형 삼각형의 질량 중심에 대한 그의 공식은 삼각형의 가장자리에 대한 합을 포함합니다. 따라서 다각형을 삼각형으로 나눌 필요없이 다각형에 적용하도록 일반화 할 수 있습니다.
cffk

베셀로 인한 면적 계산에 대한 참조도 제공 하시겠습니까? 어디에서나 찾을 수없는 것 같습니다. 빠르고 정확한 구형 다각형 영역 루틴을 작성하는 데 관심이 있습니다. 감사!
Jason Davies

나는 그것을 발견하고 그것을 당신이 영어로 번역한다는 것을 깨달았습니다. 감사합니다. :)
Jason Davies
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.