JP2022055747A - Appearance time providing method, appearance time providing device, and appearance time providing program - Google Patents

Abstract

To provide a method of providing celestial body appearance time on a skyline.SOLUTION: A method of providing celestial body appearance time on a skyline, comprises, using an appearance time providing device: an input step of receiving an input of a date, and position information on a reference point; an orbit coordinate acquisition step of acquiring an orbit coordinate of a celestial body when viewing from the reference point on the date, using the date and position information; a height angle acquisition step of, as to a plurality of candidate points existing on a prescribed azimuth when viewing from the reference point, acquiring a height angle based on a distance from the reference point and a height of the candidate point; a skyline coordinate acquisition step of acquiring a maximum value of the height angle in the azimuth to thereby acquire a skyline coordinate; and a time determination step of determining appearance time of the celestial body on the skyline on the basis of the skyline coordinate and orbit coordinate.SELECTED DRAWING: Figure 1

Description

The present invention relates to a method of providing a haunting time, a haunting time providing device, and a haunting time providing program for providing a haunting time of a celestial body in a skyline.

Conventionally, the time of appearance of celestial bodies such as the sun and the moon (in the case of the sun, the time of sunrise and / or the time of sunset (sunset time)) has been provided. For example, in Non-Patent Document 1, the horizon is obtained, and the time when the upper edge of the sun appears on the horizon is calculated as the sunrise time, or the time when the sun sets on the horizon is calculated as the sunset time. The horizon here means that the visible range changes depending on the altitude of the observation point, the effect of refraction of sunlight by the atmosphere, and the orbit that the earth revolves around the sun is an elliptical orbit, so the apparent size of the sun changes. It is decided in consideration of things to do. This horizon will be called the sun's horizon. The conventional sunrise and sunset times are calculated as the times when the ecliptic (the apparent orbit of the sun) intersects the horizon of the sun.

Further, as Patent Document 1, a device for displaying the time, direction, and scene of entering and exiting the sun and the full moon is known. In Patent Document 1, the polar coordinate system graph corresponds to the position of a specific observer, and the moon enters and exits limited to a specific moon phase from the sunrise / sunset times and directions provided by the National Astronomical Observatory. It is explained that the time and direction are plotted and created, respectively, thereby providing a device for giving the time and direction of entering and exiting the sun and the moon. In addition to entering and exiting from the horizon and horizon, the time direction of entering and exiting from the background of the object is given by specifying the position and height of a specific object from the map and calculating the elevation angle.

Further, as Patent Document 2, the mountain shape display function displays the state of the mountain seen from an arbitrary position by using the map information or the altitude information in the database, and the mountain name / altitude display function, the wide area display function, and the sky information. There is known a technique for creating a landscape by displaying a display function and a target display function respectively. In Patent Document 2, a landscape display device that creates a landscape in an arbitrary direction at an arbitrary place is used to create a mountain shape by using map data, a position input means for inputting position information, and altitude information in the map data. And display it on the scenery screen. It is also shown here that the constellations are displayed along with the scenery. There is also a description that the orbital trajectory information of the sun is also stored in a database and the orbital trajectory of the sun is displayed together with the scenery.

Ko Nagasawa "Calculation of sunrise and sunset How to find the time of appearance of celestial bodies", Chijin Shokan (1999)

Utility Model Registration No. 315471 Gazette Japanese Unexamined Patent Publication No. 8-201076

In most places in Japan, mountains exist in the direction of appearance of celestial bodies, and the appearance time assuming the horizon of the sun as in Non-Patent Document 1 is usually different from the actual sunrise and sunset times. In order to calculate the sunrise and sunset times more realistically, it is necessary to calculate the appearance time on the contour line (hereinafter referred to as the skyline) drawn by the ridgeline of the mountain with the sky as the background, not the horizon of the sun. be.

Even if the sunrise and sunset times are calculated in consideration of the terrain, as in Patent Document 1, the data provided by the National Astronomical Observatory and the topographic map of the National Land Research Institute are used to determine the appearance time and direction at a specific location in advance. However, it is not possible to provide the sunrise, sunset time and direction immediately upon request at any location and date.

Patent Document 2 provides a device for displaying a landscape using altitude information in map data in order to display a landscape from an arbitrary point. It is also stated here that the vertical cut plane between arbitrary points can be obtained from the altitude information in the map data, but the specific method and means for that purpose are not disclosed. In addition, considering the skyline, constellations, sunrises, and sunsets can be displayed together with the scenery, but it does not disclose how to specifically calculate and provide them at any place and date. Furthermore, it takes time for these calculations, and the method and means for making them practical are not disclosed.

The present invention has been made in view of the above circumstances, and it is an object to be solved to provide a method for providing the appearance time of a celestial body in a skyline.

In order to solve the above problems, the present invention is a method of providing a haunting time that provides a haunting time of a celestial body in a skyline.
An input step in which the haunting time providing device accepts input of the date and the position information of the reference point, and
An orbital coordinate acquisition step for acquiring the orbital coordinates of the celestial body viewed from the reference point on the date using the date and position information.
A height angle acquisition step for acquiring a height angle based on the distance from the reference point and the height of the candidate point for a plurality of candidate points existing in a predetermined direction as viewed from the reference point.
The skyline coordinate acquisition step for acquiring the skyline coordinates by acquiring the maximum value of the height angle in the direction, and
It has a time determination step for determining the appearance time of a celestial body in the skyline based on the skyline coordinates and the orbital coordinates.

