• Progress in Superconductivity
• /
• v.7 no.1
• /
• pp.28-35
• /
• 2005

We propose a method to measure atrial arrhythmias (AA) such as atrial fibrillation (Afb) and atrial flutter (Afl) with a SQUID magnetocardiograph (MCG) system. To detect AA is one of challenging topics in MCG. As the AA generally have irregular rhythm and atrio-ventricular conduction, the MCG signal cannot be improved by QRS averaging; therefore a SQUID MCG system having a high SNR is required to measure informative atrial excitation with a single scan. In the case of Afb, diminished f waves are much smaller than normal P waves because the sources are usually located on the posterior wall of the heart. In this study, we utilize an MCG system measuring tangential field components, which is known to be more sensitive to a deeper current source. The average noise spectral density of the whole system in a magnetic shielded room was $10\;fT/{\surd}Hz(a)\;1\;Hz\;and\;5\;fT/{\surd}Hz\;(a)\;100\;Hz$. We measured the MCG signals of patients with chronic Afb and Afl. Before the AA measurement, the comparison between the measurements in supine and prone positions for P waves has been conducted and the experiment gave a result that the supine position is more suitable to measure the atrial excitation. Therefore, the AA was measured in subject's supine position. Clinical potential of AA measurement in MCG is to find an aspect of a reentry circuit and to localize the abnormal stimulation noninvasively. To give useful information about the abnormal excitation, we have developed a method, separative synthetic aperture magnetometry (sSAM). The basic idea of sSAM is to visualize current source distribution corresponding to the atrial excitation, which are separated from the ventricular excitation and the Gaussian sensor noises. By using sSAM, we localized the source of an Afl successfully.

• Progress in Medical Physics
• /
• v.13 no.1
• /
• pp.21-26
• /
• 2002

A practical calculation algorithm which calculates the relative output factor(ROF) for irregular shaped electron field has been developed and evaluated the accuracy of the algorithm. The algorithm adapted two-source model, which assumes that the electron dose can be express as sum of the primary source component and the scattered component from the shielding block. Original two-source model has been modified in order to make the algorithm simpler and to reduce the number of parameters needed in the calculation, while the calculation error remains within clinical tolerance range. The primary source is assumed to have Gaussian distribution, while the scattered component follows the inverse square law. Depth and angular dependency of the primary and the scattered are ignored ROF can be calculated with three parameters such as, the effective source distance, the variance of primary source, and the scattering power of the block. The coefficients are obtained from the square shaped-block measurements and the algorithm is confirmed from the rectangular or irregular shaped-fields used in the clinic. The results showed less than 1.0 ％ difference between the calculation and measurements for most cases. None of cases which have bigger than 2.1 ％ have been found. By improving the algorithm for the aperture region which shows the largest error, the algorithm could be practically used in the clinic, since one can acquire the 1011 parameter's with minimum measurements(5∼6 measurements per cones) and generates accurate results within the clinically acceptable range.

• Journal of Korea Water Resources Association
• /
• v.36 no.4
• /
• pp.533-545
• /
• 2003

This study analyzed several storm events observed in the Seolma-chun basin to derive the characteristic velocity of GIUH (Geomophological Instantaneous Unit Hydrograph) as well as its variability. Especially, this study focused on the variation of characteristic velocity due to the change of rainfall characteristics. The IUH of the Seolma-chun basin was derived using the HEC-1, whose peak discharge and time were then compared with those of the GIUH to derive the characteristic velocities. The characteristics velocities were analyzed by comparing with the GcIUH (Geomorphoclimatic IUH) as well as the characteristics of rainfall. Results are summarized as follows. (1) The characteristic velocity of GIUH was estimated higher with higher variability than the GcIUH, but their trends were found similar (2) Total amount of effective rainfall (or, mean effective rainfall) well explains the characteristic velocity of GIUH. This could be assured by the regression analysis, whose coefficient of determination was estimated about 0.6. (3) The duration and the maximum intensity of rainfall were found not to affect significantly on the characteristic velocity of GIUH. The coefficients of determination were estimated less than 0.3 for all cases considered. (4) For the rainfall events used in this study, the characteristic velocities of GIUH were found to follow the Gaussian distribution with its mean and the standard deviation 0.402 m/s and 0.173 m/s, respectively. Most of the values are within the range of 0.4∼0.5 m/s, and its coefficient of variation was estimated to be 0.43, much less than that of the runoff itself (about 1.0).

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.60 no.3
• /
• pp.639-647
• /
• 2011

