WO2022013948A1

WO2022013948A1 - Trajectory calculation device, trajectory calculation method, and trajectory calculation system

Info

Publication number: WO2022013948A1
Application number: PCT/JP2020/027399
Authority: WO
Inventors: 貴雄遠藤; 俊行安藤; 隆高根澤; 多伸福田
Original assignee: 三菱電機株式会社
Priority date: 2020-07-14
Filing date: 2020-07-14
Publication date: 2022-01-20
Also published as: JPWO2022013948A1; JP7134382B2

Abstract

A trajectory calculation device comprising a first angle calculation unit for using the pointing directions at a first time and a second time of a first light condensing device for condensing light emitted from an illumination light source and reflected by a flying object and the light reflection positions from the flying object at the first time and second time to calculate first angles of the flying object as viewed from the first light condensing device at the first time and second time, a second angle calculation unit for using the pointing directions at the first time and second time of a second light condensing device for condensing light reflected by the flying object and the light reflection positions from the flying object at the first time and second time to calculate second angles of the flying object as viewed from the second light condensing device at the first time and second time, and a trajectory calculation unit for using the calculated first angles, the calculated second angles, and the positions of the first light condensing device and second light condensing device at the first time and second time to calculate the trajectory of the flying object.

Description

Orbit calculation device, orbit calculation method and orbit calculation system

The present disclosure relates to an orbit calculation device, an orbit calculation method, and an orbit calculation system for calculating the orbit of a flying object.

In the orbit calculation system that calculates the orbit of the flying object, the electromagnetic wave is transmitted toward the flying object, while the transmitting / receiving antenna that receives the electromagnetic wave reflected by the flying object and the electromagnetic wave transmitted / received by the transmitting / receiving antenna are used. , There is an orbit calculation system (hereinafter referred to as "conventional orbit calculation system") including a signal processing unit for calculating the orbit of a flying object.
In the conventional orbit calculation system, since the electromagnetic wave reciprocates between the transmitting / receiving antenna and the flying object, the propagation distance of the electromagnetic wave is twice the distance between the transmitting / receiving antenna and the flying object.

By the way, it is equipped with a plurality of receiving devices that receive electromagnetic waves transmitted from a spacecraft that is a flying object, and is used for the time change rate of the phase difference of a plurality of electromagnetic waves received by the plurality of receiving devices to determine the orbit of the spacecraft. There is an orbit calculation device for calculating (see, for example, Patent Document 1). In the orbit arithmetic unit, since the flying object transmits the electromagnetic wave, the propagation distance of the electromagnetic wave in the orbit arithmetic unit is half of the propagation distance of the electromagnetic wave in the conventional orbit calculation system.

JP-A-2007-256004

In the conventional orbit calculation system, the power of the electromagnetic wave transmitted and received by the transmission / reception antenna becomes weaker in inverse proportion to the square of the propagation distance of the electromagnetic wave. Therefore, the conventional orbit calculation system has a problem that the orbit of the flying object may not be calculated depending on the distance between the transmitting / receiving antenna and the flying object.
The propagation distance of the electromagnetic wave in the orbit calculation device disclosed in Patent Document 1 is shorter than the propagation distance of the electromagnetic wave in the conventional orbit calculation system. Therefore, even for a distant flying object whose orbit cannot be calculated by the conventional orbit calculation system, the orbit calculation device disclosed in Patent Document 1 may be able to calculate the orbit. However, in the orbit calculation device, if the flying object does not have a transmitter for transmitting electromagnetic waves, the electromagnetic waves cannot be transmitted. Therefore, the orbit cannot be calculated for the flying object without a transmitter. Therefore, if the orbit arithmetic unit is applied to a conventional orbit calculation system, the orbit of a flying object without a transmitter cannot be calculated. That is, the orbit arithmetic unit cannot be applied to an orbit calculation system that needs to calculate an orbit even for a flying object that does not have a transmitter.

This disclosure is made to solve the above-mentioned problems, and may be able to calculate the orbit of a distant projectile that cannot be calculated by the conventional orbit calculation system. Orbit calculation device, orbit calculation method, and orbit calculation. The purpose is to get the system.

The orbit calculation device according to the present disclosure includes the direction directions of the first light collector that collects the light reflected by the projectile after being emitted from the illumination light source at the first time and the second time. , The first time and the second time as the first angle when the light is viewed from the first condensing device from the respective reflection positions of the light with respect to the flying object at the first time and the second time. The direction of direction at the first time and the second time of the first angle calculation unit that calculates each first angle in the above and the second light collector that collects the light reflected by the flying object. And, from the respective reflection positions of the light with respect to the flying object at the first time and the second time, the first time and the second time are set as the second angle when the flying object is viewed from the second condensing device. A second angle calculation unit that calculates each second angle at time, a first angle calculated by the first angle calculation unit, and a second angle calculated by the second angle calculation unit. It is provided with an orbit calculation unit that calculates the orbit of the flying object from the position of the first light condensing device and the position of the second light condensing device at the first time and the second time.

According to the present disclosure, it may be possible to calculate the orbit of a distant projectile that cannot be calculated by the conventional orbit calculation system.

It is explanatory drawing which shows an example of the positional relationship of the illumination light source 1, the flying object 2, the 1st optical angle measuring station 3 and the 2nd optical angle measuring station 4 in the orbit calculation system which concerns on Embodiment 1. FIG. It is a block diagram which shows the trajectory calculation system which concerns on Embodiment 1. It is a hardware block diagram which shows the hardware of the trajectory calculation apparatus 13 which concerns on Embodiment 1. FIG. 6 is a hardware configuration diagram of a computer when the trajectory calculation device 13 is realized by software, firmware, or the like. It is a concrete block diagram which shows the trajectory calculation system which concerns on Embodiment 1. It is a flowchart which shows the trajectory calculation method which is the processing procedure of the trajectory calculation apparatus 13 shown in FIG. It is explanatory drawing which shows the calculation accuracy of the distance by the angle calculation unit 53 of the 1st optical angle measuring station 3. It is explanatory drawing which shows the calculation accuracy of each angle and the calculation accuracy of each distance at three observation time t ₀ , t ₁ , t _2. Measurement accuracy when the angle measuring station installed on the ground emits radio waves and the angle measuring station receives the radio waves reflected by the flying object 2 to calculate the distance from the angle measuring station to the flying object 2. It is explanatory drawing which shows. It is explanatory drawing which shows the calculation accuracy of each angle and the calculation accuracy of each distance at two observation times t ₁ and t _2. It is explanatory drawing which shows the calculation accuracy of an angle and the calculation accuracy of a distance in the trajectory calculation system which concerns on Embodiment 1. FIG. It is explanatory drawing which shows the calculation accuracy of the orbit in the orbit calculation system which concerns on Embodiment 1. FIG. It is explanatory drawing which shows an example of the positional relationship of the illumination light source 1, the flying object 2, the 1st optical angle measuring station 3 and the 2nd optical angle measuring station 4 in the orbit calculation system which concerns on Embodiment 2. FIG. It is explanatory drawing which shows the calculation accuracy of an angle and the calculation accuracy of a distance in the trajectory calculation system which concerns on Embodiment 2. FIG. It is explanatory drawing which shows an example of the positional relationship of the flying object 2, the 1st optical angle measuring station 3 and the 2nd optical angle measuring station 4 in the orbit calculation system which concerns on Embodiment 3. FIG. It is a concrete block diagram which shows the trajectory calculation system which concerns on Embodiment 3. It is explanatory drawing which shows an example of the positional relationship of the flying object 2, the 1st optical angle measuring station 3 and the 2nd optical angle measuring station 4 in the orbit calculation system which concerns on Embodiment 4. FIG.

Hereinafter, in order to explain the present disclosure in more detail, a mode for carrying out the present disclosure will be described in accordance with the attached drawings.

Embodiment 1.
FIG. 1 is an explanatory diagram showing an example of the positional relationship between the illumination light source 1, the flying object 2, the first optical angle measuring station 3 and the second optical measuring station 4 in the orbit calculation system according to the first embodiment. ..
The illumination light source 1 is a star such as the sun that emits light. However, the illumination light source 1 may be any light source that emits light, and is not limited to a star.
The projectile 2 is a moving body that orbits the earth and reflects the light emitted from the illumination light source 1.

The first optical angle measuring station 3 is a fixed station fixed on the ground of the earth. In the orbit calculation system shown in FIG. 1, the first optical angle measuring station 3 is a fixed station. However, this is only an example, and the first optical measuring station 3 may be a portable station whose position changes.
The first optical angle measuring station 3 includes at least the first condensing device 11, the first angle calculating unit 21, and the orbit calculating unit 23 shown in FIG. 2. The specific components mounted on the first optical angle measuring station 3 are shown in FIG.
The second optical angle measuring station 4 is a fixed station or a portable station.
The second optical angle measuring station 4 is equipped with at least the second light collecting device 12 and the second angle calculating unit 22 shown in FIG. The specific components mounted on the second optical angle measuring station 4 are shown in FIG.
The reference signal source 5 uses the reference signal used by the first optical measuring station 3 and the second optical measuring station 4 for time synchronization as the first optical measuring station 3 and the second optical measuring station 4. Send to each of.

FIG. 2 is a configuration diagram showing an orbit calculation system according to the first embodiment.
The orbit calculation system shown in FIG. 2 includes a first light condensing device 11, a second light condensing device 12, and an orbit calculation device 13.
FIG. 3 is a hardware configuration diagram showing the hardware of the trajectory calculation device 13 according to the first embodiment.
The first condensing device 11 is realized by, for example, one or more telescopes among a refraction type telescope that utilizes refraction of light and a reflection type telescope that utilizes reflection of light.
The first condensing device 11 condenses the light emitted by the illumination light source 1 and then reflected by the projectile 2.
The second condensing device 12 is realized by, for example, one or more telescopes among a refraction type telescope and a reflection type telescope.
The second condensing device 12 condenses the light emitted by the illumination light source 1 and then reflected by the projectile 2.

The orbit calculation device 13 includes a first angle calculation unit 21, a second angle calculation unit 22, and an orbit calculation unit 23.
The first angle calculation unit 21 is realized by, for example, the first angle calculation circuit 31 shown in FIG.
The first angle calculation unit 21, each of the directivity direction in the first time t ₁ and second time t ₂ of the first condenser 11, a first time t ₁ and the light with respect to the projectile 2 from the respective reflection positions in the second time t _2, the calculating the first angle is an angle viewed projectile 2 from the first condenser 11. The first angle includes an azimuth angle Az _A seen from the first light collecting device 11 and an elevation angle Elv _A seen from the first light collecting device 11. .. Incidentally, the azimuth angle at a first time _{t 1} is, _{Az A1,} the azimuth angle at a second time _{t 2} is, _{Az A2,} elevation angle at a first time _{t 1} is, _{Elv A1,} the second time t The elevation angle at ₂ _{is Elv A2} .
The first angle calculation unit 21 outputs each of the azimuth angles Az _A1 and Az _A2 and the elevation angles Elv _A1 and Elv _A2 to the orbit calculation unit 23 as the calculated first angle.

The second angle calculation unit 22 is realized by, for example, the second angle calculation circuit 32 shown in FIG.
The second angle calculation unit 22, each of the directivity direction in the first time t ₁ and second time t ₂ of the second condenser 12, a first time t ₁ and the light with respect to the projectile 2 From each reflection position at the second time t2, a _second angle, which is the angle at which the projectile 2 is viewed from the second light collecting device 12, is calculated. The second angle includes an azimuth angle Az _B as seen from the second light collector 12 and an elevation angle Elv _B as seen from the second light collector 12. .. Incidentally, the azimuth angle at a first time _{t 1} is, _{Az B1,} the azimuth angle at a second time _{t 2} is, _{Az B2,} elevation angle at a first time _{t 1} is, _{Elv B1,} the second time t The elevation angle at ₂ _{is Elv B2} .
The second angle calculation unit 22 outputs each of the azimuth angles Az _B1 and Az _B2 and the elevation angles Elv _B1 and Elv _B2 to the orbit calculation unit 23 as the calculated second angle.

The orbit calculation unit 23 is realized by, for example, the orbit calculation circuit 33 shown in FIG.
The orbit calculation unit 23 is calculated by the first angle at the first time t ₁ and the second time t ₂ calculated by the first angle calculation unit 21, and by the second angle calculation unit 22. In addition, the second angles at _{the first time t 1} and the second time t _{2 are acquired.}
Trajectory calculation unit 23, and the respective second angle and each of the first angle at a first time t ₁ and second time t _2, the at a first time t ₁ and second time t _2, the The positions of the first light collector 11 at the first time t ₁ and the second time t _2, and the positions of the second light collector 12 at the first time t ₁ and the second time t ₂ , respectively. The trajectory of the projectile 2 is calculated from the position of.

In FIG. 2, each of the first angle calculation unit 21, the second angle calculation unit 22, and the trajectory calculation unit 23, which are the components of the trajectory calculation device 13, is realized by dedicated hardware as shown in FIG. I'm assuming something. That is, it is assumed that the trajectory calculation device 13 is realized by the first angle calculation circuit 31, the second angle calculation circuit 32, and the trajectory calculation circuit 33.
Each of the first angle calculation circuit 31, the second angle calculation circuit 32, and the trajectory calculation circuit 33 is, for example, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, or an ASIC (Application Specific Integrated Circuit). ), FPGA (Field-Programmable Gate Array), or a combination thereof.

