Расчет времени истинного полдня на PHP с графиками изменения

[LeonidOS]

По поводу размещения мы посоветуемся. А обсуждать конечно можете здесь.

[Novoster]

Спасибо за ответ.

Думаю, мы тогда будем пересылать программы по почте,
чтобы вас не нагружать.

По большому счету, это мало кому интересно. Важен будет
результат, то есть готовый продукт, если он получится.
Надеюсь, что получится.

[Сергей]

Novoster
Нашел уже исходники на Паскале из книги
"Астрономия на персональном компьютере":
http://oleg.proc.ru/books/pascal/astro_pas.zip

Да, это действительно проги из книги, но у меня эти проги на С, так решил издатель.

[Александр]

Ни в коем случае!!!
С точки зрения безопасности, я категорически против!
Опыт у нас был - и очень печальный....

у нас на сервере ваши исполняемые файлы?

[LeonidOS]

Ни в коем случае!!!
С точки зрения безопасности я категорически против!
Опыт у нас был - и очень печальный....

Я о том же... А проверить на 100% мы не сможем

Не обижайтесь, но Вас мы совсем не знаем. А угрозы в наш адрес были...

[Novoster]

Не обижайтесь, но Вас мы совсем не знаем. А угрозы в наш адрес были...

Господа, все нормально. Я совсем не обижаюсь. Хочется очень быстро и толково решить вопрос закатов/восходов и предоставить рабочую программу. Реально нет смысла углубляться в обиды. Потому на все согласен.

Но может хотя бы картинки можно размещать на вашем сервере в виде Attachments, как это сделал Сергей? Все же интерфейс обсуждать было бы не лишним.

Леонид, книгу "Астрометрические алгоритмы" Меюса скачал. Спасибо.
Вчера писал об этом, но почему-то сначала отправилось два сообщения (и-нет сбоит), а когда я удалил второе, то исчезло и первое. Нашел главу в которой речь идет именно об этом:
Chapter 14
Rising, Transit and Setting

Прочел все, сключая сноски. Смутило вот это:

From an accurate ephemeris, we take the following values:

"Не в глаз, так в бровь" - если не из календаря берем данные,
так из точных эфемерид. Поэтому разбор алгоритмов Меюса продолжим после разбора SUNSET.PAS

[Novoster]

Закончил делать интерфейс и ввод данных в программу.
Перехожу к расчетам.

По ходу наткнулся на сайт расчета времени восхода/заката/зенита:
http://travel.org.ua/sunrise/index.php

А здесь координаты некоторых городов:
http://travel.org.ua/sunrise/coordinates.php

[Novoster]

Приехали. В PHP 5 есть встроенная функция определения заката/рассвета.

http://shublog.ru/php/vremya-voskhoda-v ... go-solnca/

Итак, речь пойдёт о двух функциях PHP с названием date_sunrise и date_sunset. Рассмотрим их примеры.

date_sunrise(time(),SUNFUNCS_RET_STRING,35.68, 139.6, 90+50/60, 9);
date_sunset(time(),SUNFUNCS_RET_STRING,35.68, 139.6, 90+50/60, 9);

Функции принимают 6 параметров:

1.timestamp. Время и дата, для которой определяется восход и закат в формате UNIX.

2.format. Формат возвращаемого значения. Может принимать одно из трёх значений: SUNFUNCS_RET_STRING (возвращает строку со значением от 00:00 до 23:59), SUNFUNCS_RET_DOUBLE (возвращает время в виде величины с плавающей точкой от 0 до 23.99), SUNFUNCS_RET_TIMESTAMP (возвращает время в UNIX формате).

3.latitude. Или по русски “широта”. Для северной широты используются положительные числа, для южной – отрицательные.

4.longitude. Или “долгота”. Для восточное долготы используются положительные числа, для западной – отрицательные.

5.zenith. Зенит солнца (в градусах). Выделяют четыре основных значения. Официальный (offical): 90 градусов 50’. Гражданский (civil): 96 градусов. Навигационный (nautical): 102 градуса. Астрономический (astronomical): 108 градусов.

6.gmtoffset. Смещение от времени по Гринвичу в часах. Иными словами часовой пояс по GMT.
Теперь смотрим на пример. В качестве первого параметра мы используем функцию time(). В качестве второго параметра выбираем SUNFUNCS_RET_STRING.