In this paper, we introduce a design methodology of data-centroid Radial Basis Function neural networks with extended polynomial function. The two underlying design mechanisms of such networks involve K-means clustering method and Particle Swarm Optimization(PSO). The proposed algorithm is based on K-means clustering method for efficient processing of data and the optimization of model was carried out using PSO. In this paper, as the connection weight of RBF neural networks, we are able to use four types of polynomials such as simplified, linear, quadratic, and modified quadratic. Using K-means clustering, the center values of Gaussian function as activation function are selected. And the PSO-based RBF neural networks results in a structurally optimized structure and comes with a higher level of flexibility than the one encountered in the conventional RBF neural networks. The PSO-based design procedure being applied at each node of RBF neural networks leads to the selection of preferred parameters with specific local characteristics (such as the number of input variables, a specific set of input variables, and the distribution constant value in activation function) available within the RBF neural networks. To evaluate the performance of the proposed data-centroid RBF neural network with extended polynomial function, the model is experimented with using the nonlinear process data(2-Dimensional synthetic data and Mackey-Glass time series process data) and the Machine Learning dataset(NOx emission process data in gas turbine plant, Automobile Miles per Gallon(MPG) data, and Boston housing data). For the characteristic analysis of the given entire dataset with non-linearity as well as the efficient construction and evaluation of the dynamic network model, the partition of the given entire dataset distinguishes between two cases of Division I(training dataset and testing dataset) and Division II(training dataset, validation dataset, and testing dataset). A comparative analysis shows that the proposed RBF neural networks produces model with higher accuracy as well as more superb predictive capability than other intelligent models presented previously.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.199-206
• /
• 2003

Characteristics of the static shift are discussed by comparing the three-dimensional MT inversion with/without static shift parameterization. The galvanic distortion by small-scale shallow feature often leads severe distortion in inverted resistivity structures. The new inversion algorithm is applied to four numerical data sets contaminated by different amount of static shift. In real field data interpretations, we generally do not have any a-priori information about how much the data contains the static shift. In this study, we developed an algorithm for finding both Lagrangian multiplier for smoothness and the trade-off parameter for static shift, simultaneously in 3-D MT inversion. Applications of this inversion routine for the numerical data sets showed quite reasonable estimation of static shift parameters without any a-priori information. The inversion scheme is successfully applied to all the four data sets, even when the static shift does not obey the Gaussian distribution. Allowing the static shift parameters have non-zero degree of freedom to the inversion, we could get more accurate block resistivities as well as static shifts in the data. When inversion does not consider the static shift as inversion parameters (conventional MT inversion), the block resistivities on the surface are modified considerably to match possible static shift. The inhomogeneous blocks on the surface can generate the static shift at low frequencies. By those mechanisms, the conventional 3-D MT inversion can reconstruct the resistivity structures to some extent in the deeper parts even when moderate static shifts are in the data. As frequency increased, however, the galvanic distortion is not frequency independent any more, and thus the conventional inversion failed to fit the apparent resistivity and phase, especially when strong static shift is added. Even in such case, however, reasonable estimation of block resistivity as well as static shift parameters were obtained by 3-D MT inversion with static shift parameterization.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.29 no.2
• /
• pp.139-146
• /
• 2023

Maritime object detection systems, which detects small maritime obstacles such as fish farm buoys and visualizes distance and direction, is equipped with a 3-axis gimbal to compensate for errors caused by hull motion, but there is a limit to distance error corrections necessitated by the vertical movement of the camera and the maritime object due to wave motions. Therefore, in this study, the distance error of maritime object detection systems caused by the movement of the water surface according to the external environment is analyzed and corrected using average filter and moving average filter. Random numbers following a Gaussian standard normal distribution were added to or subtracted from the image coordinates to reproduce the rise or fall of the buoy under irregular waves. The distance calculated according to the change of image coordinates, the predicted distance through the average filter and the moving average filter, and the actual distance measured by laser distance meter were compared. In phases 1 and 2, the error rate increased to a maximum of 98.5% due to the changes of image coordinates due to irregular waves, but the error rate decreased to 16.3% with the moving average filter. This error correction capability was better than with the average filter, but there was a limit due to failure to respond to the distance change. Therefore, it is considered that use of the moving average filter to correct the distance error of the maritime object detection system will enhance responses to the real-time distance change and greatly improve the error rate.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.6
• /
• pp.28-34
• /
• 2006