The components of the trajectory calculation device 13 are not limited to those realized by dedicated hardware, but the trajectory calculation device 13 is realized by software, firmware, or a combination of software and firmware. It is also good.
The software or firmware is stored as a program in the memory of the computer. A computer means hardware for executing a program, and corresponds to, for example, a CPU (Central Processing Unit), a central processing unit, a processing unit, a computing device, a microprocessor, a microcomputer, a processor, or a DSP (Digital Signal Processor). do.

FIG. 4 is a hardware configuration diagram of a computer when the trajectory calculation device 13 is realized by software, firmware, or the like.
When the trajectory calculation device 13 is realized by software, firmware, or the like, a program for causing a computer to execute each processing procedure in the first angle calculation unit 21, the second angle calculation unit 22, and the trajectory calculation unit 23 is provided. It is stored in the memory 41. Then, the processor 42 of the computer executes the program stored in the memory 41.
Further, FIG. 3 shows an example in which each of the components of the trajectory calculation device 13 is realized by dedicated hardware, and FIG. 4 shows an example in which the trajectory calculation device 13 is realized by software, firmware, or the like. .. However, this is only an example, and some components in the trajectory calculation device 13 may be realized by dedicated hardware, and the remaining components may be realized by software, firmware, or the like.

FIG. 5 is a specific configuration diagram showing the trajectory calculation system according to the first embodiment.
In the trajectory calculation system shown in FIG. 5, the first optical angle measuring station 3 has a first light collecting device 11, a light shielding device 51, a position detection unit 52, an angle calculation unit 53, a time calibration unit 54, a counter 55, and a control. It includes a device 56, a pointing device 57, a communication unit 58, a recording device 59, and an orbit calculation unit 60.
In the trajectory calculation system shown in FIG. 5, the second optical angle measuring station 4 has a second light collecting device 12, a light shielding device 71, a position detection unit 72, an angle calculation unit 73, a time calibration unit 74, a counter 75, and a control. It includes a device 76, a directional device 77, and a communication unit 78.

The shading device 51 is realized by, for example, a mechanical shutter or an electronic shutter.
The shading device 51 is arranged in an optical path between the first light collecting device 11 and the position detection unit 52.
The light-shielding device 51 controls the opening and closing of the shutter according to the exposure time controlled by the control device 56. By controlling the opening and closing of the shutter, the light-shielding device 51 alternately blocks the light collected by the first light-collecting device 11 and transmits the light.

The first angle calculation unit 21 shown in FIG. 5 is the same angle calculation unit as the first angle calculation unit 21 shown in FIG.
The first angle calculation unit 21 shown in FIG. 5 includes a position detection unit 52 and an angle calculation unit 53.
The position detection unit 52 is realized by, for example, a CCD (Charge Coupled Device) image sensor or a CMOS (Complementary Metal Oxide Sensor) image sensor.
_{When the position detection unit 52 acquires the image pickup command output from the control device 56 at the first time t 1} and the second time t ₂ , the position detection unit 52 detects the light transmitted through the light shielding device 51 by detecting the light transmitted through the light shielding device 51. An image of the light intensity showing the projectile 2 is captured.
That is, the position detection unit 52 captures a light intensity image in which the projectile 2 is reflected by detecting the light transmitted through the light shielding device 51 at _{the first time t 1, and the second time t.} By detecting the light transmitted through the shading device 51 at the time of _{2, the light intensity image in which the projectile 2 is reflected is captured.}
The position detection unit 52 outputs each of the light intensity image at _{the first time t 1 and} the light intensity image at the second time t _{2 to the angle calculation unit 53.}
In the trajectory calculation system shown in FIG. 5, the position detection unit 52 is provided inside the first angle calculation unit 21. However, this is only an example, and the position detection unit 52 may be provided outside the first angle calculation unit 21.

In the trajectory calculation system shown in FIG. 5, the position detection unit 52 is realized by a CCD image sensor or a CMOS image sensor. However, this is only an example, and the position detection unit 52 is simply a sensor that detects the presence or absence of light, and is realized by a sensor that detects light when the projectile 2 enters the field of view. You may. However, the field of view of the sensor that simply detects the presence or absence of light is narrower than the field of view of the CCD image sensor. Therefore, when the position detection unit 52 is realized by a sensor that simply detects the presence or absence of light, the directing direction of the first light collecting device 11 is finely controlled as compared with the case where the position detection unit 52 is realized by a CCD image sensor. There is a need.
When the position detection unit 52 is realized by a CCD image sensor or a CMOS image sensor, even if the directing direction of the first light collecting device 11 deviates from the center position of the light intensity image, the flying object If 2 is included in the light intensity image, the flying object 2 can be detected. Therefore, when the position detection unit 52 is realized by a CCD image sensor or a CMOS image sensor, the orientation of the first light collecting device 11 is higher than that realized by a sensor that simply detects the presence or absence of light. The direction can be easily controlled. Further, when the position detection unit 52 is realized by a CCD image sensor or a CMOS image sensor, it grasps whether the moving direction of the flying object 2 is the vertical direction of the visual field, the horizontal direction of the visual field, or the like. can do. In addition, it is possible to grasp the presence or absence of leaked light from a star or another flying object.

Angle calculation unit 53, the controller 56 obtains the orientation of the first focusing device 11 at a first time t _1, from the position detection unit 52, obtains the light intensity image at a first time t ₁ do.
The angle calculation unit 53 detects the position of the light in the light intensity image at the _first time t 1 as the reflection position of the light with respect to the projectile 2 at the first time t _1.
The angle calculation unit 53 calculates the first angle at the first time t _{1 from} the directing direction of the first light collecting device 11 at the first time t _{1 and} the reflection position at the first time t ₁ . do.
The angle calculation unit 53, the controller 56 obtains the first orientation of the condenser 11 at the second time t _2, from the position detection unit 52, the light intensity image at a second time t ₂ To get.
The angle calculation unit 53 detects the position of the light in the light intensity image at the _second time t 2 as the reflection position of the light with respect to the projectile 2 at the second time t _2.
The angle calculation unit 53 calculates the first angle at the second time t _{2 from} the directing direction of the first light collecting device 11 at the second time t _{2 and} the reflection position at the second time t ₂ . do.

The time calibration unit 54 includes, for example, a GPS (Global Positioning System) receiver, a receiving antenna, and a clock.
The receiving antenna of the time calibration unit 54 receives the reference signal transmitted from the reference signal source 5.
When the receiving antenna receives the reference signal, the GPS receiver of the time calibration unit 54 repeatedly receives the GPS signal transmitted from the GPS satellite.
The time calibration unit 54 acquires the position of the first light collecting device 11 and the current time from the GPS signal received by the GPS receiver. The position information included in the GPS signal indicates the position of the first optical angle measuring station 3 equipped with the time calibration unit 54. Since the first condensing device 11 is mounted on the first optical angle measuring station 3, the time calibration unit 54 sets the position of the first optical measuring station 3 to the first condensing. It is acquired as the position of the device 11.
The time calibration unit 54 calibrates the time of the internal clock using the acquired current time, and outputs the time after calibration to the counter 55.
In the trajectory calculation system shown in FIG. 5, the time calibration unit 54 calibrates the time of the internal clock and outputs the time after calibration to the counter 55. However, this is only an example, and the time calibration unit 54 may output the current time acquired from the GPS signal to the counter 55 as the time after calibration.
In the trajectory calculation system shown in FIG. 5, the time calibration unit 54 acquires the current time from the radio wave received by the GPS receiver. However, this is only an example, and the time calibration unit 54 may acquire the current time from a standard radio wave represented by a radio clock or the like, or a network device represented by NTP (Network Time Protocol). The current time may be acquired by using the protocol for time synchronization of.
For example, the sun moves about 15 arcseconds (= 15/3600 degrees) per second due to diurnal motion, so if accuracy of the arc arc is required to determine the position, the time calibration unit 54 may use the time of the internal clock. Is calibrated with a sub-second accuracy.
The time calibration unit 54 outputs the position of the first light collecting device 11 to the control device 56.
Further, the time calibration unit 54 acquires the position of the first _{light collecting device 11 at the first time t 1} and the position of the first light collecting device 11 at the second time t _{2 from the control device 56.} Then, the _{position of the first light collecting device 11 at the first time t 1} and the position of the first light collecting device 11 at the second time t ₂ are output to the recording device 59.

The counter 55 acquires the time after calibration by the time calibration unit 54 and measures the elapsed time from a certain typical time. A typical time is, for example, the calculation start time of the orbit of the projectile 2.
The counter 55 outputs the measured elapsed time to the control device 56.

The control device 56 is realized by, for example, a semiconductor integrated circuit on which a CPU is mounted, or a program board such as an FPGA.
The control device 56 grasps each of _{the first time t 1} and the second time t ₂ as the observation time of the projectile 2 according to the elapsed time output from the counter 55.
Further, the control device 56 acquires the past orbit information of the flying object 2 from the orbit database 80 via the communication unit 58, and obtains the approximate past orbit information of the flying object 2 at the _{first time t1 from the past orbit information.} The position and the approximate position of the projectile _{2 at} the second time t2 are grasped.

Controller 56, among the position of the first condenser 11 output from the time correcting unit 54, the position of the first focusing device 11 at a first time t _1, at a second time t ₂ The position of the first light collector 11 is specified.
Controller 56, the position of the first focusing device 11 at a first time t _1, and a rough position of the projectile 2 in the first time t _1, the projectile with respect to the first condenser 11 calculating the relative position at a first time t ₁ of 2.
Further, the control device 56, the position of the first condenser 11 at the second time t _2, the from the approximate position of the projectile 2 in the second time t _2, the relative first condenser 11 calculating the relative position at a second time t ₂ of the projectile 2.

Controller 56, from the relative position at a first time t _1, the direction in which the first condenser 11 sees the projectile 2 at the first time t _1, i.e., first at a first time t ₁ The direction of direction of the light collecting device 11 of 1 is grasped.
Further, the control device 56, from the relative position at the second time t _2, the direction in which the first condenser 11 sees the projectile 2 at the second time t _2, i.e., the second time t ₂ The direction of direction of the first light collector 11 in the above is grasped.

Controller 56, the orientation of the first focusing device 11 at a first time t ₁ and outputs to each of the directing device 57 and the angle calculation unit 53, a first condenser at a second time t ₂ The pointing direction of 11 is output to each of the pointing device 57 and the angle calculation unit 53.
The control device 56 controls the light-shielding device 51 so as to open the shutter for a predetermined exposure time at _{the first time t 1} and the second time t _2.
The control device 56 outputs an image pickup command to the position detection unit 52 at _{the first time t 1} and the second time t _2.

The directional device 57 is realized by a rotary stage and a motor having one or more rotation axes. An altazimuth mount, an equatorial mount, or the like can be considered as a rotation stage having one or more rotation axes.
A first condensing device 11 is mounted on the rotating stage of the directional device 57.
Directing device 57, the orientation direction of the first condenser 11 at a first time t ₁ is consistent with the orientation of the first focusing device 11 at a first time t ₁ the controller 56 is outputted By driving the motor, the rotary stage is rotated.
Directing device 57, the orientation direction of the first condenser 11 at a second time t ₂ is consistent with the orientation of the first focusing device 11 at a second time t ₂ the controller 56 is outputted By driving the motor, the rotary stage is rotated.

The communication unit 58 acquires the past orbit information of the flying object 2 from the orbit database 80, and outputs the orbit information to the control device 56.
The communication unit 58 is transmitted from the second communication unit 78 of the optical angle measuring station 4, receives the respective second angle and the second angle at a second time t ₂ at a first time t ₁ ..
The communication unit 58 receives transmitted from the second communication unit 78 of the optical angle measuring station 4, the first respective positions at the time t ₁ and second time t ₂ of the second condenser 12 ..
The communication unit 58 to record each of the second angle and the second angle at a second time t ₂ at a first time t ₁ to the recording device 59, a first time of the second condenser 12 Have the recording device 59 record each position at _{t 1} and the second time t _2.

The orbit calculation unit 23 shown in FIG. 5 is the same orbit calculation unit as the orbit calculation unit 23 shown in FIG.
The orbit calculation unit 23 shown in FIG. 5 includes a recording device 59 and an orbit calculation unit 60.
The recording device 59 is realized by a recording medium such as a RAM (Random Access Memory) or a hard disk.
Recording apparatus 59, each of the first angle and the first angle at a second time t ₂ at a first time t ₁ is recorded, a second angle and a second time at a first time t ₁ Record each of the second angles at t _2.
Recording apparatus 59, each of the positions at a first time t ₁ and second time t ₂ of the first condenser 11, a first time t ₁ and the second of the second condenser 12 recording the respective positions at the time t _2.