Further, the present invention is an infestation time providing device that provides an infestation time of a celestial body in a skyline.
A data input unit that accepts input of date and location information of the reference point,
An orbital coordinate acquisition unit that acquires the orbital coordinates of the celestial body viewed from the reference point on the date using the date and position information.
A height angle acquisition unit that acquires a height angle based on the distance from the reference point and the height of the candidate point for a plurality of candidate points existing in a predetermined direction as viewed from the reference point.
A skyline coordinate acquisition unit that acquires the skyline coordinates by acquiring the maximum value of the height angle in the direction, and
It has a time determination unit that determines the appearance time of a celestial body in the skyline based on the skyline coordinates and the orbital coordinates.

Further, the present invention is a haunting time providing program that provides haunting times of celestial bodies in the skyline.
A computer, a data input unit that accepts input of date and location information of a reference point,
An orbital coordinate acquisition unit that acquires the orbital coordinates of the celestial body viewed from the reference point on the date using the date and position information.
A height angle acquisition unit that acquires a height angle based on the distance from the reference point and the height of the candidate point for a plurality of candidate points existing in a predetermined direction as viewed from the reference point.
A skyline coordinate acquisition unit that acquires the skyline coordinates by acquiring the maximum value of the height angle in the direction, and
It functions as a time determination unit that determines the appearance time of a celestial body in the skyline based on the skyline coordinates and the orbital coordinates.

With such a configuration, the skyline can be calculated and the time of appearance of the celestial body in the skyline can be provided.

In a preferred embodiment of the present invention, the height angle acquisition step includes a step of acquiring the height angle for a plurality of candidate points existing at first intervals.
Around the candidate point where the height angle is maximum, the step of acquiring the height angle for a plurality of candidate points existing at a second interval is further provided.
The second interval is narrower than the first interval.

In a preferred embodiment of the present invention, the step of acquiring the height angle for the plurality of candidate points existing at the first interval is added to the plurality of candidate points existing at the first interval. The height angles are obtained for a plurality of candidate points existing at a third interval, and the height angles are obtained.
The third interval is narrower than the first interval, and all of the plurality of sets of candidate points existing across the third interval are from the plurality of sets of candidate points existing across the first interval. It exists near the reference point and
The second interval is narrower than the first interval and the third interval.

With such a configuration, the number of calculation of skyline coordinates can be reduced.

In a preferred embodiment of the present invention, the orbital coordinate acquisition step includes a step of obtaining a plurality of steps at predetermined time intervals with the direction and height angle at which the celestial body is located as calculation points.
It has a step of connecting the plurality of calculation points with an approximation line and calculating an approximation function for obtaining an arbitrary direction, height angle, and time.

With such a configuration, the number of orbital coordinate calculations can be reduced.

In a preferred embodiment of the present invention, the height angle acquisition step acquires height angles for a plurality of candidate points on the direction obtained in the orbit coordinate acquisition step.
The skyline coordinate acquisition step acquires the skyline coordinates in the direction, and further, the orbit coordinate acquisition step acquires the orbit coordinates from the height angle of the skyline coordinates acquired in the skyline coordinate acquisition step, and a new direction. By performing convergence calculation until the convergence condition is satisfied, the intersection of the skyline and the orbit is obtained.
The time determination step determines the time of appearance of the celestial body at the intersection.
As an initial value of the convergence calculation, the orbital coordinate acquisition step acquires the direction in which the celestial body appears and disappears on the horizon.

With such a configuration, the number of calculations of skyline coordinates and orbital coordinates can be reduced.

In a preferred embodiment of the present invention, a skyline data generation step for generating skyline data by connecting the skyline coordinates acquired in the skyline coordinate acquisition step and a skyline data generation step.
A skyline data storage step of associating the position information of the reference point and the skyline data when acquiring the skyline coordinates and storing the skyline data in the skyline database, and
It has a utilization step of searching skyline data stored in the skyline database based on the position information of the reference point received in the input step.

With such a configuration, it is possible to reduce the calculation of skyline coordinates.

According to the present invention, it is possible to provide a method of providing the appearance time of a celestial body in a skyline.

It is a schematic diagram of the method of providing the appearance time according to Embodiment 1 of this invention. It is a block diagram of the haunting time providing system which concerns on Embodiment 1 of this invention. It is a processing flowchart for providing the appearance time which concerns on Embodiment 1 of this invention. It is a screen display example of the operation screen which concerns on Embodiment 1 of this invention. It is a figure which shows the calculation method of the orbit with respect to Embodiment 1 of this invention. It is a schematic diagram of the height angle acquisition step which concerns on Embodiment 1 of this invention. It is a figure which shows the farthest distance of the elevation search which concerns on Embodiment 1 of this invention. It is a conceptual diagram which shows the search method of the candidate point which concerns on Embodiment 1 of this invention. It is a conceptual diagram which shows the method of the convergence calculation which concerns on Embodiment 2 of this invention.

<Embodiment 1>
Hereinafter, the haunting time providing system according to the first embodiment of the present invention will be described with reference to the drawings. The embodiments shown below are examples of the present invention, and the present invention is not limited to the following embodiments, and various configurations can be adopted. The first embodiment describes a method of obtaining the skyline data described later and using the skyline data to obtain the appearance time of the celestial body.