Laboratory column tests using $Cl^-$ tracer were conducted to study the correlation of soil particle distribution and hydrodynamic dispersion mechanism with three kinds of ununiformed soil samples, in which the ratio of gravel and sand versus silt and clay is 24.5 for S-1 soil, 4.48 for S-2 soil, and 0.4 for S-3 soil. Chloride breakthrough curves with time were fitted with gaussian functions. The relative concentrations of chloride were converged to 1.0 after 0.7 hours for S-1, 6.3 hours for S-2, and 389 hours for S-3. Average linear velocity, longitudinal dispersion coefficient, and longitudinal dispersivity were calculated by chloride breakthrough curves. Longitudinal dispersion coefficients were $1.20{\times}10^{-4}\;m^2/sec$ for S-1, $8.87{\times}10^{-7}\;m^2/sec$ for S-2, and $1.94{\times}10^{-9}\;m^2/sec$ for S-3. Peclet numbers calculated by the molecular diffusion coefficient of chloride and the mean grain diameters of soils were $2.59{\times}10^2$ for S-1, $6.27{\times}10^0$ for S-2, and $1.35{\times}10^{-4}$ for S-3. Mechanical dispersion was dominant for the hydrodynamic dispersion mechanism of S-1. Both mechanical dispersion and molecular diffusion were dominant for the hydrodynamic dispersion mechanism of S-2, but mechanical dispersion was ascendant over molecular diffusion. Hydrodynamic dispersion in S-3 was occurred mainly by molecular diffusion. When plotting three soils on the graph of $D_L/D_m$ versus Peclet number produced by Bijeljic and Blunt (2006), the values of $D_L/D_m$ for S-1 and S-2 were more than 2.0 order compared to their graph. S-3 was not plotted on their graph because the Peclet number was as small as $1.35{\times}10^{-4}$.

• Progress in Medical Physics
• /
• v.23 no.2
• /
• pp.81-90
• /
• 2012

The effect of setup uncertainties on CTV dose and the correlation between setup uncertainties and setup margin were evaluated by Monte Carlo based numerical simulation. Patient specific information of IMRT treatment plan for rectal cancer designed on the VARIAN Eclipse planning system was utilized for the Monte Carlo simulation program including the planned dose distribution and tumor volume information of a rectal cancer patient. The simulation program was developed for the purpose of the study on Linux environment using open source packages, GNU C++ and ROOT data analysis framework. All misalignments of patient setup were assumed to follow the central limit theorem. Thus systematic and random errors were generated according to the gaussian statistics with a given standard deviation as simulation input parameter. After the setup error simulations, the change of dose in CTV volume was analyzed with the simulation result. In order to verify the conventional margin recipe, the correlation between setup error and setup margin was compared with the margin formula developed on three dimensional conformal radiation therapy. The simulation was performed total 2,000 times for each simulation input of systematic and random errors independently. The size of standard deviation for generating patient setup errors was changed from 1 mm to 10 mm with 1 mm step. In case for the systematic error the minimum dose on CTV $D_{min}^{stat{\cdot}}$ was decreased from 100.4 to 72.50% and the mean dose $\bar{D}_{syst{\cdot}}$ was decreased from 100.45% to 97.88%. However the standard deviation of dose distribution in CTV volume was increased from 0.02% to 3.33%. The effect of random error gave the same result of a reduction of mean and minimum dose to CTV volume. It was found that the minimum dose on CTV volume $D_{min}^{rand{\cdot}}$ was reduced from 100.45% to 94.80% and the mean dose to CTV $\bar{D}_{rand{\cdot}}$ was decreased from 100.46% to 97.87%. Like systematic error, the standard deviation of CTV dose ${\Delta}D_{rand}$ was increased from 0.01% to 0.63%. After calculating a size of margin for each systematic and random error the "population ratio" was introduced and applied to verify margin recipe. It was found that the conventional margin formula satisfy margin object on IMRT treatment for rectal cancer. It is considered that the developed Monte-carlo based simulation program might be useful to study for patient setup error and dose coverage in CTV volume due to variations of margin size and setup error.