Trajectory calculation unit 60 from the recording device 59, respectively and each of the first angle at a first time t ₁ and second time t _2, at a first time t ₁ and second time t ₂ the Get the angle of 2.
In addition, the orbit calculation unit 60 is from the recording device 59 to the _{respective positions of the first light collecting device 11 at the first time t 1} and the second time t ₂ and the first of the second light collecting device 12. The respective positions at the time t ₁ and the second time t _{2 of the first time t 1 and the second time t 2 are acquired.}
Trajectory calculation unit 60, a first condenser respective position at time t ₁ and second time t ₂ the first 11, a first time t ₁ and the second second condenser 12 time from the respective positions at time t _2, the time of the first time t _1, a first condenser 11 and the inter-device distance is a distance between the second condenser 12, the second The distance between devices at _{t 2 is calculated.}
Trajectory calculation unit 60, each of the inter-device distances at a first time t ₁ and second time t _2, the respectively first angle at a first time t ₁ and second time t _2, the first from the first time t ₁ and second time t ₂ and each of the second angle, the respective positions at a first time t ₁ and second time t ₂ of the first condenser 11, the The first distance from the light collecting device 11 of 1 to the projectile 2 is calculated. That is, the orbit calculation unit 60 calculates the first distances at the first _{time t 1} and the second time t _2.
Trajectory calculation unit 60, each of the inter-device distances at a first time t ₁ and second time t _2, the respectively first angle at a first time t ₁ and second time t _2, the first from the first time t ₁ and second time t ₂ and each of the second angle, the respective positions at a first time t ₁ and second time t ₂ of the second condenser 12, the The second distance from the light collector 12 of 2 to the projectile 2 is calculated. That is, the orbit calculation unit 60 calculates the second distances at _{the first time t 1} and the second time t _2.
Trajectory calculation unit 60, and the respective second angle and each of the first angle at a first time t ₁ and second time t _2, the at a first time t ₁ and second time t _2, the from each of the first distance at a first time t ₁ and second time t _2, the a respective second distance at a first time t ₁ and second time t _2, the position of the projectile 2 And each of the speeds are calculated. That is, the orbit calculation unit 60 calculates the respective positions of the flying object 2 at the first time t ₁ and the second time t _2. Alternatively, trajectory calculation unit 60, and each of the first angle at a first time t ₁ and second time t _2, each of the second angle at a first time t ₁ and second time t ₂ When, from each of the first distance at a first time t ₁ and second time t _2, the a respective second distance at a first time t ₁ and second time t _2, the projectile 2 calculating a first position and velocity at time t ₁ of the. The velocity is obtained as the difference between the positions of the projectile _{2 at} _{the time from time t 1} to time t 2.
The orbit calculation unit 60 determines the flying object from _{the respective positions of the flying object 2 at the first time t 1} and the second time t ₂ , or the position and the speed of the first time t ₁ of the flying object 2. Calculate the orbit of 2.
In the trajectory calculation system shown in FIG. 5, the recording device 59 is provided inside the trajectory calculation unit 23. However, this is only an example, and the recording device 59 may be provided outside the trajectory calculation unit 23.

The shading device 71 is realized by, for example, a mechanical shutter or an electronic shutter.
The shading device 71 is arranged in an optical path between the second light collecting device 12 and the position detection unit 72.
The light-shielding device 71 controls the opening and closing of the shutter according to the exposure time controlled by the control device 76. By controlling the opening and closing of the shutter, the light-shielding device 71 alternately blocks the light collected by the second light-collecting device 12 and transmits the light.

The second angle calculation unit 22 shown in FIG. 5 is the same angle calculation unit as the second angle calculation unit 22 shown in FIG.
The second angle calculation unit 22 shown in FIG. 5 includes a position detection unit 72 and an angle calculation unit 73.
The position detection unit 72 is realized by, for example, a CCD image sensor or a CMOS image sensor.
_{When the position detection unit 72 acquires the image pickup command output from the control device 76 at the first time t 1} and the second time t ₂ , the position detection unit 72 detects the light transmitted through the light shielding device 71 by detecting the light transmitted through the light shielding device 71. An image of the light intensity showing the projectile 2 is captured.
That is, the position detection unit 72 captures a light intensity image in which the projectile 2 is reflected by detecting the light transmitted through the light shielding device 71 at _{the first time t 1, and the second time t.} By detecting the light transmitted through the shading device 71 at the time of _{2, the light intensity image in which the projectile 2 is reflected is captured.}
The position detection unit 72 outputs each of the light intensity image at _{the first time t 1 and} the light intensity image at the second time t _{2 to the angle calculation unit 73.}

Angle calculation unit 73, the controller 76 obtains the orientation of the second condenser 12 at a first time t _1, the position detection unit 72, obtains the light intensity image at a first time t ₁ do.
The angle calculation unit 73 detects the position of the light in the light intensity image at the _first time t 1 as the reflection position of the light with respect to the projectile 2 at the first time t _1.
The angle calculation unit 73 calculates _{the second angle at the first time t 1 from} the directing direction of the second light collecting device 12 at the first time t _{1 and} the reflection position at the first time t ₁ . do.
Further, the angle calculation unit 73 acquires the pointing direction of the second light collecting device 12 at the second _{time t 2} from the control device 76, and the position detection unit 72 obtains the light intensity image at the _{second time t 2.} To get.
The angle calculation unit 73 detects the position of the light in the light intensity image at the _second time t 2 as the reflection position of the light with respect to the projectile 2 at the second time t _2.
The angle calculation unit 73 calculates the second angle at the second time t _{2 from} the directing direction of the second light collecting device 12 at the second time t _{2 and} the reflection position at the second time t ₂ . do.
The angle calculation unit 73 outputs each of the second angle at _{the first time t 1} and the second angle at the second time t _{2 to the communication unit 78.}
In the trajectory calculation system shown in FIG. 5, the position detection unit 72 is provided inside the second angle calculation unit 22. However, this is only an example, and the position detection unit 72 may be provided outside the second angle calculation unit 22.

The time calibration unit 74 includes, for example, a GPS receiver, a receiving antenna, and a clock.
The receiving antenna of the time calibration unit 74 receives the reference signal transmitted from the reference signal source 5.
When the receiving antenna receives the reference signal, the GPS receiver of the time calibration unit 74 repeatedly receives the GPS signal transmitted from the GPS satellite.
The time calibration unit 74 acquires the position of the second light collecting device 12 and the current time from the GPS signal received by the GPS receiver. The position information included in the GPS signal indicates the position of the second optical angle measuring station 4 equipped with the time calibration unit 74. Since the second light collecting device 12 is mounted on the second light measuring station 4, the time calibration unit 74 sets the position of the second light measuring station 4 to the second light collecting station 4. It is acquired as the position of the device 12.
The time calibration unit 74 calibrates the time of the internal clock using the acquired current time.
In the trajectory calculation system shown in FIG. 5, the time calibration unit 74 calibrates the time of the internal clock and outputs the time after calibration to the counter 75. However, this is only an example, and the time calibration unit 74 may output the current time acquired from the GPS signal to the counter 75 as the time after calibration.
In the trajectory calculation system shown in FIG. 5, the time calibration unit 74 acquires the current time from the radio wave received by the GPS receiver. However, this is only an example, and the time calibration unit 74 may acquire the current time from a standard radio wave represented by a radio clock or the like, or for time synchronization of a network device represented by NTP. The current time may be obtained by using the protocol of.
For example, the sun moves about 15 arcseconds per second due to diurnal motion, so if accuracy of arcseconds is required to determine the position, the time calibration unit 74 calibrates the time of the internal clock with subsecond accuracy. do.
The time calibration unit 74 outputs the position of the second light collecting device 12 to the control device 76.
Further, the time calibration unit 74 acquires the position of the second light collecting device 12 at _{the first time t 1} and the position of the second light collecting device 12 at the second time t _{2 from the control device 76.} Then, the position of the second light collecting device 12 at the first time t ₁ and the position of the second light collecting device 12 at the second time t ₂ are output to the communication unit 78.

The counter 75 acquires the time after calibration by the time calibration unit 74 and measures the elapsed time from a certain typical time.
The counter 75 outputs the measured elapsed time to the control device 76.

The control device 76 is realized by, for example, a semiconductor integrated circuit on which a CPU is mounted, or a program board such as an FPGA.
The control device 76 grasps each of _{the first time t 1} and the second time t ₂ as the observation time of the projectile 2 according to the elapsed time output from the counter 75.
Further, the control device 76 acquires the past orbit information of the flying object 2 from the orbit database 80 via the communication unit 78, and from the past orbit information, the control device 76 approximates the flying object 2 at the _{first time t1.} The position and the approximate position of the projectile _{2 at} the second time t2 are grasped.

Controller 76, among the position of the second condenser 12 output from the time correcting unit 74, and the position of the second condenser 12 at a first time t _1, at a second time t ₂ The position of the second light collector 12 is specified.
Control device 76, a position of the second condenser 12 at a first time t _1, and a rough position of the projectile 2 in the first time t _1, the projectile with respect to the second condenser 12 calculating the relative position at a first time t ₁ of 2.
Further, the control device 76, the position of the second condenser 12 at the second time t _2, the from the approximate position of the projectile 2 in the second time t _2, the relative second condenser 12 calculating the relative position at a second time t ₂ of the projectile 2.

Control device 76, from the relative position at a first time t _1, the direction in which the second condenser 12 sees the projectile 2 at the first time t _1, i.e., first at a first time t ₁ Grasp the directing direction of the light collecting device 12 of 2.
Further, the control device 76, from the relative position at the second time t _2, the direction in which the second condenser 12 when the second time t ₂ to see the flying body 2, i.e., the second time t ₂ The direction of direction of the second light collecting device 12 in the above is grasped.

Controller 76, the orientation of the second condenser 12 at a first time t ₁ and outputs to each of the directing device 77 and the angle calculation unit 73, the second condensing device at a second time t ₂ The direction of the 12 is output to the pointing device 77 and the angle calculation unit 73, respectively.
The control device 76 controls the light-shielding device 71 so that the shutter is opened for a predetermined exposure time at _{the first time t 1} and the second time t _2.
The control device 76 outputs an image pickup command to the position detection unit 72 at _{the first time t 1} and the second time t _2.

The directional device 77 is realized by a rotary stage and a motor having one or more rotation axes.
A second condensing device 12 is mounted on the rotating stage of the directional device 77.
Directing device 77, the directivity direction of the second condenser 12 at a first time t ₁ is consistent with the orientation of the second condenser 12 at a first time t ₁ the controller 76 is output By driving the motor, the rotary stage is rotated.
Directing device 77, the directivity direction of the second condenser 12 at a second time t ₂ is consistent with the orientation of the second condenser 12 at a second time t ₂ the controller 76 is output By driving the motor, the rotary stage is rotated.

The communication unit 78 acquires the past orbit information of the flying object 2 from the orbit database 80, and outputs the orbit information to the control device 76.
The communication unit 78 is outputted from the angle calculating unit 73, and transmits the respective second angle and the second angle at a second time t ₂ at a first time t ₁ to the communication unit 58.
The communication unit 78 is output from the time correcting unit 74, and transmits to the second first time point t ₁ and the second communication unit 58 to respective positions at the time t ₂ of the condenser 12.

The orbital database 80 is realized by, for example, a RAM or a recording medium such as a hard disk.
The orbit database 80 stores past orbit information in the projectile 2.
In the trajectory calculation system shown in FIG. 5, the trajectory database 80 is provided outside the trajectory calculation system. However, this is only an example, and the trajectory database 80 may be provided inside the trajectory calculation system.
Further, the orbital database 80 may be connected to a network (not shown), and the orbital database 80 may be connected to the orbital calculation system via the network.

Next, the operation of the trajectory calculation system shown in FIG. 5 will be described.
FIG. 6 is a flowchart showing a trajectory calculation method which is a processing procedure of the trajectory calculation device 13 shown in FIG.
First, the reference signal source 5 uses the reference signal used by the first optical measuring station 3 and the second optical measuring station 4 for time synchronization as the first optical measuring station 3 and the second optical measuring station 3. Send to each of the stations 4.
The receiving antenna of the time calibration unit 54 in the first optical angle measuring station 3 receives the reference signal transmitted from the reference signal source 5.
When the receiving antenna receives the reference signal, the GPS receiver of the time calibration unit 54 repeatedly receives the GPS signal transmitted from the GPS satellite.
The time calibration unit 54 acquires the position of the first light collecting device 11 and the current time from the GPS signal received by the GPS receiver.
The time calibration unit 54 calibrates the time of the internal clock using the acquired current time.
The time calibration unit 54 outputs the position of the first light collecting device 11 to the control device 56.

When the receiving antenna receives the reference signal, the GPS receiver of the time calibration unit 74 repeatedly receives the GPS signal transmitted from the GPS satellite.
The time calibration unit 74 acquires the position of the second light collecting device 12 and the current time from the GPS signal received by the GPS receiver.
The time calibration unit 74 calibrates the time of the internal clock using the acquired current time.
The time calibration unit 74 outputs the position of the second light collecting device 12 to the control device 76.
The time calibration unit 54 of the first optical angle measuring station 3 calibrates the time of the internal clock, and the time calibration unit 74 of the second optical measuring station 4 calibrates the time of the internal clock. The time is synchronized between the angle measuring station 3 and the second optical angle measuring station 4.
The accuracy of time synchronization is limited by the accuracy of the shutter opening / closing time of the

shading devices

51 and 71. Since the shutter opening / closing time is about millisecond, time synchronization may be performed with an accuracy of about microsecond. The time used by the internal clocks of the