For example, in the present embodiment, the configuration, operation, and the like of the haunting time providing system will be described, but a method, an apparatus, a computer program, a recording medium, and the like having the same configuration can also exert the same effect. Further, the program may be stored in a recording medium. Using this recording medium, for example, the program can be installed in a computer. Here, the recording medium in which the program is stored may be a non-transient recording medium such as a CD-ROM.

FIG. 1 is a schematic diagram of a method of providing haunting time according to an embodiment of the present invention. The code MO indicates a mountain. The contour line (skyline SL) drawn by the ridgeline of the mountain with the sky as the background is determined by the position of the user and the positional relationship of one or more mountain MOs. The code SN is the sun, and the code SO indicates the orbit (ecliptic) of the sun SN. Here, the skyline SL in the direction of sunrise is referred to as the skyline SR in the direction of sunrise, and the skyline SL in the direction of sunset is referred to as the skyline SS in the sunset direction.

Generally, when obtaining the sunrise (sunset) time, the time TR (TS) of the sunrise (sunset) on the horizon HO of the sun is obtained. The user operates the user terminal device CP to obtain the time STR (STS) of the sunrise (sunset) in the sunrise skyline SR (sunset skyline SS). Here, when the sun SN apparently moves along the orbit SO, the direction of sunrise (sunset) on the horizon HO is the direction DR (direction DS), and the direction of sunrise (sunset) on the sunrise skyline SR (sunset skyline SS) is the direction. Let it be SDR (direction SDS). Further, the height angle of the sunrise (sunset) on the horizon HO is defined as the height angle HR (height angle HS), and the height angle of the sunrise (sunset) on the skyline SL is defined as the height angle SHR (height angle SHS).

In this embodiment, the method of providing the sunrise time STR and the sunset time STS will be described. Similarly, the appearance time of the moon and other celestial bodies can be obtained by obtaining the apparent orbit of the moon and other celestial bodies. Can be done. Further, in the present embodiment, the case where the skyline SL is determined by the mountain MO will be described, but the skyline SL is not limited to the terrain, and may be determined in consideration of, for example, a building. When considering a building, the latitude, longitude and height of the building may be registered in the map information database. In the present invention, the "appearance" may refer to only one of the apparent appearance and the apparent disappearance of the celestial body. For example, the infestation time providing system has a configuration in which only the sunrise time is obtained. May be good.

The haunting time providing system 1 includes a haunting time providing device 100 and a user terminal device CP operated by the user. The haunting time providing device 100 and the user terminal device CP are configured to be communicable by the network NW. The wireless base station BS is a device that is connected to the network NW, establishes a wireless communication connection with the user terminal device CP, and provides data communication or the like by the user terminal device CP. In this embodiment, an embodiment of a client-server model in which the haunting time providing device 100 for executing the haunting time providing program is a server and the user terminal device CP is a client will be described, but the functions described later of the haunting time providing device 100 will be described. A part or all of the components may be configured to be included in the user terminal device CP.

FIG. 2 is a block diagram showing functional components of the haunting time providing system according to the present embodiment. The haunting time providing device 100 includes a control unit 101, a storage unit 102, and a communication unit 103.

The control unit 101 includes a processor such as a CPU (Central Processing Unit), and controls the entire haunting time providing device 100 by executing an OS, a web browser application, other applications, and the like. The control unit 101 can use an arbitrary processor such as a semiconductor device equipped with a CPU, a hardware circuit using a DSP (Digital Signal Processor), an ASIC, an FPGA, or the like.

The storage unit 102 includes a storage device such as a RAM or a ROM, and the RAM is used as a working memory or the like of the CPU (control unit 101) and temporarily stores the calculation result by the CPU (control unit 101). A program for the CPU (control unit 101) is recorded in the ROM, and program execution instructions are sequentially output in response to a request from the CPU (control unit 101). In the present embodiment, the storage unit 102 stores the operating system (OS), the appearance time providing application, and the like. The communication unit 103 connects the user terminal device CP to the network NW and controls communication with the user terminal device CP and the like.

The control unit 101 includes an application execution unit 111. By executing the haunting time providing program of the present invention, the application executing unit 111 virtually constructs each functional module of the haunting time providing system 1 on the control unit 101.

The user terminal device CP is a terminal device operated by the user, and is a smartphone in the illustrated example, but for example, a PC, a tablet terminal, or the like can also be used. The user terminal device CP is an application for using the haunting time to request the haunting time providing device 100 to provide the haunting time, display an operation screen, or the like, or a web browser for realizing these by internet communication or the like. Remembers applications etc.

The user terminal device CP includes a control unit 901, a storage unit 902, a communication unit 903, an input unit 904, an output unit 905, an internal clock 906, and a GNSS (Global Navigation Satellite System) position information interface 907.

The control unit 901 includes a processor such as a CPU (Central Processing Unit), and controls the entire user terminal device CP by executing an OS, a web browser application, other applications, and the like. The control unit 901 can use any processor such as a semiconductor device equipped with a CPU, a hardware circuit using a DSP (Digital Signal Processor), an ASIC, an FPGA, or the like.

