Processing pipeline for data mining of the primary orbits of the Solar System objects using the Väisälä method Volodymyr Troianskyi1,2,3, , Olexander Bazyey2,4, , Sergii Khlamov5, , Tetiana Trunova5 and Stanislav Orlov6 1 Astronomical Institute of the Slovak Academy of Sciences, Vysoké Tatry, Tatranska Lomnica, 05960, Slovak Republic 2 Odesa I. I. Mechnikov National University, Pastera street 42, Odesa, 65082, Ukraine 3 Vasyl Stefanyk Precarpathian National University, Shevchenko Street 57, Ivano-Frankivsk, 76000, Ukraine 4 Odesa National Maritime University, Mechnykov Street 34, Odesa, 65029, Ukraine 5 Kharkiv National University of Radio Electronics, Nauki avenue 14, Kharkiv, 61166, Ukraine 6 National Aerospace University Kharkiv Aviation Institute, 17, Kharkiv, 61000, Ukraine Abstract In this paper we presented the processing pipeline for data mining of the primary orbits and discovery of the new Solar System objects using the CoLiTec software and classical Väisälä method. The 1st step of a such method is performing the sequential astronomical observations and the classical astrometric data reduction. It is implemented as a combination of the general mathematical algorithms for processing of the astronomical frames incapsulated in the CoLiTec software. It contains the inverse median filter, astronomical calibration with usage of the master-frames, detection of the astronomical objects and its trajectories in series of frames and others features. After that the especial parameters orbit (right ascension and declination) are calculated at the certain time. When we receive at least two observations of the same object under study at different moments, the classic Väisälä algorithm can be applied. For this, the primary orbit from the two nearest observations is calculated. Then the geocentric rectangular coordinates and the appropriate geocentric velocity components are computed. This give us an opportunity for determining the Keplerian elements of orbit of the investigated object at any time we are interested in. The developed algorithm is implemented as a processing pipeline that includes a combination of the CoLiTec software and created tool with encapsulation of the Väisälä method. The created processing pipeline for data mining of the primary orbits and discovery of the new Solar System objects using the CoLiTec software and classical Väisälä method is practically verified. Using it a few new asteroids are first reported and a few lost small bodies of the Solar System are found in the Kyiv comet station and the Odesa - Mayaky observatory. Keywords Data mining, processing pipeline, astronomical observations, CoLiTec software, Väisälä method, image processing, image calibration, astrometric reduction, object detection, positional coordinates, orbit determination, Solar System objects1 1. Introduction The asteroid-comet hazard is a big potential problem in the XXI century. It leads to the developing astronomical scientific direction like an image processing, which includes the following algorithms for astrometric reduction [1], photometric reduction [2], detection or even discovery of the Solar System objects [3], etc. A common goal of all scientific and technological algorithms is to automate as much as available processes without any human actions. In general cases it can be done by the different astronomical scientific processing pipelines. In these pipelines the various data mining [4] and knowledge discovery in databases (KDD) [5] tasks are used for speeding up and optimizing the astronomical ICST-2024: Information Control Systems & Technologies, September 23-25, 2023, Odesa, Ukraine v.troianskyi@onu.edu.ua (V. Troianskyi); o.bazyey@onu.edu.ua (O. Bazyey); sergii.khlamov@gmail.com (S. Khlamov); tetiana.trunova@nure.ua (T. Trunova); s.v.orlov@student.khai.edu (S. Orlov) 0000-0002-5899-2300 (V. Troianskyi); 0000-0001-6886-4575 (O. Bazyey); 0000-0001-9434-1081 (S. Khlamov); 0000-0003- 2689-2679 (T. Trunova); 0009-0008-0680-206X (S. Orlov) © 2024 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR ceur-ws.org Workshop ISSN 1613-0073 Proceedings data processing. The astronomical scientific information is collected from the different historical archives, clusters, Virtual Observatories [6], clouds, astronomical astrometric and photometric catalogues [7], different servers and other storages. Almost all astronomical scientific observations are created by the charge-coupled device (CCD) [8] and performed during the specified observational period of the investigated small celestial objects of the Solar System (SSOs) [9] (like, comets or asteroids [10]), as well as the artificial satellites [11]. After performing the series of observations of the investigated SSOs it is required to analyze the results of observation, which can include the period and shape of rotations determining of such investigated SSOs. That is why it is required to create a mathematical method to determine a primary orbit and discovery of the new Solar System objects using the CoLiTec software and the classical Väisälä method. Such a method gives a possibility for estimating the geocentric rectangular positional coordinates of a SSO under study with and calculating the different Keplerian elements of the SSO's orbit. At the same time this algorithm can improve the conditional expectation of the correct discovery [12] of the real SSOs. The Väisälä algorithm is a very useful algorithm when used in the different situations [13]. Its importance becomes noticeable with the particularly short observation arcs, which are insufficient to confirm an exact orbit and predict a position of the investigated small astronomical Solar System object during the next week. 2. Related Works Each Solar System object presented in a CCD-frame has a typical form of its image [14]. The common methods for the image processing [15] and machine vision [16] are developed for detection/recognition such images of SSOs and an estimation of their positional and motion parameters [17]. Such methods are based on the analysis of only those pixels that potentially belong to the investigated object. The disadvantages of such methods are very low accuracies when the typical form of object has a different shape [18]. The methods for pixelization and segmentation by the aperture brightness [19] will work only with a single image of each SSO. A method for the matched filtration [20] uses a model image, which is analytically determined. The disadvantages of the methods above are the big complexity and low accuracy during the processing, when an image has a several peaks of magnitude. Methods for the Wavelet analysis [21] or even time series analysis [22] are not so effective, because we do not have a big volume of the input data to be analyzed. Also, the disadvantage of such algorithms is the corrupting of the general statistics and possibility to process only clear measurements without any deviations in the typical form of image. 3. Processing pipeline for data mining of the primary orbits of the Solar System objects 3.1. Astronomical calibration Astronomical image calibration is a critical process in astrophotography and observational astronomy [23]. These calibration frames are used to correct different imperfections and artifacts introduced during the image acquisition process, ensuring the accuracy and quality of the final images: • Bias Frames. Bias frames capture the electronic signal present in the camera sensor when no light is allowed to enter the telescope. They record the baseline electronic offset or zero signal level of the sensor. Bias frames help correct for the sensor's inherent electronic noise, such as readout noise and amplifier offsets; • Dark Frames. Dark frames are images taken with the telescope or camera covered, capturing the thermal signal produced by the sensor itself. This signal is a result of heat generating random electronic noise in the detector. Dark frames help in correcting thermal noise, hot pixels, and other sensor-specific imperfections; • Flat Frames. Flat frames are images of a uniformly illuminated surface or a blank, evenly illuminated part of the sky. They are used to correct for variations in sensitivity across the camera sensor, as well as for dust specks or imperfections in the optical system. Flat frames help to normalize the pixel-to-pixel sensitivity differences in the sensor and remove any vignetting or dust shadows. The calibration process involves the following steps: 1. Acquisition: Bias, dark, and flat frames are captured under specific conditions matching the settings used for the actual astrophotography session. It's important to ensure that the exposure times and temperatures match those of the light frames (images of celestial objects). 2. Subtraction: The bias frame is subtracted from all the frames (light, dark, and flat) to remove the electronic offset and noise present in the images. 3. Dark Subtraction: Dark frames are subtracted from the light frames to eliminate the thermal signal caused by the camera sensor's temperature during image acquisition. 4. Flat Fielding: Flat frames are used to normalize the response of the sensor. Light frames are divided by the normalized flat frames to correct variations in pixel sensitivity and eliminate artifacts caused by dust, vignetting, or optical imperfections. 5. Post-Processing: Once the calibration process is completed, further image processing techniques, such as stacking, alignment, and color calibration, are often applied to enhance the final image quality. Calibrating images with bias, dark, and flat frames is crucial in producing high-quality astrophotographs and ensuring accurate scientific analysis by removing unwanted noise and artifacts introduced during the imaging process. In the developed algorithm we performed the different actions for the astronomical data calibration, like calibration dark frames subtracting, calibration bias frames subtracting and calibration flat frames fielding/dividing [24]. 3.2. Orbit determination Every astronomical observation of the investigated astronomical object provides the topocentric directions to it:  X i = −G1i cos i cos i ,  (1) Yi = −G1i cos i sin i ,  Z = −G sin  ,  i 2i i where a coefficient G and a local sidereal time ( ) are calculated by the following equations: ae G1i = + Hi , 1 − ( 2 f − f 2 ) sin 2 i (2) (1 − f ) ae 2 G2i = + Hi , 1 − ( 2 f − f 2 ) sin 2 i d (3) i =  g 0 + ( ti − t0 ) + Ei , dt where i, Ei are the geodetic latitude and longitude of the observer; Hi is the altitude of the observer; Earth ellipsoid model recommended by the IERS [25]; semi-major axis (aE) of 6378171.364331512 m; inverse of polar flattening (1/f) of 297.7736994668283; g0 is Greenwich Mean Sidereal Time (GMST). After astrometric reduction to the Earth center, we received the right ascension ( i) (RA) and declination ( i) (DE) for the current time (ti). This gives us a possibility to create a system of connection equations:  i cos  i cos  i = xi + X i ,  (4)  i sin  i cos  i = yi + Yi ,   sin  = z + Z ,  i i i i where Xi, Yi, Zi are rectangular geocentric equatorial coordinates of the Sun; xi, yi, zi are rectangular heliocentric equatorial coordinates of an investigated astronomical object. Geocentric coordinates of the Sun can be taken from SS model. And geocentric distances ( i) and heliocentric coordinates (xi, yi, zi) of an investigated astronomical are unknown in system 4. 3.3. Väisälä algorithm General representation of the Väisälä algorithm requires the 2 nearest astronomical observations of the same investigated object at times t1 and t2. So, we can calculate the primary orbit from such two observations using a set of the following parameters (t1, 1, 1, t2, 2, 2). To calculate the geocentric rectangular coordinates of an investigated object we can use the following system of equations:  x2 =  2 cos  2 + X 2 ,  (5)  y2 =  2 sin  2 + Y2 ,  z =  tan  + Z ,  2 2 2 2 where  2 =  2 cos  2 . To calculate the geocentric velocity components of an object we can use the following system:  1 cos 1 − F1 x2 + X 1  x2 = ,  G 1  1 sin 1 − F1 y2 + Y1 (6)  y2 = ,  G 1   tan 1 − F1 z2 + Z1  z2 = 1 ,  G1 F1r2 + X 1 x1 + Y1 y1 + Z1 z1 1 = , x2 cos 1 + y2 sin 1 + z2 tan 1 F1 = 1 − A 2 , G1 =  − B 3 , (7) 3 r A A= ,B = , 2 3 where  = k ( t1 − t2 ) is a time interval between two nearest observations; k is a gravitation constant. The determined parameters ( x2 , y2 , z2 , x2 , y2 , z2 ) will be used for calculation the Keplerian elements of the orbit of an investigated object at a time t2. 3.4. Processing pipeline for data mining of the primary orbits of the Solar System objects An architecture of a processing pipeline for data mining of the primary orbits and discovery of the new Solar System objects using the CoLiTec software and classical Väisälä method is presented in Figure 1. Such method contains the next steps. 1. Perform the observations of the investigated SSOs. 2. Perform the classic data reduction process using bias, flat, and dark calibration frames. 3. Perform the astrometric reduction of the investigated SSOs by the CoLiTec software. 4. Calculation of a primary orbit of the investigated SSOs by the Väisälä algorithm. 5. Repeat the astronomical observations of the investigated SSOs to confirm a discovery of the new SSOs. Figure 1: Architecture of a processing pipeline for data mining of the primary orbits and discovery of the new Solar System objects using the CoLiTec software and classical Väisälä method 4. Results 4.1. CoLiTec software The created algorithm to determine a primary orbit and discoveries of the SSOs using the Väisälä algorithm works in combination with the CoLiTec software (https://colitec.space) [26]. The specific functional features related to the mathematical method for determining the primary orbits and discovery of asteroids in the CoLiTec software are: • automated frame calibration [9]; • cosmetic frame correction; • track-and-stack feature; • brightness equalization [2]; • background alignment; • astronomical image filtering [27]; • determining the contours of objects [28]; • fully automated robust method of the astrometric reduction [1]; • fully automated robust method of the photometric reduction [23]; • support of the multi-threaded processing; • On-Line Data Analysis System (OLDAS) for managing the processing pipeline; • transferring of astronomical data with intermediate storage [29]; More extended details about the CoLiTec software are presented in these works and research [30, 31, 32]. The high-level processing pipeline with the developed modules and implemented methods of the CoLiTec software is presented in the Figure 2. Figure 2: The high-level processing pipeline of the CoLiTec software The CoLiTec software is installed at the different telescopes at the various observatories in Ukraine and around the world: • OMT-800 and AZT-3 telescopes installed at the Odesa-Mayaky observatory [33]; • SANTEL-400AN telescope installed at the ISON-NM observatory; • ISON-Uzhgorod [34]; • VNT and Celestron C11 telescopes installed at the Vihorlat Observatory [23]; • PROMPT-8 telescope installed at Cerro Tololo observatory [35]; • NARIT (National Astronomical Research Institute of Thailand) [36]; • AZT8 and Takahashi BRC-250M telescopes. All listed above telescopes installed at the observatories have an official identifiers received from the Minor Planet Center (MPC) (https://minorplanetcenter.net) MPC code from the International Astronomical Union (IAU) (https://iau.org). 4.2. Odesa-Mayaky (observatory) [code: 583] Astronomical observations were made on the OMT-800 telescope installed at the Mayaky observational station of the Astronomical Observatory of the Odesa I. I. Mechnykov National University. The official MPC code of the Odesa-Mayaky station is 583 (see Figure 3) [37]: Figure 3: The official MPC codes of the Odesa-Mayaky observatory and the Kyiv comet station The 0.8-m main hyperbolic mirror is installed at the OMT-800 telescope and has a useful focus ratio as f = 1/2.7. The FLI ML09000 CCD-camera is used as an image CCD-detector along with a 4- lens field corrector for a primary focus, which allows receiving the field of view (FOV) as 58.6'×58.6' and the scale of image as 1.15 arcseconds per pixel. During the astronomical observations at the OMT-800 telescope and research after processing the results, the next new asteroids were firstly discovered: • 2017 BC94 (MPS 766258); • 2017 QX33 (MPS 813479); • 2017 RV12 (MPS 817137); • 2017 QJ36 (MPS 817122). During the astronomical observations at the OMT-800 telescope and research after processing the results, the next lost SSOs were found again: • 2012 FN61 = 2010 VP219 (MPS 441842); • 2017 AB8 = 2014 OD380 (MPS 525811). 4.3. Kyiv comet station (observatory) [code: 585] Astronomical observations were made on the AZT-8 telescope installed at the Kyiv comet station of the Lisnyky observational station of the Taras Shevchenko National University of Kyiv. The official MPC code of the Kyiv comet station is 585 (see Figure 3). The installed AZT-8 telescope is a telescope-reflector. The FLI PL4710 CCD-camera (63.5 mm) was used as an image CCD-detector, which allows receiving the FOV as 16.2'×16.7' and the scale of image as 0.948 arcseconds per pixel. During the astronomical observations on the Kyiv comet station and research after processing the results, the new asteroid 2017 SV39 (MPS 828365) was firstly discovered and the next lost SSOs were found again: • (536266) = 2007 HU101 = 2015 CX48 = 2017 ST39 (MPS 891037); • (540584) = 2000 WN134 = 2007 XY37 = 2010 RD162 = 2015 DH168 = 2017 TS7 (MPS 23003). 5. Conclusions We present a processing pipeline for data mining of the primary orbits and discovery of the new Solar System objects using the CoLiTec software and classical Väisälä method. The literature review showed the benefits and disadvantages of the existing methods and algorithms. So, we made a decision to improve such algorithms for the better accuracy and quality of detection of the new or lost SSOs. The first step of such computational algorithm is to perform the observations and the standard data reduction. It was performed in scope of the common mathematical methods for the astronomical image processing and numeric analysis [38] implemented in the CoLiTec software. It includes inverse median filtration, calibration using the master frames, object detection, etc. Then the orbit parameters, like the right ascension and declination, were determined for a certain moment. After receiving the two observations of an investigated object at different times, the classic Väisälä method is applied. The created algorithm is implemented as a processing pipeline that contains a cooperation between the CoLiTec software pipeline and a special implemented plugin based on the general Väisälä algorithm. Such a plugin is realized as a console application based on the programming language Delphi. The CoLiTec software was selected among another astronomical software tools for astrometric reduction like Astrometrica [39] as the more accurate software with more precise measuring of positional coordinates of celestial objects. The created mathematical method was verified in practice based on the classical mathematical models [40] and fuzzy clustering data arrays with omitted observations [41] as well as statistical modeling [42]. With help of the developed processing pipeline, the several new SSOs (Main-belt asteroid 2017 AB8 (2014 OD380), Mars-crossing asteroid 2017 QX33, Main-belt asteroid 2017 RV12) were firstly discovered and some of the lost SSOs were found again during the astronomical observations at the Kyiv comet station and the Odesa-Mayaky observatory. Acknowledgements The research was supported by the National Scholarship Programme of the Slovak Republic for 2023- 2024 years and the Ukrainian project of fundamental scientific research Development of computational methods for detecting objects with near-zero and locally constant motion by optical- electronic 0124U000259 in 2024-2026 years. References [1] S. Khlamov, I. Tabakova, T. Trunova, and Z. Deineko, Machine Vision for Astronomical Images Using the Canny Edge Detector, CEUR Workshop Proceedings, vol. 3384, pp. 1 10, 2022. [2] A new tool for an automated reduction of photometric observations, Contributions of the Astronomical Observatory Skalnate Pleso, 49 2 (2019) 151- 153. [3] V. Troianskyi, V. Godunova, A. Serebryanskiy, et al., Optical observations of the potentially hazardous asteroid (4660) Nereus at opposition 2021, Icarus 420 (2024) 116146. doi: 10.1016/j.icarus.2024.116146. [4] S. Khlamov, et al., Automated data mining of the reference stars from astronomical CCD frames, CEUR Workshop Proceedings, vol. 3668, pp. 83 97, 2024. [5] K. D. Borne, Scientific data mining in astronomy, Next Generation of Data Mining, Chapman and Hall/CRC, 2008. [6] I. B. Vavilova, Y. S. Yatskiv, L. K. Pakuliak, et al., UkrVO astroinformatics software and web- services, in: Proceedings of the International Astronomical Union, vol. 12, issue S325, 2016, pp. 361-366. doi: 10.1017/S1743921317001661. [7] V. Akhmetov, et al., Fast coordinate cross-match tool for large astronomical catalogue, Advances in Intelligent Systems and Computing 871 (2019) 3 16. doi: 10.1007/978-3-030-01069-0_1. [8] G. Adam, P. Kontaxis, L. Doulos, E.-N. Madias, C. Bouroussis, F. Topalis, Embedded Microcontroller with a CCD Camera as a Digital Lighting Control System, Electronics 8 1 (2019). doi: 10.3390/electronics8010033. [9] D. Oszkiewicz, et al., Spins and shapes of basaltic asteroids and the missing mantle problem, Icarus, 397 (2023) 115520. doi: 10.1016/j.icarus.2023.115520. [10] V. Troianskyi, P. Kankiewicz, D. Oszkiewicz, Dynamical evolution of basaltic asteroids outside the Vesta family in the inner main belt, Astronomy and Astrophysics, 672 (2023). doi: 10.1051/0004-6361/202245678. [11] V. Akhmetov, et al., Cloud computing analysis of Indian ASAT test on March 27, 2019, in: Proceedings of the 2019 IEEE International Scientific-Practical Conference: Problems of Infocommunications Science and Technology, PIC S and T 2019, 2019, pp. 315 318,. doi: 10.1109/PICST47496.2019.9061243. [12] E. L. Lehmann, J. P. Romano, G. Casella, Testing statistical hypotheses, New York: Springer, 2005. [13] V. Troianskyi, V. Kashuba, O. Bazyey, et al., First reported observation of asteroids 2017 AB8, 2017 QX33, and 2017 RV12, Contributions of the Astronomical Observatory Skalnaté Pleso 53 (2023) 5-15. doi: 10.31577/caosp.2023.53.2.5. [14] V. Savanevych, et al., Formation of a typical form of an object image in a series of digital frames, Eastern-European Journal of Enterprise Technologies, 6 2 (2022) 51 59. doi: 10.15587/1729- 4061.2022.266988. [15] W. Burger, M. Burge, Principles of digital image processing: fundamental techniques, New York, NY: Springer, 2009. [16] R. Klette, Concise computer vision. An Introduction into Theory and Algorithms, London: Springer, 2014. [17] V. Savanevych, et al., Mathematical methods for an accurate navigation of the robotic telescopes, Mathematics 11 10 (2023) 2246. doi: 10.3390/math11102246. [18] T. Ando, Bayesian model selection and statistical modeling, CRC Press, 2010. [19] V. Vlasenko, et al., Devising a procedure for the brightness alignment of astronomical frames background by a high frequency filtration to improve accuracy of the brightness estimation of objects, Eastern-European Journal of Enterprise Technologies 2 2 (2024) 31 38. doi: 10.15587/1729-4061.2024.301327. [20] S. Khlamov, et al., Development of computational method for matched filtration with analytic profile of the blurred digital image, EEJET 5 4 (2022) 24 32. doi: 10.15587/1729-4061.2022.265309. [21] M. Dadkhah, et al., Methodology of wavelet analysis in research of dynamics of phishing attacks, International Journal of Advanced Intelligence Paradigms 12 (2019) 220-238. doi: 10.1504/IJAIP.2019.098561. [22] L. Kirichenko, A.S.A. Alghawli, T. Radivilova, Generalized approach to analysis of multifractal properties from short time series, International Journal of Advanced Computer Science and Applications, 11 5 (2020) 183 198. doi: 10.14569/IJACSA.2020.0110527. [23] I. Kudzej, et al., CoLiTecVS A new tool for the automated reduction of photometric observations, Astronomische Nachrichten 340 (2019) 68 70. doi: 10.1002/asna.201913562. [24] D. Oszkiewicz, et al., Spin rates of V-type asteroids, Astronomy & Astrophysics 643, (2020) A117. doi: 10.1051/0004-6361/202038062. [25] Hu Xuanyu, C. K. Shum, M. Bevis, A triaxial reference ellipsoid for the Earth, Journal of Geodesy 97 29 (2023). doi: 10.1007/s00190-023-01717-1. [26] S. Khlamov, et al., Machine vision for astronomical images using the modern image processing algorithms implemented in the CoLiTec software, Measurements and Instrumentation for Machine Vision, Chapter 12: CRC Press, Taylor & Francis Group, (2024) 269-310. doi: 10.1201/9781003343783-12. [27] S. Khlamov, I. Tabakova, T. Trunova, Recognition of the astronomical images using the Sobel filter, in: Proceedings of the 29th IEEE International Conference on Systems, Signals, and Image Processing, IWSSIP 2022, Sofia, Bulgaria, June 1st 3rd, 2022. doi: 10.1109/IWSSIP55020.2022.9854425. [28] H. Khudov, I. Ruban, O. Makoveichuk, et al., Development of methods for determining the contours of objects for a complex structured color image based on the ant colony optimization algorithm, EUREKA, Physics and Engineering 2020 1 (2020) 34 47. doi: 10.21303/2461- 4262.2020.001108. [29] V. Tkachov, et al., Method for transfer of data with intermediate storage, in: Proceedings of the 1st International Scientific-Practical Conference Problems of Infocommunications Science and Technology, PIC S and T 2014, Kharkiv, Ukraine, October 14th 17th, 2014, pp. 105 106. doi: 10.1109/INFOCOMMST.2014.6992315. [30] S. Khlamov, et al., Development of the matched filtration of a blurred digital image using its typical form, Eastern-European Journal of Enterprise Technologies 1 9 (2023) 62 71. doi: 10.15587/1729-4061.2023.273674. [31] S. Khlamov, et al., The astronomical object recognition and its near-zero motion detection in series of images by in situ modeling, in:Proceedings of the 29th IEEE International Conference on Systems, Signals, and Image Processing, IWSSIP 2022, Sofia, Bulgaria, June 1st 3rd, 2022. doi: 10.1109/IWSSIP55020.2022.9854475. [32] S. Khlamov, et al., Astronomical knowledge discovery in databases by the CoLiTec software, in: Proceedings of the 12th International Conference on Advanced Computer Information Technologies, ACIT 2022, Ruzomberok, Slovakia, September 26th 28th, 2022, pp. 583 586. doi: 10.1109/ACIT54803.2022.9913188. [33] S. Kashuba, et al., The Simeiz plate collection of the ODESSA astronomical observatory, in: Proceedings of the 11th Bulgarian-Serbian Astronomical Conference, 2018, pp. 207-216. [34] V. Kudak, V. Epishev, V. Perig, I. Neybauer, Determining the orientation and spin period of TOPEX/Poseidon satellite by a photometric method, Astrophysical Bulletin 72 3 (2017) 340-348. doi: 10.1134/S1990341317030233. [35] T. Li, D. DePoy, J. Marshall, D. Nagasawa, D. Carona, S. Boada, Monitoring the atmospheric throughput at Cerro Tololo Inter-American Observatory with aTmCam, Ground-based and Airborne Instrumentation for Astronomy V 9147 (2014) 2194 2205. [36] S. Rattanasoon, E. Semenko, D. Mkrtichian, S. Poshyachinda, Spectroscopic Devices for Asteroseismology With Small Telescopes in NARIT, arXiv -> astro-ph, vol. 2307.03985, 2023. doi: 10.48550/arXiv.2307.03985. [37] Minor Planet Center, List Of Observatory Codes. URL: https://www.minorplanetcenter.net/iau/lists/ObsCodesF.html. [38] R. Burden, J. Faires, A. Reynolds, Numerical Analysis, Brooks/Cole: Boston, 2010. [39] V. E. Savanevych, A. B. Briukhovetskyi, Yu. N. Ivashchenko, et al., Comparative analysis of the positional accuracy of CCD measurements of small bodies in the solar system software CoLiTec and Astrometrica, Kinematics and Physics of Celestial Bodies 31 6 (2015) 302-313. doi: 10.3103/S0884591315060045. [40] I. Grebennik, O. Chorna, I. Urniaieva, Distribution of Permutations with Different Cyclic Structure in Mathematical Models of Transportation Problems, in: Proceedings of the 12th International Conference on Advanced Computer Information Technologies, ACIT 2022, Ruzomberok, Slovakia, September 26th 28th, 2022, pp. 18 21. doi: 10.1109/ACIT54803.2022.9913183. [41] Z. Hu, Y. Bodyanskiy, O. Tyshchenko, V. Tkachov, Fuzzy clustering data arrays with omitted observations, International Journal of Intelligent Systems and Applications 9 6 (2017) 24 32. doi: 10.5815/ijisa.2017.06.03. [42] V. Shvedun, et al., Statistical modelling for determination of perspective number of advertising legislation violations, Actual Problems of Economics 184 10 (2016) 389-396.