time calibration units

54 and 74 is, for example, Coordinated Universal Time.

The counter 55 of the first optical angle measuring station 3 acquires the time after calibration by the time calibration unit 54 and measures the elapsed time t from a certain representative time. A typical time is, for example, the calculation start time of the orbit of the projectile 2.
The counter 55 outputs the measured elapsed time t to the control device 56.
The counter 75 of the second optical angle measuring station 4 acquires the time after calibration by the time calibration unit 74 and measures the elapsed time t from a certain representative time.
The counter 75 outputs the measured elapsed time t to the control device 76.
The communication unit 58 of the first optical angle measuring station 3 acquires the past orbit information of the flying object 2 from the orbit database 80, and outputs the orbit information to the control device 56.
The communication unit 78 of the second optical angle measuring station 4 acquires the past orbit information of the flying object 2 from the orbit database 80, and outputs the orbit information to the control device 76.

The control device 56 of the first optical angle measuring station 3 acquires the elapsed time t output from the counter 55.
The control device 56 grasps each of _{the first time t 1} and the second time t ₂ as the observation time of the flying object 2 according to the elapsed time t.
Further, the control unit 56 acquires orbit information output from the communication unit 58, from the track information, and approximate location of the projectile 2 in the first time t _1, the projectile at the second time t ₂ 2 Grasp the approximate position.

Controller 56, among the position of the first condenser 11 output from the time correcting unit 54, the position of the first focusing device 11 at a first time t _1, at a second time t ₂ The position of the first light collector 11 is specified.
Controller 56, the position of the first focusing device 11 at a first time t _1, and a rough position of the projectile 2 in the first time t _1, the projectile with respect to the first condenser 11 calculating the relative position at a first time t ₁ of 2. Since the process of calculating the relative position of the projectile 2 with respect to the first condensing device 11 is a known technique, detailed description thereof will be omitted.
Further, the control device 56, the position of the first condenser 11 at the second time t _2, the from the approximate position of the projectile 2 in the second time t _2, the relative first condenser 11 calculating the relative position at a second time t ₂ of the projectile 2.

Controller 56, the orientation of the first focusing device 11 at a first time t ₁ and outputs to each of the directing device 57 and the angle calculation unit 53, a first condenser at a second time t ₂ The pointing direction of 11 is output to each of the pointing device 57 and the angle calculation unit 53.
The control device 56 controls the light-shielding device 51 so as to open the shutter for a predetermined exposure time at _{the first time t 1} and the second time t _2.
The control device 56 outputs an image pickup command to the position detection unit 52 at _{the first time t 1} and the second time t _2.
The control device 56 outputs _{the position of the first light collecting device 11 at the first time t 1} and the position of the first light collecting device 11 at the second time t ₂ to the time calibration unit 54.

Directing device 57 of the first optical angle measuring station 3, the orientation direction of the first condenser 11 at a first time t ₁ is, controller 56 first time t ₁ in the first output The rotary stage is rotated by driving the motor so as to coincide with the directing direction of the light collecting device 11.
Directing device 57, the orientation direction of the first condenser 11 at a second time t ₂ is consistent with the orientation of the first focusing device 11 at a second time t ₂ the controller 56 is outputted By driving the motor, the rotary stage is rotated.
Since the first condensing device 11 has a function of projecting an object surface onto an image plane and has a reference optical axis, the directing device 57 controls the optical axis of the first condensing device 11. The directing direction of the first condensing device 11 is controlled by adjusting to the directing direction output from the device 56.
Incidentally, the orientation of the control of the first condenser 11 by directing device 57, it takes a certain amount of time, the control device 56, before the observation time is a first time t _1, a first time the orientation of the first focusing device 11 outputs to the directing device 57 in the t _1. Further, the control device 56 outputs the pointing direction of the first condensing device 11 at the second time t ₂ to the pointing device 57 before the observation time reaches the second time t _2.

In the trajectory calculation system shown in FIG. 5, the control device 56 directs the directing direction of the first condensing device 11 to the directing device 57 so that the projectile 2 is within the field of view of the first condensing device 11. I'm instructing. When the field of view of the first light collecting device 11 and the field of view of the position detection unit 52 are different, the control device 56 uses the narrower of the field of view of the first light collecting device 11 and the field of view of the position detection unit 52. The directing device 57 is instructed to direct the directing direction of the first condensing device 11 so that the projectile 2 is within the range of the field of view of.
If the flying object 2 is a geostationary satellite and the first optical measuring station 3 is a fixed station, the relative positions of the first optical measuring station 3 and the flying object 2 do not change. When the relative position does not change, the control device 56 once instructs the directing device 57 in the direction of the first condensing device 11 so that the projectile 2 is within the field of view of the first condensing device 11. Then, after that, the directing direction of the first light collecting device 11 may be fixed. Strictly speaking, the relative position changes slightly due to the difference in the orbital inclination angle, but if the observation time is about several seconds, there is no problem even if the directing direction of the first condensing device 11 is fixed. ..

The first light collecting device 11 of the first optical angle measuring station 3 collects the light emitted from the illumination light source 1 and then reflected by the projectile 2.
The first light-shielding device 51 of the optical angle measuring station 3, when the first time t ₁ and second time t _2, the by the controller 56, a predetermined exposure time, are controlled to open the shutters ..
When the first time t _1, a predetermined exposure time, by opening the shutter of the shading device 51, is condensed light by the first condenser apparatus 11 through the shading device 51, the light It reaches the position detection unit 52.
Further, when the second time t _2, the predetermined exposure time, by opening the shutter of the shading device 51, light condensed by the first condensing device 11 is transmitted through the shading device 51, the The light reaches the position detection unit 52.

Position detection unit 52 of the first optical angle measuring station 3, when the first time t _1, when acquiring the imaging command output from the controller 56, the detection processing of the light transmitted through the shading device 51 Start and take a light intensity image showing the projectile 2.
Position detector 52, the light intensity image captured, as a light intensity image at a first time t _1, and outputs the angle calculation unit 53.
The position detection unit 52, when the second time t _2, the acquiring the imaging command output from the control unit 56 starts the detection process of light transmitted through the shading device 51, the projectile 2 Take an image of the reflected light intensity.
Position detector 52, the light intensity image captured, as a light intensity image at a second time t _2, the output to the angle calculation unit 53.

Angle calculating portion 53 of the first optical angle measuring station 3, the controller 56 obtains the orientation of the first focusing device 11 at a first time t _1, from the position detection unit 52, a first for obtaining an optical intensity image at time t _1.
The angle calculation unit 53 detects the position of the light in the light intensity image at the _first time t 1 as the reflection position of the light with respect to the projectile 2 at the first time t _1.
The angle calculation unit 53 calculates the first angle at the first time t _{1 from} the directing direction of the first light collecting device 11 at the first time t _{1 and} the reflection position at the first time t ₁ . (Step ST1 in FIG. 6).
Angle calculation unit 53, a first removal process of blackout flow noise contained in the light intensity image at time t ₁ of, removal processing of the background light noise included in the light intensity image, or included in the light intensity image It is also possible to carry out attenuation processing of the peripheral light amount and the like, and calculate the first angle from the light intensity image after these processings.
Hereinafter, the first angle calculation process by the angle calculation unit 53 will be specifically described.

The angle calculation unit 53 obtains the viewing angle of the light intensity image on the assumption that the directing direction of the first light collecting device 11 is the center of the field of view of the position detecting unit 52. The viewing angle per pixel can be obtained from the size of one pixel included in the light intensity image and the focal length of the first condensing device 11. The size of one pixel and the focal length of the first light collecting device 11 are already values in the angle calculation unit 53.
The viewing angle of the light intensity image is obtained by multiplying the viewing angle per pixel by the number of pixels from the pixel at the center of the light intensity image to the pixel at the center of the field of view of the light intensity image.
The angle calculation unit 53 can obtain the direction of the position detection unit 52 and the apparent size of the projectile 2 from the center of the field of view of the light intensity image and the viewing angle of the light intensity image. If the direction of the position detection unit 52 and the apparent size of the projectile 2 are obtained, the angle calculation unit 53 determines the positions of the coordinates on the light intensity image and the celestial sphere coordinates indicating the position of the illumination light source 1 in space. You can ask for a relationship.
The coordinate system of celestial sphere coordinates may be an equator coordinate system based on the celestial equator, a horizontal coordinate system based on the horizon seen by the observer, or a galaxy coordinate system based on the galaxy surface. good. The position on the light intensity image can be represented by the position of the celestial sphere coordinates, right ascension, declination, or an angle such as an azimuth angle and an elevation angle.

The angle calculation unit 53 determines that the point image, which is a group of pixels in which pixels having a pixel value equal to or larger than a threshold value among a plurality of pixels included in the light intensity image, represents the projectile 2. Then, the position of the center of gravity of the point image is calculated. The threshold value is a value larger than 0 and smaller than the maximum value of the pixel value, and is stored in the internal memory of the angle calculation unit 53, for example.
When the size of the first light collector 11 is finite, the point image has a spread beyond the diffraction limit. In reality, the point image captured when the projectile 2 is small is further expanded because it is affected by the fluctuation of the atmosphere or seeing.
The angle calculation unit 53 calculates the position of the center of gravity of the point image in order to obtain the position of the point image having a spread as the position of light reflection with respect to the projectile 2.
As described above, the position on the light intensity image can be represented by an angle such as an azimuth angle and an elevation angle. Therefore, the angle calculation unit 53 expresses the position of the center of gravity of the calculated point image by an angle. Find the angle of 1.
The first angle is the angle at which the flying object 2 is viewed from the first condensing device 11, the azimuth angle Az _A when the flying object 2 is viewed from the first condensing device 11, and the first condensing device. _{It includes the elevation angle Elv A} , which is the view of the projectile 2 from 11.
Angle calculation unit 53, a first angle at a first time _{t 1,} the azimuth angle _{Az A1} at a first time _{t 1,} and the elevation angle _{Elv A1} at a first time _{t 1} to the recording device 59 Have them record.

Angle calculation unit 53, the controller 56 obtains the orientation of the first focusing device 11 at the second time t _2, from the position detection unit 52, obtains the light intensity image at a second time t ₂ do.
The angle calculation unit 53 detects the position of the light in the light intensity image at the _second time t 2 as the reflection position of the light with respect to the projectile 2 at the second time t _2.
The angle calculation unit 53 calculates the first angle at the second time t _{2 from} the directing direction of the first light collecting device 11 at the second time t _{2 and} the reflection position at the second time t ₂ . (Step ST2 in FIG. 6).
Angle calculation unit 53, the removal process of the blackout flow noise contained in the light intensity image at a second time t _2, the process of removing the background light noise included in the light intensity image, or included in the light intensity image It is also possible to carry out attenuation processing of the peripheral light amount and the like, and calculate the first angle from the light intensity image after these processings.
Since the calculation process of the first angle at the second time t ₂ is the same as the calculation process of the first angle at the first time t _1, detailed description thereof will be omitted.
Angle calculation unit 53, a first angle at the second time _{t 2,} the azimuth angle _{Az A2} at the second time _{t 2,} the elevation angle _{Elv A2} at a second time _{t 2} to the recording device 59 Have them record.

Time correcting unit 54, output from the controller 56, the position of the first focusing device 11 at a first time t _1, and the position of the first focusing device 11 at a second time t ₂ get.
Time correcting unit 54, the position of the first focusing device 11 at a first time t _1, and records the position of the first focusing device 11 at a second time t ₂ to the recording device 59.
Here, time correcting unit 54, the recording position of the first focusing device 11 at a first time t _1, and the position of the first focusing device 11 at a second time t ₂ to the recording device 59 I'm letting you. If the first optical angle measuring station 3 is a fixed station, the position of the first light collecting device 11 does not change, so that the position of the first light collecting device 11 at any time is recorded in the recording device 59. Just let me do it.

The control device 76 of the second optical angle measuring station 4 acquires the elapsed time t output from the counter 75.
The control device 76 grasps each of _{the first time t 1} and the second time t ₂ as the observation time of the flying object 2 according to the elapsed time t.
Further, the control unit 76 acquires orbit information output from the communication unit 78, from the track information, and approximate location of the projectile 2 in the first time t _1, the projectile at the second time t ₂ 2 Grasp the approximate position.

Controller 76, the orientation of the second condenser 12 at a first time t ₁ and outputs to each of the directing device 77 and the angle calculation unit 73, the second condensing device at a second time t ₂ The direction of the 12 is output to the pointing device 77 and the angle calculation unit 73, respectively.
The control device 76 controls the light-shielding device 71 so that the shutter is opened for a predetermined exposure time at _{the first time t 1} and the second time t _2.
The control device 76 outputs an image pickup command to the position detection unit 72 at _{the first time t 1} and the second time t _2.
The control device 76 outputs _{the position of the second light collecting device 12 at the first time t 1} and the position of the second light collecting device 12 at the second time t ₂ to the time calibration unit 74.