The storage unit 902 includes a storage device such as a RAM or a ROM, and the RAM is used as a working memory or the like of the CPU (control unit 901) and temporarily stores the calculation result by the CPU (control unit 901). A program for the CPU (control unit 901) is recorded in the ROM, and program execution instructions are sequentially output in response to a request from the CPU (control unit 901). In the present embodiment, the storage unit 902 stores the operating system (OS), the haunting time utilization application, the web browser application, and the like. The communication unit 903 connects the user terminal device CP to the network NW and controls communication with the haunting time providing device 100 and the like. There may be various means adopted as a network NW, but it is common to realize communication by mobile phone communication or Internet connection from a wireless LAN.

The input unit 904 receives an operation signal from an operation device such as a touch panel, a button, a mouse, and a keyboard, and transmits the received operation signal to the control unit 901, whereby the user can perform an operation on the OS or an application. The output unit 905 transmits a video signal or an audio signal in order to output a video or audio from an output device such as a display or a speaker. An operation device such as a touch panel is connected to the input unit 904 and the output unit 905, and the display data of the operation screen is displayed on the output device such as a display, and the operation signal is output by performing a touch operation or the like accordingly. The GUI to acquire is provided. The internal clock 906 provides date and time data to the control unit 901. The GNSS position information interface 907 receives the GNSS radio wave, acquires the position information at the place, and gives it to the control unit 901.

As shown in FIG. 2, in the haunting time providing device 100 in which the haunting time providing program is executed, the application execution unit 111 of the control unit 101 includes a data input unit 2, an orbital coordinate acquisition unit 3, a height angle acquisition unit 4, and a skyline. It has a coordinate acquisition unit 5, a skyline data generation unit 6, a skyline data storage unit 7, a sky data utilization unit 8, and a time determination unit 9. Further, the haunting time providing device 100 has a calculation parameter file 21, a map database 22, and a skyline database 23 in the storage unit 102.

The storage unit 102 records predetermined parameters used when calculating the orbits and skylines of celestial bodies as the calculation parameter file 21. The map database 22 stores topographic map information associated with the latitude, longitude, and altitude of each point. Further, the skyline database 23 stores the skyline data described later. After that, the intervention of these data will be omitted as long as it is not unknown.

In the present embodiment, the topographic map information is stored in the storage unit 102 as the map database 22, but it may be acquired from an external service by, for example, API.

Next, a method of providing the haunting time will be described with reference to FIGS. 3 to 8. FIG. 3 is a processing flowchart for providing the appearance time.

<Input of reference point>
FIG. 4 is a screen display example of the operation screen. The user terminal device CP displays the operation screen W90 on the output device via the output unit 905. The operation screen W90 has a date input portion W91, today's date button W92, a position input portion W93, a current position acquisition button W94, a sunrise / entry button W95, a sunrise / entry display portion W96, a sunrise terrain image display portion W97, and a sunset. A terrain image display portion W98 is provided. The user operates the input device and the input unit 904 to write the date data for which the appearance time is to be acquired in the date input portion W91. When today's date button W92 is selected and operated, today's date is acquired from the internal clock 906 and written in the date input portion W91.

Further, the user operates the input device and the input unit 904 to write the latitude (north latitude), longitude (east longitude), and altitude as the position information of the reference point in the position input portion W93. When the current position acquisition button W94 is selected and operated, the latitude / longitude / elevation data of the current location acquired from the GNSS position information interface 907 is input to the position input portion W93 as position information. After inputting the date and location information, the user selects and operates the sunrise / entry button W95. When the sunrise / entry button W95 is selected and operated, the date and position information are transmitted to the haunting time providing device 100. The data input unit 2 of the haunting time providing device 100 receives the transmitted date and position information in step S31.

Further, when the haunting time providing device 100 returns the haunting calculation result data calculated by the method described below to the user terminal device CP, the user terminal device CP displays the sunrise / sunset time on the sunrise / sunset display portion W96. It also displays the direction (horizontal angle measured around the east with 0 ° to the north). Further, the user terminal device CP displays the orbit SO, the sunrise skyline SR, and the sunset skyline SS on the sunrise terrain image display portion W97 and the sunset terrain image display portion W98, and presents the directions of sunrise and sunset by images. ..

<Determining the orbital coordinates and the angle of appearance and height on the horizon>
When the data input unit 2 receives the position information of the date and the reference point, the orbital coordinate acquisition unit 3 uses the date and the position information to set a plurality of apparent orbital coordinates of the celestial body seen from this reference point on this date. At the same time as acquiring (step S32), the direction DR (direction DS) and the height angle HR (height angle HS) of the sunrise (sunset) at the horizon HO are calculated (step S33). The orbital coordinates P have an apparent direction X of the celestial body as seen from the reference point, a height angle Y of the celestial body, and a time T.

In the present embodiment, since sunrise and sunset are targeted, the orbital coordinate acquisition unit 3 calculates the approximately assumed time range at predetermined time intervals. For the sunrise, the time T ₀ of the sunrise on the horizon HO is first calculated as the orbital coordinate P ₀ of the sunrise on the horizon HO, and the direction X ₀ is also calculated (conventional technique). Assuming that the calculation accuracy (calculation interval) is 1 minute, the direction X _i (i = 0, 1, 2, ... N) and the height angle Y _i (indicating the altitude of the sun) indicate the altitude of the sun every minute from P ₀ . i = 0, 1, 2, ... N) is calculated sequentially. The number of times N as the calculation range is, for example, N = 2 × 60 = 120 when the time accuracy is 1 minute until 2 hours after sunrise from the horizon HO. Therefore, the orbit (ecliptic) is created on the storage unit 102 (RAM) as the sequence data (i = 0, 1, 2, ... N) of Pi (Ti _, X _{i ,} Y _i ). The sunset is the same as the above sunrise, but it is calculated sequentially in the direction opposite to the sunrise (sunset direction: from west to south) with respect to the direction of sunset on the horizon. Since such a calculation range changes depending on the request, it is recorded in the storage unit 102 as a calculation parameter file, and the orbital coordinate acquisition unit 3 takes in appropriate calculation parameters based on the date and time and / or the position information of the reference point. ..

