[PHP] class of lunar, 달의 위치, 음양력변환, 별자리, 일월식 계산 [author] - Chilbong Kim, - http://linuxchannel.net/ [changes] - see http://ftp.linuxchannel.net/devel/php_lunar/ [conjunction error by approximative method] - avg: -6538 seconds(+598 seconds, -11787 seconds) [conjunction error by approximative method day patch] - avg: -368 seconds(+592 seconds, -1081 seconds) - but, wrong day not exists :) [valid date] - solar : 1902-01-10 - 2038-01-18 - lunar : 1901-12-01 - 2037-12-14 [download & online source view] - http://ftp.linuxchannel.net/devel/php_calendar/ - http://ftp.linuxchannel.net/devel/php_lunar/ [demo] - http://linuxchannel.net/gaggle/lunar.php [docs] - http://linuxchannel.net/docs/lunar.txt [study] - http://www.kao.re.kr/html/study/qna/index.html - http://sunearth.gsfc.nasa.gov/eclipse/SEsaros/SEsaros.html - Synodic Month (new moon to new moon) 29.53059 days = 29d 12h 44m (A) - Draconic Month (node to node) 27.21222 days = 27d 05h 06m (B) - Anomalistic Month (perigee to perigee) 27.55455 days = 27d 13h 19m - (223 * A) = 6585.36 days, (242 * B) = 6585.32 days - 1 Saros = 6585.32 days = 18 years 11 days 8 hours [eclipses] - Solar Eclipses: T(Total), A(Annular), H(Hybrid(Annular/Total)), P(Partial) - Lunar Eclipses: t(Total), p(Umbral(Partial)), n(Penumbral) [eclipses(e) by this program] - S A 104 0.010582273616323 1.06751983653440(1.067743) - S T 91 0.014674627401046 1.04166992331860 - S H 6 0.210833791736160 0.96121766199154(0.964180) - S P 106 0.829402137331320 1.59530999819400 - L t 114 0.010630181434784 0.52538271749109(0.527908) - L p 85 0.294681791935650 1.14183526374750 - L n 111 0.832900661071600 1.57436507046800(1.578244) [reference] - http://user.chollian.net/~kimdbin/re/moonpos.html - http://zimmer.csufresno.edu/~fringwal/skycalc.c - http://williams.best.vwh.net/sunrise_sunset_example.htm - http://aa.usno.navy.mil/faq/docs/SunApprox.html - ftp://ssd.jpl.nasa.gov/pub/eph/export/C-versions/hoffman/ - http://www.linuxchannel.net/docs/solar-24terms.txt - http://www.stjarnhimlen.se/english.html // good - http://www.stargazing.net/kepler/ // good - http://sunearth.gsfc.nasa.gov/eclipse/eclipse.html // good - http://sunearth.gsfc.nasa.gov/eclipse/phase/phasecat.html - http://sunearth.gsfc.nasa.gov/eclipse/phase/phase2001gmt.html - http://sunearth.gsfc.nasa.gov/eclipse/LEvis/LEaltitude.html - http://sunearth.gsfc.nasa.gov/eclipse/resource.html - http://www.mreclipse.com/Special/SEprimer.html // Solar Eclipses for Beginners - http://www.mreclipse.com/Special/LEprimer.html // Lunar Eclipses for Beginners - http://www.mreclipse.com/MrEclipse.html // for Beginners - http://astronote.org/ - http://myhome.naver.com/dudwn1109/appearance/solor_eclipse.htm - http://www.kao.re.kr/html/study/ ------------------------------------------- *주1) 이 문서에 대한 최신 내용은 아래 URL에서 확인할 수 있습니다. http://linuxchannel.net/docs/lunar.txt ------------------------------------------- 0. 배경 및 개요 1. 기본지식(알고 있거나 개념정도) 2. 어떻게 계산했나(계산방법)? 2-0. 음력과 윤달 2-1. 양력에서 음력으로 2-2. 음력에서 양력으로 2-3. 별자리(12궁도) 2-4. 일식 2-5. 월식 3. PHP 클래스 사용법 3-0. 함수(메소드) 요약 설명 3-1. 특정시각에 대한 달(moon)의 위치를 구하려면? 3-2. 양력을 음력으로 변환하려면? 3-3. 음력을 양력으로 변환(어르신 생신에 자주 사용되죠)하려면? 3-4. 특정 날짜에 대한 일식여부 확인해 보려면? 3-5. 특정 날짜에 대한 월식여부 확인해 보려면? 3-6. 특정 날짜에 대한 12궁도 별자리지를 알아보려면? 4. 후기 -------------------------------------------- 0. 배경 및 개요 생략 (굳지 찾자면 달의 위치를 구하고 싶어서.. 어렸을 때에는 별헤는 소년이었답니다. :) 1. 기본지식(알고 있거나 개념정도) 율리우스적일(Julian day), 황도, 백도, 황도좌표, 24절기, 절기와중기, 황경, 황위, 천문단위(AU), 적경, 적위, 윤달, 무중치윤법, 12궁도, 일식, 월식, 합삭, 망, 월령, 개기월식, 부분월식, 반영식, 금환식, 개기일식, 부분월식, Hybird, 사로스(Saros)주기 그 외 다수의 천문용어 2. 어떻게 계산했나(계산방법)? 2-0. 음력과 윤달 음력 한달(Synodic Month)은 어느 합삭시간부터 다음 합삭시간까지의 평균 간격 29.53059 일을 사용합니다. 이것을 `삭망월'이라고 합니다. 29일 일때에는 `소월', 30일 일때는 `대월'이라고도 하죠. 물론 천문학에서는 이 시간 외에 Draconic Month(27.21222 일), Anomalistic Month(27.55455 일)등이 있습니다. 참고로) - 태양력(양력) 현재 우리가 사용하고 있는 `양력'은 그레고리역으로 태양력입니다. 그냥 줄여서 `양력'이라고 하죠. 그외 음력의 종류(?)가 몇가지 있는데, - 태음력(순태음력) 삭망월 평균 29.5306 일 기준으로 함(현재 사용하고 있지 않음) - 태음태양력(음력) 현재 우리가 사용하고 있는 음력으로, 1삭망월(29.5306일)과 1회귀년(태양년), 24절기, 그리고 무중치윤법을 적용한 윤달이 있는 달력 체계로, 줄여서 그냥 `음력'이라고 합니다. 양력에서 양(陽)은 `태양(sun)'을 의미하고, 음력에서 음(陰)은 `달(moon)'을 의미 합니다. 모두 천체의 운동(운행)에 기반한 역(曆)입니다. (천문분야에서 사용하는 `천체력'이라는 달력도 있지요.) [합삭 - 초하루] 음력 초하루는 태양과 달이 합삭(conjunction)되는 날을 음력 초하루로 잡습니다. `합삭'은 천체의 운행에서 `태양-달-지구'의 순으로, 일직선 일때를 말하며, 합삭시각은 태양과 달의 황경값이 같은 때의 시각을 말합니다. 따라서 태양과 달의 황경을 구하면 초하루는 쉽게(?) 계산할 수 있죠. 문제는 얼마나 정확하게 태양과 달의 (지심시)황경(황도에서 겉보기 천체의 위치)을 구하느냐가 과제입니다. [망 - 보름] 반면 망(보통 `보름'이라고 함)은 `태양-지구-달'의 순으로 일직선 일때를 말하며, 태양과 달의 황경 차이가 180도 일 때를 말합니다. 음력으로 15일이나 16일에 해당 됩니다. 역시 황경값이 관건이네요. [윤달] 윤달은 어떻게 계산되나? 개념은 간단한데 수식이나 논리적으로 풀려면 상당히 어렵습니다(프로그램시 상당히 고전했던 부분이라서..) 윤달을 설정하는 이유는, (태)음력으로만 계속사용하면 계절의 차이가 많기 때문에 이 차이(1년에 약 11일 차이로 태음력이 작음)를 없애고자 윤달을 둡니다. 음력은 달의 이름을 짓기 위해서 최소한 한번의 중기(24절기에서 중기)가 들어가야 정상인데, 태음력이 짧기 때문에 어떤 달은 중기가 아예 없습니다. 이것을 `무중달' 이라고 하고 윤달의 기본이 됩니다. 물론 중기가 두번 들어가는 달도 있는데 이것은 윤달이 아닙니다. 대략 3년에 1번 윤달이 들어가게 됩니다. 그러나 무중달이라고 해서 무조건 모두 윤달로 설정하지 않습니다. 19년 사이에 7번의 윤달을 설정하고, 한해에 두개의 윤달이 있으면 먼저 온 녀석(?)을 윤달로 잡습니다. 윤달의 이름은 전달의 이름을 그대로 사용합니다. 이것을 `무중치윤법(無中 置閏法)'이라고 하고, 현재 우리나라 달력 자료를 제공하는 `천문연구원'에서는 이방법을 사용합니다. 음... 윤달을 설정하는데 24절기(태양의 황경으로 구분)를 사용한다는 점에서 무지 복잡하네요.. 음력에도 양력성분이 들어가기 때문에 `태음태양력'이라고 하는 이유가 여기에 있죠. (태음태양력의 준말이 `음력'임) 2-1. 양력에서 음력으로 양력에서 음력으로 계산하는 방법은 생각보다 어렵지(?) 않습니다. 즉 태양과 달의 시황경이 같은 시각을 초하루로 잡고, 특정 날에 대한 시간을 일 수로 계산하면 됩니다(말은 싶지만 구현하려면 머리털 다 뽑힙니다..크크). 예제) 양력 2005년 3월 1일을 음력으로 계산하려면? 1) 2005.03.01 의 태양의 시황경(Ls)과 달의 시황경(Ll)을 구합니다. 2) Ls 와 Ll 이 같도록(합삭) 시간을 뒤로 빼가면서 시각을 계산합니다. 3) 2)의 합삭시각이 구해지면(적분 비슷한 공식을 구현해야함) 이것이 초하루가 됩니다. 4) 초하루 시각부터 주어진 날에 대한 일 수를 계산합니다(일(day) 계산 완료). 5) 다음 달의 합삭시각(초하루)을 1),2),3)의 과정으로 구합니다. 6) 그해 몇번째 달인지 결정하기 위해서 합삭시각에 대한 각각의 절기, 즉 해당 월(양력)에 대한 절기와 중기, 그리고 다음 월(양력)에 대한 절기와 중기를 구합니다. 7) 달을 결정하기 위해서 이번달 초하루와 다음달의 초하루, 그리고 2개의 중기를 놓고 윤달을 판별합니다. 대개 이번달 초하루가 해당 달(양력)의 중기 이전에 있으면 (해당달-1)이 달의 이름이 됩니다. 또한 중기가 하나도 안 들어 있으면 윤달이 됩니다(실제 음력에 대해서 검산이 필요함). 윤달이면 해당 달에서 1을 빼주면 되죠. 8) 년도를 구합니다(12월이나 1월일 경우에는 잘못된 년도가 나오지 않도록 판별) 본 class.lunar.php 파일에서 tolunar() 함수를 분석해 보면 위와 같습니다... :) (분석하지 마세요.. 크크) 각 항목에 대한 설명은 다음 절의 함수 사용법에서 아주 간단하게 요약해 두었습니다. :) 태양과 달의 시황경 구하는 공식은 모두 근사식입니다. 그러나 오차가 그리 크지 않기 때문에 본 프로그램에서 채택했습니다(정확한 값을 구하려면 PHP 에서는 무리가 있음). 2-2. 음력에서 양력으로 음력에서 양력계산은 좀더 복잡합니다(음력생일을 양력으로 계산할 때 많이 사용하죠). 양력에 대한 특정 시각이 주어지지 않았기 때문에 곧바로 태양과 달의 시황경을 구하지 못합니다. 때문에 최적의 특정 시각을 정해서 그 시각을 기준으로 합삭시간을 찾아가는 알고리즘을 구현해야 합니다(역시 말을 싶지만 적분 비슷한 알고리즘이..) 이것 역시 합삭시간을 구하는 것이 우선입니다. 예제) 음력 2005년 1월 1일(설날)을 양력으로 계산하려면? 1) 양력으로 임시시각(년-월-일 시분초)을 잡습니다(주어진 음력을 기준). 2) 1)의 시각에 대한 합삭시각을 구하고 해당 정보(윤달) 등등을 판별합니다. 3) 판별후 정보가 맞지 않으면, 1),2)의 과정을 되풀이 합니다. 또한 주어진 정보가 정확한지도 역시 판별해야 합니다(아주복잡). 4) 3)에서 합삭시각(unix timestamp)이 구해지면 그 시각에 대한 양력을 구합니다 (이건 쉽죠). 5) 4)의 양력 일수에 주어진 음력 일수를 더해서 최종적으로 양력 날짜를 구합니다. 이 계산 역시, `양력에서 음력' 계산이 틀리면 덩달아(?) 음력에서 양력 계산도 함께 틀려지기 때문에 `양력에서 음력' 계산이 아주 중요합니다. 2-3. 별자리(12궁도) 원래 class.lunar.php 에서는 예정에 없는 함수인데, 태양의 황경, 음양력 변환을 계산해 놓으니 별자리 계산은 간단하더군요. 그래서 그냥 재미로 한번 넣어 두었습니다. 별자리는 점성술(?)에서 사용되곤 하지요. 모두 태어난 날을 기준으로 자신의 별자리를 가지게 되죠(제건 `게자리'입니다. 히히) 우선, 별자리를 계산하는 시간체계는 양력입니다. 음력으로 계산하지 않습니다. 이유는 자신의 별자리는 점성술(Horoscope)에서 언급되어 있는데, 이때 계산의 기준이 되는 것이 황도(먹는것 아닙니다. T.T)입니다. 황도는 태양의 겉보기 운동에 따른 가상의 원이고, 별자리는 이 황도를 따라 12 개의 일정 간격으로 나누는데 이것을 `12궁도'라고 합니다. 궁 하나하나에 별자리 이름을 붙이죠. 궁의 의미는 `궁합을 보다'에서 그 `宮(집)'입니다. 12궁도를 나누는 기준은 역시 황도에서 태양이 위치한 지점, 즉 태양의 시황경 값이 그 기준(30도 간격)입니다(24절기도 태양의 황경값으로 계산하는데 15도가 기준이죠). 실제 천문계산이 아니고서는 이 `황경값' 구하기가 상당히 어렵습니다. 대부분 싸이트는 년도별로 DB를 구축해 놓고 그 값을 검색해서 이용할 뿐이죠. 별자리 정보를 찾다보면 특이한 점을 발견할 수 있는데, 어떤 싸이트는 `천칭자리'를 9/23 - 10/21 까지라고 하고, 또 다른 싸이트는 9/23 - 10/23 까지라고 못박고 있습니다. 둘다 맞을 수 있고 틀릴 수 도 있는데, 둘다 틀릴 확률이 높습니다(이 무슨 괴변 ..) 이유는 해마다 같은 황경값에 대한 날짜는 조금씩 상이하기 때문이죠. 용어) - Horoscope 12 궁도 - House 궁(宮) - constellation 별자리 - asterism 성군(星群), 성좌 계산) 별자리(양력으로 계산) i = floor(태양의황경/30) 황경 i 구간(대략) 이름 영문 0 0 0321 - 0420 산양 Aries 30 1 0421 - 0520 황소 Taurus 60 2 0521 - 0621 쌍둥이 Gemini 90 3 0622 - 0722 게 Cancer 120 4 0723 - 0822 사자 Leo 150 5 0823 - 0922 처녀 Virgo 180 6 0923 - 1021 천칭 Libra 210 7 1022 - 1121 전갈 Scorpius 240 8 1122 - 1221 궁수 Sagittarius 270 9 1222 - 0119 염소 Capricon 300 10 0120 - 0218 물병 Aquarius 330 11 0219 - 0320 물고기 Pisces 본 class.lunar.php 에서는 태양의 황경값은 직접 구했기 때문에 어느 싸이트보다도 더 신뢰가 갈겁니다 ^.^(에~~). 참고로, 별자리에 정보에 대한 설명(번역 문서인듯)은 http://mizian.com.ne.kr/englishwiz/library/names/zodiac/signs_of_zodiac.htm 에 있으니 참고하세요. 2-4. 일식 (문서 작성 시작후 4시간째..... 쩝 힘드네요...) 일식(일본식 음식이 절대로 아닙니다... 크크)이란? 일식은 합삭일(음력초하루)에 일어납니다. 즉 `태양-달-지구'가 일직선일때 일어나고 지구에서 볼때 달이 태양을 가리는 현상을 말합니다.(아주 옛날에는 일식이 일어나는 날을 아주 불길한 증조로 여겼다고 하죠). 그러나 모든 합삭일에 일식이 일어나질 않습니다. 그 이유는 태양과 달의 황위나 적위가 유효범위안(?)에 들어가지 않는 경우가 있기 때문이죠. (천문용어를 빌리자면 황도와 백도가 약 5도 정도 기울어져 있기 때문이죠) 대략 1년에 1-3회 정도 일식이 일어납니다. 올해(2005년) 첫번째 일식(금환식)은 한국시각으로 4월 9일 새벽에 있습니다. 그러나 불행이도 우리나라에서는 관측이 불가능하다고 하네요. 쩌~업 본론으로 들어가서, 일식계산이나 월식계산은 아주 복잡하다고 합니다. 즉 태양과 달의 황경, 적경, 황위, 적위, 지구와의 거리및 시직경(눈에 보이는 천체의 구경), 관측지의 고도와 방위각, 에고... 이런 것들을 정확하게 구해야 합니다. 본 class.lunar.php 에서는 시도는 해볼려고 했는데 그 범위를 넘어설것 같아서 처음부터 편법을 사용했습니다. 즉 2000년을 기준으로, 전후 100년씩 약 200년간의 데이터를 수치화 했습니다. 수치화의 기준의 물론 합삭일에 대한 `달의 황위'입니다. Solar Eclipses 수치화 결과) 달의 황위 절대값 범위 - T(Total,개기일식): 91번, 범위: 0.014674627401046 1.04166992331860 - A(Annular,금환식): 104번, 범위: 0.010582273616323 1.06751983653440(1.067743) - H(Hybrid(Annular/Total),두가지): 6번, 범위: 0.210833791736160 0.96121766199154(0.964180) - P(Partial,부분일식): 106번, 범위: 0.829402137331320 1.59530999819400 즉 편법으로 합삭시각에 대한 달의 황위 절대값을 구해서 판별했습니다. 참고로 http://sunearth.gsfc.nasa.gov/eclipse/eclipse.html 에 가시면 일식/월식에 대한 자료를 구할 수 있습니다. 위의 자료에서 시간 단위는 KST, GMT 도 아닌 TT(지구시) 임. - 한국표준시(KST), 세계시(UT), 지구시(TT) - KST = UT + 9시 - UT = TT - 델타T 2-5. 월식 월식은 망(`망'은 천문용어, `보름'은 일반용어)에 일어납니다. 즉 보름날에 일어나죠 그러나 정확히 말하면 보름이 아니고 `망시각'에 (전후로) 일어납니다. 일식과 반대로 `태양-지구-달' 순이고, 나머지 조건은 일식과 같거나 비슷한 원리입니다. 망시각은 일반적으로 음력 15일으로 알려져 있지만, 정확히 계산하면 음력 15일이나 16일에 있습니다. 즉 합삭에서 망으로 가는 시간이 모두 틀리기 때문이죠. (통계에 대한 그래프를 그려보면 sin() 곡선을 그립니다.) 월식 역시, 일식과 같은 편법을 사용했습니다. Lunar Eclipses 수치화 결과) 달의 황위 절대값 범위 - t(Total,개기월식): 114번, 범위: 0.010630181434784 0.52538271749109(0.527908) - p(Umbral(Partial),부분월식): 85번, 범위: 0.294681791935650 1.14183526374750 - n(Penumbral,반영식): 111번, 범위: 0.832900661071600 1.57436507046800(1.578244) 만약, 어떤 날 저녁에 빨간 보름달을 보셨다면 아마 반영식일 가능성이 높습니다. http://www.mreclipse.com/Special/photo.html 에 가면 월식,일식에 관한 사진이 있습니다(외국에는 이런 사진도 판매하네요. 부럽습니다.) http://www.mreclipse.com/MrEclipse.html 의 메인에 있는 달을 한번 감상해 보세죠(무지 이쁩니다.ㅋㅋ) 참고로, 일식과 월식에서 중요한 주기가 있는데 `사로스(Saros)주기' 입니다. (제가 천문전공이 아니라서, 설명하기가 벅찹니다. 검색엔진의 도움을 받으세요. :) 3. PHP 클래스 사용법 3-0. 함수(메소드) 요약 설명 `_'으로 시작한 클래스나 함수는 모두 내부적으로 구현하기 위해서 만든 것들입니다. - class.lunar.php 파일구조) class _solar { } // 기본유틸리티나 solar 에 관한 object class lunar extends _solar { } // lunar 에 관한 object - array _solar::_getjd($utime) JD(율리우스적일, Junlian day) 구하는 메소드 - float _solar::_deg2valid($deg) 360 도각 안으로 변환해주는 유틸리티 - integer _solar::_moon2valid($moon) 12개월 안으로 변환해주는 간단한 유틸리티 - float _solar::_deg2daytime($deg) 시간각을 시간초로 바꾸어주는 유틸리티 - string _solar::_deg2angle($deg, $singed=FALSE) 시간각 단위를 `도분초'로 변환해주는 유틸리티 - float _solar::_deg2solartime($deg) 시간각을 1태양년 시간으로 바꾸어주는 유틸리티 - array _solar::_sunl($D) 특정일($D)에 대한 태양의 시황경과 보정값(시간초)을 구하는 메소드 - array _solar::_terms($year=0, $smoon=1, $length=12) 특정 년-월에 대한 24절기 구하는 메소드 - array _solar::_constellation($y, $m, $d) 특정 년-월-일에 대한 12궁도 별자리 구하는 메소드 - string lunar::_time2mtime($time, $singed=FALSE) 특정 시각(unix timestamp)에 대한 '시분초'로 변환해 주는 유틸리티 - array lunar::_moonl($T) 특정 세기($T, Julian century)에 대한 달의 시황경 구하는 메소드 - float lunar::_moone($T) 특정 세기에 대한 달위 황위를 구하는 메소드 - array lunar::moon($utime=0) 특정 시각(unix timestamp)에 대한 달의 자세한 정보를 구하는 메소드 - float lunar::_gettd($utime) 특정 시각에 대한 태양-달의 시황경 차이(0도 기준)를 구하는 메소드 - integer lunar::_newmoon($_y, $_m=0, $_d=0, $ctime=0) 특정 년-월-일에 대한 합삭시각을 구하는 메소드(적분 비슷한 알고리즘 사용) - float lunar::_gettm($utime) 특정 시각에 대한 태양-달의 시황경 차이(180도 기준)를 구하는 메소드 - integer lunar::_fullmoon($_y, $_m=0, $_d=0) 특정 년-월-일에 대한 망시각을 구하는 메소드(적분 비슷한 알고리즘 사용) - string lunar::_solareclipse($utime) 특정 시각에 대한 일식을 판별하여 문자열로 리턴하는 메소드(내부적으로 사용됨) - string lunar::solareclipse($y, $m, $d) 특정 년-월-일에 대한 일식 여부를 판별하여 결과를 문자열로 리턴하는 메소드 - string lunar::_lunareclipse($utime) 특정 시각에 대한 월식을 판별하여 문자열로 리턴하는 메소드(내부적으로 사용됨) - string lunar::_lunareclipse($y, $m, $d) 특정 년-월-일에 대한 월식 여부를 판별하여 결과를 문자열로 리턴하는 메소드 - array lunar::tolunar($_y, $_m, $_d) 특정 년-월-일(양력)을 음력으로 계산해주는 메소드 - array lunar::_newmooninfo($_y, $_o, $_m, $_leap, $_stop=0) 양력을 음력으로 변환하기 위해서 내부적으로 사용되는 메소드 - array lunar::tosolar($_y, $_m, $_d, $_leap=0) 특정 년-월-일(음력)을 양력으로 계산해주는 메소드 - array lunar::zodiac($y, $m, $d, $lunar=0, $leap=0) 특정 년-월-일(음력,윤달 체크 가능)에 대한 12궁도 별자리를 구하는 메소드 실제 사용자(개발자)가 사용될 메소드(함수)는 5-6개 정도 밖에 되질 않습니다. 다음이 이에 대한 예제입니다. 3-1. 특정시각에 대한 달(moon)의 위치를 구하려면? $utime = time(); // 현재시각 $info = linar::moon($utime); print_r($info); 결과) Array ( [JD] => 2453431.0722106481 // 율리우스 적일 [J] => J2005.1638 // J2000.0 부터 흐른 율리우스 년도 표기 [L] => 223.800612462 // 달의 시황경 각도 [e] => -1.85002040264 // 달의 황위 각도 [RA] => 220.763546765 // 달의 적경 [d] => -17.7469781192 // 달위 적위 [mtime] => 53712.1470 // 시황경 시간각에 대한 시간초 [dtime] => 52983.2512 // 적경에 대한 시간초 [utime] => 1109684575 // 주어진 unix timestamp [date] => Tue, 01 Mar 2005 22:42:55 KST // 날짜 시간정보 KST [gmdate] => Tue, 01 Mar 2005 13:42:55 GMT // 날짜 시간정보 GMT [_L] => 223°48′2.2049″ // 달의 시황경(도분초) [_e] => -1°51′0.0734″ // 달의 황위(도분초) [_RA] => 220°45′48.7684″ // 달위 적경(도분초) [_d] => -17°44′49.1212″ // 달위 적위(도분초) [_mtime] => 14 55 12 // mtime 에 대한 시분초 [_dtime] => 14 43 03 // dtime 에 대한 시분초 ) 3-2. 양력을 음력으로 변환하려면? $lunar = lunar::tolunar(1945,8,15); // 양력 1945년 8월 15일 print_r($lunar); 결과) Array ( [0] => 1945-07-08 // 음력 1945년 7월 8일이군요/ [1] => Array // 해당 음력달에 대한 정보 ( [0] => 1945 // 음력 년도 [1] => 7 // 음력 월 [2] => 8 // 음력 일 [3] => 0 // 윤달 여부, 0(평달), 1(윤달) [4] => 29 // 29일까지 있나 보네요.. ㅎㅎ [5] => 7.10 // 주어진 날 정오에 대한 달의 월령(달의 나이?) [6] => N // 일식/월식 여부, `N' 은 아무 식도 안 일어남 ) [2] => Array // 해당 달의 합삭시각에 대한 정보 ( [0] => 1945-08-08 // 합삭 년-월-일 [1] => 09:42:10 // 합삭 시:분:초 [2] => -769994270 // 합삭시각(unix timestamp) [3] => N // 일식 여부 ) [3] => Array // 해당 달의 망시각에 대한 정보 ( [0] => 1945-08-23 // 망 년-월-일 [1] => 21:09:01 // 망 시:분:초 [2] => -768657059 // 망시각(unix timestamp) [3] => N // 월식 여부 ) [4] => Array // 다음 달의 합삭시각에 대한 정보 ( [0] => 1945-09-06 // 다음 합삭 년-월-일 [1] => 22:51:09 // 다음 합삭 시:분:초 [2] => -767441331 // 다음 합삭시각(unix timestamp) [3] => N // 다음 일식 여부 ) ) 3-3. 음력을 양력으로 변환(어르신 생신에 자주 사용되죠)하려면? $year = 1969; $month = 8; $day = 17; $leap = FALSE; // 윤달 여부 $solar = lunar::tosolar($year,$month,$day,$leap); print_r($solar); 결과) Array ( [0] => 1969-09-28 // 양력 1969년 9월 28일 이군요. [1] => 1 // 1이면 주어진 정보가 맞음, 0이면 틀린 정보(예, 윤달) [2] => 1969 // 양력 년도 [3] => 9 // 양력 월 [4] => 28 // 양력 일 [5] => 0 // 양력 요일 인덱스, 0(일), 1(월), 2(화), ... 6(토) [6] => Array // 해당 음력에 대한 합삭시각 정보 ( [0] => 1969-09-12 // 해당 음력에 대한 합삭날짜 년-월-일(초하루) [1] => 04:38:52 // 해당 음력에 대한 합삭시각 시:분:초 [2] => -9606068 // 해당 음력에 대한 합삭시각(unix timestamp) [3] => ATH // 일식여부 (금환식이거나 개기 또는 Hybird 가 있었군요) [4] => 29 // 음력 29일까지 있군요. 소월 ) ) 위의 음력은 제가 아는 어느 분의 생일이랍니다(좋은 예제일 것 같아서). ㅎㅎ 3-4. 특정 날짜에 대한 일식여부 확인해 보려면? $eclipse = lunar::solareclipse(2005,4,9); // 양력 2005년 4월 9일 print_r($eclipse); 결과) ATH // 금환식(A)이거나(확률적으로 제일 높음) 개기일식(T) 또는 Hybird(H) *부연설명) T(Total,개기일식), A(Annular,금환식) H(Hybrid(Annular/Total),두가지) P(Partial,부분일식) N(식없음). 본 프로그램에서는 정확하게 계산하기 어렵기 때문에 먼저 나열된 문자가 될 확률이 높습니다. 그러나 나열된 일식중에서 적어도 하나의 식은 반드시 일어남. :) *주의) 실제 2005년 4월 9일(KST)에는 일식이 있습니다. 식의 종류는 `Hybrid' 이고, 한국에서는 관측이 불가능하다고 합니다. T.T http://sunearth.gsfc.nasa.gov/eclipse/SEplot/SEplot2001/SE2005Apr08H.GIF 3-5. 특정 날짜에 대한 월식여부 확인해 보려면? $eclipse = lunar::lunareclipse(2005,4,24); // 양력 2005년 4월 24일 print_r($eclipse); 결과) np // 반영식(n)이거나 부분월식(p) *부연설명) t(Total,개기월식), p(Umbral(Partial),부분월식), p(Penumbral,반영식) *주의) 실제 2005년 4월 24일(KST)에는 월식이 있습니다. 식의 종류는 `반영식'이고, 한국에서도 관측이 가능하다고 합니다. :) http://sunearth.gsfc.nasa.gov/eclipse/LEplot/LEplot2001/LE2005Apr24N.GIF 3-6. 특정날짜에 대한 12궁도 별자리지를 알아보려면? $year = 1970; $month = 5; $day = 22; $lunar = TRUE; // 위의 정보(1970,5,22)이 음력이라는 의미 $leap = FALSE; // 음력이고 윤달이 아닌 평달이라는 의미 $zodiac = lunar::zodiac($year,$month,$day,$lunar,$leap); print_r($zodiac); 결과) Array ( [0] => 게 [1] => Cancer ) 위의 음력날짜는 제 생일입니다.. 히히 제 별자리는 `게자리' 이네요. :) 이상 예제를 마칩니다. 4. 후기 아직 문서가 `알파' 수준이라서 나중에 시간나면 탈고해야할 것 같습니다. (말만... ㅋㅋ) 본 PHP class 작성하는데 한 달 정도 걸린 것 같습니다. 물론 한달 내내는 아니고, 저녁에 시간나는 대로 틈틈히 하나씩 원리를 깨우쳐(?)가면서, 그리고 알고리즘을 구현하다보니 시간이 이렇게 오래 걸린것 같습니다. 아마 대부분 오차와의 싸움인것 같더군요(작성해 놓고 오류 검증...검증...) 또한, 다 계산해 놓고 PHP mktime()인지 date() 함수 버그인지는 모르겠지만 검증시 엉뚱한 결과가 나와서 첨부터 다시 손본적이 두번 있었던 것 같습니다. (중간에 포기할까 하다가도... 집념이 무섭긴 무섭네요.. 크크) PHP 버그에 대한 결과) http://linuxchannel.net/board/read.php?table=alpha&no=159 특히 수학 `적분'에 대한 알고리즘을 PHP 수식으로 풀어갈려고 하니 머리 아프더군요.. ㅠ.ㅠ (암튼 어거지로 끝냈음) ... 천문계산에서 한번 더 욕심이 가는 부분이 있는데 바로 `행성의 궤도와 좌표' 계산 입니다. 점성술과 약간 관련이 있는데, 즉 태어난 날 행성의 위치를 알고 싶은 게 제 욕심(목적)입니다. 약간의 그래프 기능이 들어 가야겠죠. 물론 시중에 이와 비슷한 프로그램(대부분 윈도우용)이 있지만 덩치가 크고, 자료도 구하기 어렵습니다. 정확하게 계산하려면 PHP 로는 도저히 무리고, 그냥 웹에서 간단하게 자신이 태어날 날 행성의 위치를 개략적으로 계산 및 그려 낼 수 있으면 좋을텐데...요... 이제보니 일월출시각, 시민/천문박명도 계산해 보고 싶네요... 욕심이.. ㅠㅠ 아... sysinfo, class.mysql.status.php 도 업데이트해야 하는데, 게으른 탓인지..ㅋㅋ EOF