Directing device 77 of the second optical angle measuring station 4, the directivity direction of the second condenser 12 at a first time t ₁ is the control unit 76 is a first time t ₁ in the second output The rotary stage is rotated by driving the motor so as to match the direction of the light collector 12.
Directing device 77, the directivity direction of the second condenser 12 at a second time t ₂ is consistent with the orientation of the second condenser 12 at a second time t ₂ the controller 76 is output By driving the motor, the rotary stage is rotated.
Since the second condensing device 12 has a function of projecting an object surface onto the image plane and has a reference optical axis, the directing device 77 controls the optical axis of the second condensing device 12. The directing direction of the second light collecting device 12 is controlled by adjusting to the directing direction output from the device 76.
Incidentally, the control of the orientation of the second condenser 12 by directing device 77, it takes a certain amount of time, the control device 76, before the observation time is a first time t _1, a first time the orientation of the second condenser 12 and outputs to the directing device 77 in the t _1. Further, the control device 76 outputs the pointing direction of the second condensing device 12 at the second time t ₂ to the pointing device 77 before the observation time reaches the second time t _2.

In the orbit calculation system shown in FIG. 5, the control device 76 directs the directing direction of the second condensing device 12 to the directing device 77 so that the projectile 2 is within the field of view of the second condensing device 12. I'm instructing. When the field of view of the second light collecting device 12 and the field of view of the position detection unit 72 are different, the control device 76 is the narrower of the field of view of the second light collecting device 12 and the field of view of the position detection unit 72. The directing device 77 is instructed to direct the direction of the second condensing device 12 so that the projectile 2 is within the range of the field of view of.
If the flying object 2 is a geostationary satellite and the second optical measuring station 4 is a fixed station, the relative positions of the second optical measuring station 4 and the flying object 2 do not change. If the relative position does not change, the control device 76 once instructs the directing device 77 in the direction of the second condensing device 12 so that the projectile 2 is within the field of view of the second condensing device 12. Then, after that, the directing direction of the second light collecting device 12 may be fixed. Strictly speaking, the relative position changes slightly due to the difference in the orbital inclination angle, but if the observation time is about several seconds, there is no problem even if the direction of the second light collector 12 is fixed. ..

The second light collecting device 12 of the second optical measuring station 4 collects the light emitted from the illumination light source 1 and then reflected by the projectile 2.
The shading device 71 of the second optical measuring station 4 is controlled by the control device 76 to open the shutter for a predetermined exposure time at _{the first time t 1} and the second time t _2. ..
When the first time t _1, a predetermined exposure time, by opening the shutter of the shading device 71, is condensed light by the second focusing device 12 through the shading device 71, the light It reaches the position detection unit 72.
Further, when the second time t _2, the predetermined exposure time, by opening the shutter of the shading device 71, is condensed light passes through the shading device 71 by the second condenser 12, the The light reaches the position detection unit 72.

Position detector 72 of the second optical angle measuring station 4, when the first time t _1, when acquiring the imaging command output from the controller 76, the detection processing of the light transmitted through the shading device 71 Start and take a light intensity image showing the projectile 2.
Position detector 72, the light intensity image captured, as a light intensity image at a first time t _1, and outputs the angle calculation unit 73.
The position detection unit 72, when the second time t _2, the acquiring the imaging command outputted from the control unit 76 starts the detection process of light transmitted through the shading device 71, the projectile 2 Take an image of the reflected light intensity.
Position detector 72, the light intensity image captured, as a light intensity image at a second time t _2, the output to the angle calculation unit 73.

Angle calculating portion 73 of the second optical angle measuring station 4, the controller 76 obtains the orientation of the second condenser 12 at a first time t _1, the position detector 72, a first for obtaining an optical intensity image at time t _1.
The angle calculation unit 73 detects the position of the light in the light intensity image at the _first time t 1 as the reflection position of the light with respect to the projectile 2 at the first time t _1.
The angle calculation unit 73 calculates _{the second angle at the first time t 1 from} the directing direction of the second light collecting device 12 at the first time t _{1 and} the reflection position at the first time t ₁ . (Step ST3 in FIG. 6).
Angle calculation unit 73, a first removal process of blackout flow noise contained in the light intensity image at time t ₁ of, removal processing of the background light noise included in the light intensity image, or included in the light intensity image A second angle may be calculated from the light intensity images after these treatments by performing attenuation processing of the peripheral light amount and the like.
Since the second angle calculation process by the angle calculation unit 73 is the same as the first angle calculation process by the angle calculation unit 53, detailed description thereof will be omitted.
Angle calculation unit 73, a second angle at a first time _{t 1,} the azimuth angle _{Az B1} at a first time _{t 1,} the communication unit 78 and the elevation angle _{Elv B1} at a first time _{t 1} Output.

The angle calculation unit 73 acquires the pointing direction of the second light collecting device 12 at the second _{time t 2} from the control device 76, and acquires the light intensity image at the _{second time t 2 from the position detection unit 72.} do.
The angle calculation unit 73 detects the position of the light in the light intensity image at the _second time t 2 as the reflection position of the light with respect to the projectile 2 at the second time t _2.
The angle calculation unit 73 determines the second angle at the second time t _{2 from} the directing direction of the second light collecting device 12 at the second time t _{2 and} the reflection position at the second time t ₂ . Calculate (step ST4 in FIG. 6).
Angle calculation unit 73, the removal process of the blackout flow noise contained in the light intensity image at a second time t _2, the process of removing the background light noise included in the light intensity image, or included in the light intensity image A second angle may be calculated from the light intensity images after these treatments by performing attenuation processing of the peripheral light amount and the like.
Angle calculation unit 73, a second angle at the second time _{t 2,} the azimuth angle _{Az B2} at the second time _{t 2,} the communication unit 78 and the elevation angle _{Elv B2} at a second time _{t 2} Output.

Time correcting unit 74, output from the controller 76, and the position of the second condenser 12 at a first time t _1, and the position of the second condenser 12 at a second time t ₂ get.
The time calibration unit 74 outputs _{the position of the second light collecting device 12 at the first time t 1} and the position of the second light collecting device 12 at the second time t ₂ to the communication unit 78.

The communication unit 78 has a _{second angle at the first} time t 1, a second angle at the second time t ₂ , a position of the second light collector 12 at the first time t _{1, and a second time.} transmitting the respective positions of the second focusing device 12 at t ₂ to the communication unit 58.
The communication unit 58, from the communication unit 78, has a _{second angle at the first} time t 1, a second angle at the second time t ₂ , and a position of the second light collector 12 at the first time t _1. and receiving the respective positions of the second condenser 12 at the second time t _2.
The communication unit 58 has a _{second angle at the first} time t 1, a second angle at the second time t ₂ , a position of the second light collector 12 at the first time t _{1, and a second time.} the respective positions of the second condenser 12 is recorded in the recording device 59 in the t _2.
Here, the communication unit 58 causes the recording device 59 to record the position of the second light collecting device 12 at _{the first time t 1} and the position of the second light collecting device 12 at the second time t _2. ing. If the second optical measuring station 4 is a fixed station, the position of the second light collecting device 12 does not change, so that the recording device 59 records the position of the second light collecting device 12 at any time. Just do it.

Trajectory calculation unit 60 from the recording device 59, respectively and each of the first angle at a first time t ₁ and second time t _2, at a first time t ₁ and second time t ₂ the Get the angle of 2.
From the recording device 59, the orbit calculation unit 60 determines _{the position of the first light collecting device 11 at the first time t 1} and the position of the first light collecting device 11 at the second time t ₂ . The position of the second light collecting device 12 at the time t ₁ and the position of the second light collecting device 12 at the second time t ₂ are acquired.

Trajectory calculation unit 60, the position of the first focusing device 11 at a first time _{_{_{t 1 (x A1, y A1}}} , z A1), the position of the second condenser 12 at a first time _{t 1} since _{_{_{(x B1, y B1, z}}} B1) and, as shown in FIG. 12, at a first time _{t 1,} between the devices between the first condenser 11 and the second condenser apparatus 12 a distance _{L AB1} Is calculated.
FIG. 12 is an explanatory diagram showing the accuracy of orbit calculation in the orbit calculation system according to the first embodiment.
In a first time t _1, the first condenser apparatus 11 and the projectile 2 the second condenser 12 is disposed in each of the vertices in the triangle. Therefore, the distance between the devices L _AB1 , the first angle (Az _A1 , Elv _A1 ), the second angle (Az _B1 , Elv _B1 ), and the position of the first light collector 11 (x _A1 , y _A1). , Z _A1 ) or the position of the second light collector 12 (x _B1 , y _B1 , z _B1 ), by using the triangular cosine theorem, the first light collector 11 and the second collection. The distance from each of the optical devices 12 to the projectile 2 can be calculated geometrically.
The orbit calculation unit 60 has a distance L _AB1 between devices, a first angle (Az _A1 , Elv _A1 ), a second angle (Az _B1 , Elv _B1 ), and a position (x) of the first condensing device 11. since _{_A1,} _y _A1, _z _A1) and, using the cosine theorem, as the distance to the projectile 2 from the first condenser 11, calculates a first distance Renge _A1 at a first time _{t 1.}
Further, the orbit calculation unit 60 has a distance L _AB1 between devices, a first angle (Az _A1 , Elv _A1 ), a second angle (Az _B1 , Elv _B1 ), and a position of the second light collector 12. since the _{_{_{(x B1, y B1, z}}} B1), using the cosine theorem, as the distance to the projectile 2 from the second condenser 12, calculates a second distance Renge _B1 at a first time _{t 1} do.

Trajectory calculation unit 60, the position of the first focusing device 11 at a second time _{_{_{t 2 (x A2, y A2}}} , z A2), the position of the second condenser 12 at a second time _{t 2} since _{_{_{(x B2, y B2, z}}} B2) and, at the second time _{t 2, the} calculated first condenser 11 inter-device distance _{L AB2} between the second condenser 12.
If each of the first optical angle measuring station 3 and the second optical angle measuring station 4 is a fixed station, the inter-device distance L _AB2 at the second time t _2, the inter-device distances at a first time t ₁ is the same as L _AB1, it may be omitted for calculating the inter-device distance _{L AB2} at the second time _{t 2.}
The orbit calculation unit 60 includes the inter-device distance L _AB2 , the first angle (Az _A2 , Elv _A2 ), the second angle (Az _B2 , Elv _B2 ), and the position (x) of the first condensing device 11. since _{_A2,} _y _A2, _z _A2) and, using the cosine theorem, as the distance to the projectile 2 from the first condenser 11, calculates a first distance Renge _A2 at the second time _{t 2.}
Further, the orbit calculation unit 60 has a distance between the devices L _AB2 , a first angle (Az _A2 , Elv _A2 ), a second angle (Az _B2 , Elv _B2 ), and a position of the second light collecting device 12. since the _{_{_{(x B2, y B2, z}}} B2), by using the cosine theorem, as the distance to the projectile 2 from the second condenser 12, calculates a second distance Renge _B2 at a second time _{t 2} do.