Depending on the area, sunrise and sunset may occur at short intervals or multiple times. Therefore, the calculation range of sunrise and / or sunset at a certain point in a certain season does not necessarily have to be one, and the calculation range may be used for both calculation of sunrise and sunset.

FIG. 5 is a diagram showing a method of calculating the orbit. The orbital coordinate acquisition unit 3 calculates the orbital coordinates at predetermined time intervals for an approximately assumed time range. In the present embodiment, the orbital coordinate acquisition unit 3 obtains a plurality of calculation points at predetermined time intervals with the direction and height angle at which the celestial body is located as calculation points, and connects the plurality of calculation points with an approximate line to obtain any arbitrary value. Calculate an approximate function to find the direction, height angle and time. Here, as a predetermined time interval, for example, every 30 minutes, the orbital coordinates Pi (T _i , X _i , Y _i ) ( _i = 1, 2, ..., N) are calculated, and continuous coordinate points are obtained. Functions connected by straight lines are approximated by linear functions using equations (1) to (3).

According to the equation (1), the direction X is given and the height angle Y at that time is given, and according to the equation (2), the height angle Y is given and the direction X at that time is given. The time T at that time can be obtained by giving the direction X. For example, when it is desired to know the height angle Y in the direction X (X ₄ <X <X ₅ ) as shown in the illustrated example, the equation (1) is used. By connecting the orbital coordinates calculated in this way with an approximate line and calculating the direction X, the height angle Y, and the time T with an approximate function, the calculation of the orbital coordinates with a high calculation load is simplified and each is obtained immediately. be able to.

It is not always necessary to use the approximation function, and the orbital coordinate acquisition unit 3 may calculate the orbital coordinates at time intervals such as 1 minute. Further, for example, the orbital coordinates are calculated at time intervals such as 30 minutes, the orbital coordinates between them are calculated at time intervals such as 1 minute, and the array data of _Pi ( _Ti, Xi _{, Yi} ) ( It may be stored in the storage unit 102 (RAM) as i = 0, 1, 2, ... N).

<Skyline data creation, sunrise direction>
In step S34, the skyline SR in the direction of sunrise is calculated. FIG. 6 is a schematic diagram of the height angle acquisition step in the direction of sunrise. The height angle acquisition unit 4 acquires height angles based on the distance from the reference point and the height of the candidate points for a plurality of candidate points existing in a predetermined direction as seen from the reference point.

First, the height angle acquisition unit 4 refers to the map database 22 for the sunrise direction DR calculated by the orbital coordinate acquisition unit 3, and searches for a candidate point on the direction d ₀ (DR). The details of the farthest distance of the search and the method of determining the candidate point will be described later. The height angle acquisition unit 4 searches for a plurality of candidate points far from the reference point (elevation _h0 of the reference point is given as position information), and each candidate point searched from the map database 22. The sinking distance ΔR due to the roundness of the earth (ΔR = 0.5L ² / R with an approximation in which L is sufficiently smaller than the earth radius R) is calculated from the distance L between the reference point and the candidate point with respect to the altitude h of. Let the elevation correction data be h-h ₀ -ΔR, and calculate the height angle φ of the candidate point seen from the reference point by the formula (4) based on the distance L to that point.

<Determining the farthest distance>
The altitude search interval of the candidate point is calculated, for example, every 10 m, and the distance is calculated up to the farthest distance shown below. If the farthest distance is determined uniformly, it is useless when it makes no sense to search for altitude, for example, following the plain to the top of the sea, or when there is clearly no maximum viewing angle at a distance. Elevation search and calculation will be performed. Therefore, as shown in FIG. 7, the farthest distance in the area divided by the solid line 31 is 100 km on the east side and 150 km on the west side (within the range of the alternate long and short dash line 32). A distance that does not have a maximum looking-up angle even if the altitude of the above is searched is defined in advance and stored as a calculation parameter file 21. In the example of FIG. 7, the search distance in the east-west direction of the elevation search at the points within the rectangular range formed by the two points of the northwest point and the southeast point is defined as follows. By dividing Japan into such a rectangular range and defining the farthest distance, it is possible to shorten the time required for elevation search.

Central 1: Northwest point latitude = 37 ° Northwest point longitude = 137 ° Southeast point latitude = 34.5 ° Southeast point longitude = 138 °
The farthest distance on the east side = 100km The farthest distance on the west side = 150km

The skyline coordinate acquisition unit 5 acquires the maximum value of the height angle in a certain direction and acquires the skyline coordinates. The skyline coordinate acquisition unit 5 determines a candidate point for which the height angle has the maximum value among the plurality of candidate points calculated by the height angle acquisition unit 4, and sets this maximum angle as one point on the skyline SR in this direction. The skyline coordinates are.

