• Title/Summary/Keyword: 계산 복잡도

Search Result 2,579, Processing Time 0.038 seconds

Analysis of Uncertainty in Ocean Color Products by Water Vapor Vertical Profile (수증기 연직 분포에 의한 GOCI-II 해색 산출물 오차 분석)

  • Kyeong-Sang Lee;Sujung Bae;Eunkyung Lee;Jae-Hyun Ahn
    • Korean Journal of Remote Sensing
    • /
    • v.39 no.6_2
    • /
    • pp.1591-1604
    • /
    • 2023
  • In ocean color remote sensing, atmospheric correction is a vital process for ensuring the accuracy and reliability of ocean color products. Furthermore, in recent years, the remote sensing community has intensified its requirements for understanding errors in satellite data. Accordingly, research is currently addressing errors in remote sensing reflectance (Rrs) resulting from inaccuracies in meteorological variables (total ozone, pressure, wind field, and total precipitable water) used as auxiliary data for atmospheric correction. However, there has been no investigation into the error in Rrs caused by the variability of the water vapor profile, despite it being a recognized error source. In this study, we used the Second Simulation of a Satellite Signal Vector version 2.1 simulation to compute errors in water vapor transmittance arising from variations in the water vapor profile within the GOCI-II observation area. Subsequently, we conducted an analysis of the associated errors in ocean color products. The observed water vapor profile not only exhibited a complex shape but also showed significant variations near the surface, leading to differences of up to 0.007 compared to the US standard 62 water vapor profile used in the GOCI-II atmospheric correction. The resulting variation in water vapor transmittance led to a difference in aerosol reflectance estimation, consequently introducing errors in Rrs across all GOCI-II bands. However, the error of Rrs in the 412-555 nm due to the difference in the water vapor profile band was found to be below 2%, which is lower than the required accuracy. Also, similar errors were shown in other ocean color products such as chlorophyll-a concentration, colored dissolved organic matter, and total suspended matter concentration. The results of this study indicate that the variability in water vapor profiles has minimal impact on the accuracy of atmospheric correction and ocean color products. Therefore, improving the accuracy of the input data related to the water vapor column concentration is even more critical for enhancing the accuracy of ocean color products in terms of water vapor absorption correction.

Estimation of Fractional Urban Tree Canopy Cover through Machine Learning Using Optical Satellite Images (기계학습을 이용한 광학 위성 영상 기반의 도시 내 수목 피복률 추정)

  • Sejeong Bae ;Bokyung Son ;Taejun Sung ;Yeonsu Lee ;Jungho Im ;Yoojin Kang
    • Korean Journal of Remote Sensing
    • /
    • v.39 no.5_3
    • /
    • pp.1009-1029
    • /
    • 2023
  • Urban trees play a vital role in urban ecosystems,significantly reducing impervious surfaces and impacting carbon cycling within the city. Although previous research has demonstrated the efficacy of employing artificial intelligence in conjunction with airborne light detection and ranging (LiDAR) data to generate urban tree information, the availability and cost constraints associated with LiDAR data pose limitations. Consequently, this study employed freely accessible, high-resolution multispectral satellite imagery (i.e., Sentinel-2 data) to estimate fractional tree canopy cover (FTC) within the urban confines of Suwon, South Korea, employing machine learning techniques. This study leveraged a median composite image derived from a time series of Sentinel-2 images. In order to account for the diverse land cover found in urban areas, the model incorporated three types of input variables: average (mean) and standard deviation (std) values within a 30-meter grid from 10 m resolution of optical indices from Sentinel-2, and fractional coverage for distinct land cover classes within 30 m grids from the existing level 3 land cover map. Four schemes with different combinations of input variables were compared. Notably, when all three factors (i.e., mean, std, and fractional cover) were used to consider the variation of landcover in urban areas(Scheme 4, S4), the machine learning model exhibited improved performance compared to using only the mean of optical indices (Scheme 1). Of the various models proposed, the random forest (RF) model with S4 demonstrated the most remarkable performance, achieving R2 of 0.8196, and mean absolute error (MAE) of 0.0749, and a root mean squared error (RMSE) of 0.1022. The std variable exhibited the highest impact on model outputs within the heterogeneous land covers based on the variable importance analysis. This trained RF model with S4 was then applied to the entire Suwon region, consistently delivering robust results with an R2 of 0.8702, MAE of 0.0873, and RMSE of 0.1335. The FTC estimation method developed in this study is expected to offer advantages for application in various regions, providing fundamental data for a better understanding of carbon dynamics in urban ecosystems in the future.

