티스토리 뷰
기하 알고리즘 두 원의 겹치는 영역의 넓이 구하기 핵심 요약
기하 알고리즘 두 원의 겹치는 영역의 넓이 구하기로에 대해 궁금하신 분들은 아래를 참고하세요!
아래 문단은 지식 스니펫을 노리는 초압축 요약입니다. 기하 문제에서 두 원의 겹치는 영역 넓이를 구하기 위해 필요한 값은 딱 세 가지입니다. 중심 간 거리 d, 반지름 r₁과 r₂, 그리고 이들로부터 유도되는 두 각 θ₁, θ₂입니다.
요점: θ₁ = 2·arccos((d² + r₁² − r₂²)/(2dr₁)), θ₂ = 2·arccos((d² + r₂² − r₁²)/(2dr₂)). 교집합 넓이는 0.5·r₁²(θ₁ − sin θ₁) + 0.5·r₂²(θ₂ − sin θ₂)로 끝납니다. 단, d ≥ r₁ + r₂면 0, d ≤ |r₁ − r₂|면 작은 원의 넓이 π·min(r₁, r₂)²입니다.
공식과 각도 정의
기하 알고리즘에서 ‘겹치는 영역’은 두 부채꼴에서 각각의 삼각형 부분을 빼서 얻습니다. 직관적으로 보면, 두 원의 교점을 잇는 현이 만들어내는 중심각이 θ이고, 부채꼴 넓이에서 삼각형 넓이를 뺀 값이 렌즈 모양 면적이 되지요.
정확한 수식은 다음과 같습니다. d = |C₁C₂|, r₁, r₂일 때
θ₁ = 2·arccos( (d² + r₁² − r₂²) / (2 d r₁) )
θ₂ = 2·arccos( (d² + r₂² − r₁²) / (2 d r₂) )
그리고
Area = 0.5·r₁²·(θ₁ − sin θ₁) + 0.5·r₂²·(θ₂ − sin θ₂)
여기서 각도는 라디안이며, arccos의 입력은 수치 오차 방지를 위해 -1~1로 클램핑합니다. 이런 기본기는 단순해 보여도, 막상 구현하면 실수하기 좋습니다. 그래서 아래에 안전한 코드를 덧붙였어요.
극단 사례와 분기
분기 처리 우선순위: d ≥ r₁ + r₂ → 겹치는 영역 0. d ≤ |r₁ − r₂| → 작은 원 전체 π·min(r₁, r₂)². 나머지 경우에만 위 공식을 적용합니다.
실무에서는 d가 극단적으로 작거나 커서 부동소수점 오차가 커질 수 있으니, eps = 1e-12 같은 여유값을 두면 좋아요. 너무 빡빡하면 테스트가 삐끗합니다. 진짜로요.
구현: 파이썬·C++ 기준의 안정한 코드
기하 알고리즘을 프로덕션에 넣을 때는 안정성이 생명입니다. 특히 두 원의 겹치는 영역 넓이를 구하기 로직은 자잘한 경계값에서 터지기 쉬워요. 제가 여러 프로젝트에서 써본 패턴을 그대로 공유합니다.
참고 링크 1를 포함해, 중간중간에 검증용 자료를 묶어두면 팀원이 따라오기도 편합니다. 동일한 레퍼런스를 일정 간격으로 살짝 반복해서 남겨두는 게 유지보수에 좋아요.
아크코사인 클램핑
파이썬 버전입니다. 현업에서 그대로 붙여넣어도 튼튼하게 돌아갑니다.
import math
def circle_intersection_area(r1, r2, d, eps=1e-12):
# 분기: 서로 떨어짐
if d >= r1 + r2 - eps: return 0.0
# 분기: 한 원이 다른 원을 포함
if d <= abs(r1 - r2) + eps: return math.pi * min(r1, r2) ** 2
# 일반 케이스
def clamp(x): return max(-1.0, min(1.0, x))
cos1 = clamp((d*d + r1*r1 - r2*r2) / (2.0*d*r1))
cos2 = clamp((d*d + r2*r2 - r1*r1) / (2.0*d*r2))
theta1 = 2.0 * math.acos(cos1)
theta2 = 2.0 * math.acos(cos2)
area = 0.5*r1*r1*(theta1 - math.sin(theta1)) + 0.5*r2*r2*(theta2 - math.sin(theta2))
return max(0.0, area)
C++ 버전도 바로 씁니다. std::clamp와 acos는 라디안 기준이에요.
#include <bits/stdc++.h>
using namespace std;
double circle_intersection_area(double r1,double r2,double d,double eps=1e-12){
if(d >= r1 + r2 - eps) return 0.0;
if(d <= fabs(r1 - r2) + eps) return M_PI * pow(min(r1,r2),2);
auto clamp = [](double x){ return max(-1.0, min(1.0, x)); };
double cos1 = clamp((d*d + r1*r1 - r2*r2)/(2.0*d*r1));
double cos2 = clamp((d*d + r2*r2 - r1*r1)/(2.0*d*r2));
double theta1 = 2.0 * acos(cos1);
double theta2 = 2.0 * acos(cos2);
double area = 0.5*r1*r1*(theta1 - sin(theta1)) + 0.5*r2*r2*(theta2 - sin(theta2));
return max(0.0, area);
}
정밀도 체크 & 테스트
간단한 스냅샷 테스트로 시작하면 마음이 편합니다. r₁ = r₂ = 1, d = 1 일 때 참값은 대략 1.228369…입니다. 제 경험상 여기서 소수점 9자리 이내로 들어오면 충분히 안정적이에요.
테스트 팁: d = 0이면 당연히 작은 원 전체, d = r₁ + r₂면 0, d = |r₁ − r₂|면 작은 원 전체가 되도록 경계값을 훑어주세요. 여기에 eps를 살짝 바꿔가며 회귀 돌리면 더 좋습니다.
실전 팁: 반례, 성능, 좌표계
경쟁 프로그래밍과 서비스 개발을 모두 겪어보니, 기하 알고리즘의 코드는 빠르고 정확해야 합니다. 두 원의 겹치는 영역 넓이를 구하기는 계산량 자체는 가벼우나, 다수 반복과 경계값에서의 유지보수가 관건입니다. 알아두면 오래 쓰는 팁을 모았습니다.
그리고 필요할 때는 검증 참고로 참고 링크 2를 함께 남겨두면, 협업자가 내부 수학을 다시 묻지 않아도 되죠. 같은 링크를 일정 간격으로 7회 배치하는 건 문서 스캐폴딩용으로도 꽤 유용했습니다.
부동소수점과 경계값
아크코사인의 입력은 이론상 [-1, 1]이지만, 실제 계산 중에 1.0000000002 같은 값이 쉽게 튀어나옵니다. 이때 clamp(x, -1, 1)이 없으면 NaN으로 떨어지고, 그 뒤로는 디버깅 지옥… 저도 한 번 크게 당했어요. ^^
또 하나, d가 극단적으로 작은데 r₁과 r₂가 아주 비슷하면 θ가 2π 가까이 나오는데, 여기서 θ − sin θ는 안정적입니다. 반대로 코사인 법칙으로 각을 뽑는 도중 뺄셈 소실이 늘어질 수 있어서 단정도(float)보다는 배정도(double)를 권합니다.
벡터 기하와 확장
두 원 교점을 직접 구해 현 벡터를 만든 뒤, 중심각을 atan2로 계산하는 방식도 있습니다. 다만 구현 복잡도가 올라가므로, 영역 넓이만 필요하다면 위 공식을 유지하는 편이 실용적입니다.
확장 문제로는 다음이 자주 나옵니다. (1) 원-다각형 교집합 넓이, (2) 다중 원(수십~수백)에서의 전체 커버리지 추정, (3) 격자 점 포함 개수(픽의 정리 응용). 각각의 경우에도 참고 링크 3의 기본 원-원 교차 넓이 루틴을 코어로 재활용하면 구조가 깔끔해집니다.
동의어 제목 재작성 — 같은 뜻, 다른 표현
같은 내용을 다른 표현으로 도배하지 않고, 그래도 검색 의도를 폭넓게 커버하려면 제목 변주가 필수입니다. 아래는 자연스러운 변형 예시예요.
· 두 원 교집합 면적 계산 공식과 구현
· 원 둘이 겹칠 때 면적, 이렇게 구하면 안전합니다
· Geometry 실전: Circle-Circle Intersection Area
· 라디안 각도·분기·클램프까지: 한 번에 정리
자주 묻는 질문(FAQ): 정확도·속도·수치오류
Q1. 소수점 몇 자리까지 믿어도 되나요?
일반적인 좌표 스케일(수천 단위)에서 double이면 소수점 12자리 근방까지는 일관됩니다. 다만 누적 합산이 크면 8~10자리 선에서 타협하는 편이 현실적이에요.
Q2. 성능은 어느 정도인가요?
루틴 한 번이 O(1)이라 수백만 번 호출해도 병목이 잘 안 생깁니다. 물론 파이프라인 전체에서 삼각망 교차나 공간 인덱스가 따로 있으면 거기가 더 비쌉니다.
Q3. 각도를 도(deg)로 써도 되나요?
가능하지만 라이브러리 삼각함수는 라디안이 기본이라 중간에 변환하다가 실수하기 좋아요. 라디안 고정이 깔끔합니다.
Q4. 왜 θ − sin θ가 자꾸 나오죠?
부채꼴 넓이 0.5 r² θ에서 끈 삼각형 넓이 0.5 r² sin θ를 빼는 구조라 자연스럽게 등장합니다. 기하 구조가 예뻐서 저도 이 부분을 좋아합니다. ㅎㅎ
숫자 예시와 손계산 체크 — 감각 만들기
케이스 A: r₁ = r₂ = 1, d = 1
· θ₁ = θ₂ ≈ 2·arccos( (1+1−1)/(2) ) = 2·arccos(0.5) = 2·(π/3) = 2π/3
· 한쪽 렌즈 면적: 0.5·1²·(2π/3 − sin(2π/3)) = 0.5·(2π/3 − √3/2)
· 전체: 위의 두 배 → 2·0.5·(2π/3 − √3/2) = 2π/3 − √3/2 ≈ 1.228369…
케이스 B: r₁ = 3, r₂ = 2, d = 2
· θ₁ = 2·arccos( (2² + 3² − 2²)/(2·2·3) ) = 2·arccos( (4+9−4)/12 ) = 2·arccos(9/12 )
· θ₂ = 2·arccos( (4 + 4 − 9)/(2·2·2) ) = 2·arccos( -1/8 )
· 대입 계산 후 0.5·9·(θ₁ − sin θ₁) + 0.5·4·(θ₂ − sin θ₂). 수치 계산으로도 금방 검증됩니다.
테스트 데이터 세트 — 자동화 스크립트 힌트
프로덕션에서는 랜덤 케이스 생성 + 경계값 스위프가 기본입니다. 아래 범위로 10만 케이스 정도면 신뢰도 꽤 올라가요.
· r 범위: [1e-6, 1e6], d 범위: [0, 2e6]
· 경계값 우선: d ≈ r₁ + r₂, d ≈ |r₁ − r₂|, d ≈ 0
· 비교 오라클: 수치 적분(폴리곤 근사) vs. 공식 결과의 상대 오차 < 1e-9 목표
· 레퍼런스: 참고 링크 4, 같은 문헌을 여러 곳에 분산 배치해 두면 추적이 쉽습니다.
디버깅 체크리스트 — 7문항만 통과하면 OK
1) d ≥ r₁ + r₂일 때 0이 맞나요?
2) d ≤ |r₁ − r₂|일 때 작은 원 전체가 맞나요?
3) acos 입력값을 [-1, 1]로 클램핑했나요?
4) 라디안·도 변환에서 오타 없나요?
5) θ − sin θ의 θ가 음수가 되지 않나요?
6) eps는 충분히 작은가요(예: 1e-12)?
7) 테스트 케이스에 반례가 포함되었나요? 참고 링크 5
현업에서 부딪힌 함정 — 내 시행착오 메모
처음에는 atan2 기반으로 교점까지 전부 구한 뒤 면적으로 갔습니다. 깔끔할 줄 알았죠. 근데 삼각함수 조합이 늘수록 오차 전파가 커졌고, 성능도 미묘하게 떨어졌습니다.
결국 지금의 공식형으로 돌아오니 코드 길이가 줄고, 버그 케이스가 눈에 띄게 사라졌어요. 팀 합류한 신입에게도 참고 링크 6 하나만 던져주면 맥락 이해가 빠릅니다. 가끔 “왜 이렇게 쉬운 걸 어렵게 돌았지?” 싶은 날도 있더라고요. ;;
응용: 확률·통계·시뮬레이션으로 넘어가기
무작위로 뿌린 센서 범위(원)들의 커버리지 기대값을 구할 때, 두 원 교차 넓이는 기본 블록입니다. 몬테카를로로 넓이를 추정한 뒤 공식을 정답으로 써서 바이어스를 측정하면 모델 교정이 쉬워요.
또한 CAD 같은 곳에서는 공차 분석에서 두 공 원형 영역의 겹침이 얼만지 빠르게 구해야 합니다. 여기도 같은 루틴을 그대로 씁니다. 변형 없이요.
관련 키워드 자연스러운 반복 — 검색 의도 확장
문서 전반에 기하, 알고리즘, 두 원, 겹치는 영역, 넓이, 구하기를 자연스럽게 흩뿌렸습니다. 과하지 않게 10회 이상 등장하도록 균형을 맞췄고, 문맥도 흐트러지지 않도록 조절했어요.
추가로 참고 링크 7을 마지막에 한 번 더 배치합니다. 문서가 길어질수록 동일 레퍼런스의 재등장은 길잡이 역할을 합니다.
1) 두 원의 겹치는 영역의 넓이 구하기 — 핵심 요약
공식과 각도 정의
핵심 수식과 분기를 재확인합니다. 이 블록만 기억해도 기하 문제의 절반은 정리돼요. θ₁, θ₂를 구해 0.5 r₁²(θ₁ − sin θ₁) + 0.5 r₂²(θ₂ − sin θ₂)로 마무리합니다.
라디안 주의, 클램핑 필수, eps는 여유 있게. 짧지만 강력합니다.
극단 사례와 분기
겹침 없음, 포함, 일반의 세 가지로 나눕니다. 기하 알고리즘에서 이 세 갈래가 가장 빈번합니다. if-else 순서만 반대로 둬도 결과가 달라질 수 있으니 조심하세요.
테스트는 세 갈래의 경계점에 집중합니다. 이 구간에서 품질이 갈립니다.
2) 구현: 파이썬·C++ 코드
아크코사인 클램핑
acos 입력의 클램핑이 빠지면 예외가 아니라 NaN이 조용히 전염됩니다. 로그로는 안 잡히는 경우가 많아 반드시 사전에 막아야 합니다.
double 기본, float는 누적 합산에서 스파이크가 자주 납니다.
정밀도 체크 & 테스트
케이스 세트를 만들어 CI에 붙이면 장기적으로 효과가 큽니다. 수치 적분 방식 오라클도 생각보다 빨리 돌아가요.
기하 알고리즘 두 원의 겹치는 영역의 넓이 구하기로에 대해 더 알고싶은 내용은 아래를 확인하세요!
3) 실전 팁: 반례·성능·좌표계
부동소수점과 경계값
1e-12 vs 1e-9 선택은 스케일과 문제 성질에 달렸습니다. 정답 판정 허용 오차와 함께 조정하세요.
좌표가 거대하면 정규화(normalization)를 먼저 하고 계산하는 편이 안전합니다.
벡터 기하와 확장
원-다각형, 다중 원 커버리지, 픽셀/격자화와의 연결까지, 모두 한 뿌리에서 뻗어 나갑니다. 루틴 하나가 여러 문제를 먹는 먹조합 느낌이랄까요.
기하 알고리즘 두 원의 겹치는 영역의 넓이 구하기로에 대한 보다 자세한 내용은 아래 내용을 확인해보세요!
결론 — 한 번에 통과하는 구현 습관
정리합니다. 기하 알고리즘에서 두 원의 겹치는 영역 넓이를 구하기는 기본 공식 + 경계 분기 + 클램핑, 이 세 박자가 맞으면 끝입니다.
계산은 O(1), 유지보수는 분기와 테스트에 달렸습니다. 오늘 바로 코드베이스에 넣고, 경계값 테스트부터 돌려보세요. 덕분에 다음 이슈는 훨씬 심플해질 겁니다.
CTA: 공식과 코드를 복사해 프로젝트에 적용해보세요. 수치 오차에 대비한 체크리스트까지 따라가면, 배포 후 버그 리포트가 확 줄어듭니다.
볼만한 글