In FIG. 6A, if the elevation of each candidate point is searched with the reference point O as the starting point, the cross-sectional shape of the terrain in the direction d can be obtained as the elevation search result. In the example of FIG. 6A, it is assumed that the points A and B are selected as candidate points in the direction d ₀ (DR). At this time, since the height angle φ _0A looking up at the point A from the reference point O is larger than the height angle φ _{B 0} looking up at the point B, the maximum height angle φ 0 _MAX in the direction d ₀ (DR) is obtained. The skyline coordinate acquisition unit 5 determines the skyline coordinates c ₀ having a direction = d ₀ and a height angle φ 0 _MAX (= height angle φ ₀ A). The skyline is created on the storage unit 102 (RAM) as sequence data (i = 0, 1, 2, ... N) of C _i ( _Ti, X _i, Z _i ).

The height angle acquisition unit 4 and the skyline coordinate acquisition unit 5 perform the above calculation for each predetermined angle from the sunrise direction DR (direction d ₀ ) to the south, and determine the skyline coordinates c _{0 to N.} In the illustrated example, the height angle acquisition unit 4 and the skyline coordinate acquisition unit 5 perform the above calculation for each of the directions d ₀ (DR), d _1, d ₂ , ..., D _N , and the skyline coordinates c. Calculate _{0 to N.} FIG. 6B is a diagram showing the sunrise skyline SR generated by arranging the calculated skyline coordinates c _{0 to N.} The maximum value of this predetermined direction (direction d _N of the southernmost point of the skyline coordinates) is, for example, the direction X _N of the orbital coordinates _PN .

If the directions of the skyline coordinates c _{0 to N} match the direction Xi of the orbital coordinates, the calculation of the haunting time, which will be described later, will be simplified.

<Skyline data creation, sunset direction>
Further, in step S35, the skyline SS in the sunset direction is similarly calculated. The height angle acquisition unit 4 refers to the map database 22 for the sunset direction DS calculated by the orbital coordinate acquisition unit 3, and searches for a candidate point on the direction d ₀ (DS). Further, the height angle acquisition unit 4 and the skyline coordinate acquisition unit 5 perform the above calculation for each predetermined angle from the sunset direction DS (direction d ₀ ) to the south, and determine the skyline coordinates c _{0 to N.}

<Reuse of skyline data>
The skyline data generation unit 6 generates skyline data by connecting the skyline coordinates acquired by the skyline coordinate acquisition unit 5. In the present embodiment, the skyline data storage unit 7 associates the position information of the reference point when acquiring the skyline coordinates, the skyline data generated by the skyline data generation unit 6, and the generation date and time of the skyline data in the skyline database 23. Store.

The sky data utilization unit 8 searches for and uses the skyline data stored in the skyline database 23 based on the position information of the reference point received by the data input unit 2. In the present embodiment, in step S34 (step S35), the skyline data in the skyline database 23 is searched based on the position information received in step S31. Then, when it is determined that the reference point recorded in the skyline database 23 is close to the position information received in step S31, the skyline coordinates are calculated by the height angle acquisition unit 4 and the skyline coordinate acquisition unit 5. Instead, get the skyline data recorded in the skyline database 23 and proceed to the next step.

<Calculation of haunting time>
In step S36, the time determination unit 9 calculates the appearance time of the celestial body in the skyline based on the skyline coordinates and the orbital coordinates. The orbital coordinate acquisition unit 3 extracts the direction _Xi from the array of orbital coordinates ( _Ti , Xi, _Yi ) calculated in step S32 and the array of _skylines ( _Ti , Xi, _{Zi )} , and the direction is that direction. Compare the maximum height angle Z _i calculated for each X _i with Y _i . The height angle Z _i and the height angle _Y i do not always match, but in the case of sunset, if Z i < ₁ > = Y _{i +} 1 immediately after Z _i <Y _i , then Z in the case of sunrise. Immediately after _i > Y _i , Z _{i + 1} = <Y _{i + 1} ), and that is regarded as a coincidence point.

Alternatively, an orbital function defined by the orbital coordinate acquisition unit 3 complementing the orbital coordinates calculated in step S32, and a skyline function complementing the skyline coordinates calculated by the skyline coordinate acquisition unit 5 in step S34 (step S35). The sunrise (sunset) time may be calculated by finding the intersection of. The interpolation method is arbitrary, and the calculation is performed by, for example, linear interpolation.

The time determination unit 9 sets this coincidence point (or intersection) as the sunrise direction SDR (sunset direction SDS) and the sunrise height angle SHR (sunset height angle SHS) at the intersection of the sunrise skyline SR (sunset skyline SS) and the orbit SO. ) Is calculated.

In step S37, the time determination unit 9 calculates a time that coincides with the sunrise height angle SHR (sunset height angle SHS), that is, the skyline sunrise time STR (sunset sunset time STS). For example, the time STR (time STS) is calculated using the above-mentioned equation (3).

In step S38, the sunrise time STR, the direction SDR, the sunset time STS, and the direction SDS obtained by the time determination unit 9 are transmitted (returned) to the user terminal device CP as haunting calculation result data. The user terminal device CP displays and processes the received haunting calculation result data, and displays the sunrise / sunrise display portion W96, the sunrise terrain image display portion W97, and the sunset terrain image display portion W98 as illustrated in FIG. To.