IR Study on the Adsorption of Carbon Monoxide on Silica Supported Ruthenium-Nickel Alloy (실리카 지지 루테늄-니켈 합금에 있어서 일산화탄소의 흡착에 관한 IR 연구)

  • Park, Sang-Youn;Yoon, Dong-Wook
    • Applied Chemistry for Engineering
    • /
    • v.17 no.4
    • /
    • pp.349-356
    • /
    • 2006
  • We have investigated adsorption and desorption properties of CO adsorption on silica supported Ru/Ni alloys at various Ru/Ni mole content ratio as well as CO partial pressures using Fourier transform infrared spectrometer (FT-IR). For Ru-$SiO_{2}$ sample, four bands were observed at $2080.0cm^{-1}$, $2021.0{\sim}2030.7cm^{-1}$, $1778.9{\sim}1799.3cm^{-1}$, $1623.8cm^{-1}$ on adsorption and three bands were observed at $2138.7cm^{-1}$, $2069.3cm^{-1}$, $1988.3{\sim}2030.7cm^{-1}$ on vacumn desorption. For Ni-$SiO_{2}$ sample, four bands were observed at $2057.7cm^{-1}$, $2019.1{\sim}2040.3cm^{-1}$, $1862.9{\sim}1868.7cm^{-1}$, $1625.7cm^{-1}$ on adsorption and two bands were observed at $2009.5{\sim}2040.3cm^{-1}$, $1828.4{\sim}1868.7cm^{-1}$ on vacumn desorption. These absorption bands correspond with those of the previous reports approximately. For Ru/Ni(9/1, 8/2, 7/3, 6/4, 5/5; mole content ratio)-$SiO_{2}$ samples, three bands were observed at $2001.8{\sim}2057.7cm^{-1}$, $1812.8{\sim}1926.5cm^{-1}$, $1623.8{\sim}1625.7cm^{-1}$ on adsorption and three bands were observed at $2140.6cm^{-1}$, $2073.1cm^{-1}$, $1969.0{\sim}2057.7cm^{-1}$ on vacumn desorption. The spectrum pattern observed for Ru/Ni-$SiO_{2}$ sample at 9/1 Ru/Ni mole content ratio on CO adsorption and on vacumn desorption is almost like the spectrum pattern observed for Ru-$SiO_{2}$ sample. But the spectrum patterns observed for Ru/Ni-$SiO_{2}$ samples under 8/2 Ru/Ni mole content ratio on CO adsorption and vacumn desorption are almost like the pattern observed for $Ni-SiO_{2}$ sample. It may be suggested surfaces of alloy clusters on the Ru/Ni-$SiO_{2}$ samples contain more Ni components than the mole content ratio of the sample considering the above phenomena. With Ru/Ni-$SiO_{2}$ samples the absorption band shifts may be ascribed to variations of surface concentration, strain variation due to atomic size difference, variation of bonding energy and electronic densities, and changes of surface geometries according to surface concentration variation. Studies for CO adsorption on Ru/Ni alloy cluster surface by LEED and Auger spectroscopy, interation between Ru/Ni alloy cluster and $SiO_{2}$, and MO calculation for the system would be needed to look into the phenomena.