Trajectory calculation unit 60 includes a first angle at a first time _{_{_{t 1 (Az A1, Elv A1}}} ), the position of the first focusing device 11 at a first time _{_{_{t 1 (x A1, y A1}}} , z and _A1), the first distance Renge _A1 Metropolitan at a first time _{t 1,} the position of the projectile 2 in the first time _{_{_{t 1 (x 1 ', y}}} 1', z 1 ') is calculated.
Further, the trajectory calculating section 23 includes a second angle at a first time _{_{_{t 1 (Az B1, Elv B1}}} ), the position of the second condenser 12 at a first time _{t 1} _(x _{B1, y} B1 calculates a _{z B1),} the second distance Renge _B1 Metropolitan at a first time _{t 1,} the position of the projectile 2 in the first time _{_{_{t 1 (x 1 ", y}}} 1", the _z 1 ") ..
If the position (x ₁ ', y ₁ ', z ₁ ') and the position (x ₁ ", y ₁ ", z ₁ ") of the orbit calculation unit 60 are different, for example, the following equation (1) As shown in, as the _{position (x 1} , y ₁ , z ₁ ) of the projectile 2 at the _first time t 1, the position (x ₁ ', y ₁ ', z ₁ ') and the position (x ₁ ', Calculate the average value with y ₁ ", z _1").
x ₁ = (x ₁ '+ x ₁ ") / 2
y ₁ = (y ₁ '+ y ₁ ") / 2
z ₁ = (z ₁ '+ z ₁ ") / 2 (1)
If the position (x ₁ ', y ₁ ', z ₁ ') and the position (x ₁ ", y ₁ ", z ₁ ") are the same, the _{orbit calculation unit 60 will perform the position (x 1} ', y _1". ', Z ₁ '), or the position (x ₁ ", y ₁ ", z ₁ ") is _{defined as the position (x 1} , y ₁ , z ₁ ) of the projectile 2 at the _{first time t 1.}

Trajectory calculation unit 60 includes a first angle at a second time _{_{_{t 2 (Az A2, Elv A2}}} ), the position of the first focusing device 11 at a second time _{_{_{t 2 (x A2, y A2}}} , z and _A2), the first distance Renge _A2 Metropolitan at the second time _{t 2, the} position of the projectile 2 at the second time _{_{_{t 2 (x 2 ', y}}} 2', z 2 ') is calculated.
Further, the trajectory calculating section 23 includes a second angle at a second time _{_{_{t 2 (Az B2, Elv B2}}} ), the position of the second condenser 12 at a second time _{t 2} _(x _{B2, y} B2 calculates a _{z B2),} the second distance Renge _B2 Metropolitan at the second time _{t 2, the} position of the projectile 2 at the second time _{_{_{t 2 (x 2 ", y}}} 2", the _z 2 ") ..
If the position (x ₂ ', y ₂ ', z ₂ ') and the position (x ₂ ", y ₂ ", z ₂ ") of the orbit calculation unit 60 are different, for example, the following equation (2) As shown in, as the position (x ₂ , y ₂ , z ₂ _{) of the projectile 2 at} the second time t 2, the position (x ₂ ', y ₂ ', z ₂ ') and the position (x ₂ ', Calculate the average value with y ₂ ", z _2").
x ₂ = (x ₂ '+ x ₂ ") / 2
y ₂ = (y ₂ '+ y ₂ ") / 2
z ₂ = (z ₂ '+ z ₂ ") / 2 (2)
If the position (x ₂ ', y ₂ ', z ₂ ') and the position (x ₂ ", y ₂ ", z ₂ ") are the same, the _{orbit calculation unit 60 will perform the position (x 2} ', y _2". ', Z ₂ '), or the position (x ₂ ", y ₂ ", z ₂ ") is the position (x ₂ , y ₂ , z ₂ ) _{of the projectile 2 at} the second time t 2.

The orbit calculation unit 60 determines the position of the projectile 2 at the _first _{time t 1 (x 1} , y ₁ , z ₁ ) and the position of the projectile _{2 at} the second time t _{2 (x 2} , y ₂ , z). ₂ ) and is calculated. Alternatively, the orbit calculation unit 60 calculates the position and velocity of the projectile 2 at the _{first time t1.}
Therefore, the position of the projectile 2 at the _first _{time t 1 (x 1} , y ₁ , z ₁ ) and the position of the projectile _{2 at} the second time t _{2 (x 2} , y ₂ , z ₂ ), or. If the _{position (x 1} , y ₁ , z ₁ ) and velocity ( _{v x 1} , vy ₁ , vz ₁ ) of the projectile 2 at the _first time t 1 are known, the trajectory of the projectile 2 can be calculated. .. The orbit represents the time change of the position of the projectile and the trajectory, and assuming that the projectile 2 is moving mainly according to the gravity of the earth, the motion is Kepler motion. , The locus can be approximated by a quadratic curve. Quadratic curve is eccentricity, long radius, defined by track 6 elements represented by the orbit inclination angle, etc., the parameters of the track 6 elements, the position at time t _1, t ₂ with a projectile (x _1, y _{It is determined from 1} , z ₁ ), (x ₂ , y ₂ , z ₂ ), and the orbit is obtained. However, the orbit of a moving object near the earth changes from moment to moment, and the six elements of the orbit are insufficient to represent the orbit, so the position (x ₁ , y ₁ , z ₁ ) and velocity ( _{v x 1} , vy ₁ , Sequentially finding vz ₁ ) may be regarded as finding the orbit.
The orbit calculation unit 60 determines the position of the projectile 2 at the _first _{time t 1 (x 1} , y ₁ , z ₁ ) and the position of the projectile _{2 at} the second time t _{2 (x 2} , y ₂ , z). ₂ ) Or, the orbit of the projectile 2 is calculated from _{the position (x 1} , y ₁ , z ₁ ) and the velocity ( _{v x 1} , v y ₁ , vz ₁ ) of the projectile 2 at the _{first time t 1.} (Step ST5 in FIG. 6).

Hereinafter, in the orbit calculation system shown in FIG. 5, the reason why the time is synchronized between the first optical angle measuring station 3 and the second optical angle measuring station 4 will be described.
The angle calculation unit 53 of the first optical angle measuring station 3 can calculate the first angle from the directing direction of the first light collecting device 11 and the light intensity image, but the first light collecting device It is not possible to calculate the distance from the first condensing device 11 to the projectile 2 with high accuracy only by the information of the directing direction of 11 and the light intensity image.
The angle calculation unit 73 of the second optical angle measuring station 4 can calculate the second angle from the directing direction of the second condensing device 12 and the light intensity image, but the second condensing device It is not possible to calculate the distance from the second condensing device 12 to the projectile 2 with high accuracy only by the information of the directing direction of 12 and the light intensity image.
FIG. 7 is an explanatory diagram showing the accuracy of distance calculation by the angle calculation unit 53 of the first optical angle measuring station 3.
In FIG. 7, the shaded ellipse region shows the distance calculation accuracy and the angle calculation accuracy, respectively. Since the shaded elliptical region has an elliptical shape that is long in the range direction, it indicates that the calculation accuracy of the distance is poor.

As shown in FIG. 7, even when the distance from the first condensing device 11 to the projectile 2 cannot be calculated with high accuracy, as shown in FIG. 8, at three observation times t ₀ , t ₁ , and t ₂ . If the angle calculation unit 53 calculates the first angle, the trajectory of the projectile 2 can be calculated.
FIG. 8 is an explanatory diagram showing the calculation accuracy of each angle and the calculation accuracy of each distance at the _{three observation times t 0} , t ₁ , and t _2.
With respect to the first optical angle measuring station 3, since the flying body 2 is relatively moved, the angle calculation unit 53, the first angle at observation time _{_{_{t 0 (Az A0, Elv A0}}} ), observation time can be calculated respective first angle in _{_{_{t 1 (Az A1, Elv A1}}} ) and a first angle at observation time _{_{_{t 2 (Az A2, Elv A2}}} ).
Assuming that the projectile 2 is moving mainly according to the gravity of the earth, the orbit of the projectile 2 becomes a quadratic curve due to the Kepler motion. Even if the distance to 2 is not known, the distance to the projectile 2 can be estimated from the first angle of each of the _{observation times t 0} , t ₁ , and t _{2 and the time when the projectile 2 moves at each point.} Therefore, the orbit of the projectile 2 is uniquely determined. Actually, it may contain an error because it is affected by the gravity anomaly of the earth or radiation pressure, but in principle, if there are three observations of the projectile 2, the projectile 2 The orbit can be calculated.

In FIG. 9, the angle measuring station installed on the ground emits radio waves, and the angle measuring station receives the radio waves reflected by the flying object 2, thereby calculating the distance from the angle measuring station to the flying object 2. It is explanatory drawing which shows the measurement accuracy of a case.
When the angle measuring station emits radio waves and the angle measuring station receives the radio waves reflected by the projectile 2, the distance from the angle measuring station to the projectile 2 is calculated. The distance is calculated from the time difference from the time or the phase difference between the phase of the emitted radio wave and the phase of the received radio wave. Therefore, the angle measuring station can calculate the angle (Az, Elv) of the projectile 2 with respect to the angle measuring station in addition to the distance.
In FIG. 9, the shaded ellipse region shows the distance calculation accuracy and the angle calculation accuracy, respectively.
Since the shaded ellipse region shown in FIG. 9 has a shorter length in the range direction than the shaded ellipse region shown in FIG. 7, the calculation accuracy of the distance in the angle measuring station shown in FIG. 9 is the calculation of the distance in FIG. It is higher than the accuracy. However, since the wavelength (cm) of the radio wave is longer than the wavelength (micron) of visible light, for example, the angle calculation accuracy in the angle measuring station shown in FIG. 9 is worse than the angle calculation accuracy in FIG. There is.

When both the distance and the angle can be obtained, as shown in FIG. 10, if both the distance and the angle are calculated at the _{two observation times t 1} and t _{2, the trajectory of the projectile 2 can be calculated.}
FIG. 10 is an explanatory diagram showing the calculation accuracy of each angle and the calculation accuracy of each distance at the _{two observation times t 1} and t _2.
From the respective distances and angles at the _two observation times t ₁ and t 2, the position of the projectile 2 at the _{observation time t 1} _{(x 1} , y ₁ , z ₁ ) and the position of the projectile _{2 at} the observation time t 2. (X ₂ , y ₂ , z ₂ ) can be obtained.
Alternatively, the speed of the projectile 2 from the position of the projectile 2 at the observation time _{_{_{t 1 (x 1, y 1}}} , z 1) at the observation time _{_{_{t 1 (vx 1, vy 1}}} , vz 1) is obtained. Therefore, if both the distance and the angle can be calculated at the two observation times t ₁ and t _{2, the orbit of the projectile 2 can be calculated.}

FIG. 11 is an explanatory diagram showing the accuracy of angle calculation and the accuracy of distance calculation in the trajectory calculation system according to the first embodiment.
In the orbit calculation system according to the first embodiment, the first optical angle measuring station 3 and the second optical measuring station 4, which are two optical measuring stations having different installation positions, fly at the same time. The angle with respect to the body 2 is calculated.
Therefore, the first condensing device 11 in the first optical angle measuring station 3, the second condensing device 12 in the second optical measuring station 4, and the projectile 2 have their respective vertices in the triangle. The triangular cosine theorem can be applied as if it were placed in. By using the triangular cosine theorem, the distance from each of the first condensing device 11 and the second condensing device 12 to the projectile 2 can be calculated geometrically.
Since the triangular cosine theorem cannot be applied unless the first optical measuring station 3 and the second optical measuring station 4 calculate the angle with respect to the same projectile 2 at the same time, FIG. 7 As shown by the shaded elliptical region of, the calculation accuracy of the distance deteriorates.
In FIG. 11, the diamond-shaped region including the projectile 2 shows the accuracy of calculating the distance and the accuracy of calculating the angle, respectively. The diamond-shaped region has a shorter length in the range direction than the shaded elliptical region shown in FIG. Therefore, it can be seen that the distance calculation accuracy in the trajectory calculation system shown in FIG. 5 is higher than that shown in FIG. 7.
From the above, the reason why the time is synchronized between the first optical measuring station 3 and the second optical measuring station 4 is to improve the accuracy of distance calculation. This is also to enable the orbit of the projectile 2 to be calculated even if the observation time is two.

In the above embodiment 1, the directing at the first time and the second time of the first condensing device 11 that condenses the light emitted by the illumination light source 1 and then reflected by the projectile 2. The first time as the first angle when the light collector 2 is viewed from the first condensing device 11 from the direction and the respective reflection positions of the light with respect to the projectile 2 at the first time and the second time. The first time and the second of the first angle calculation unit 21 that calculates the respective first angles at the second time and the second light collector 12 that collects the light reflected by the projectile 2. A second view of the flying object 2 from the second condensing device 12 from the respective pointing directions at the time of 2 and the respective reflection positions of the light with respect to the flying object 2 at the first time and the second time. As angles, a second angle calculation unit 22 that calculates each second angle at the first time and a second time, a first angle calculated by the first angle calculation unit 21, and a second angle. From the second angle calculated by the angle calculation unit 22 of the above, and the position of the first light collector 11 and the position of the second light collector 12 at the first time and the second time, the flying object. The orbit calculation device 13 is configured to include the orbit calculation unit 23 for calculating the orbit of 2. Therefore, it may be possible to calculate the orbit of the distant projectile 2 which cannot be calculated by the conventional orbit calculation system.

The orbit calculation system shown in FIG. 5 includes a first optical angle measuring station 3 and a second optical measuring station 4 as two optical measuring stations. However, depending on the weather, for example, even if the first optical measuring station 3 can collect the light reflected by the projectile 2, the second optical measuring station 4 is reflected by the flying object 2. It may cause a situation where light cannot be focused. In such a situation, the orbit calculation system shown in FIG. 5 cannot calculate the orbit of the projectile 2. Even in such a situation, in order to be able to calculate the orbit of the projectile 2, a third optical angle measuring station is provided in addition to the first optical measuring station 3 and the second optical measuring station 4. You may. The configuration of the third optical angle measuring station is the same as the configuration of the second optical measuring station 4, and the installation position of the third optical measuring station is the first optical measuring station 3 and the second optical measuring station. It is different from each installation position in the optical angle measuring station 4.

Embodiment 2.
In the second embodiment, an orbit calculation system in which the first optical angle measuring station 3 is a fixed station and the second optical angle measuring station 4 is a portable station will be described.

FIG. 13 is an explanatory diagram showing an example of the positional relationship between the illumination light source 1, the projectile 2, the first optical angle measuring station 3 and the second optical measuring station 4 in the orbit calculation system according to the second embodiment. ..
In the orbit calculation system shown in FIG. 13, the first optical angle measuring station 3 is a fixed station, and the second optical angle measuring station 4 is a portable station.
FIG. 14 is an explanatory diagram showing the accuracy of angle calculation and the accuracy of distance calculation in the trajectory calculation system according to the second embodiment.

Even in the orbit calculation system shown in FIG. 13, if the position of the first optical angle measuring station 3 and the position of the second optical measuring station 4 are different, the same principle as that of the orbit calculation system shown in FIG. 5 can be used. The trajectory of the projectile 2 can be calculated.
However, when the position of the sun, which is the illumination light source 1, or the position of the second optical angle measuring station 4, which is a portable station, changes, the first optical measuring station 3 and the second optical measuring station 3 are changed. 4 may be an observation condition in which light from the same projectile 2 cannot be received at the same time. For example, when the second optical measuring station 4 is located on the back side of the earth when viewed from the first optical measuring station 3, the first optical measuring station 3 and the second optical measuring station 3 are measured. The corner station 4 cannot receive the light from the same projectile 2 at the same time.
As long as the first optical angle measuring station 3 and the second optical measuring station 4 can receive the light from the same projectile 2 at the same time, the orbit calculation system shown in FIG. 13 can be shown in FIG. Similar to the trajectory calculation system shown, the trajectory of the projectile 2 can be calculated.

In the orbit calculation system shown in FIG. 13, similarly to the orbit calculation system shown in FIG. 5, the communication unit 78 of the second optical angle measuring station 4 determines the second angle and the position of the second condensing device 12. , Is transmitted to the communication unit 58 of the first optical angle measuring station 3.
However, since the second optical measuring station 4 is moving in orbit, for example, the communication unit 78 may be used depending on the positional relationship between the first optical measuring station 3 and the second optical measuring station 4. , The second angle and the position of the second light collecting device 12 may not be directly transmitted to the communication unit 58 of the first optical measuring station 3.
In such a case, the communication unit 78 of the second optical angle measuring station 4 transmits the second angle and the position of the second light collecting device 12 to the ground station 6 shown in FIGS. 13 and 14. The ground station 6 may transfer the second angle and the position of the second light collecting device 12 to the communication unit 58 of the first optical measuring station 3.
The ground station 6 is a ground station that operates the second optical angle measuring station 4, and controls, for example, the posture and position of the second optical measuring station 4.

Embodiment 3.
In the third embodiment, an orbit calculation system in which the first optical angle measuring station 3 is provided with the ranging device 91 and the second optical angle measuring station 4 is provided with the retroreflective unit 81 will be described.

FIG. 15 is an explanatory diagram showing an example of the positional relationship between the flying object 2, the first optical angle measuring station 3, and the second optical angle measuring station 4 in the orbit calculation system according to the third embodiment.
The second optical angle measuring station 4 includes a retroreflective unit 81.
The retroreflective unit 81 is realized by an optical member having a function of emitting light in the incident direction. As the optical member, a corner cube reflector or the like can be considered in addition to a member in which beads are spread.
The retroreflective unit 81 retroreflects the light radiated from the ranging device 91 of the first optical angle measuring station 3.

FIG. 16 is a specific configuration diagram showing the trajectory calculation system according to the third embodiment. In FIG. 16, the same reference numerals as those in FIG. 5 indicate the same or corresponding parts, and thus the description thereof will be omitted.
The distance measuring device 91 is mounted on the first optical angle measuring station 3.
The distance measuring device 91 includes a light collecting device 92, an optical transmitting unit 93, an optical receiving unit 94, a distance calculation unit 95, and a directing device 96.
The distance measuring device 91 receives the light reflected by the retroreflective unit 81 after radiating the light toward the retroreflective unit 81, and from the emitted light and the received light, the first optical angle measuring station 3 And the second optical measuring station 4 are calculated.

The condensing device 92 is realized by, for example, one or more telescopes among a refraction type telescope and a reflection type telescope.
The light collecting device 92 radiates the laser light output from the light transmitting unit 93 toward the retroreflective unit 81, and collects the laser light reflected by the retroreflective unit 81.
The optical transmission unit 93 is realized by, for example, a laser light source, an optical amplifier, and an optical modulator.
The light transmitting unit 93 outputs the laser light to the condensing device 92, so that the laser light is radiated from the condensing device 92 toward the retroreflective unit 81.
The optical receiver 94 is realized by, for example, a semiconductor optical receiver. As the semiconductor optical receiver, a photodiode, a quadrant photodiode, or the like can be considered.
The light receiving unit 94 receives the laser light collected by the condensing device 92.

The distance calculation unit 95 is based on the time difference between the time when the radio wave is radiated from the optical transmission unit 93 and the time when the radio wave is received by the optical reception unit 94, or the phase of the radio wave output from the optical transmission unit 93 and the optical reception unit 94. From the phase difference with the phase of the received radio wave, the distance between the devices, which is the distance between the first optical measuring station 3 and the second optical measuring station 4, is calculated.
That is, the distance calculation unit 95, the inter-device distance _{L AB1} at a first time _{t 1} is calculated and the inter-device distance _{L AB2} at the second time _{t 2, the} a unit distance _{L AB1} and apparatus distance _{L AB2} Is output to the orbit calculation unit 61 of the orbit calculation unit 23.

The directional device 96 is realized by a rotary stage and a motor having one or more rotation axes. An altazimuth mount, an equatorial mount, or the like can be considered as a rotation stage having one or more rotation axes.
A condensing device 92 is mounted on the rotating stage of the directional device 96.
Directing device 96, so that the directivity direction of the light collection device 92 when the first time t ₁ coincides with the orientation direction of the condenser 92 when the first time t ₁ the controller 56 is outputted In addition, the rotary stage is rotated by driving the motor.
Directing device 96, so that the directivity direction of the light collection device 92 when the second time t ₂ coincides with the orientation direction of the condenser 92 when the second time t ₂ the controller 56 is outputted In addition, the rotary stage is rotated by driving the motor.

The orbit calculation unit 23 includes a recording device 59 and an orbit calculation unit 61.
The orbit calculation unit 61, from the distance calculation unit 95, indicates _{the direction directions of the light collector 92 at the first time t 1} and the second time t ₂ , and the first time t ₁ and the second time t _2. _{The distances L AB1} and L _AB2 between the devices in the above are acquired.
Trajectory calculation unit 61 from the recording apparatus 59 acquires the first position at time t ₁ and second time t ₂ of the first focusing device 11.
The orbit calculation unit 61 has the direction directions of the light collector 92 at the first time t ₁ and the second time t ₂ _{, and the first time t 1} and the second time of the first light collector 11. Using _{the respective positions at t 2} and the distances between the devices _{LA AB1} and _{LA AB2 at} _{the first time t 1} and the second time t ₂ , the first time t of the second light collector 12 is used. and it calculates the respective positions in the _first and second time t _2.

Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 has a first angle and a first _{time t at the first} time t 1 and the second time t _{2 from the recording device 59.} obtaining a respective second angle in the _first and second time t _2.
Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 has a first angle, a second angle, and the first condensing device 11 and the second condensing device 12. from the respective positions, as the distance to the projectile 2 from the first condenser 11, to calculate a respective first distance at a first time t ₁ and second time t _2.
Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 is a flying object from the second light collector 12 from each first angle, each second angle, and each position. as the distance to 2, and calculates the respective second distance at a first time t ₁ and second time t _2.
Trajectory calculation unit 61, like the trajectory calculation unit 60 shown in FIG. 5, first, respectively at a first time t ₁ and second time t ₂ the angle a, the first time t ₁ and the second a respective second angle at time t _2, the respective a first distance at a first time t ₁ and second time t _2, the respective at a first time t ₁ and second time t ₂ From the second distance, each of the position and speed of the projectile 2 is calculated. That is, the orbit calculation unit 61 calculates the respective positions _{of the projectile 2 at the first time t 1} and the second time t ₂ , or the position and velocity of the projectile 2 at the _{first time t 1.} ..
Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 is _{located at each position of the projectile 2 at the first time t 1} and the second time t ₂ , or the first time of the projectile 2. The orbit of the projectile 2 is calculated from the position and velocity of t _1.

Next, the operation of the trajectory calculation system shown in FIG. 16 will be described.
In the orbit calculation system shown in FIG. 5, the orbit calculation unit 60 uses the first light collector 11 and the second light collector 12 from their respective positions in the first light collector 11 and the second light collector 12. _{The distances L AB1} and L _AB2 between the devices with 12 are calculated.
The trajectory calculation system shown in FIG. 16, the distance measuring device 91, as the distance between the distance _L _{AB1, L AB2} between the first condenser 11 and the second condenser 12, a first optical angle measuring station 3 It differs from the orbit calculation system shown in FIG. 5 in that the distance between the light and the second optical measuring station 4 is calculated.
The accuracy of distance calculation by the distance measuring device 91 is higher than the accuracy of distance calculation by the trajectory calculation system shown in FIG. The accuracy of distance calculation by the trajectory calculation system shown in FIG. 5 is the length in the range direction in the diamond-shaped region shown in FIG.
However, since the intensity of the laser beam becomes weaker in inverse proportion to the square of the propagation distance of the laser beam, when the distance between the first optical angle measuring station 3 and the second optical measuring station 4 is long, the intensity of the laser light becomes weaker. The distance measuring device 91 may not be able to calculate the _{inter-device distances LA AB1} and _{LA AB2.} In the orbit calculation system shown in FIG. 16, since the second optical measuring station 4 includes the retroreflecting unit 81, the attenuation of the laser beam intensity due to the reflection of the second optical measuring station 4 is suppressed. ing. Therefore, in the orbit calculation system shown in FIG. 16, the distance that can be calculated is longer than that in the second optical measuring station 4 that does not have the retroreflection unit 81.

The control device 56 _{outputs the pointing direction of the condensing device 92 at the first} time t1 to the directional device 96, and outputs the directional direction of the condensing device 92 at the second time t ₂ to the directional device 96.
The control device 56 outputs a laser beam radiation command to each of the light transmission unit 93 and the distance calculation unit 95 at _{the first time t 1} and the second time t _2.
For example, the orbit information of the second optical measuring station 4 is recorded in the internal memory of the control device 56. Controller 56, based on the orbital information to identify the orientation of the condenser 92 at a first time t _1, to identify the orientation of the condenser 92 at the second time t _2.

Directing device 96, so that the directivity direction of the light collection device 92 when the first time t ₁ coincides with the orientation direction of the condenser 92 when the first time t ₁ the controller 56 is outputted In addition, the rotary stage is rotated by driving the motor.
Directing device 96, so that the directivity direction of the light collection device 92 when the second time t ₂ coincides with the orientation direction of the condenser 92 when the second time t ₂ the controller 56 is outputted In addition, the rotary stage is rotated by driving the motor.
Since it takes a certain amount of time to control the directivity direction of the condensing device 92 by the directivity device 96, the control device 56 sets the observation time at the first time t ₁ before the observation time reaches the first time t ₁ . The directivity direction of the light collector 92 is output to the directivity device 96. Further, the control device 56 outputs the directivity direction of the condensing device 92 at the _second time t 2 to the directivity device 96 before the observation time reaches the second time t _2.

When the optical transmission unit 93 _{receives a laser light emission command from the control device 56 at the first time t 1} and the second time t ₂ , the light transmission unit 93 outputs the laser light to the condensing device 92 by outputting the laser light to the condensing device 92. Laser light is radiated from the light collector 92 toward the retroreflecting unit 81.
The light collecting device 92 radiates the laser light output from the light transmitting unit 93 toward the retroreflecting unit 81.
When the first time t _1, the laser light emitted from the light collector 92 is reflected by the retroreflective portion 81 attached to the second optical angle measuring station 4.
The laser light reflected by the retroreflective unit 81 is collected by the condensing device 92.
The light receiving unit 94 receives the laser light focused by the condensing device 92, and outputs the received signal of the laser light to the distance calculation unit 95.

The distance calculation unit 95 is based on the time difference between the time when the radio wave is radiated from the optical transmission unit 93 and the time when the radio wave is received by the optical reception unit 94, or the phase of the radio wave output from the optical transmission unit 93 and the optical reception unit 94. The distance between the first optical measuring station 3 and the second optical measuring station 4 is calculated from the phase difference with the phase of the received radio wave.
That is, the distance calculation unit 95, from the first time t ₁ is the time when the emission command is output from the controller 56, based on the time until receiving the reception signal of the laser beam from the optical receiving unit 94, the first calculating the inter-device distance _{L AB1} at time _{t 1} of the.
The distance calculation unit 95, second from the time t ₂ is the time at which the emission command is output from the controller 56, based on the time until receiving the reception signal of the laser beam from the optical receiving unit 94, the second calculating the inter-device distance _{L AB2} at time _{t 2} of the.
Distance calculation unit 95, the pointing direction of the condenser 92 at a first time _{t 1,} the pointing direction of the condenser 92 at the second time _{t 2, the} inter-device distance _{L AB1} at a first time _{t 1} and the Each of the inter-device distances _LAB2 at time t2 of ₂ is output to the orbit calculation unit 61 of the orbit calculation unit 23.

The orbit calculation unit 61, from the distance calculation unit 95, indicates _{the direction directions of the light collector 92 at the first time t 1} and the second time t ₂ , and the first time t ₁ and the second time t _2. _{The distances L AB1} and L _AB2 between the devices in the above are acquired.
Trajectory calculation unit 61 from the recording apparatus 59 acquires the first respective positions at the time t ₁ and second time t ₂ of the first focusing device 11.
Trajectory calculation unit 61, the orientation of the condenser 92 at a first time t _1, the position of the first focusing device 11 at a first time t _1, the inter-device distances at a first time t ₁ by using the L _AB1, it calculates the position of the second condenser 12 at a first time _{t 1.}
Trajectory calculation unit 61, the orientation of the condenser 92 at the second time t _2, the position of the first focusing device 11 at the second time t _2, the inter-device distance at a second time t ₂ by using the L _AB2, it calculates the position of the second condenser 12 at the second time _{t 2.}

Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 has a _{first angle at the first time t 1} and a second time t ₂ , respectively, and the first time t ₁ and the second time t 1. obtaining a respective second angle at time t _2.
Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 has a first angle, a second angle, a position in the first condensing device 11, and a second angle. from the respective positions of the light collection device 12, as the distance to the projectile 2 from the first condenser 11, to calculate a respective first distance at a first time t ₁ and second time t ₂ ..
Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 has a first angle, a second angle, a position in the first condensing device 11, and a second angle. from the respective positions of the light collection device 12, as the distance to the projectile 2 from the second condenser 12, to calculate the respective second distance at a first time t ₁ and second time t ₂ ..

Trajectory calculation unit 61, like the trajectory calculation unit 60 shown in FIG. 5, first, respectively at a first time t ₁ and second time t ₂ the angle a, the first time t ₁ and the second a respective second angle at time t _2, the respective a first distance at a first time t ₁ and second time t _2, the respective at a first time t ₁ and second time t ₂ From the second distance, each of the position and speed of the projectile 2 is calculated. That is, the orbit calculation unit 61 calculates the respective positions _{of the projectile 2 at the first time t 1} and the second time t ₂ , or the position and velocity of the projectile 2 at the _{first time t 1.} ..
Similar to the orbit calculation unit 60 shown in FIG. 5, the orbit calculation unit 61 is _{located at each position of the projectile 2 at the first time t 1} and the second time t ₂ , or the first time of the projectile 2. from the position and velocity of the t _1, to calculate the trajectory of the projectile 2.

In the above embodiment 3, the second optical angle measuring station 4 includes a retroreflecting unit 81 that retroreflects light, and the first optical angle measuring station 3 radiates light toward the retroreflecting unit 81. After that, the light reflected by the retroreflection unit 81 is received, and the distance between the first optical angle measuring station 3 and the second optical measuring station 4 is calculated from the emitted light and the received light. The orbit calculation system shown in FIG. 16 is provided so that the orbit calculation unit 23 calculates the position of the second condensing device 12 by using the distance calculated by the distance measuring device 91. Configured. Therefore, the orbit calculation system shown in FIG. 16 may be able to calculate the orbit of the distant projectile 2 which cannot be calculated by the conventional orbit calculation system, similar to the orbit calculation system shown in FIG. Further, the trajectory calculation system shown in FIG. 16 has improved accuracy of calculating the trajectory of the projectile 2 as compared with the trajectory calculation system shown in FIG.

Embodiment 4.
In the fourth embodiment, the orbit calculation system in which the second optical angle measuring station 4 is arranged in an orbit having an altitude different from that of the orbit of the projectile 2 will be described.

FIG. 17 is an explanatory diagram showing an example of the positional relationship between the flying object 2, the first optical angle measuring station 3, and the second optical angle measuring station 4 in the orbit calculation system according to the fourth embodiment.
17, 101 shows a trajectory of the second optical angle measuring station 4, _{h 1} represents the height of the track 101.
102 shows the orbit of the flying object 2, and h ₂ shows the altitude of the orbit 102.
103 shows a trajectory of the second optical angle measuring station 4, _{h 3} shows a high degree of track 103.
The second optical angle measuring station 4 may be arranged in an orbit 101 having an altitude lower than the orbit 102 of the projectile 2, and may be arranged in an orbit 103 having an altitude higher than the orbit 102 of the projectile 2. Sometimes.
Further, the orbit of the second optical measuring station 4 may be switched from the orbit 101 to the orbit 103, while the orbit 103 may be switched to the orbit 101.

The orbit around the earth is determined by the gravity and altitude of the earth. Speed of the second optical angle measuring station 4 which moves the trajectory v _{i (i} = 1, 3) is expressed by the following equation (3). v ₁ is the speed of the second optical measuring station 4 moving in the orbit 101 at the altitude h ₁ _{, and v 3} is the speed of the second optical measuring station 4 moving in the orbit 103 at the altitude h _3. v ₂ is the velocity of the projectile 2.

In equation (3), r _E is the radius of the earth, G is the gravitational constant, and M is the mass of the earth.

Second optical angle measuring station 4 speed _v i (i = 1, 3), as is apparent from equation (3), as the altitude _{h i} is high, slower.
Therefore, when the second optical measuring station 4 is arranged in the orbit 101 having an altitude higher than the orbit 102 of the flying object 2, the second optical measuring station 4 becomes the flying object 2 with the passage of time. You may be overtaken. On the other hand, when the second optical measuring station 4 is arranged in the orbit 101 whose altitude is lower than the orbit 102 of the flying object 2, the second optical measuring station 4 sets the flying object 2 with the passage of time. You may overtake.
The temporal change in the geometrical arrangement of the first optical measuring station 3, the second optical measuring station 4, and the projectile 2 is the altitude of the orbit where the second optical measuring station 4 is arranged. Depends on.

The orbit calculation system according to the fourth embodiment has, for example, a second optical angle measuring station 4 arranged in an orbit 101 having an altitude lower than the orbit 102 of the projectile 2, and an altitude higher than the orbit 102 of the projectile 2. May include both with a second optical measuring station 4 located in a high orbit 101.
If the orbit calculation system includes two second optical measuring stations 4, the geometric arrangement related to the second optical measuring station 4 arranged in the orbit 101 and the arrangement in the orbit 103 are arranged. The orbit 102 of the projectile 2 can be calculated based on any of the geometrical arrangements related to the second optical angle measuring station 4. In one geometric arrangement, even if the first optical measuring station 3 and the second optical measuring station 4 are under observation conditions that the light from the same projectile 2 cannot be received at the same time, the other. In the geometric arrangement of, there is a possibility of observation conditions that the first optical measuring station 3 and the second optical measuring station 4 can receive the light from the same projectile 2 at the same time. Therefore, in the orbit calculation system according to the fourth embodiment, it is possible to shorten the time during which the orbit of the projectile 2 cannot be calculated.

In the present disclosure, any combination of the embodiments can be freely combined, any component of the embodiment can be modified, or any component can be omitted in each embodiment.

This disclosure is suitable for an orbit calculation device, an orbit calculation method, and an orbit calculation system for calculating the orbit of a flying object.

1 Illumination light source, 2 Flying object, 3 1st light measuring station, 4 2nd light measuring station, 5 Reference signal source, 6 Ground station, 11 1st light collecting device, 12 2nd light collecting device , 13 orbit calculation device, 21 first angle calculation unit, 22 second angle calculation unit, 23 orbit calculation unit, 31 first angle calculation circuit, 32 second angle calculation circuit, 33 orbit calculation circuit, 41 memory , 42 processor, 51 shading device, 52 position detection unit, 53 angle calculation unit, 54 time calibration unit, 55 counter, 56 control device, 57 pointing device, 58 communication unit, 59 recording device, 60 orbit calculation unit, 61 orbit calculation Unit, 71 shading device, 72 position detection unit, 73 angle calculation unit, 74 time calibration unit, 75 counter, 76 control device, 77 directional device, 78 communication unit, 80 orbit database, 81 retroreflection unit, 91 distance measuring device, 92 condensing device, 93 optical transmitter, 94 optical receiver, 95 distance calculation unit, 96 directional device, 101, 102, 103 orbits.

Claims

The direction of each of the first time and the second time of the first condensing device that condenses the light reflected by the projectile after being emitted from the illumination light source, and the first of the light to the projectile. From the time of 1 and the respective reflection positions at the time of the second, the first angle of the flying object seen from the first condensing device is the first time and the second time, respectively. The first angle calculation unit that calculates the first angle of
The direction of each of the first time and the second time of the second light collecting device that collects the light reflected by the flying object, the first time of the light with respect to the flying object, and the above. From the respective reflection positions at the second time, as the second angle when the flying object is viewed from the second condensing device, the second angle at the first time and the second angle at the second time. The second angle calculation unit that calculates
The first angle calculated by the first angle calculation unit, the second angle calculated by the second angle calculation unit, and the first time at the first time and the second time. An orbit calculation device including an orbit calculation unit that calculates the orbit of the flying object from the position of the light condensing device and the position of the second light condensing device.
The first angle calculation unit is
As the reflection position of the light with respect to the flying object at the first time, the position of the light in the image obtained by capturing the light collected by the first condensing device at the first time is detected. death,
As the reflection position of the light with respect to the flying object at the second time, the position of the light in the image obtained by capturing the light collected by the first condensing device at the second time is detected. death,
The first angle at the first time is calculated from the directivity direction of the first light collector at the first time and the reflection position at the first time.
The claim is characterized in that the first angle at the second time is calculated from the directing direction of the first light collecting device at the second time and the reflection position at the second time. 1. The orbit calculation device according to 1.
The second angle calculation unit is
As the reflection position of the light with respect to the flying object at the first time, the position of the light in the image obtained by capturing the light collected by the second condensing device at the first time is detected. death,
As the reflection position of the light with respect to the flying object at the second time, the position of the light in the image obtained by capturing the light collected by the second condensing device at the second time is detected. death,
The second angle at the first time is calculated from the directivity direction of the second light collector at the first time and the reflection position at the first time.
The claim is characterized in that the second angle at the second time is calculated from the directing direction of the second light collecting device at the second time and the reflection position at the second time. 1. The orbit calculation device according to 1.
The orbit calculation unit
From the respective positions of the first light collector at the first time and the second time, and from the respective positions of the second light collector at the first time and the second time. The distance between the devices, which is the distance between the first light collecting device and the second light collecting device at the first time, and the distance between the devices at the second time are calculated. death,
The distance between the devices at the first time and the second time, the first angle at the first time and the second time, the first time and the second time. The first from the first condensing device to the projectile from the respective second angles in the above and the respective positions of the first condensing device at the first time and the second time. As the distance of, the first distances at the first time and the second time are calculated.
The distance between the devices at the first time and the second time, the first angle at the first time and the second time, the first time and the second time. From the respective second angle in the above and the respective positions of the second light collector at the first time and the second time, the second light collector to the projectile. As the distance of, the second distances at the first time and the second time are calculated.
The first angle at the first time and the second time, the second angle at the first time and the second time, the first time and the second time. From the respective first distances at the time and the respective second distances at the first time and the second time, the respective positions of the projectile at the first time and the second time. Is calculated,
The orbit calculation device according to claim 1, wherein the orbit of the flying object is calculated from the respective positions of the flying object at the first time and the second time.
The first angle calculation unit collects the light reflected by the projectile after being radiated from the illumination light source. The first time as the first angle when the flying object is viewed from the first light collecting device from the respective reflection positions of the light with respect to the flying object at the first time and the second time. And each first angle at the second time is calculated.
The second angle calculation unit collects the light reflected by the flying object, the respective directing directions at the first time and the second time of the second condensing device, and the light with respect to the flying object. From the first time and the reflection position at the second time, the first time and the second time are used as the second angle when the flying object is viewed from the second light collecting device. Calculate each second angle at time and
The orbit calculation unit has a first angle calculated by the first angle calculation unit, a second angle calculated by the second angle calculation unit, the first time, and the second time. A method for calculating an orbit in which the orbit of the flying object is calculated from the position of the first condensing device and the position of the second condensing device.
A first condensing device that condenses the light that is radiated from the illumination light source and then reflected by the projectile,
From the respective pointing directions of the first light collector at the first time and the second time, and the reflection positions of the light with respect to the flying object at the first time and the second time, the said As a first angle when the flying object is viewed from the first condensing device, a first angle calculation unit for calculating the first angle at the first time and the second time, respectively.
A second condensing device that condenses the light that is radiated from the illumination light source and then reflected by the projectile,
From the respective pointing directions of the second light collector at the first time and the second time, and the reflection positions of the light with respect to the flying object at the first time and the second time. As a second angle when the flying object is viewed from the second condensing device, a second angle calculation unit for calculating the second angles at the first time and the second time, respectively.
The first angle calculated by the first angle calculation unit, the second angle calculated by the second angle calculation unit, and the first time at the first time and the second time. An orbit calculation system including an orbit calculation unit that calculates the orbit of the flying object from the position of the light collector and the position of the second light collector.
The first condensing device, the first angle calculation unit, and the orbit calculation unit are mounted on the first optical angle measuring station.
The orbital calculation system according to claim 6, wherein the second light collector and the second angle calculation unit are mounted on a second optical angle measuring station.
The first optical angle measuring station is a fixed station whose position is fixed.
The orbit calculation system according to claim 7, wherein the second optical angle measuring station is a portable station whose position changes.
The second optical measuring station has a second angle at the first time and the second time, and at the first time and the second time of the second condensing device. The orbit calculation system according to claim 7, wherein each position is transmitted to the first optical angle measuring station.
The second optical measuring station has a second angle at the first time and the second time, and at the first time and the second time of the second condensing device. The orbit calculation system according to claim 8, wherein each position is transmitted to the first optical angle measuring station via a ground station.
The second optical measuring station includes a retroreflective unit that retroreflects light.
The first optical angle measuring station emits light toward the retroreflective unit, then receives the light reflected by the retroreflective unit, and the emitted light and the received light are used to obtain the first light. It is equipped with a distance measuring device that calculates the distance between the optical angle measuring station and the second optical angle measuring station.
The orbit calculation system according to claim 8, wherein the orbit calculation unit calculates the position of the second condensing device by using the distance calculated by the distance measuring device.
The orbit calculation system according to claim 8, wherein the second optical angle measuring station is arranged in an orbit having an altitude different from that of the orbit of the flying object.