<How to determine candidate points>
Here, in step S34 (step S35), the search method of the candidate point performed by the height angle acquisition unit 4 will be described in detail. FIG. 8 is a conceptual diagram showing a search method for candidate points.

In the calculation method shown in FIG. 8A, the height angle acquisition unit 4 first searches for a wide search interval in the first interval D1, searches for the first candidate point, and then searches for the first candidate point, and then the height angle is obtained from the search interval. Find the candidate point A, which is the maximum point of. Next, with respect to the interval before and after the candidate point A, the search interval is narrowed to the second interval D2 (first interval> second interval), and the altitude search is repeated again. In the illustrated example, new candidate points a1, a2, a3, a4 ... Are obtained by searching at the second interval D2 between the candidate points A1 and A2 as the interval before and after the candidate point A. .. Then, the height angle acquisition unit 4 sets the candidate point having the maximum angle as the height angle of the skyline coordinates.

In the present embodiment, the height angle acquisition unit 4 determines that the candidate point a4 that has become the maximum angle in the second search is determined to be the maximum angle, but for example, the previous search (here, the first interval of the first time). The candidate points obtained in the search in D1) may be added to the comparison target to determine the maximum angle. The repeated search that narrows the search interval is not limited to twice, but may be performed three or four times. Further, in the present embodiment, when determining the distance before and after the candidate point A, the distance to the candidate points before and after the candidate point A is set (twice the first interval), but this distance may be any distance.

As for the candidate point near the reference point O, as shown in equation (1), even if the altitude h is low, it is inevitable that the distance L is small and the height angle φ is large. Therefore, in the calculation method shown in FIG. 8B, the height angle acquisition unit 4 searches at the third interval D3 by narrowing the search interval in a range close to the reference point O (for example, L <10 km). , The search interval in a distant range (10 km ≦ L) is widened (first interval), and the altitude is searched (first interval> third interval). As a result, by reducing the number of elevation searches and finding the maximum point of the looking-up angle, it is possible to avoid overlooking the point that becomes the skyline in the vicinity. It should be noted that this search method may be combined at the time of the first search of the search method shown in FIG. 8A, or when the search distance after the second search is close to the reference point O or includes the reference point O. good. By doing so, the number of searches for elevation points can be significantly reduced, which has a great effect on shortening the calculation time.

<Embodiment 2>
Next, the haunting time providing system according to the second embodiment of the present invention will be described with reference to FIG. The same components as those in the first embodiment are designated by the same reference numerals and the description thereof will be omitted.

In the second embodiment, instead of step S34 to step S36 in FIG. 3, the sunrise direction SDR (sunset direction SDS) and the sunrise height angle SHR (sunset height angle SHS) are calculated by the convergence calculation without obtaining the skyline data. do.

As a convergence calculation, the height angle acquisition unit 4 acquires height angles for a plurality of candidate points on the direction obtained by the orbital coordinate acquisition unit 3, and the skyline coordinate acquisition unit 5 acquires the skyline coordinates in the direction. Further, the orbital coordinate acquisition unit 3 acquires the orbital coordinates from the height angle of the skyline coordinates acquired by the skyline coordinate acquisition unit 5, and performs a convergence calculation until the convergence condition is satisfied. Then, the intersection of the skyline and the orbit is obtained, the time determination unit 9 determines the appearance time of the celestial body at this intersection, and the orbital coordinate acquisition unit 3 sets the direction in which the celestial body appears on the horizon as the initial value of the convergence calculation. ..

A more specific description will be given with reference to FIG. FIG. 9 is a conceptual diagram showing a method of convergence calculation in the direction of sunrise. First, the orbital coordinate acquisition unit 3 obtains the sunrise point _p0 on the horizon HO. Then, the height angle acquisition unit 4 calculates the height angle for the direction X ₀ (direction DR) of the point p ₀ , and the skyline coordinate acquisition unit 5 calculates the point p ₁ (skyline) on the skyline SR at the direction X 0. Coordinates). Next, the point p ₂ on the orbit at which the height angle of p ₁ is φ1 is calculated using the equation (1).

Next, the height angle acquisition unit 4 calculates the height angle for the direction X ₂ of the point p ₂ , and the skyline coordinate acquisition unit 5 obtains the point p ₃ (skyline coordinates) on the skyline SR in the direction X ₂ . ..

Similarly, the direction X ₄ of the point p ₄ on the orbit having the same height angle as this p ₃ is calculated by the equation (1). By performing such a convergence calculation until the determination condition is reached, it is possible to converge to the intersection of the skyline and the ecliptic. The determination of whether or not the product has converged is arbitrary, but it is determined that the product has converged, for example, when the amount of change in the direction X and / or the height angle is below the threshold value. For example, when the difference between the points _pi and _{pi + 1} in the direction X is within 0.1 °, it is determined to be convergent.

When forming a skyline with buildings, etc., it is understood that the sun cannot be seen from the reference point (the reference point is not exposed to direct sunlight) unless the sunrise time STR of a certain day is obtained in areas such as Japan. .. Further, when a building or the like exists at a position that divides the orbit of the sun, a plurality of sets of sunrise time STR and sunset time STS may be obtained. Based on this, for example, it is possible to perform a simulation regarding sunshine such as the presence or absence of direct sunlight and the sunshine duration.