Dosimetric Effect on Selectable Optimization Parameters of Volumatric Modulated Arc Therapy (선택적 최적화 변수(Selectable Optimization Parameters)에 따른 부피적조절회전방사선치료(VMAT)의 선량학적 영향)

  • Jung, Jae-Yong;Shin, Yong-Joo;Sohn, Seung-Chang;Kim, Yeon-Rae;Min, Jung-Wan;Suh, Tae-Suk
    • Progress in Medical Physics
    • /
    • v.23 no.1
    • /
    • pp.15-25
    • /
    • 2012
  • The aim of this study is to evaluate plan quality and dose accuracy for Volumetric Modulated Arc Therapy (VMAT) on the TG-119 and is to investigate the effects on variation of the selectable optimization parameters of VMAT. VMAT treatment planning was implemented on a Varian iX linear accelerator with ARIA record and verify system (Varian Mecical System Palo Alto, CA) and Oncentra MasterPlan treatment planning system (Nucletron BV, Veenendaal, Netherlands). Plan quality and dosimetric accuracy were evaluated by effect of varying a number of arc, gantry spacing and delivery time for the test geometries provided in TG-119. Plan quality for the target and OAR was evaluated by the mean value and the standard deviation of the Dose Volume Histograms (DVHs). The ionization chamber and $Delta^{4PT}$ bi-planar diode array were used for the dose evaluation. For treatment planning evaluation, all structure sets closed to the goals in the case of single arc, except for the C-shape (hard), and all structure sets achieved the goals in the case of dual arc, except for C-shape (hard). For the variation of a number of arc, the simple structure such as a prostate did not have the difference between single arc and dual arc, whereas the complex structure such as a head and neck showed a superior result in the case of dual arc. The dose distribution with gantry spacing of $4^{\circ}$ was shown better plan quality than the gantry spacing of $6^{\circ}$, but was similar results compared with gantry spacing of $2^{\circ}$. For the verification of dose accuracy with single arc and dual arc, the mean value of a relative error between measured and calculated value were within 3% and 4% for point dose and confidence limit values, respectively. For the verification on dose accuracy with the gantry intervals of $2^{\circ}$, $4^{\circ}$ and $6^{\circ}$, the mean values of relative error were within 3% and 5% for point dose and confidence limit values, respectively. In the verification of dose distribution with $Delta^{4PT}$ bi-planar diode array, gamma passing rate was $98.72{\pm}1.52%$ and $98.3{\pm}1.5%$ for single arc and dual arc, respectively. The confidence limit values were within 4%. The smaller the gantry spacing, the more accuracy results were shown. In this study, we performed the VMAT QA based on TG-119 procedure, and demonstrated that all structure sets were satisfied with acceptance criteria. And also, the results for the selective optimization variables informed the importance of selection for the suitable variables according to the clinical cases.

