Click here to Skip to main content
11,483,003 members (72,742 online)
Click here to Skip to main content

Calculate and Draw Moon Phase

, 12 Aug 2010 CPOL 99.4K 2K 37
Rate this:
Please Sign up or sign in to vote.
Calculate Moon's age and draw Moon Phase on any given day

Introduction

Many communities pay close attention to the lunar calendar next to the solar calendar, so I looked at many sites on the internet to know how to calculate the age of the moon on any given day. I found many sites offering different ways, I took what I found to give results closer to the truth.

I've noticed that most sites agree on the expense of the Julian date but do not agree to calculate the age of the moon, and found the difference between these sites up to one day, and when the moon's age is 30 days, the result is zero in some sites.

In this program, I calculated the approximate age of the moon in days and did not give attention to the parts of the day of the hours and minutes.

In order for the program would be more useful, I add PictureBox control to display the lighted part of the moon and darkness part of the moon commensurate with the age of the moon.

img039.JPG

Background

I created two projects, I wrote code of one in C# (2003) and wrote another in VB.NET (2003).

The MoonPhase project has one form (frmMoon) with the following controls:

  • MonthCalendar control (MyCalendar)
  • Button control (btnToDay)
  • Button control (btnClose)
  • PictureBox control (PicMoon)
  • Label control (lblAge)

About the Code

Convert date to Julian date:

private int JulianDate(int d, int m, int y)
{ 
    int mm, yy;
    int k1, k2, k3;
    int j;

    yy = y - (int)((12 - m) / 10);
    mm = m + 9;
    if (mm >= 12)
    {
        mm = mm - 12;
    }
    k1 = (int)(365.25 * (yy + 4712));
    k2 = (int)(30.6001 * mm + 0.5);
    k3 = (int)((int)((yy / 100) + 49) * 0.75) - 38;
    // 'j' for dates in Julian calendar:
    j = k1 + k2 + d + 59;
    if (j > 2299160)
    {
        // For Gregorian calendar:
        j = j - k3; // 'j' is the Julian date at 12h UT (Universal Time)
    }
    return j;
}

Calculate the approximate moon's age in days:

private double MoonAge(int d, int m, int y)
{ 
    int j = JulianDate(d, m, y);
    //Calculate the approximate phase of the moon
    ip = (j + 4.867) / 29.53059;
    ip = ip - Math.Floor(ip); 
    //After several trials I've seen to add the following lines, 
    //which gave the result was not bad 
    if(ip < 0.5)
        ag = ip * 29.53059 + 29.53059 / 2;
    else
        ag = ip * 29.53059 - 29.53059 / 2;
    // Moon's age in days
    ag = Math.Floor(ag) + 1;
    return ag;
}

Draw moon:

private void DrawMoon()
{
    int Xpos, Ypos, Rpos;
    int Xpos1, Xpos2;
    double Phase;

    Phase = ip;
    // Width of 'ImageToDraw' Object = Width of 'PicMoon' control    
    int PageWidth = PicMoon.Width; 
    // Height of 'ImageToDraw' Object = Height of 'PicMoon' control    
    int PageHeight = PicMoon.Height; 
    // Initiate 'ImageToDraw' Object with size = size of control 'PicMoon' control    
    Bitmap ImageToDraw = new Bitmap(PageWidth, PageHeight); 
    // Create graphics object for alteration.    
    Graphics newGraphics = Graphics.FromImage(ImageToDraw);

    Pen PenB = new Pen(Color.Black); // For darkness part of the moon
    Pen PenW = new Pen(Color.White); // For the lighted part of the moon

    for (Ypos=0; Ypos<= 45; Ypos++)
    {
        Xpos = (int)(Math.Sqrt(45*45 - Ypos*Ypos));
        // Draw darkness part of the moon        
        Point pB1 = new Point(90-Xpos, Ypos+90);
        Point pB2 = new Point(Xpos+90, Ypos+90);
        Point pB3 = new Point(90-Xpos, 90-Ypos);
        Point pB4 = new Point(Xpos+90, 90-Ypos);
        newGraphics.DrawLine(PenB, pB1, pB2); 
        newGraphics.DrawLine(PenB, pB3, pB4);        
       // Determine the edges of the lighted part of the moon       
       Rpos = 2 * Xpos;
       if (Phase < 0.5)
       {
           Xpos1 = - Xpos;
           Xpos2 = (int)(Rpos - 2*Phase*Rpos - Xpos);
       }
       else
       {
          Xpos1 = Xpos;
          Xpos2 = (int)(Xpos - 2*Phase*Rpos + Rpos);
       }       
       // Draw the lighted part of the moon       
       Point pW1 = new Point(Xpos1+90, 90-Ypos);
       Point pW2 = new Point(Xpos2+90, 90-Ypos);
       Point pW3 = new Point(Xpos1+90, Ypos+90);
       Point pW4 = new Point(Xpos2+90, Ypos+90);
       newGraphics.DrawLine(PenW, pW1, pW2);
       newGraphics.DrawLine(PenW, pW3, pW4);
    }    
    // Display the bitmap in the picture box.    
    PicMoon.Image = ImageToDraw;    
    // Release graphics object    
    PenB.Dispose();
    PenW.Dispose();
    newGraphics.Dispose();
    ImageToDraw = null;
} 