1: Appearance time providing system 2: Data input unit 3: Orbit coordinate acquisition unit 4: Height angle acquisition unit 5: Skyline coordinate acquisition unit 6: Skyline data generation unit 7: Skyline data storage unit 8: Sky data utilization unit 9: Time determination unit 21: Calculation parameter file 22: Map database 23: Skyline database 100: Appearance time providing device 101: Control unit 102: Storage unit 103: Communication unit 111: Application execution unit 901: Control unit 902: Storage unit 903: Communication Section 904: Input section 905: Output section 906: Internal clock 907: GNSS position information interface BS: Radio base station CP: User terminal device HO: Horizon NW: Network SL, SR, SS: Skyline SN: Sun SO: Orbit W90 : Operation screen W91: Date input part W92: Date button W93: Position input part W94: Current position acquisition button W95: Sunrise / entry button W96: Sunrise / entry display portion W97: Sunrise terrain image display portion W98: Sunset terrain image display part

Claims (8)

It is a method of providing the infestation time that provides the infestation time of the celestial body in the skyline.
An input step in which the haunting time providing device accepts input of the date and the position information of the reference point, and
An orbital coordinate acquisition step for acquiring the orbital coordinates of the celestial body viewed from the reference point on the date using the date and position information.
A height angle acquisition step for acquiring a height angle based on the distance from the reference point and the height of the candidate point for a plurality of candidate points existing in a predetermined direction as viewed from the reference point.
The skyline coordinate acquisition step for acquiring the skyline coordinates by acquiring the maximum value of the height angle in the direction, and
A method for providing an infestation time, comprising a time determination step for determining an infestation time of a celestial body in the skyline based on the skyline coordinates and the orbital coordinates. The height angle acquisition step includes a step of acquiring the height angle for a plurality of candidate points existing at the first interval and a step of acquiring the height angle.
Around the candidate point where the height angle is maximum, the step of acquiring the height angle for a plurality of candidate points existing at a second interval is further provided.
The method for providing a haunting time according to claim 1, wherein the second interval is an interval narrower than the first interval. The step of acquiring the height angle for the plurality of candidate points existing at the first interval is separated by a third interval in addition to the plurality of candidate points existing at the first interval. Obtain the height angle for a plurality of existing candidate points, and obtain the height angle.
The third interval is narrower than the first interval, and all of the plurality of sets of candidate points existing across the third interval are from the plurality of sets of candidate points existing across the first interval. It exists near the reference point and
The method for providing haunting time according to claim 2, wherein the second interval is an interval narrower than the first interval and the third interval. The orbital coordinate acquisition step includes a step of obtaining a plurality of orbital coordinates at predetermined time intervals using the direction and height angle at which the celestial body is located as calculation points.
The method for providing haunting time according to any one of claims 1 to 3, further comprising a step of connecting the plurality of calculation points with an approximate line to calculate an approximate function for obtaining an arbitrary direction, height angle, and time. The process of acquiring the height angle for a plurality of candidate points on the direction obtained by the height angle acquisition step in the orbital coordinate acquisition step, and
The skyline coordinate acquisition step acquires the skyline coordinates in the direction, and further, the orbit coordinate acquisition step acquires the orbit coordinates from the height angle of the skyline coordinates acquired in the skyline coordinate acquisition step, and a new direction. By performing convergence calculation until the convergence condition is satisfied, the intersection of the skyline and the orbit is obtained.
The time determination step determines the time of appearance of the celestial body at the intersection.
The haunting time providing method according to any one of claims 1 to 4, wherein as the initial value of the convergence calculation, the orbital coordinate acquisition step is to acquire the direction in which the celestial body appears on the horizon. A skyline data generation step that generates skyline data by connecting the skyline coordinates acquired in the skyline coordinate acquisition step,
A skyline data storage step of associating the position information of the reference point and the skyline data when acquiring the skyline coordinates and storing the skyline data in the skyline database, and
The haunting time providing method according to any one of claims 1 to 4, further comprising a usage step of searching skyline data stored in the skyline database based on the position information of the reference point received in the input step. It is a haunting time providing device that provides the haunting time of celestial bodies in the skyline.
A data input unit that accepts input of date and location information of the reference point,
An orbital coordinate acquisition unit that acquires the orbital coordinates of the celestial body viewed from the reference point on the date using the date and position information.
A height angle acquisition unit that acquires a height angle based on the distance from the reference point and the height of the candidate point for a plurality of candidate points existing in a predetermined direction as viewed from the reference point.
A skyline coordinate acquisition unit that acquires the skyline coordinates by acquiring the maximum value of the height angle in the direction, and
A haunting time providing device having a time determining unit for determining the haunting time of a celestial body in the skyline based on the skyline coordinates and the orbital coordinates. It is a program that provides the time of appearance of celestial bodies in the skyline.
A computer, a data input unit that accepts input of date and location information of a reference point,
An orbital coordinate acquisition unit that acquires the orbital coordinates of the celestial body viewed from the reference point on the date using the date and position information.
A height angle acquisition unit that acquires a height angle based on the distance from the reference point and the height of the candidate point for a plurality of candidate points existing in a predetermined direction as viewed from the reference point.
A skyline coordinate acquisition unit that acquires the skyline coordinates by acquiring the maximum value of the height angle in the direction, and
A program for providing an infestation time that functions as a time determination unit that determines the infestation time of a celestial body in the skyline based on the skyline coordinates and the orbital coordinates.

Images