|
|
using System; |
|
|
|
|
|
namespace UnityEngine.AzureSky |
|
|
{ |
|
|
[Serializable] |
|
|
public class AzureSkyTimeOfDayComponent |
|
|
{ |
|
|
//Used by Calendar. |
|
|
public int selectableDayInt = 15; |
|
|
public int dayOfYear = 1; |
|
|
public string[] selectableDayList = new string[] |
|
|
{ |
|
|
"0", "1", "2", "3", "4" , "5" , "6", "7" , "8", "9", "10", |
|
|
"11", "12", "13", "14" , "15" , "16", "17" , "18", "19", "20", |
|
|
"21", "22", "23", "24" , "25" , "26", "27" , "28", "29", "30", |
|
|
"31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41" |
|
|
}; |
|
|
|
|
|
//Used by Calendar and Time of Day tab. |
|
|
public float hour = 6.0f; |
|
|
public int day = 15; |
|
|
public int month = 9; |
|
|
public int year = 1903; |
|
|
public int dayInMonth = 30; |
|
|
public float utc = 0; |
|
|
public float latitude = 0; |
|
|
public float longitude = 0; |
|
|
|
|
|
//Used by Time of Day tab. |
|
|
public int timeMode = 0; |
|
|
public float dayCycle = 10.0f;//In minutes |
|
|
public bool setTimeByCurve = false; |
|
|
public float hourByCurve = 6.0f; |
|
|
public AnimationCurve dayCycleCurve = AnimationCurve.Linear(0, 0, 24, 24); |
|
|
public float lst; |
|
|
private float m_radians, m_radLatitude, m_sinLatitude, m_cosLatitude; |
|
|
|
|
|
//Internal uses. |
|
|
private Vector2 m_hourAndMinutes; |
|
|
private Vector3 m_date; |
|
|
private DateTime m_dateTime; |
|
|
private int m_daysInMonth = 30; |
|
|
private int m_dayOfWeek = 0; |
|
|
private int m_fixLeapYear = 0; |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController and AzureSkyControllerEditor scripts to correctly adjusts the day buttons and day profiles on the calendar. |
|
|
/// </summary> |
|
|
/// <param name="dayList">List containing 42 strings with the day numbers to rename the calendar buttons.</param> |
|
|
public void UpdateCalendar (string[] dayList) |
|
|
{ |
|
|
m_daysInMonth = DateTime.DaysInMonth(year, month); |
|
|
if (day > m_daysInMonth) { day = m_daysInMonth; } |
|
|
if (day < 1) { day = 1; } |
|
|
|
|
|
m_dateTime = new DateTime(year, month, 1); |
|
|
m_dayOfWeek = (int)m_dateTime.DayOfWeek; |
|
|
selectableDayInt = day - 1 + m_dayOfWeek; |
|
|
|
|
|
m_dateTime = new DateTime(year, month, day); |
|
|
for (int i = 0; i < dayList.Length; i++) |
|
|
{ |
|
|
if (i < m_dayOfWeek || i >= (m_dayOfWeek + m_daysInMonth)) |
|
|
{ |
|
|
dayList[i] = ""; |
|
|
continue; |
|
|
} |
|
|
m_dateTime = new DateTime(year, month, (i - m_dayOfWeek) + 1); |
|
|
dayList[i] = m_dateTime.Day.ToString(); |
|
|
} |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Get the current day of the week and return an integer between 0 and 6. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public int GetDayOfWeek () |
|
|
{ |
|
|
m_dateTime = new DateTime(year, month, day); |
|
|
return (int)m_dateTime.DayOfWeek; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyControllerEditor script to get the day of the week from a specific day of the month. Returns an integer between 0 and 6. |
|
|
/// </summary> |
|
|
/// <param name="dayNumber">Number of the day in the month.</param> |
|
|
/// <returns></returns> |
|
|
public int GetDayOfWeek (int dayNumber) |
|
|
{ |
|
|
m_dateTime = new DateTime(year, month, dayNumber); |
|
|
return (int)m_dateTime.DayOfWeek; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Returns the current day number of a 366-day year. |
|
|
/// The return value is fix when it is not leap year, this method is ideal for accessing, adding, or removing day profiles from the calendar. |
|
|
/// If you want to get the correct number for the day in the year, then use System.DateTime.DayOfYear instead. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public int GetDayOfYear () |
|
|
{ |
|
|
m_dateTime = new DateTime(year, month, day); |
|
|
dayOfYear = m_dateTime.DayOfYear; |
|
|
if ((dayOfYear > 59 && !IsLeapYear())) |
|
|
{ |
|
|
m_fixLeapYear = 1; |
|
|
} |
|
|
else { m_fixLeapYear = 0; } |
|
|
return m_dateTime.DayOfYear - 1 + m_fixLeapYear; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Return true if the current year is a leap year. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public bool IsLeapYear () |
|
|
{ |
|
|
return DateTime.IsLeapYear(year); |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Converts the timeline in hours and minutes and returns as a Vector2(hours, minutes). |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public Vector2 GetTime () |
|
|
{ |
|
|
if (setTimeByCurve) |
|
|
{ |
|
|
m_hourAndMinutes.x = Mathf.Floor(hourByCurve); |
|
|
m_hourAndMinutes.y = 60.0f * (hourByCurve - m_hourAndMinutes.x); |
|
|
m_hourAndMinutes.y = Mathf.Floor(m_hourAndMinutes.y); |
|
|
} |
|
|
else |
|
|
{ |
|
|
m_hourAndMinutes.x = Mathf.Floor(hour); |
|
|
m_hourAndMinutes.y = 60.0f * (hour - m_hourAndMinutes.x); |
|
|
m_hourAndMinutes.y = Mathf.Floor(m_hourAndMinutes.y); |
|
|
} |
|
|
return m_hourAndMinutes; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Get the current date and returns as a Vector3(month, day, year) |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public Vector3 GetDate () |
|
|
{ |
|
|
m_date.x = month; |
|
|
m_date.y = day; |
|
|
m_date.z = year; |
|
|
return m_date; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Get the current system time and returns as a Float. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public float GetSystemTime () |
|
|
{ |
|
|
return DateTime.Now.Hour + ((1.0f / 60.0f) * DateTime.Now.Minute); |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to convert current system time into a float and apply to timeline. |
|
|
/// </summary> |
|
|
public void ApplySystemTime () |
|
|
{ |
|
|
hour = DateTime.Now.Hour + ((1.0f / 60.0f) * DateTime.Now.Minute); |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to apply the current system date. |
|
|
/// </summary> |
|
|
public void ApplySystemDate () |
|
|
{ |
|
|
month = DateTime.Now.Month; |
|
|
day = DateTime.Now.Day; |
|
|
year = DateTime.Now.Year; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Starts the next day and resets the time based on the Repeat Mode set in the Inspector. This method is commonly called at midnight by AzureSkyController script. |
|
|
/// </summary> |
|
|
/// <param name="repeatMode">Place here the Repeat Mode setted on the Options tab.</param> |
|
|
public void StartNextDay (int repeatMode) |
|
|
{ |
|
|
if(repeatMode != 1)day += 1; |
|
|
hour = 0; |
|
|
if(day > m_daysInMonth) |
|
|
{ |
|
|
day = 1; |
|
|
if (repeatMode != 2) |
|
|
{ |
|
|
month += 1; |
|
|
if (month > 12) |
|
|
{ |
|
|
month = 1; |
|
|
if (repeatMode != 3) |
|
|
{ |
|
|
year += 1; |
|
|
if (year > 9999) |
|
|
{ |
|
|
year = 1; |
|
|
} |
|
|
} |
|
|
} |
|
|
} |
|
|
} |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Skips to the next day and ignores the Repeat Mode set in the Inspector. |
|
|
/// </summary> |
|
|
/// <param name="keepTime">Keep the current time?</param> |
|
|
public void JumpToNextDay ( bool keepTime = true) |
|
|
{ |
|
|
day += 1; |
|
|
if (!keepTime) |
|
|
hour = 0; |
|
|
if (day > m_daysInMonth) |
|
|
{ |
|
|
day = 1; |
|
|
month += 1; |
|
|
if (month > 12) |
|
|
{ |
|
|
month = 1; |
|
|
year += 1; |
|
|
if (year > 9999) |
|
|
{ |
|
|
year = 1; |
|
|
} |
|
|
} |
|
|
} |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Go to a custom date. |
|
|
/// </summary> |
|
|
/// <param name="month">The number of the month you want to go.</param> |
|
|
/// <param name="day">The number of the day you want to go.</param> |
|
|
/// <param name="year">The number of the year you want to go.</param> |
|
|
public void GotoDate (int month, int day, int year) |
|
|
{ |
|
|
this.month = month; |
|
|
this.day = day; |
|
|
this.year = year; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Go to a custom time. |
|
|
/// </summary> |
|
|
/// <param name="hours">The number of the hour you want to go.</param> |
|
|
/// <param name="minutes">The number of the minute you want to go.</param> |
|
|
public void GotoTime (int hours, int minutes) |
|
|
{ |
|
|
this.hour = hours + ((1.0f / 60.0f) * minutes); |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to apply the time progression. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public float GetDayLength () |
|
|
{ |
|
|
if (dayCycle > 0.0f) |
|
|
return (24.0f / 60.0f) / dayCycle; |
|
|
else |
|
|
return 0.0f; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to calculate the time of day based on the curve time. |
|
|
/// </summary> |
|
|
public void CalculateTimeByCurve () |
|
|
{ |
|
|
hourByCurve = dayCycleCurve.Evaluate(hour); |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to set the simple sun position. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public float SetSimpleSunPosition () |
|
|
{ |
|
|
if (setTimeByCurve) |
|
|
{ |
|
|
return ((hourByCurve + utc) * 360.0f / 24.0f) - 90.0f; |
|
|
} |
|
|
else |
|
|
{ |
|
|
return ((hour + utc) * 360.0f / 24.0f) - 90.0f; |
|
|
} |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to set the simple moon position. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public Quaternion SetSimpleMoonPosition ( Transform sun) |
|
|
{ |
|
|
return sun.transform.rotation * Quaternion.Euler(0, -180, 0); |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to set the realistic sun position based on Time, Date and Location. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public Vector3 SetRealisticSunPosition () |
|
|
{ |
|
|
m_radians = (Mathf.PI * 2.0f) / 360.0f;//Used to convert degress to radians. |
|
|
m_radLatitude = m_radians * latitude; |
|
|
m_sinLatitude = Mathf.Sin(m_radLatitude); |
|
|
m_cosLatitude = Mathf.Cos(m_radLatitude); |
|
|
|
|
|
float hour = this.hour - utc; |
|
|
//Time Scale. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//d = 367*y - 7 * ( y + (m+9)/12 ) / 4 + 275*m/9 + D - 730530 |
|
|
//d = d + UT/24.0 |
|
|
float d = 367 * year - 7 * (year + (month + 9) / 12) / 4 + 275 * month / 9 + day - 730530; |
|
|
d = d + hour / 24.0f; |
|
|
|
|
|
//Tilt of earth's axis. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//obliquity of the ecliptic. |
|
|
float ecliptic = 23.4393f - 3.563E-7f * d; |
|
|
//Need convert to radians before apply sine and cosine. |
|
|
float radEcliptic = m_radians * ecliptic; |
|
|
float sinEcliptic = Mathf.Sin(radEcliptic); |
|
|
float cosEcliptic = Mathf.Cos(radEcliptic); |
|
|
|
|
|
//Orbital elements of the Sun. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//float N = 0.0; |
|
|
//float i = 0.0; |
|
|
float w = 282.9404f + 4.70935E-5f * d; |
|
|
//float a = 1.000000f; |
|
|
float e = 0.016709f - 1.151E-9f * d; |
|
|
float M = 356.0470f + 0.9856002585f * d; |
|
|
|
|
|
//Eccentric anomaly. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//E = M + e*(180/pi) * sin(M) * ( 1.0 + e * cos(M) ) in degress. |
|
|
//E = M + e * sin(M) * ( 1.0 + e * cos(M) ) in radians. |
|
|
//Need convert to radians before apply sine and cosine. |
|
|
float radM = m_radians * M; |
|
|
float sinM = Mathf.Sin(radM); |
|
|
float cosM = Mathf.Cos(radM); |
|
|
|
|
|
//Need convert to radians before apply sine and cosine. |
|
|
float radE = radM + e * sinM * (1.0f + e * cosM); |
|
|
float sinE = Mathf.Sin(radE); |
|
|
float cosE = Mathf.Cos(radE); |
|
|
|
|
|
//Sun's distance (r) and its true anomaly (v). |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//Xv = r * cos (v) = cos (E) - e |
|
|
//Yv = r * sen (v) = sqrt (1,0 - e * e) * sen (E) |
|
|
float xv = cosE - e; |
|
|
float yv = Mathf.Sqrt(1.0f - e * e) * sinE; |
|
|
|
|
|
//V = atan2 (yv, xv) |
|
|
//R = sqrt (xv * xv + yv * yv) |
|
|
float v = Mathf.Rad2Deg * Mathf.Atan2(yv, xv); |
|
|
float r = Mathf.Sqrt(xv * xv + yv * yv); |
|
|
|
|
|
//Sun's true longitude. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float radLongitude = m_radians * (v + w); |
|
|
float sinLongitude = Mathf.Sin(radLongitude); |
|
|
float cosLongitude = Mathf.Cos(radLongitude); |
|
|
|
|
|
float xs = r * cosLongitude; |
|
|
float ys = r * sinLongitude; |
|
|
|
|
|
//Equatorial coordinates. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float xe = xs; |
|
|
float ye = ys * cosEcliptic; |
|
|
float ze = ys * sinEcliptic; |
|
|
|
|
|
//Sun's Right Ascension(RA) and Declination(Dec). |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float RA = Mathf.Atan2(ye, xe); |
|
|
float Dec = Mathf.Atan2(ze, Mathf.Sqrt(xe * xe + ye * ye)); |
|
|
float sinDec = Mathf.Sin(Dec); |
|
|
float cosDec = Mathf.Cos(Dec); |
|
|
|
|
|
//The Sidereal Time. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float Ls = v + w; |
|
|
|
|
|
float GMST0 = Ls + 180.0f; |
|
|
float UT = 15.0f * hour;//Universal Time. |
|
|
float GMST = GMST0 + UT; |
|
|
float LST = m_radians * (GMST + longitude); |
|
|
|
|
|
//Store local sideral time. |
|
|
lst = LST; |
|
|
|
|
|
//Azimuthal coordinates. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float HA = LST - RA; |
|
|
float sinHA = Mathf.Sin(HA); |
|
|
float cosHA = Mathf.Cos(HA); |
|
|
|
|
|
float x = cosHA * cosDec; |
|
|
float y = sinHA * cosDec; |
|
|
float z = sinDec; |
|
|
|
|
|
float xhor = x * m_sinLatitude - z * m_cosLatitude; |
|
|
float yhor = y; |
|
|
float zhor = x * m_cosLatitude + z * m_sinLatitude; |
|
|
|
|
|
//az = atan2( yhor, xhor ) + 180_degrees |
|
|
//alt = asin( zhor ) = atan2( zhor, sqrt(xhor*xhor+yhor*yhor) ) |
|
|
float azimuth = Mathf.Atan2(yhor, xhor) + m_radians * 180.0f; |
|
|
float altitude = Mathf.Asin(zhor); |
|
|
|
|
|
//Zenith angle. |
|
|
//Zenith=90°−α Where α is the elevation angle. |
|
|
float zenith = 90.0f * m_radians - altitude; |
|
|
|
|
|
//Converts from Spherical(radius r, zenith-inclination θ, azimuth φ) to Cartesian(x,y,z) coordinates. |
|
|
//https://en.wikipedia.org/wiki/Spherical_coordinate_system |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//x = r sin(θ)cos(φ) |
|
|
//y = r sin(θ)sin(φ) |
|
|
//z = r cos(θ) |
|
|
Vector3 ret; |
|
|
|
|
|
//radius = 1 |
|
|
ret.z = Mathf.Sin(zenith) * Mathf.Cos(azimuth); |
|
|
ret.x = Mathf.Sin(zenith) * Mathf.Sin(azimuth); |
|
|
ret.y = Mathf.Cos(zenith); |
|
|
|
|
|
return ret * -1.0f; |
|
|
} |
|
|
|
|
|
/// <summary> |
|
|
/// Used by AzureSkyController script to set the realistic moon position based on Time, Date and Location. |
|
|
/// </summary> |
|
|
/// <returns></returns> |
|
|
public Vector3 SetRealisticMoonPosition () |
|
|
{ |
|
|
float hour = this.hour - utc; |
|
|
|
|
|
//Time Scale. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//d = 367*y - 7 * ( y + (m+9)/12 ) / 4 + 275*m/9 + D - 730530 |
|
|
//d = d + UT/24.0 |
|
|
float d = 367 * year - 7 * (year + (month + 9) / 12) / 4 + 275 * month / 9 + day - 730530; |
|
|
d = d + hour / 24.0f; |
|
|
|
|
|
//Tilt of earth's axis. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//obliquity of the ecliptic. |
|
|
float ecliptic = 23.4393f - 3.563E-7f * d; |
|
|
//Need convert to radians before apply sine and cosine. |
|
|
float radEcliptic = m_radians * ecliptic; |
|
|
float sinEcliptic = Mathf.Sin(radEcliptic); |
|
|
float cosEcliptic = Mathf.Cos(radEcliptic); |
|
|
|
|
|
//Orbital elements of the Moon. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float N = 125.1228f - 0.0529538083f * d; |
|
|
float i = 5.1454f; |
|
|
float w = 318.0634f + 0.1643573223f * d; |
|
|
float a = 60.2666f; |
|
|
float e = 0.054900f; |
|
|
float M = 115.3654f + 13.0649929509f * d; |
|
|
|
|
|
//Eccentric anomaly. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//E = M + e*(180/pi) * sin(M) * ( 1.0 + e * cos(M) ) |
|
|
float radM = m_radians * M; |
|
|
float E = radM + e * Mathf.Sin(radM) * (1f + e * Mathf.Cos(radM)); |
|
|
|
|
|
//Planet's distance and true anomaly. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//xv = r * cos(v) = a * ( cos(E) - e ) |
|
|
//yv = r * sin(v) = a * ( sqrt(1.0 - e*e) * sin(E) ) |
|
|
float xv = a * (Mathf.Cos(E) - e); |
|
|
float yv = a * (Mathf.Sqrt(1f - e * e) * Mathf.Sin(E)); |
|
|
//V = atan2 (yv, xv) |
|
|
//R = sqrt (xv * xv + yv * yv) |
|
|
float v = Mathf.Rad2Deg * Mathf.Atan2(yv, xv); |
|
|
float r = Mathf.Sqrt(xv * xv + yv * yv); |
|
|
|
|
|
//Moon position in 3D space. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float radLongitude = m_radians * (v + w); |
|
|
float sinLongitude = Mathf.Sin(radLongitude); |
|
|
float cosLongitude = Mathf.Cos(radLongitude); |
|
|
|
|
|
//Geocentric (Earth-centered) coordinates. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//xh = r * ( cos(N) * cos(v+w) - sin(N) * sin(v+w) * cos(i) ) |
|
|
//yh = r * ( sin(N) * cos(v+w) + cos(N) * sin(v+w) * cos(i) ) |
|
|
//zh = r * ( sin(v+w) * sin(i) ) |
|
|
float radN = m_radians * N; |
|
|
float radI = m_radians * i; |
|
|
|
|
|
float xh = r * (Mathf.Cos(radN) * cosLongitude - Mathf.Sin(radN) * sinLongitude * Mathf.Cos(radI)); |
|
|
float yh = r * (Mathf.Sin(radN) * cosLongitude + Mathf.Cos(radN) * sinLongitude * Mathf.Cos(radI)); |
|
|
float zh = r * (sinLongitude * Mathf.Sin(radI)); |
|
|
|
|
|
//float xg = xh; //No needed to the moon. |
|
|
//float yg = yh; |
|
|
//float zg = zh; |
|
|
|
|
|
//Equatorial coordinates. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float xe = xh; |
|
|
float ye = yh * cosEcliptic - zh * sinEcliptic; |
|
|
float ze = yh * sinEcliptic + zh * cosEcliptic; |
|
|
|
|
|
//Planet's Right Ascension (RA) and Declination (Dec). |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float RA = Mathf.Atan2(ye, xe); |
|
|
float Dec = Mathf.Atan2(ze, Mathf.Sqrt(xe * xe + ye * ye)); |
|
|
|
|
|
//The Sidereal Time. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//It is already calculated for the sun and stored in the lst, it is not necessary to calculate again for the moon. |
|
|
//float Ls = ls; |
|
|
|
|
|
//float GMST0 = Ls + 180.0f; |
|
|
//float UT = 15.0f * hour; |
|
|
//float GMST = GMST0 + UT; |
|
|
//float LST = radians * (GMST + Azure_Longitude); |
|
|
|
|
|
//Azimuthal coordinates. |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
float HA = lst - RA; |
|
|
|
|
|
float x = Mathf.Cos(HA) * Mathf.Cos(Dec); |
|
|
float y = Mathf.Sin(HA) * Mathf.Cos(Dec); |
|
|
float z = Mathf.Sin(Dec); |
|
|
|
|
|
float xhor = x * m_sinLatitude - z * m_cosLatitude; |
|
|
float yhor = y; |
|
|
float zhor = x * m_cosLatitude + z * m_sinLatitude; |
|
|
|
|
|
//az = atan2( yhor, xhor ) + 180_degrees |
|
|
//alt = asin( zhor ) = atan2( zhor, sqrt(xhor*xhor+yhor*yhor) ) |
|
|
float azimuth = Mathf.Atan2(yhor, xhor) + m_radians * 180.0f; |
|
|
float altitude = Mathf.Asin(zhor); |
|
|
|
|
|
//Zenith angle. |
|
|
//Zenith = 90°−α where α is the elevation angle. |
|
|
float zenith = 90.0f * m_radians - altitude; |
|
|
|
|
|
//Converts from Spherical(radius r, zenith-inclination θ, azimuth φ) to Cartesian(x,y,z) coordinates. |
|
|
//https://en.wikipedia.org/wiki/Spherical_coordinate_system |
|
|
//--------------------------------------------------------------------------------------------------- |
|
|
//x = r sin(θ)cos(φ) |
|
|
//y = r sin(θ)sin(φ) |
|
|
//z = r cos(θ) |
|
|
Vector3 ret; |
|
|
|
|
|
//radius = 1 |
|
|
ret.z = Mathf.Sin(zenith) * Mathf.Cos(azimuth); |
|
|
ret.x = Mathf.Sin(zenith) * Mathf.Sin(azimuth); |
|
|
ret.y = Mathf.Cos(zenith); |
|
|
|
|
|
return ret * -1.0f; |
|
|
} |
|
|
} |
|
|
} |