Третий и четвёртый параметр – координаты нужного нам города. В данном случае Токио. Однако не забываем, что 35°41’ будет записано не как 35.41, а как 35.68. Думаю не надо объяснять почему. В интернете я нашёл интересный сервис, который даёт и обычные координаты, и десятичные. Например, для Токио я нашёл координаты здесь.

Пятый параметр указываем по официальному – 90+50/60 (90°50’). Шестой параметр – для Токио часовой пояс +9. Ставим цифру 9.


Теперь мы знаем дату восхода и заката для Токио. По ссылке вы найдёте работающий пример, где написан его код, и указано время для Токио и для Санкт-Петербурга.

Сегодня в Санкт-Петербурге восход в 10:00, а закат в 16:09
В то время как в Токио восход в 06:50, а закат в 16:41


А вот и код.

< ?
$spb_sunrise = date_sunrise(time(), SUNFUNCS_RET_STRING, 60.25, 29.6, 90+50/60, 3);
$spb_sunset = date_sunset(time(), SUNFUNCS_RET_STRING, 60.25, 29.6, 90+50/60, 3);

$tokyo_sunrise = date_sunrise (time(), SUNFUNCS_RET_STRING, 35.68, 139.6, 90+50/60, 9);
$tokyo_sunset = date_sunset (time(), SUNFUNCS_RET_STRING, 35.68, 139.6, 90+50/60, 9);

echo "Сегодня в Санкт-Петербурге восход в <b>".$spb_sunrise."</b>, а закат в <b>".$spb_sunset."</b><br>";
echo "В то время как в Токио восход в <b>".$tokyo_sunrise."</b>, а закат в <b>".$tokyo_sunset."</b>";
? >

Но на Delphi уже самому интересно доделать.

[Novoster]

Готово

Отправляю Сергею рабочую программу.

[Типичный марсианин]

Если интересно, могу выложить книгу "Астрометрические алгоритмы" Меюса

А еще есть книга Астрономические формулы для калькуляторов, того же автора. У меня такая имеется дома (честно прихватизировал из библиотеки-- в книжных было просто не найти). Хотел тут недавно ее просканировать-- типа иметь в эл. виде, порылся в сети и нашел готовую, ну и надобность в сканировании отпала.
Интересно, что некоторые проги (типа HNSKY) пишут с открытыми кодами и на сайте лежит листинг проги на Паскале. Думаю, что авторы не обидятся, если кто-то будет использовать готовые блоки из их работы (ну хотя бы как образец или справочник). Формулы из программы очень похожи на формулы из книжки Мееса.

[LeonidOS]

А еще есть книга Астрономические формулы для калькуляторов, того же автора.

Да, но Алгоритмы полнее. Хотя опять же, все это поверхностно. Алгоритмы тоже на месте не стоят

[Novoster]

Господа, алгоритмы может быть для кого-то прекрасны. Но я честно
признаюсь, что без исходников

O. Montenbruck, T. Pfleger
ASTRONOMY ON THE PERSONAL
COMPUTER
(3rd edition 1997)

меня вряд ли спасли бы алгоритмы. Я наивно полагал, что смогу в
них разобраться. Возможно за пару недель что-то и вышло бы, но
это слишком сложно и роскошно две недели считать все эти циферки.

Чтобы было понятно, о чем я говорю, хочу показать процедуры
(фрагменты кода), которые позволяют делать расчеты без всяких
дополнительных календарей, альманахов и таблиц.

Грубые Солнечные координаты с точностью 1':

Код: Выделить всё

