[misc] 항성 적경값으로 개략적인 자정 남중날짜 계산
적경 날짜 남중날짜
0 3.21 9.21 // 추분
6 6.22 12.22 // 동지
12 9.21 3.21 // 춘분
18 12.22 6.22 // 하지
예) 적경 약 19h 항성의 자정 남중날짜는?
- 적경에서 날짜 offset = 적경/2 + 3개월 21일
- 남중날짜 offset = offset + 6개월
따라서
- 적경/2 + 9개월 + 21일 // 13보다 크면 -12개월
또는
- 적경/2 + 10월 -7일 // 13보다 크면 -12, 1보다 작으면 +12 개월 더해줌
또는
- 적경/2 - 2월 -7일 // 1보다 작으면 +12 개월 더해줌
만약 계산된 값이 또는 계산된 값이
19/2 - 2 = 7.5 (7월 중순)
따라서 7월 중순에서 7일을 빼주면 7월 8일 자정즈음에 이 별은 남중함
<?php
function dms2deg($dms)
{
if(!is_array($dms) && preg_match('/^[\d.+-]+$/',$dms)) return sprintf('%.4f',$dms);
if(is_array($dms)) $dms = implode(' ',$dms); // to string
$dms = preg_replace('/[^\d.+-]+/',' ',$dms);
list($d,$m,$s) = preg_split('/\s+/',$dms);
if(preg_match('/^-/',$d)) $deg = sprintf('%.4f',$d-$m/60-$s/3600); // -0 to < 0
else $deg = sprintf('%.4f',$d+$m/60+$s/3600);
return $deg;
}
function ra2middate($ra)
{
$ra = dms2deg($ra); // match h to d
$sec = $ra/24.0 * 365.24160 * 86400; // to seconds
$sec += 22723200; // gmmktime(0,0,0,9,21,1970);
return gmdate('n.d',$sec);
}
$f = file('ra.txt');
foreach($f AS $v)
{
if(!$v = trim($v)) continue;
list($m,$ra) = preg_split('/\s+/',$v);
if($ra == '') contine;
$m = preg_replace('/[^\d]+/','',$m);
$r[$m] = ra2middate($ra);
}
ksort($r);
foreach($r AS $k=>$v)
{
printf("%d\t%s\n",$k,$v);
}
?> |
내일은 대한입니다.
/board/delete.php:소스보기 