• The Korean Journal of Nuclear Medicine
• /
• v.38 no.6
• /
• pp.486-491
• /
• 2004

Purpose: factor analysis and independent component analysis (ICA) has been used for handling dynamic image sequences. Theoretical advantages of a newly suggested ICA method, ensemble ICA, leaded us to consider applying this method to the analysis of dynamic myocardial $H_2^{15}O$ PET data. In this study, we quantified patients' blood flow using the ensemble ICA method. Materials and Methods: Twenty subjects underwent $H_2^{15}O$ PET scans using ECAT EXACT 47 scanner and myocardial perfusion SPECT using Vertex scanner. After transmission scanning, dynamic emission scans were initiated simultaneously with the injection of $555{\sim}740$ MBq $H_2^{15}O$. Hidden independent components can be extracted from the observed mixed data (PET image) by means of ICA algorithms. Ensemble learning is a variational Bayesian method that provides an analytical approximation to the parameter posterior using a tractable distribution. Variational approximation forms a lower bound on the ensemble likelihood and the maximization of the lower bound is achieved through minimizing the Kullback-Leibler divergence between the true posterior and the variational posterior. In this study, posterior pdf was approximated by a rectified Gaussian distribution to incorporate non-negativity constraint, which is suitable to dynamic images in nuclear medicine. Blood flow was measured in 9 regions - apex, four areas in mid wall, and four areas in base wall. Myocardial perfusion SPECT score and angiography results were compared with the regional blood flow. Results: Major cardiac components were separated successfully by the ensemble ICA method and blood flow could be estimated in 15 among 20 patients. Mean myocardial blood flow was $1.2{\pm}0.40$ ml/min/g in rest, $1.85{\pm}1.12$ ml/min/g in stress state. Blood flow values obtained by an operator in two different occasion were highly correlated (r=0.99). In myocardium component image, the image contrast between left ventricle and myocardium was 1:2.7 in average. Perfusion reserve was significantly different between the regions with and without stenosis detected by the coronary angiography (P<0.01). In 66 segment with stenosis confirmed by angiography, the segments with reversible perfusion decrease in perfusion SPECT showed lower perfusion reserve values in $H_2^{15}O$ PET. Conclusions: Myocardial blood flow could be estimated using an ICA method with ensemble learning. We suggest that the ensemble ICA incorporating non-negative constraint is a feasible method to handle dynamic image sequence obtained by the nuclear medicine techniques.

• Investigative Magnetic Resonance Imaging
• /
• v.18 no.4
• /
• pp.303-313
• /
• 2014

Purpose : In-vivo conductivity reconstruction using transmit field ($B_1{^+}$) information of MRI was proposed. We assessed the accuracy of conductivity reconstruction in the presence of statistical noise in complex $B_1{^+}$ map and provided a parametric model of the conductivity-to-noise ratio value. Materials and Methods: The $B_1{^+}$ distribution was simulated for a cylindrical phantom model. By adding complex Gaussian noise to the simulated $B_1{^+}$ map, quantitative conductivity estimation error was evaluated. The quantitative evaluation process was repeated over several different parameters such as Larmor frequency, object radius and SNR of $B_1{^+}$ map. A parametric model for the conductivity-to-noise ratio was developed according to these various parameters. Results: According to the simulation results, conductivity estimation is more sensitive to statistical noise in $B_1{^+}$ phase than to noise in $B_1{^+}$ magnitude. The conductivity estimate of the object of interest does not depend on the external object surrounding it. The conductivity-to-noise ratio is proportional to the signal-to-noise ratio of the $B_1{^+}$ map, Larmor frequency, the conductivity value itself and the number of averaged pixels. To estimate accurate conductivity value of the targeted tissue, SNR of $B_1{^+}$ map and adequate filtering size have to be taken into account for conductivity reconstruction process. In addition, the simulation result was verified at 3T conventional MRI scanner. Conclusion: Through all these relationships, quantitative conductivity estimation error due to statistical noise in $B_1{^+}$ map is modeled. By using this model, further issues regarding filtering and reconstruction algorithms can be investigated for MREPT.