You can go back to the source file to read the VB.NET code after you extract Moon_VB.zip file.

Final Words

If you have any idea or find another code to calculate the age of the moon, please tell me. Thanks to Code Project and thanks to all.

Mostafa Kaisoun
m_kaisoun@hotmail.com

License

This article, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)

Share

About the Author

No Biography provided

Comments and Discussions

 
QuestionAdd more feature Pin
Christianirwan26-Jun-13 1:39
memberChristianirwan26-Jun-13 1:39 
AnswerRe: Add more feature Pin
Mostafa Kaisoun26-Jun-13 4:02
memberMostafa Kaisoun26-Jun-13 4:02 
SuggestionMoon phase Pin
Edward Omowa16-Apr-13 12:29
memberEdward Omowa16-Apr-13 12:29 
GeneralRe: Moon phase Pin
Mostafa Kaisoun16-Apr-13 15:43
memberMostafa Kaisoun16-Apr-13 15:43 
GeneralRe: Moon phase Pin
t9mike28-Mar-15 16:44
membert9mike28-Mar-15 16:44 
Questionfrmoon Pin
Edward Omowa16-Apr-13 12:23
memberEdward Omowa16-Apr-13 12:23 
GeneralMy vote of 5 Pin
DanielSheets5-Mar-13 8:31
memberDanielSheets5-Mar-13 8:31 
GeneralRe: My vote of 5 Pin
Mostafa Kaisoun6-Mar-13 6:40
memberMostafa Kaisoun6-Mar-13 6:40 
GeneralMy vote of 5 Pin
manoj kumar choubey1-Mar-12 23:18
membermanoj kumar choubey1-Mar-12 23:18 
QuestionError, and another algorithm Pin
Robert Baker II7-Sep-11 6:50
memberRobert Baker II7-Sep-11 6:50 
AnswerRe: Error, and another algorithm Pin
Robert Baker II7-Sep-11 6:55
memberRobert Baker II7-Sep-11 6:55 
AnswerRe: Error, and another algorithm Pin
Mostafa Kaisoun7-Sep-11 10:03
memberMostafa Kaisoun7-Sep-11 10:03 
GeneralDraw Moon Phase Pin
Member 799874712-Jun-11 11:02
memberMember 799874712-Jun-11 11:02 
GeneralRe: Draw Moon Phase Pin
Mostafa Kaisoun13-Jun-11 15:24
memberMostafa Kaisoun13-Jun-11 15:24 
GeneralMy vote of 5 Pin
dim13123-Dec-10 5:38
memberdim13123-Dec-10 5:38 
GeneralMy vote of 5 Pin
Paul Conrad15-Sep-10 14:50
memberPaul Conrad15-Sep-10 14:50 
GeneralGood job - Thanks for sharing... Pin
Steve Hageman16-Aug-10 20:11
memberSteve Hageman16-Aug-10 20:11 
GeneralRe: Good job - Thanks for sharing... Pin
Mostafa Kaisoun17-Aug-10 9:34
memberMostafa Kaisoun17-Aug-10 9:34 
GeneralGood article Pin
dwilliss10-Aug-10 16:50
memberdwilliss10-Aug-10 16:50 
GeneralRe: Good article Pin
Mostafa Kaisoun11-Aug-10 4:58
memberMostafa Kaisoun11-Aug-10 4:58 
GeneralMy vote of 5 Pin
myker10-Aug-10 5:53
membermyker10-Aug-10 5:53 
GeneralRe: My vote of 5 Pin
Mostafa Kaisoun11-Aug-10 4:55
memberMostafa Kaisoun11-Aug-10 4:55 
GeneralDownload MoonPhase C# code Pin
Mostafa Kaisoun9-Aug-10 7:43
memberMostafa Kaisoun9-Aug-10 7:43 
Generalnot Bad, Keep going Mostafa... Pin
Omar Amin Ibrahim9-Aug-10 5:11
memberOmar Amin Ibrahim9-Aug-10 5:11 
GeneralRe: not Bad, Keep going Mostafa... Pin
Mostafa Kaisoun9-Aug-10 5:48
memberMostafa Kaisoun9-Aug-10 5:48 

General General    News News    Suggestion Suggestion    Question Question    Bug Bug    Answer Answer    Joke Joke    Rant Rant    Admin Admin   

Use Ctrl+Left/Right to switch messages, Ctrl+Up/Down to switch threads, Ctrl+Shift+Left/Right to switch pages.

| Advertise | Privacy | Terms of Use | Mobile
Web02 | 2.8.150520.1 | Last Updated 12 Aug 2010
Article Copyright 2010 by Mostafa Kaisoun
Everything else Copyright © CodeProject, 1999-2015
Layout: fixed | fluid