On the vibration influence to the running power plant facilities when the foundation excavated of the cautious blasting works. (노천굴착에서 발파진동의 크기를 감량 시키기 위한 정밀파실험식)

  • Huh Ginn
    • Explosives and Blasting
    • /
    • v.9 no.1
    • /
    • pp.3-13
    • /
    • 1991
  • The cautious blasting works had been used with emulsion explosion electric M/S delay caps. Drill depth was from 3m to 6m with Crawler Drill ${\phi}70mm$ on the calcalious sand stone (soft -modelate -semi hard Rock). The total numbers of test blast were 88. Scale distance were induced 15.52-60.32. It was applied to propagation Law in blasting vibration as follows. Propagtion Law in Blasting Vibration $V=K(\frac{D}{W^b})^n$ were V : Peak partical velocity(cm/sec) D : Distance between explosion and recording sites(m) W : Maximum charge per delay-period of eight milliseconds or more (kg) K : Ground transmission constant, empirically determind on the Rocks, Explosive and drilling pattern ets. b : Charge exponents n : Reduced exponents where the quantity $\frac{D}{W^b}$ is known as the scale distance. Above equation is worked by the U.S Bureau of Mines to determine peak particle velocity. The propagation Law can be catagorized in three groups. Cubic root Scaling charge per delay Square root Scaling of charge per delay Site-specific Scaling of charge Per delay Plots of peak particle velocity versus distoance were made on log-log coordinates. The data are grouped by test and P.P.V. The linear grouping of the data permits their representation by an equation of the form ; $V=K(\frac{D}{W^{\frac{1}{3}})^{-n}$ The value of K(41 or 124) and n(1.41 or 1.66) were determined for each set of data by the method of least squores. Statistical tests showed that a common slope, n, could be used for all data of a given components. Charge and reduction exponents carried out by multiple regressional analysis. It's divided into under loom over loom distance because the frequency is verified by the distance from blast site. Empirical equation of cautious blasting vibration is as follows. Over 30m ------- under l00m ${\cdots\cdots\cdots}{\;}41(D/sqrt[2]{W})^{-1.41}{\;}{\cdots\cdots\cdots\cdots\cdots}{\;}A$ Over 100m ${\cdots\cdots\cdots\cdots\cdots}{\;}121(D/sqrt[3]{W})^{-1.66}{\;}{\cdots\cdots\cdots\cdots\cdots}{\;}B$ where ; V is peak particle velocity In cm / sec D is distance in m and W, maximLlm charge weight per day in kg K value on the above equation has to be more specified for further understaring about the effect of explosives, Rock strength. And Drilling pattern on the vibration levels, it is necessary to carry out more tests.

  • PDF

Association Between Temporomandibular Disorders and Cervical Muscle Pressure Pain (측두하악장애와 경부근육 압통 간의 상관성)

  • Im, Yeong-Gwan;Kim, Jae-Hyeong;Kim, Byung-Gook
    • Journal of Oral Medicine and Pain
    • /
    • v.33 no.4
    • /
    • pp.339-352
    • /
    • 2008
  • Aims: The aims of this study were to identify the association between cervical muscle pain and TMD by pressure pain response, and to find cervical muscles showing moderate to severe pressure pain that are correlated with masticatory muscle pain. Methods: Patients(n=129, female 65.9%, mean age 28.8 years) answered a TMD questionnaire asking about headache, neck pain, emotional stress, sleep disturbance, parafunction habits, and pain intensity. A clinical examination of the masticatory system was performed. Of the neck muscles, (1) the upper sternocleidomastoid, (2) the middle sternocleidomastoid, (3) the upper trapezius, (4) the splenius capitis, (5) the semispinalis capitis, (6) the scalene medius, and (7) the levator scapulae muscles were examined by palpation. Pressure pain or tenderness of all palpation sites was scored from 0 to 3 according to the pain response. The variables of sum of pressure pain scores were calculated from pressure pain scores and were used for statistical analyses. Results: Eighty patients(62.0%) answered that they suffer from neck pain in the TMD questionnaire. More than 40% of sternocleidomastoid and upper trapezius examination sites showed moderate to severe tenderness in the cervical muscles, and 36% of middle masseter in the masticatory muscles. For the 129 patients, the sum of cervical muscle pain scores(mean=12.88, SD=8.06) and the sum of TMD pain scores(mean=5.36, SD=5.10) were moderately correlated($\rho$ = 0.502, P < 0.001). The sum of TMD pain scores tends to increase as the sum of cervical muscle pain scores increases(Y = 0.395${\cdot}$X, $R^2$ = 0.659, P < 0.001). In the patients with masticatory muscle disorders, the sum of sternocleidomastoid and upper trapezius pain scores(mean = 8.67, SD = 4.95) and the sum of temporalis and masseter pain scores(mean = 3.37, SD = 3.56) showed moderate correlation($\rho$ = 0.375, P < 0.001). Those two variables were in a proportionate relationship(Y = 0.359${\cdot}$X, $R^2$ = 0.538, P < 0.001). In a partial correlation analysis of the sum of unilateral pain scores, the sum of right cervical muscle pain scores and the sum of left cervical muscle pain scores showed the highest correlation(r = 0.802, P < 0.001). The sum of right TMD pain scores and the sum of left TMD pain scores were moderately correlated(r = 0.481, P < 0.001). For the twenty patients with unilateral TMD pain, the partial correlation coefficient between the sum of ipsilateral cervical muscle pain scores and the sum of contralateral cervical muscle pain scores was the largest(r = 0.597, P = 0.009). A partial correlation between the sum of primary TMD side pain scores and the sum of ipsilateral cervical muscle pain scores was 0.564(P = 0.015). Conclusions: TMD is associated with cervical muscle pain on condition of pressure pain response to palpation. Of the cervical muscles, sternocleidomastoid and upper trapezius frequently exhibit moderate to severe pressure pain, and they are closely related to the masticatory muscle pain. The characteristic of symmetric involvement of pain is prominent in cervical muscles; however, TMD can affect the level of cervical muscle pain to modify its symmetric nature.

Comparative analysis of Glomerular Filtration Rate measurement and estimated glomerular filtration rate using 99mTc-DTPA in kidney transplant donors. (신장이식 공여자에서 99mTc-DTPA를 이용한 Glomerular Filtration Rate 측정과 추정사구체여과율의 비교분석)

  • Cheon, Jun Hong;Yoo, Nam Ho;Lee, Sun Ho
    • The Korean Journal of Nuclear Medicine Technology
    • /
    • v.25 no.2
    • /
    • pp.35-40
    • /
    • 2021
  • Purpose Glomerular filtration rate(GFR) is an important indicator for the diagnosis, treatment, and follow-up of kidney disease and is also used by healthy individuals for drug use and evaluating kidney function in donors. The gold standard method of the GFR test is to measure by continuously injecting the inulin which is extrinsic marker, but it takes a long time and the test method is complicated. so, the method of measuring the serum concentration of creatinine is used. Estimated glomerular filtration rate (eGFR) is used instead. However, creatinine is known to be affected by age, gender, muscle mass, etc. eGFR formulas that are currently used include the Cockroft-Gault formula, the modification of diet in renal disease (MDRD) formula, and the chronic kidney disease epidemilogy collaboration (CKD-EPI) formula for adults. For children, the Schwartz formula is used. Measurement of GFR using 51Cr-EDTA (diethylenetriamine tetraacetic acid), 99mTc-DTPA (diethylenetriamine pentaacetic acid) can replace inulin and is currently in use. Therefore, We compared the GFR measured using 99mTc-DTPA with the eGFR using CKD-EPI formula. Materials and Methods For 200 kidney transplant donors who visited Asan medical center.(96 males, 104 females, 47.3 years ± 12.7 years old) GFR was measured using plasma(Two-plasma-sample-method, TPSM) obtained by intravenous administration of 99mTc-DTPA(0.5mCi, 18.5 MBq). eGFR was derived using CKD-EPI formula based on serum creatinine concentration. Results GFR average measured using 99mTc-DTPA for 200 kidney transplant donors is 97.27±19.46(ml/min/1.73m2), and the eGFR average value using the CKD-EPI formula is 96.84±17.74(ml/min/1.73m2), The concentration of serum creatinine is 0.84±0.39(mg/dL). Regression formula of 99mTc-DTPA GFR for serum creatinine-based eGFR was Y = 0.5073X + 48.186, and the correlation coefficient was 0.698 (P<0.01). Difference (%) was 1.52±18.28. Conclusion The correlation coefficient between the 99mTc-DTPA and the eGFR derived on serum creatinine concentration was confirmed to be moderate. This is estimated that eGFR is affected by external factors such as age, gender, and muscle mass and use of formulas made for kidney disease patients. By using 99mTc-DTPA, we can provide reliable GFR results, which is used for diagnosis, treatment and observation of kidney disease, and kidney evaluation of kidney transplant patients.

Estimation of GARCH Models and Performance Analysis of Volatility Trading System using Support Vector Regression (Support Vector Regression을 이용한 GARCH 모형의 추정과 투자전략의 성과분석)

  • Kim, Sun Woong;Choi, Heung Sik
    • Journal of Intelligence and Information Systems
    • /
    • v.23 no.2
    • /
    • pp.107-122
    • /
    • 2017
  • Volatility in the stock market returns is a measure of investment risk. It plays a central role in portfolio optimization, asset pricing and risk management as well as most theoretical financial models. Engle(1982) presented a pioneering paper on the stock market volatility that explains the time-variant characteristics embedded in the stock market return volatility. His model, Autoregressive Conditional Heteroscedasticity (ARCH), was generalized by Bollerslev(1986) as GARCH models. Empirical studies have shown that GARCH models describes well the fat-tailed return distributions and volatility clustering phenomenon appearing in stock prices. The parameters of the GARCH models are generally estimated by the maximum likelihood estimation (MLE) based on the standard normal density. But, since 1987 Black Monday, the stock market prices have become very complex and shown a lot of noisy terms. Recent studies start to apply artificial intelligent approach in estimating the GARCH parameters as a substitute for the MLE. The paper presents SVR-based GARCH process and compares with MLE-based GARCH process to estimate the parameters of GARCH models which are known to well forecast stock market volatility. Kernel functions used in SVR estimation process are linear, polynomial and radial. We analyzed the suggested models with KOSPI 200 Index. This index is constituted by 200 blue chip stocks listed in the Korea Exchange. We sampled KOSPI 200 daily closing values from 2010 to 2015. Sample observations are 1487 days. We used 1187 days to train the suggested GARCH models and the remaining 300 days were used as testing data. First, symmetric and asymmetric GARCH models are estimated by MLE. We forecasted KOSPI 200 Index return volatility and the statistical metric MSE shows better results for the asymmetric GARCH models such as E-GARCH or GJR-GARCH. This is consistent with the documented non-normal return distribution characteristics with fat-tail and leptokurtosis. Compared with MLE estimation process, SVR-based GARCH models outperform the MLE methodology in KOSPI 200 Index return volatility forecasting. Polynomial kernel function shows exceptionally lower forecasting accuracy. We suggested Intelligent Volatility Trading System (IVTS) that utilizes the forecasted volatility results. IVTS entry rules are as follows. If forecasted tomorrow volatility will increase then buy volatility today. If forecasted tomorrow volatility will decrease then sell volatility today. If forecasted volatility direction does not change we hold the existing buy or sell positions. IVTS is assumed to buy and sell historical volatility values. This is somewhat unreal because we cannot trade historical volatility values themselves. But our simulation results are meaningful since the Korea Exchange introduced volatility futures contract that traders can trade since November 2014. The trading systems with SVR-based GARCH models show higher returns than MLE-based GARCH in the testing period. And trading profitable percentages of MLE-based GARCH IVTS models range from 47.5% to 50.0%, trading profitable percentages of SVR-based GARCH IVTS models range from 51.8% to 59.7%. MLE-based symmetric S-GARCH shows +150.2% return and SVR-based symmetric S-GARCH shows +526.4% return. MLE-based asymmetric E-GARCH shows -72% return and SVR-based asymmetric E-GARCH shows +245.6% return. MLE-based asymmetric GJR-GARCH shows -98.7% return and SVR-based asymmetric GJR-GARCH shows +126.3% return. Linear kernel function shows higher trading returns than radial kernel function. Best performance of SVR-based IVTS is +526.4% and that of MLE-based IVTS is +150.2%. SVR-based GARCH IVTS shows higher trading frequency. This study has some limitations. Our models are solely based on SVR. Other artificial intelligence models are needed to search for better performance. We do not consider costs incurred in the trading process including brokerage commissions and slippage costs. IVTS trading performance is unreal since we use historical volatility values as trading objects. The exact forecasting of stock market volatility is essential in the real trading as well as asset pricing models. Further studies on other machine learning-based GARCH models can give better information for the stock market investors.

A Study on the Forest Yield Regulation by Systems Analysis (시스템분석(分析)에 의(依)한 삼림수확조절(森林收穫調節)에 관(關)한 연구(硏究))

  • Cho, Eung-hyouk
    • Korean Journal of Agricultural Science
    • /
    • v.4 no.2
    • /
    • pp.344-390
    • /
    • 1977
  • The purpose of this paper was to schedule optimum cutting strategy which could maximize the total yield under certain restrictions on periodic timber removals and harvest areas from an industrial forest, based on a linear programming technique. Sensitivity of the regulation model to variations in restrictions has also been analyzed to get information on the changes of total yield in the planning period. The regulation procedure has been made on the experimental forest of the Agricultural College of Seoul National University. The forest is composed of 219 cutting units, and characterized by younger age group which is very common in Korea. The planning period is devided into 10 cutting periods of five years each, and cutting is permissible only on the stands of age groups 5-9. It is also assumed in the study that the subsequent forests are established immediately after cutting existing forests, non-stocked forest lands are planted in first cutting period, and established forests are fully stocked until next harvest. All feasible cutting regimes have been defined to each unit depending on their age groups. Total yield (Vi, k) of each regime expected in the planning period has been projected using stand yield tables and forest inventory data, and the regime which gives highest Vi, k has been selected as a optimum cutting regime. After calculating periodic yields and cutting areas, and total yield from the optimum regimes selected without any restrictions, the upper and lower limits of periodic yields(Vj-max, Vj-min) and those of periodic cutting areas (Aj-max, Aj-min) have been decided. The optimum regimes under such restrictions have been selected by linear programming. The results of the study may be summarized as follows:- 1. The fluctuations of periodic harvest yields and areas under cutting regimes selected without restrictions were very great, because of irregular composition of age classes and growing stocks of existing stands. About 68.8 percent of total yield is expected in period 10, while none of yield in periods 6 and 7. 2. After inspection of the above solution, restricted optimum cutting regimes were obtained under the restrictions of Amin=150 ha, Amax=400ha, $Vmin=5,000m^3$ and $Vmax=50,000m^3$, using LP regulation model. As a result, about $50,000m^3$ of stable harvest yield per period and a relatively balanced age group distribution is expected from period 5. In this case, the loss in total yield was about 29 percent of that of unrestricted regimes. 3. Thinning schedule could be easily treated by the model presented in the study, and the thinnings made it possible to select optimum regimes which might be effective for smoothing the wood flows, not to speak of increasing total yield in the planning period. 4. It was known that the stronger the restrictions becomes in the optimum solution the earlier the period comes in which balanced harvest yields and age group distribution can be formed. There was also a tendency in this particular case that the periodic yields were strongly affected by constraints, and the fluctuations of harvest areas depended upon the amount of periodic yields. 5. Because the total yield was decreased at the increasing rate with imposing stronger restrictions, the Joss would be very great where strict sustained yield and normal age group distribution are required in the earlier periods. 6. Total yield under the same restrictions in a period was increased by lowering the felling age and extending the range of cutting age groups. Therefore, it seemed to be advantageous for producing maximum timber yield to adopt wider range of cutting age groups with the lower limit at which the smallest utilization size of timber could be produced. 7. The LP regulation model presented in the study seemed to be useful in the Korean situation from the following point of view: (1) The model can provide forest managers with the solution of where, when, and how much to cut in order to best fulfill the owners objective. (2) Planning is visualized as a continuous process where new strateges are automatically evolved as changes in the forest environment are recognized. (3) The cost (measured as decrease in total yield) of imposing restrictions can be easily evaluated. (4) Thinning schedule can be treated without difficulty. (5) The model can be applied to irregular forests. (6) Traditional regulation methods can be rainforced by the model.

  • PDF