1. Introduction
Conventional agriculture practices adopt homogenization-farming processes to enhance production, with lower costs, standardization, mass production, and production efficiency. For example, waterfowl egg hatching management keeps all waterfowl eggs in the same nature by controlling temperature and humidity in the incubator, periodically inspecting waterfowl egg characteristics, and waterfowl egg status in the hatcher. However, due to agricultural knowledge and the advancement of science and technology, it is found that the above-mentioned practices eliminate the waterfowl eggs—slightly different from the expected status in each phase of the lifecycle due to the changing natural environment and climate. However, the individual variations of waterfowl eggs is tolerable if we provide customized care, i.e., doing the right thing, in the right place, in the right way, and at the right time. That is—this elimination of waterfowl eggs can be avoided when farmers provide customized care to bring these waterfowl eggs back to the expected developing path by a universal prediction model, for the quantitative description of biomass changing with time and environment, and the corresponding automatic environment control.
The previous studies mainly explore the relationship between the environmental factors and the hatching rate of fertilized eggs. Lourens et al. [
1,
2] studied the effects of eggshell temperature and oxygen concentration on embryo development. Joseph et al. [
3] indicate the impact of the incubation environment on broiler chick development. In studies [
4,
5,
6], the authors investigate the relationship between humidity/eggshell temperature and incubation. However, most waterfowl farms in Taiwan control temperature and humidity based on experience and do not adopt a scientific hatching model. In order to dissipate the heat generated by embryo development and promote the embrittlement of the eggshell to help the young birds to peck the eggshells, it is necessary to cool the egg and spray water during the incubation process. Taiwan is humid and has high temperatures. Therefore, it is necessary to rely on human experiences to control the egg cooling of waterfowl breeding eggs and the time, frequency, and microenvironment of water spraying during the hatching process. The operation and management during the incubation process of waterfowl eggs has a significant impact on the hatching results, so it is necessary to find out customized parameters suitable for the local environment.
This study therefore applies artificial intelligence (AI) technology to build a prediction model to quantitatively describe the biomass development path along with each phase of the waterfowl-like production lifecycle and the corresponding environment settings. First, we use the concept of precision agriculture to cut into the field of waterfowl hatching and find major hatching condition parameters out of the parameters explored by the previous works. In Taiwan, the industry’s requirements for the hatching environment of waterfowl are still in the traditional agricultural stage, pursuing the uniformity of the hatching environment, and relying heavily on human experience to make parameter corrections for the individual differences of waterfowl. We intended to carry out precise sensing of each waterfowl egg, which could provide customized parameter adjustments for waterfowl eggs with large individual differences.
A thermal imager was used in our experiment to measure the hatching environment and the eggshell temperature of each fertilized duck egg. The eggshell temperature data were collected to identify the suitable parameter settings for the incubation of fertilized duck eggs. Accordingly, we employed air quality sensors to gather O2, CO2, RH, and temperature inside the incubator.
In the second step, we built a sufficient hatching prediction model to predict the hatching rate of fertilized eggs based on the explored major parameters in the first step. This model is expected to help the duck farmer, to provide customized care of each egg with minimized parameters instead of all parameters, which incurs high observation costs and, thus, it can reduce unnecessary elimination of eggs and increase hatchability.
In order to build the prediction model, this study uses the black Muscovy duck hatching as a case study (
Figure 1); we explored key characteristics to quantitatively predict the survivability of the black Muscovy duck egg. Data from the experiment were remotely collected by thermal imaging cameras and air sensors to lower the labor costs and build automatic environment control. The black Muscovy duck was used as a case study to build the prediction model [
7]. These ducks have been bred in Taiwan for more than 300 years. Their feather colors could be white, blue, black, a white–blue combination, a white–black combination, black with white rings, blue with white rings, or chocolate brown [
8]. Initially, people in Taiwan mainly bred black Muscovy ducks. Male Muscovy ducks were a source of food (ginger duck stew) for people to consume during the winter to warm their bodies [
9]. In 1962, large-sized, white Muscovy ducks were introduced to Taiwan from Australia, the United States, and the Netherlands. From 1975, they were widely used due to their economic value; thus the number of original black Muscovy ducks in Taiwan reduced.
In summary, this study provides a prediction model to quantitatively predict the hatching rate of each egg by observing parameters and the waterfowl-like production lifecycle. Two experiments, without any human interference from experienced farmers, were conducted: in the first one, we built the prediction model; the second one was to validate the effectiveness of the prediction model. The proposed prediction model will contribute to the duck industry by minimizing the training costs in the duck industry and maximizing the product qualification rate for individual agricultural products.
2. Literature Review
Many studies have focused on enhancing the hatching rate of fertilized eggs, including exploring the factors that affect the hatching rate and how to inspect the egg effectively by thermal imaging technology.
First, regarding temperature, Lourens et al. [
1] attempted to control machine temperature to maintain the eggshell temperature during incubation; they found that the eggshell temperature could vary independently from the machine temperature. Joseph et al. [
2] further studied the effects of different eggshell temperature settings in one incubation (e.g., the low eggshell temperature at the start of incubation, and high eggshell temperature at the end of incubation) on hatchability. Joseph also suggested that maintaining the eggshell temperature at 37.8 °C from days 0 to 10 leads to the best hatchability. Ipek [
6] noted that embryo development is affected by small changes in the eggshell temperature.
Many scholars have attempted to explore other possible factors on hatchability, i.e., other than temperature. Lourens et al. [
3] suggested that oxygen determines the embryo development (via heat production in the egg) rather than eggshell temperature. Molenaar and van der Pol [
4,
5] found that low relative humidity or low O
2 concentration leads to the highest embryo mortality. In [
10,
11,
12,
13], the researchers designed an automatic sprinkler cooling egg system to maintain a specific temperature, humidity, and concentration of CO
2 settings. In [
14], the authors discuss the possibility of controlling the velocity and temperature of the airflow field in an incubator to increase hatchability.
Many studies apply thermal imaging as a non-invasive and non-contact method of measuring the surface temperature of eggshells in order to eliminate the effect of human factors on egg incubation. McCafferty [
15] demonstrated several ways to use thermal imagers in poultry science, including physiological status monitoring, behavior research, and field survey on poultry. Lin et al. [
16] employed thermal imagery technology to filter fertilized chicken eggs. Sunardi et al. [
17] used thermal imaging to label the eggs and trace their development cycle.
The works discussed above provide valuable knowledge on how to evaluate the egg, with what tools or technology, and explore factors that may affect the hatching rate. However, they do not provide a sufficient hatching prediction model to demonstrate what key factors affect the hatching rate and how to apply these factors to predict the hatching rate of fertilized eggs.
3. The Proposed Method and Environmental Settings of the Experiment
3.1. Proposed Method
This study attempts to build a model based on Extreme Gradient Boosting to quantitatively predict the survivability of the black Muscovy duck egg. Extreme gradient boosting (XGBoost), a type of gradient boosted tree developed by Chen and Guestrin [
18], was used for optimization based on tree ensemble models.
Tree boosting is used to integrate multiple weak learners to form a strong learner, whereas a tree ensemble model considers the parameter optimization problems of multiple trees. In tree ensemble models, the original model is retained through additive training, and a new function is added to the model. Training is conducted according to the previous tree; it involves the addition of a new tree to the model to mitigate the deficiencies of the previous tree and further optimize the model.
where
l is a loss function. The penalty term at the back prevents model overfitting and ensures a more satisfactory model.
XGBoost, mainly used to explain supervised learning, can be applied for classification and to solve regression problems. In this study, we employed XGBoost to identify the successfully hatched fertilized eggs and those that stopped developing. XGBoost—a decision-tree-based, efficient, and expandable machine learning system—is widely used in various fields and has been frequently adopted by top winners at Kaggle competitions [
19,
20,
21].
3.2. Environmental Settings of the Experiment
The environmental settings of the experiment was as follows. First, the air quality detector was built outside the incubator (as shown on the right side of
Figure 2a) and its sensor was placed inside the incubator (as shown on the top side of
Figure 2b). All of the data are stored in the laptop placed outside the incubator, as shown in the right side of
Figure 2a. The thermal camera was placed inside the incubator (as shown in the center of
Figure 2b) to collect the thermal image of the egg tray (as shown at the bottom side of
Figure 2c); the sample of the thermal image is shown in
Figure 2d. The thermal imager was covered and protected by plastic wrap, according to the manufacturer’s instructions, in order to avoid unnecessary damage from the incubator and hatcher. There were 66 eggs randomly placed in the egg tray; the hatching period was from 11 September to 15 October, 2020, as shown in
Figure 2c–e. The hatched ducks in this experiment are shown in
Figure 2e. Climatic data of the experimental site are shown in
Table 1.
Two experiments were conducted for this study. The first experiment was applied to (1) explore the major parameters out of all the parameters to affect the hatching rate; and (2) to build a prediction model. Furthermore, the first experiment was executed in the unfamiliar environment to avoid human interference from experienced farmers. Based on the experience, the hatching rate was around 50% in the first experiment. The second experiment validated the prediction model by controlling parameters explored in the first experiment, without any human interference from experienced farmers.
The thermal imager manufacture periodically calibrated the imager based on the contract, and the accuracy was ±2% of the detected temperature value. The precisions of the O2 and CO2 sensors were:
3.3. Data Collection and Preprocessing Procedure
The sample consisted of 66 Muscovy duck fertilized eggs. We conducted analyses according to the data obtained from the hatching experiment by the candling procedure. Egg candling was implemented on day 6 (E6), day 13 (E13), and day 31 (E31). Unfertilized duck eggs and those with terminated development were removed from the egg tray to prevent the healthy duck eggs from spoiling or rupturing. The egg placement on the egg tray and status of the fertilized eggs are displayed, collected in E13. in
Figure 3.
In addition to the thermal imagers, we installed O
2 and CO
2 sensors as well as temperature and humidity sensors on both the interior and exterior of the incubator, to more precisely control the influence of weather on the hatching result. The data gathered during the incubation process were continuously transmitted to the Taiwan Waterfowl Hatching Database. After the completion of the incubation, we marked the successful hatches among all incubation data and imported the results to machine learning for classification training (
Figure 4).
Before the experiment, we recorded the weight, the long and short diameters, and the shape coefficient of each fertilized egg. The eggshell temperatures of the fertilized eggs were automatically obtained by the thermal images captured (per minute, every day) in the incubator. The experiment recorded the mean, min, and max temperature of each egg of each day. During candling, the eggs were weighted and weight loss of the eggs were calculated. Surface temperatures (eggshell temperatures) were recorded using the highest temperature, the lowest temperature, the maximum temperature difference, and the average temperature. We employed the average eggshell temperature rather than the median temperature to compensate for human error in the manual selection of the surface temperature range. This study’s objective was to determine the key characteristics and build a model to quantitatively describe the biomass development path, along with each phase of the waterfowl production lifecycle. Therefore, we interviewed duck farmers with numerous site visits and collected features that they observed based on long time experience, as shown in
Figure 5b. The interview procedure was conducted and regulated by the “Duck Production System Manual” [
22] of the Council of Agriculture (CoA), Executive Yuan in Taiwan. One of the research team members and authors, Chih-Hsiang Cheng, was the coordinator and representative of CoA, project ID 109M000001, project name “Development of sensing modules suitable for waterfowl hatching equipment”. The incubation features used in the machine learning classification models are displayed in
Figure 5a.
The lifecycle of an egg consists of two phases: (1) incubation and (2) hatching. Each phase has its own model based on XGBoost. We determined how incubation features in various stages affected egg development termination and the final hatching result. Hence, we established two machine learning classification models based on the data gathered during the incubation. Data collected before the third candling examination were used to construct Model 1, and those after the third examination, Model 2.
A FLIR A300 thermal imager was employed and the accuracy is ±2 °C with ±2% of accuracy periodically calibrated by the manufacturer.
Figure 6a is a thermal image of an egg tray with a 6 × 5 space arrangement. The thermal imager detected the changes in the surface temperature as well as the highest, lowest, and average eggshell temperatures of the fertilized eggs.
Figure 6b demonstrates the heating information of the no. EL1 fertilized egg in the incubator and a chart of temperature changes over a course of 25 min). Each fertilized egg had a distinct maximum temperature difference. We recorded the hatching result of each fertilized egg after the incubation ended. The detailed temperature data were then adopted for machine learning for establishing hatching prediction models. After the models were calibrated and optimized, the temperature parameters could serve as references for improving incubators.
4. Results and Discussion
The target of this study was to build a prediction model to quantitatively describe the biomass development path along with each phase of the waterfowl-like production lifecycle. That is, we attempted to build a prediction function:
where Y is the predicted hatching rate. The duck farmer could use this function to quantitatively describe the biomass development path, i.e., the hatching rate, and provide sufficient customized care along with each phase of the waterfowl-like production lifecycle and the corresponding environment settings.
Figure 7 is the classification tree of Model 1. We applied 18 features to construct this model, including the egg size, weight, surface temperature at different stages, and weight loss. In the first stage of the incubation, and with the temperature, humidity, and air quality held constant, three conditions were discovered to have key influences on the hatching results: a long diameter less than 6.65 cm, a weight loss less than 1.45 g between days 0 and 7, and an average eggshell temperature higher than 41.36 °C between days 7 and 14.
Figure 8 is the classification tree result of Model 2. We applied 11 features to establish this model, including the egg size, weight, surface temperature at different stages, and weight loss. In the second stage of the incubation, and with the temperature, humidity, and air quality held constant, two conditions were found to influence the hatching result: a shape coefficient greater than 0.6715 and an average eggshell temperature higher than 41.07 °C between days 14 and 31.
In addition to the analysis results of the classification trees, some crucial data were obtained from the final state of the egg tray (
Figure 9). Nine of the eighteen eggs stopped developing during the last 3 days of the incubation. Judging from the positional relationship of these eggs, some mutual influences might exist among them. The strikethroughs in
Figure 8 indicate abnormal data during the three egg candling examinations, which we regarded as manual marking errors. Hence, the data of B15, B32, B36, B60, and B64 were excluded.
Area A, random duck eggs placed to fill the remaining spaces; area B, 66 Muscovy duck fertilized eggs; strikethrough, data error; baseline, unfertilized duck eggs; grayed out boxes, duck eggs that have stopped growing; blacked out boxes, duck eggs that died in the tray.
Based on the above results, we explored four characteristics out of 25 characteristics and built a model to quantitatively predict the survivability of black Muscovy duck egg. The key characteristics include the weight loss percent of an egg, the average eggshell temperature in different phases of the product lifecycle, and the shape coefficient of an egg. Therefore, the prediction function based on the XGBoost algorithm is:
where Y is the predicted hatching rate,
x1 is the weight loss percent of an egg,
x2 is the average eggshell temperature in different phases of the product lifecycle,
x3 is the shape coefficient of an egg,
x4 is the long diameter of an egg. The duck farmer can use the thermal camera to automatically measure the four characteristics to quantitatively describe the biomass development path, i.e., the hatching rate, and provide sufficient customized care along with each phase of the waterfowl-like production lifecycle and the corresponding environment settings.
To validate the prediction model that we built in the first experiment, we applied this prediction model and conducted the second experiment by controlling parameters explored in the first experiment, without any human interference from experienced farmers.
Figure 10 shows the hatching rate rose from 47% in the first experiment to 62% in the second experiment, and this result validates the effectiveness of our prediction model.
In addition, we furthermore evaluated the effects of O
2 and CO
2 values inside the incubator and hatcher toward the survivability of the black Muscovy duck egg between two experiments.
Figure 10 illustrates the four types of air quality values inside the incubator. The data from the first and second experiments are displayed using blue and red lines, respectively. The first experiment shows a 47% hatchability, and the second experiment demonstrates a 62% hatchability. The different ‘hatchability’ between the two experiments is highly related to the difference in CO
2 values in the two experiments. This means that more eggs survived in the second experiment than the first experiment; therefore, the oxygen consumption and CO
2 value of the second experiment were higher than that of the first experiment. However, the CO
2 and O
2 values are the results of the survivability of the black Muscovy duck egg rather than the cause of the survivability of the black Muscovy duck egg. Therefore, our model automatically excludes the CO
2 and O
2 values as being the key characteristics that predict survivability.
5. Conclusions
In this study, we employed XGBoost to determine four essential incubation features out of 25 features and to build a prediction model. The effectiveness of our prediction model is validated by the experiment results, i.e., the hatching rate rose from 47% to 62% without any human interference from experienced farmers. The results show that our proposed prediction model is capable of decreasing the training costs and enhancing the product qualification rate for individual agricultural products.
This study has several implications for precision agriculture researchers and practitioners. First, this study offers a new system development path in order for farmers to adopt AI technology to predict the hatching rate. This study demonstrates how the prediction model is built and how to use this function to enhance agriculture production.
Second, farmers could apply the proposed model to develop a hatching and business strategy, rather than merely rely on experience. For example, farmers could adjust their future hatching equipment investments by this model, according to the breed of the eggshell and the corresponding nature characteristics.
Third, agriculture wholesalers could use this model to understand the productivity of the breed of the duck egg and predict the business strategy of the farmer, including the pros and cons of their short-/long-term production strategies.
In the future, we will extend our proposed model and explore the possible ensemble learning architecture to minimize false positives and false negatives on classification performance.
Author Contributions
Conceptualization, S.-C.C.; methodology, S.-C.C. and C.-H.C.; software, S.-C.C.; validation, S.-C.C. and C.-H.C.; formal analysis, S.-C.C. and Y.-H.C.; investigation, C.-H.C.; resources, C.-H.C. and Y.-H.C.; data curation, S.-C.C.; writing original draft preparation, S.-C.C.; writing review and editing, S.-C.C. and Y.-H.C.; visualization, S.-C.C.; supervision, S.-C.C.; project administration, S.-C.C. All authors have read and agreed to the published version of the manuscript.
Funding
This study is funded by the Ministry of Science and Technology of Taiwan under contract no. MOST 110-2221-E-227-001, and by the Council of Agriculture of Taiwan under project ID 109M000001 with the project name “Development of sensing modules suitable for waterfowl hatching equipment”.
Institutional Review Board Statement
The study was approved (No. LRIIL IACUC 109015) by the Institutional Animal Care and Use Committee, Council of Agriculture, Executive House, Taiwan.
Acknowledgments
The authors would like to thank the Ministry of Science and Technology of Taiwan for financially supporting this research under contract no. MOST 110-2221-E-227-001. The authors furthermore thank the Council of Agriculture of Taiwan, for financially supporting this research under project ID 109M000001 with the project name “Development of sensing modules suitable for waterfowl hatching equipment”.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Lourens, A.; van den Brand, H.; Meijerhof, R.; Kemp, B. Effect of eggshell temperature during incubation on embryo development, hatchability, and posthatch development. Poult. Sci. 2005, 84, 914–920. [Google Scholar] [CrossRef] [PubMed]
- Lourens, A.; van den Brand, H.; Heetkamp, M.J.W.; Meijerhof, R.; Kemp, B. Effects of Eggshell Temperature and Oxygen Concentration on Embryo Growth and Metabolism During Incubation. Poult. Sci. 2007, 86, 2194–2199. [Google Scholar] [CrossRef] [PubMed]
- Joseph, N.S.; Lourens, A.; Moran, E.T., Jr.; Meijerhof, R.; Kemp, B. The effects of suboptimal eggshell temperature during incubation on broiler chick quality, live performance, and further processing yield. Poult. Sci. 2006, 85, 932–938. [Google Scholar] [CrossRef] [PubMed]
- Molenaar, R.; van den Anker, I.; Meijerhof, R.; Kemp, B.; van den Brand, H. Effect of eggshell temperature and oxygen concentration during incubation on the developmental and physiological status of broiler hatchlings in the perinatal period. Poult. Sci. 2011, 90, 1257–1266. [Google Scholar] [CrossRef] [PubMed]
- van der Pol, C.W.; van Roovert-Reijrink, I.A.; Maatjens, C.M.; van den Brand, H.; Molenaar, R. Effect of relative humidity during incubation at a set eggshell temperature and brooding temperature posthatch on embryonic mortality and chick quality. Poult. Sci. 2013, 92, 2145–2155. [Google Scholar] [CrossRef] [PubMed]
- Ipek, A.; Sahan, U.; Baycan, S.C.; Sozcu, A. The effects of different eggshell temperatures on embryonic development, hatchability, chick quality, and first-week broiler performance. Poult. Sci. 2014, 93, 464–472. [Google Scholar] [CrossRef] [PubMed]
- “Wujie Black Muscovy Duck Will Become the Only Purebred Black Muscovy Duck in Taiwan”. Available online: https://www.google.com/url?sa=i&url=https%3A%2F%2Fnews.ltn.com.tw%2Fnews%2Flife%2Fbreakingnews%2F1034652&psig=AOvVaw3bsef1xkPrI9g18PPVjmVZ&ust=1633385411623000&source=images&cd=vfe&ved=2ahUKEwiPs5n7oK_zAhVENqYKHYsWDB8QjRx6BAgAEAk (accessed on 10 October 2021).
- Pingel, H.; Tieu, H.V. Duck Production; The Agricultural Publishing House: Hanoi, Vietnam, 2005. [Google Scholar]
- Zhou, G.Y.; Huang, H.H. The Duck Industry in Taiwan; Joint Commission on Rural Reconstruction: Taiwan, 1970. Available online: https://library.coa.gov.tw/lb/QueryLibrary.asp?FCode=900&NO=000005967 (accessed on 10 October 2021).
- Lan, W.-T. The Development of Automatic Waterfowl Incubator with Pre-Controlled Thermal System. Master’s Thesis, Department of Biomechatronic Engineering, National Chiayi University, Chiayi City, Taiwan, 2013. Available online: https://hdl.handle.net/11296/8re84f (accessed on 10 October 2021).
- Kuo, T.-H. Study on The Automatic Waterfowl Incubator Control Method. Master’s Thesis, Department of Biomechatronic Engineering, National Chiayi University, Chiayi City, Taiwan, 2014. Available online: https://hdl.handle.net/11296/j9mp2w (accessed on 10 October 2021).
- Yang, H.C. Study on Automatic Waterfowl Hatching and Egg-Cooling Control System Intergrated with Human-Machine Interface. Master’s Thesis, Department of Biomechatronic Engineering, National Chiayi University, Chiayi City, Taiwan, 2015. Available online: https://hdl.handle.net/11296/n29e8c (accessed on 10 October 2021).
- Shen, Y.-J. Combining CFD software of Response Surface Methodology to Simulate the Optimal Design of the Incubate Using for Waterfowl Eggs Hatching. Master’s Thesis, Department of Biomechatronic Engineering, National Chiayi University, Chiayi City, Taiwan, 2016. Available online: https://hdl.handle.net/11296/urrgsp (accessed on 10 October 2021).
- Ho, H.-Y. Flow Field Simulation and Analysis of Automatic-Cooling Incubator. Master’s Thesis, Department of Biomechatronic Engineering, National Chiayi University, Chiayi City, Taiwan, 2013. Available online: https://hdl.handle.net/11296/k33rzq (accessed on 10 October 2021).
- McCafferty, D.J. Applications of thermal imaging in avian science. Int. J. Avian Sci. 2013, 155, 4–15. [Google Scholar] [CrossRef]
- Lin, C.S.; Yeh, P.T.; Chen, D.C.; Chiou, Y.C.; Lee, C.H. The identification and filtering of fertilized eggs with a thermal imaging system. Comput. Electron. Agric. 2013, 91, 94–105. [Google Scholar] [CrossRef]
- Sunardi, S.; Yudhana, A.; Saifullah, S. Identity analysis of egg based on digital and thermal imaging: Image processing and counting object concept. Int. J. Electr. Comput. Eng. 2017, 7, 200–208. [Google Scholar] [CrossRef]
- Chenk, T.; Guestrin, C. Xgboost: A scalable tree boosting system. In Proceedings of the 22Nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; ACM: New York, NY, USA, 2016; pp. 785–794. [Google Scholar]
- Mustapha, I.B.; Saeed, F. Bioactive molecule prediction using extreme gradient boosting. Molecules 2016, 21, 983. [Google Scholar] [CrossRef] [PubMed]
- Xia, Y.; Liu, C.; Li, Y.; Liu, N. A boosted decision tree approach using Bayesian hyper-parameter optimization for credit scoring. Expert Syst. Appl. 2017, 78, 225–241. [Google Scholar] [CrossRef]
- Luckner, M.; Topolski, B.; Mazurek, M. Application of XGBoost algorithm in fingerprinting localisation task. In IFIP International Conference on Computer Information Systems and Industrial Management; Springer: Cham, Switzerland, 2017; pp. 661–671. [Google Scholar]
- Duck Production System Manual, Published by the Council of Agriculture (CoA), Executive Yuan in Taiwan. Available online: https://kmweb.coa.gov.tw/subject/subject.php?id=14027 (accessed on 1 April 2020).
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).