(*-----------------------------------------------------------------------*)
(* MINI_SUN: low precision solar coordinates (approx. 1')                *)
(*           T  : time in Julian centuries since J2000                   *)
(*                ( T=(JD-2451545)/36525 )                               *)
(*           RA : right ascension (in h; equinox of date)                *)
(*           DEC: declination (in deg; equinox of date)                  *)
(*-----------------------------------------------------------------------*)
procedure TFormMain.MINI_SUN(T: REAL; var RA, DEC: REAL);
  CONST P2  = 6.283185307; COSEPS=0.91748; SINEPS=0.39778;
  VAR   L,M,DL,SL,X,Y,Z,RHO: REAL;
  FUNCTION FRAC(X:REAL):REAL;
    BEGIN  X:=X-TRUNC(X); IF (X<0) THEN X:=X+1; FRAC:=X  END;
  BEGIN
    M  := P2*FRAC(0.993133+99.997361*T);
    DL := 6893.0*SIN(M)+72.0*SIN(2*M);
    L  := P2*FRAC(0.7859453 + M/P2 + (6191.2*T+DL)/1296E3);
    SL := SIN(L);
    X:=COS(L); Y:=COSEPS*SL; Z:=SINEPS*SL; RHO:=SQRT(1.0-Z*Z);
    DEC := (360.0/P2)*ARCTAN(Z/RHO);
    RA  := ( 48.0/P2)*ARCTAN(Y/(X+RHO)); IF (RA<0) THEN RA:=RA+24.0;
  END;


Числа, которые здесь используются, это не просто какие-то частно
взятые данные, а константы, о которых не было ни с лова в
литературе, которую удалось изучить.

С Луной еще веселее:

Код: Выделить всё

(*-----------------------------------------------------------------------*)
(* MINI_MOON: low precision lunar coordinates (approx. 5'/1')            *)
(*            T  : time in Julian centuries since J2000                  *)
(*                 ( T=(JD-2451545)/36525 )                              *)
(*            RA : right ascension (in h; equinox of date)               *)
(*            DEC: declination (in deg; equinox of date)                 *)
(*-----------------------------------------------------------------------*)
procedure TFormMain.MINI_MOON(T: REAL; var RA, DEC: REAL);
  CONST P2 =6.283185307; ARC=206264.8062;
        COSEPS=0.91748; SINEPS=0.39778;  (* cos/sin(obliquity ecliptic)  *)
  VAR   L0,L,LS,F,D,H,S,N,DL,CB    : REAL;
        L_MOON,B_MOON,V,W,X,Y,Z,RHO: REAL;
  FUNCTION FRAC(X:REAL):REAL;
    (* with some compilers it may be necessary to replace *)
    (* TRUNC by LONG_TRUNC oder INT if T<-24!             *)
    BEGIN  X:=X-TRUNC(X); IF (X<0) THEN X:=X+1; FRAC:=X  END;
  BEGIN
    (* mean elements of lunar orbit *)
    L0:=   FRAC(0.606433+1336.855225*T); (* mean longitude Moon (in rev) *)   
    L :=P2*FRAC(0.374897+1325.552410*T); (* mean anomaly of the Moon     *)
    LS:=P2*FRAC(0.993133+  99.997361*T); (* mean anomaly of the Sun      *)
    D :=P2*FRAC(0.827361+1236.853086*T); (* diff. longitude Moon-Sun     *)
    F :=P2*FRAC(0.259086+1342.227825*T); (* mean argument of latitude    *)
    DL := +22640*SIN(L) - 4586*SIN(L-2*D) + 2370*SIN(2*D) +  769*SIN(2*L)
          -668*SIN(LS)- 412*SIN(2*F) - 212*SIN(2*L-2*D) - 206*SIN(L+LS-2*D)
          +192*SIN(L+2*D) - 165*SIN(LS-2*D) - 125*SIN(D) - 110*SIN(L+LS)
          +148*SIN(L-LS) - 55*SIN(2*F-2*D);
    S := F + (DL+412*SIN(2*F)+541*SIN(LS)) / ARC;
    H := F-2*D;
    N := -526*SIN(H) + 44*SIN(L+H) - 31*SIN(-L+H) - 23*SIN(LS+H)
         + 11*SIN(-LS+H) -25*SIN(-2*L+F) + 21*SIN(-L+F);
    L_MOON := P2 * FRAC ( L0 + DL/1296E3 ); (* in rad *)
    B_MOON := ( 18520.0*SIN(S) + N ) / ARC; (* in rad *)
    (* equatorial coordinates *)
    CB:=COS(B_MOON);
    X:=CB*COS(L_MOON); V:=CB*SIN(L_MOON); W:=SIN(B_MOON);
    Y:=COSEPS*V-SINEPS*W; Z:=SINEPS*V+COSEPS*W; RHO:=SQRT(1.0-Z*Z);
    DEC := (360.0/P2)*ARCTAN(Z/RHO);
    RA  := ( 48.0/P2)*ARCTAN(Y/(X+RHO)); IF RA<0 THEN RA:=RA+24.0;
  END;


Как вам это:

Код: Выделить всё

    L0:=   FRAC(0.606433+1336.855225*T); (* mean longitude Moon (in rev) *)   
    L :=P2*FRAC(0.374897+1325.552410*T); (* mean anomaly of the Moon     *)
    LS:=P2*FRAC(0.993133+  99.997361*T); (* mean anomaly of the Sun      *)
    D :=P2*FRAC(0.827361+1236.853086*T); (* diff. longitude Moon-Sun     *)
    F :=P2*FRAC(0.259086+1342.227825*T); (* mean argument of latitude    *)


Абсолютно нереальные вычисления - нигде ничего подобного не описано!

И все это на базе "Юлианских столетий", которые посчитать тоже
совсем непросто.

Вы не представляете, как высоко я прыгал, когда весь код, взятый
у O. Montenbruck и T. Pfleger заработал и выдал правильные результаты!

Там такие дебри... На их написание в моем исполнении ушла бы не одна неделя. Другими словами, с задачей я бы не справился - бросил
бы.

[Novoster]

Спасибо Сергею, который направил меня в верном направлении. Сегодня постараюсь доделать новую версию программы и ответить на два его письма, на форуме. Он прислал их мне на e-mail.

[Александр]

Удачи вам в этом нелегком труде!
На голову шлем наденьте на всякий

[Novoster]

Удачи вам в этом нелегком труде!
На голову шлем наденьте на всякий

Спасибо за рекомендацию. Молодость - время сил. Можно
и шлемом погриметь. Как раз сегодня потратил море времени
на разгром местного провайдера - страницы с 10 раза грузятся.
Потому не успел закончить новую версию и ответить Сергею.
Завтра уже постараюсь справиться.

[Novoster]

Здравствуйте, Сергей и все, кто читает эту тему.

По ходу работы над программой нашел интересные сайты.

Анимация солнечного затемения
Animation of Annular Solar Eclipse, January 15, 2010
http://www.timeanddate.com/eclipse/anim ... -2010.html

Затмения Луны и Солнца на весь год:
http://www.timeanddate.com/eclipse/

На этом же сайте море информации:
часовые пояса, смена времени Daylight Saving Time,
фазы Луны и кое-что еще по астрономии.

Обои голубой Луны:

[Novoster]

Благодарю Вас за предоставленную мне возможность ПЕРВЫМ посмотреть результат Вашей работы, но мне придется, наверное, Вас разочаровать, так как я не смогу достоверно

Сергей, вашей проверки вполне достаточно.
Особенно на этом первом этапе.

Мне всегда импонировали люди, которые в короткий срок могли разобраться с сущностью какой либо проблемы, в нашем случае – это Юлианский период, который является основным аргументом для всех вычислений, и я рад за Вас, что Вы откопали в ИНЕТЕ программы, написанные на Паскале, мне пришлось немного труднее. Я думаю, что теперь мои исходники Вам не потребуются, раз найдены оригиналы.

Спасибо за добрые слова. Благодаря тому, что вы
направили меня на верный путь, действительно программа
по расчету восхода/заката оказалась
простым делом. Конечно, только благодаря исходникам!

По Луне.
Еще внимательно почитайте в книге о восходе и заходе Луны, так как, в некоторые дни, Луна либо не восходит, либо не заходит, как это получилось 01/01/2010. В этот день на 00 часов Луна была уже на небе, то есть не всходила, а только заходила.

Сергей, спасибо, что обратили внимание. Нужно будет еще почитать
и разобраться. Действительно бывают два дня, когда программа не
выдает никакие данные по заходу или восходу. Эти дни идут с
периодичностью 12..16 дней.

Но вот с датами у меня ситуация не совпадает.

Код: Выделить всё

   Rise   Set   
31.12.2009         0
01.01.2010         0
02.01.2010         0
03.01.2010         0
04.01.2010         0
05.01.2010   -----      0
06.01.2010   ===      1
07.01.2010         2
08.01.2010         3
09.01.2010         4
10.01.2010         5
11.01.2010         6
12.01.2010         7
13.01.2010         8
14.01.2010         9
15.01.2010         10
16.01.2010         11
17.01.2010         12
18.01.2010         13
19.01.2010         14
20.01.2010         15
21.01.2010      -----   16
22.01.2010      ===   1
23.01.2010         2
24.01.2010         3
25.01.2010         4
26.01.2010         5
27.01.2010         6
28.01.2010         7
29.01.2010         8
30.01.2010         9
31.01.2010         10
01.02.2010         11
02.02.2010         12
03.02.2010   -----      13
04.02.2010   ===      1
05.02.2010         2
06.02.2010         3
07.02.2010         4
08.02.2010         5
09.02.2010         6
10.02.2010         7
11.02.2010         8
12.02.2010         9
13.02.2010         10
14.02.2010         11
15.02.2010         12
16.02.2010         13
17.02.2010         14
18.02.2010         15
19.02.2010      -----   16
20.02.2010      ===   1
21.02.2010         2
22.02.2010         3
23.02.2010         4
24.02.2010         5
25.02.2010         6
26.02.2010         7
27.02.2010         8
28.02.2010         9
01.03.2010         10
02.03.2010         11
03.03.2010         12
04.03.2010   -----      13
05.03.2010   ===      1
06.03.2010         2
07.03.2010         3
08.03.2010         4
09.03.2010         5
10.03.2010         6
11.03.2010         7


К слову сказать мне не удалось связать дни, в которых получается
пустота по закату или восходу, ни с фазами луны, ни с чем-то еще.
Смущает так же и то, что дни без восхода и заката идут один за
другим. Если вы сможете пояснить, то будет здорово!

О временах.
Возможно Вы обратили внимание и учли это, что в астрономии ВСЕ просчитывается по Гринвичскому времени (UTC или GMT), а в месте наблюдения существуют
Среднее местное время,
Истинное местное время
Поясное время (зимнее/летнее)
Поэтому при вычислении времен восхода, заката, полдня следует учитывать эти нюансы.

Да, к счастью авторы программы это учли, а я заметил. И все
работает правильно.

Чтобы точно проверить работу Вашей программы, я посоветую найти в ИНЕТЕ, хотя бы даже и на ASTROGALAXY, программу планетария, а в этих программах обязательно должны вычисляться все времена для любой точки Земли, по крайней мере я так думаю.

Спасибо за рекомендацию. Вы правы. Нашел несколько сайтов и
проверил. Отличие не более 1 минуты по восходу/закату. Так что
можно полагать, что расчеты ведутся верно.

Код: Выделить всё

В настоящий момент я немного поостыл к астрономии, меня сейчас интересует создание устройства на базе КОМПА, контролирующее линейные сигналы низкой частоты, а более проще – переговорное устройство типа телефона, а заодно подучиться работать со звуковой картой (почти что на уровне оборудования).

Хочу обратить ваше внимание на сайт программиста
Бояна Митова. Боян специализируется в обработке и отображении
аудио/видео данных. Так же он прекрасный специалист по 2D
графике, разработал супер библиотеку по iGDI+. Если будут
вопросы конкретно по графике - обращайтесь. Я его уже замучил
вопросами и кое-что изучил. Так что буду рад помочь. Тем более,
что Боян любит именно Delphi больше, чем C++, что меня очень
радует. И любит он Delphi не зря, так как на нем делать 2D графику
намного проще.

Обращу Ваше внимание на мою ошибку по поводу восхода и захода Луны.
Недорозумение возникло из-за ассоциаций - подсознательно предполагается, что заход должен быть позже восхода, а, следовательно, расположение клетки со временем захода ожидается после клетки со временем восхода. Это меня и смутило.

Честно сказать, меня больше смущает тот факт, что исходники
не считают целых 2 дня в том случае, когда восхода/захода нет.
Этого пока понять не могу.

По ходу работы над программой случилось то, чего боялся больше
всего. Я увлекся. И теперь появилось желание развивать программу
до полноценного астрологического календаря. Дело это будет
сложное но идея зреет давно и, пожалуй, можно начинать потихоньку,
чтобы к лету, не спеша кое-что уже сделать. Но об этом чуть позже.

Пока что посылаю вам новую версию 2.11 программы Zenith (Зенит).
Добавлены почти все часовые пояса, коих 34. Не учитывал пару
поясов с 15-минутным смещением. Сделал возможность добавлять
города в базу данных. А вот удалять нужно вручную из INI-файла
zencit.ini

Так же сделал запись положения окна и всех настроек, а так же
введенные данные по широте/долготе.

Недавно мне рассказали о новом бесплатном фаловом архиве.
EX.UA
Первое, что решил посмотреть - документальные фильмы.
Там есть так же немного астрономии:
http://www.ex.ua/view/1987?r=23775

Можете скачать, пока не удалили море интересных материалов от
Discovery и National Geographic.

Вчера посмотрел Астрономия между А и Я
Вам тоже может понравиться:
http://www.ex.ua/view/475290?r=1987,23775

[Novoster]

Алексей, здравствуйте!
Посмотрел Ваше решение, и вот какие у меня возникли вопросы и пожелания:
Минеральные Воды по латинице так и пишутся Mineralnye Vody, но это так ...
Опция TIME - мне понравилась, дает широкую перспективу.
Выбор ЗИМНЕГО и ЛЕТНЕГО времени желательно сделать в автомате,
Кнопка Add не учитывакт уже введенные города и нет возможности вводить земные координаты места,
Так же, желательно, уточнить, какое время отображается - ПОЯСНОЕ или СРЕДНЕЕ МЕСТНОЕ.
А так, достаточно прилично, хорошая работа!
С уважением, Сергей.

Здравствуйте, Сергей.

Минеральные Воды по латинице так и пишутся Mineralnye Vody,
но это так ...

Идея была в том, чтобы писать именно на
английском Moscow, Mineral Waters и т.д. Но идея ваша тоже хороша.
Нужно просто решить, на сколько нужно будет ориентироваться на
зарубежного пользователя. Пока просто по инерции сделал на
английском. Или даже не по инерции, а от того, что сильно замаяли
меня проблемы с русским на всех компьютерах, кроме родного
русского Windows. Много вопросов от зарубежных пользователей
получаю по этому случаю, хотя и на тему другой программы.

Опция TIME - мне понравилась, дает широкую перспективу.

Если вы о Time Zone, то это должно быть полезно. Рад, что
вам понравилось. Снова-таки думал о том, чтобы охватить разных
пользователей, даже если будем пользоваться только мы с вами.

Выбор ЗИМНЕГО и ЛЕТНЕГО времени желательно сделать в автомате

Снова та же проблема. Далеко не во всех странах
сдвиг идет тогда, когда в России, Украине и т.д. В США, например,
это время на 3 недели отличается. А иногда сдвиг вообще составляет
2 часа вместо 1. Одним словом, сдвиг если и делать автоматом, то только
буквально для пары стран. Не более. Иначе неизбежны ошибки.
Более того, время смены ЗИМНЕГО и ЛЕТНЕГО времени может меняться.
Это как раз хорошо описано на сайте
http://www.timeanddate.com/

Кнопка Add не учитывакт уже введенные города и

Это так. Как-то не подумал об этом. Спасибо.

нет возможности вводить земные координаты места,

Как это нет? Нужно ввести координаты,
задать смещение по GMT и нажать Add.
Введенные координаты введутся и нужно лишь
добавить название города. Такая была задумка.

Так же, желательно, уточнить, какое время отображается - ПОЯСНОЕ или СРЕДНЕЕ МЕСТНОЕ.

Не понял. По идее местное время. А какое еще может быть?

А так, достаточно прилично, хорошая работа!

Спасибо. Буду благодарен, если сможете пояснить проблему
с пустыми днями восхода/заката.

[Novoster]

Пример ввода координат нового города:

[Сергей]

Я знаю следующую классификацию времен:
Звездное гринвичское время
Звездное местное время
Эти времена определяются часовым углом (экваториальная система координат 1-ого типа) точки весеннего равноденствия.

Среднее солнечное время в Гринвиче (UTC)
Истинное солнечное время в Гринвиче (надо учитывать уравнение времени)

Местное среднее время (часовой угол среднего Солнца в месте наблюдения)
Местное истинное время (часовой угол истинного Солнца в месте наблюдения)
Поясное (летнее или зимнее) время в точке наблюдения.

Почему упоминается Гринвич, потому что принято вести расчеты небесных явлений относительно Гринвича, а затем пересчитывать полученные результаты к месту наблюдения.