Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.129
Sung Ha Park, Chang-Ock Lee, J. Hahn
We propose a variational segmentation model based on statistical information of intensities in an image. The model consists of both a local region-based energy and a global region-based energy in order to handle misclassification which happens in a typical statistical variational model with an assumption that an image is a mixture of two Gaussian distributions. We find local ambiguous regions where misclassification might happen due to a small difference between two Gaussian distributions. Based on statistical information restricted to the local ambiguous regions, we design a local region-based energy in order to reduce the misclassification. We suggest an algorithm to avoid the difficulty of the Euler-Lagrange equations of the proposed variational model.
{"title":"IMAGE SEGMENTATION BASED ON THE STATISTICAL VARIATIONAL FORMULATION USING THE LOCAL REGION INFORMATION","authors":"Sung Ha Park, Chang-Ock Lee, J. Hahn","doi":"10.12941/JKSIAM.2014.18.129","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.129","url":null,"abstract":"We propose a variational segmentation model based on statistical information of intensities in an image. The model consists of both a local region-based energy and a global region-based energy in order to handle misclassification which happens in a typical statistical variational model with an assumption that an image is a mixture of two Gaussian distributions. We find local ambiguous regions where misclassification might happen due to a small difference between two Gaussian distributions. Based on statistical information restricted to the local ambiguous regions, we design a local region-based energy in order to reduce the misclassification. We suggest an algorithm to avoid the difficulty of the Euler-Lagrange equations of the proposed variational model.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"18 1","pages":"129-142"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"66365406","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.167
Munkh-Erdene Ts, E. Lee
In this paper we consider a novel reconstruction method in electrical impedance tomography (EIT) and its application for monitoring and detecting a hydrargyrum (mercury) polluted soil near to the surface of underground. We use electrodes placed on the surface of land to collect the data which provides the relations of voltage and current map and to produce a projected image of interior conductivity distribution onto the surface of land. Here the projected image reconstruction method is used to monitor the pollution in soil underneath the ground without any destruction and any digging into a land.
{"title":"LONG TERM MONITORING OF HYDRARGYRUM POLLUTED SOIL USING PROJECTED IMAGE RECONSTRUCTION IN ELECTRICAL IMPEDANCE TOMOGRAPHY","authors":"Munkh-Erdene Ts, E. Lee","doi":"10.12941/JKSIAM.2014.18.167","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.167","url":null,"abstract":"In this paper we consider a novel reconstruction method in electrical impedance tomography (EIT) and its application for monitoring and detecting a hydrargyrum (mercury) polluted soil near to the surface of underground. We use electrodes placed on the surface of land to collect the data which provides the relations of voltage and current map and to produce a projected image of interior conductivity distribution onto the surface of land. Here the projected image reconstruction method is used to monitor the pollution in soil underneath the ground without any destruction and any digging into a land.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"82 1","pages":"167-180"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87103369","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.181
C. Ahn
In this paper, we propose a robust real-time myocardial border tracking algorithm for echocardiography. Commonly, after an initial contour of LV border is traced at one or two frame from the entire cardiac cycle, LV contour tracking is performed over the remaining frames. Among a variety of tracking techniques, optical flow method is the most widely used for motion estimation of moving objects. However, when echocardiography data is heavily corrupted in some local regions, the errors bring the tracking point out of the endocardial border, resulting in distorted LV contours. This shape distortion often occurs in practice since the data acquisition is affected by ultrasound artifacts, dropout or shadowing phenomena of cardiac walls. The proposed method deals with this shape distortion problem and reflects the motion realistic LV shape by applying global deformation modeled as affine transform partitively to the contour. We partition the tracking points on the contour into a few groups and determine each affine transform governing the motion of the partitioned contour points. To compute the coefficients of each affine transform, we use the least squares method with equality constraints that are given by the relationship between the coefficients and a few contour points showing good tracking results. Many real experiments show that the proposed method supports better performance than existing methods.
{"title":"A LOCALIZED GLOBAL DEFORMATION MODEL TO TRACK MYOCARDIAL MOTION USING ECHOCARDIOGRAPHY","authors":"C. Ahn","doi":"10.12941/JKSIAM.2014.18.181","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.181","url":null,"abstract":"In this paper, we propose a robust real-time myocardial border tracking algorithm for echocardiography. Commonly, after an initial contour of LV border is traced at one or two frame from the entire cardiac cycle, LV contour tracking is performed over the remaining frames. Among a variety of tracking techniques, optical flow method is the most widely used for motion estimation of moving objects. However, when echocardiography data is heavily corrupted in some local regions, the errors bring the tracking point out of the endocardial border, resulting in distorted LV contours. This shape distortion often occurs in practice since the data acquisition is affected by ultrasound artifacts, dropout or shadowing phenomena of cardiac walls. The proposed method deals with this shape distortion problem and reflects the motion realistic LV shape by applying global deformation modeled as affine transform partitively to the contour. We partition the tracking points on the contour into a few groups and determine each affine transform governing the motion of the partitioned contour points. To compute the coefficients of each affine transform, we use the least squares method with equality constraints that are given by the relationship between the coefficients and a few contour points showing good tracking results. Many real experiments show that the proposed method supports better performance than existing methods.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"51 1","pages":"181-192"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"84535117","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.143
Y. Jung
This paper aims to give a quick view on denoising without comprehensive details. Denoising can be understood as removing unwanted parts in signals and images. Noise incorporates intrinsic random fluctuations in the data. Since noise is ubiquitous, denoising methods and models are diverse. Starting from what noise means, we briefly discuss a denoising model as maximum a posteriori estimation and relate it with a variational form or energy model. After that we present a few major branches in image and signal processing; filtering, shrinkage or thresholding, regularization and data adapted methods, although it may not be a general way of classifying denoising methods.
{"title":"A REVIEW ON DENOISING","authors":"Y. Jung","doi":"10.12941/JKSIAM.2014.18.143","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.143","url":null,"abstract":"This paper aims to give a quick view on denoising without comprehensive details. Denoising can be understood as removing unwanted parts in signals and images. Noise incorporates intrinsic random fluctuations in the data. Since noise is ubiquitous, denoising methods and models are diverse. Starting from what noise means, we briefly discuss a denoising model as maximum a posteriori estimation and relate it with a variational form or energy model. After that we present a few major branches in image and signal processing; filtering, shrinkage or thresholding, regularization and data adapted methods, although it may not be a general way of classifying denoising methods.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"43 1","pages":"143-156"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80756032","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.077
Tingting Zhang, T. Bera, E. Woo, J. Seo
Medical imaging techniques have evolved to expand our ability to visualize new contrast information of electrical, optical, and mechanical properties of tissues in the human body using noninvasive measurement methods. In particular, electrical tissue property imaging techniques have received considerable attention for the last few decades since electrical properties of biological tissues and organs change with their physiological functions and pathological states. We can express the electrical tissue properties as the frequency-dependent admittivity, which can be measured in a macroscopic scale by assessing the relation between the timeharmonic electric field and current density. The main issue is to reconstruct spectroscopic admittivity images from 10 Hz to 1 MHz, for example, with reasonably high spatial and temporal resolutions. It requires a solution of a nonlinear inverse problem involving Maxwell’s equations. To solve the inverse problem with practical significance, we need deep knowledge on its mathematical formulation of underlying physical phenomena, implementation of image reconstruction algorithms, and practical limitations associated with the measurement sensitivity, specificity, noise, and data acquisition time. This paper discusses a number of issues in electrical tissue property imaging modalities and their future directions.
{"title":"SPECTROSCOPIC ADMITTIVITY IMAGING OF BIOLOGICAL TISSUES","authors":"Tingting Zhang, T. Bera, E. Woo, J. Seo","doi":"10.12941/JKSIAM.2014.18.077","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.077","url":null,"abstract":"Medical imaging techniques have evolved to expand our ability to visualize new contrast information of electrical, optical, and mechanical properties of tissues in the human body using noninvasive measurement methods. In particular, electrical tissue property imaging techniques have received considerable attention for the last few decades since electrical properties of biological tissues and organs change with their physiological functions and pathological states. We can express the electrical tissue properties as the frequency-dependent admittivity, which can be measured in a macroscopic scale by assessing the relation between the timeharmonic electric field and current density. The main issue is to reconstruct spectroscopic admittivity images from 10 Hz to 1 MHz, for example, with reasonably high spatial and temporal resolutions. It requires a solution of a nonlinear inverse problem involving Maxwell’s equations. To solve the inverse problem with practical significance, we need deep knowledge on its mathematical formulation of underlying physical phenomena, implementation of image reconstruction algorithms, and practical limitations associated with the measurement sensitivity, specificity, noise, and data acquisition time. This paper discusses a number of issues in electrical tissue property imaging modalities and their future directions.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"1 1","pages":"77-105"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90872889","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.157
Kiwan Jeon, S. Kang, C. Ahn, Sungwhan Kim
If there are metals located in the X-ray scanned object, a point outside the metals gas its range of projection angle at which projections passing through the point are disturbed by the metals. Roughly speaking, this implies that attenuation information at the point is missing in the blocked projection range. So conventional projection completion MAR algorithms to use the undisturbed projection data on the boundary of the metaltrace is less efficient in reconstructing the attenuation coefficient in detailed parts, in particular, near the metal region. In order to overcome this problem, we propose the algebraic correction technique (ACT) to utilize a prereconstructed interim image of the attenuation coefficient outside the metal region which is obtained by solving a linear system designed to reduce computational costs. The reconstructed interim image of the attenuation coefficient is used as prior information for MAR. Numerical simulations support that the proposed correction technique shows better performance than conventional inpainting techniques such as the total variation and the harmonic inpainting.
{"title":"ALGEBRAIC CORRECTION FOR METAL ARTIFACT REDUCTION IN COMPUTED TOMOGRAPHY","authors":"Kiwan Jeon, S. Kang, C. Ahn, Sungwhan Kim","doi":"10.12941/JKSIAM.2014.18.157","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.157","url":null,"abstract":"If there are metals located in the X-ray scanned object, a point outside the metals gas its range of projection angle at which projections passing through the point are disturbed by the metals. Roughly speaking, this implies that attenuation information at the point is missing in the blocked projection range. So conventional projection completion MAR algorithms to use the undisturbed projection data on the boundary of the metaltrace is less efficient in reconstructing the attenuation coefficient in detailed parts, in particular, near the metal region. In order to overcome this problem, we propose the algebraic correction technique (ACT) to utilize a prereconstructed interim image of the attenuation coefficient outside the metal region which is obtained by solving a linear system designed to reduce computational costs. The reconstructed interim image of the attenuation coefficient is used as prior information for MAR. Numerical simulations support that the proposed correction technique shows better performance than conventional inpainting techniques such as the total variation and the harmonic inpainting.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"42 1","pages":"157-166"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77626588","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.107
M. Chipot
A BSTRACT . The goal of this note is to describe the asymptotic behaviour of problems set in cylinders when the size of them is becoming infinite. This leads to consider problems in unbounded domains as well as new singular perturbations issues.
{"title":"ℓ GOES TO PLUS INFINITY : AN UPDATE","authors":"M. Chipot","doi":"10.12941/JKSIAM.2014.18.107","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.107","url":null,"abstract":"A BSTRACT . The goal of this note is to describe the asymptotic behaviour of problems set in cylinders when the size of them is becoming infinite. This leads to consider problems in unbounded domains as well as new singular perturbations issues.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"35 1","pages":"107-127"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87508731","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-06-25DOI: 10.12941/JKSIAM.2014.18.193
J. Seo, H. Zorgati
In this paper we explore the emergence of the notion of compactness within its historical beginning through rigor versus intuition modes in the treatment of Dirichlet‘s principle. We emphasize on the intuition in Riemann‘s statement on the principle criticized byWeierstrass‘ requirement of rigor followed by Hilbert‘s restatement again criticized by Hadamard, which pushed the ascension of the notion of compactness in the analysis of PDEs. A brief overview of some techniques and problems involving compactness is presented illustrating the importance of this notion. Compactness is discussed here to raise educational issues regarding rigor vs intuition in mathematical studies. The concept of compactness advanced rapidly afterWeierstrass’s famous criticism of Riemann’s use of the Dirichlet principle. The rigor of Weierstrass contributed to establishment of the concept of compactness, but such a focus on rigor blinded mathematicians to big pictures. Fortunately, Poincar´e and Hilbert defended Riemann’s use of the Dirichlet principle and found a balance between rigor and intuition. There is no theorem without rigor, but we should not be a slave of rigor. Rigor (highly detailed examination with toy models) and intuition (broader view with real models) are essentially complementary to each other.
{"title":"COMPACTNESS AND DIRICHLET’S PRINCIPLE","authors":"J. Seo, H. Zorgati","doi":"10.12941/JKSIAM.2014.18.193","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.193","url":null,"abstract":"In this paper we explore the emergence of the notion of compactness within its historical beginning through rigor versus intuition modes in the treatment of Dirichlet‘s principle. We emphasize on the intuition in Riemann‘s statement on the principle criticized byWeierstrass‘ requirement of rigor followed by Hilbert‘s restatement again criticized by Hadamard, which pushed the ascension of the notion of compactness in the analysis of PDEs. A brief overview of some techniques and problems involving compactness is presented illustrating the importance of this notion. Compactness is discussed here to raise educational issues regarding rigor vs intuition in mathematical studies. The concept of compactness advanced rapidly afterWeierstrass’s famous criticism of Riemann’s use of the Dirichlet principle. The rigor of Weierstrass contributed to establishment of the concept of compactness, but such a focus on rigor blinded mathematicians to big pictures. Fortunately, Poincar´e and Hilbert defended Riemann’s use of the Dirichlet principle and found a balance between rigor and intuition. There is no theorem without rigor, but we should not be a slave of rigor. Rigor (highly detailed examination with toy models) and intuition (broader view with real models) are essentially complementary to each other.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"38 1","pages":"193-207"},"PeriodicalIF":0.6,"publicationDate":"2014-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86091015","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-03-25DOI: 10.12941/JKSIAM.2014.18.061
Seunggyu Lee, Yibao Li, Yongho Choi, H. Hwang, Junseok Kim
This paper presents accurate and efficient numerical methods for calculating the sensitivities of two-asset European options, the Greeks. The Greeks are important financial instruments in management of economic value at risk due to changing market conditions. The option pricing model is based on the Black-Scholes partial differential equation. The model is discretized by using a finite difference method and resulting discrete equations are solved by means of an operator splitting method. For Delta, Gamma, and Theta, we investigate the effect of high-order discretizations. For Rho and Vega, we develop an accurate and robust automatic algorithm for finding an optimal value. A cash-or-nothing option is taken to demonstrate the performance of the proposed algorithm for calculating the Greeks. The results show that the new treatment gives automatic and robust calculations for the Greeks.
{"title":"ACCURATE AND EFFICIENT COMPUTATIONS FOR THE GREEKS OF EUROPEAN MULTI-ASSET OPTIONS","authors":"Seunggyu Lee, Yibao Li, Yongho Choi, H. Hwang, Junseok Kim","doi":"10.12941/JKSIAM.2014.18.061","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.061","url":null,"abstract":"This paper presents accurate and efficient numerical methods for calculating the sensitivities of two-asset European options, the Greeks. The Greeks are important financial instruments in management of economic value at risk due to changing market conditions. The option pricing model is based on the Black-Scholes partial differential equation. The model is discretized by using a finite difference method and resulting discrete equations are solved by means of an operator splitting method. For Delta, Gamma, and Theta, we investigate the effect of high-order discretizations. For Rho and Vega, we develop an accurate and robust automatic algorithm for finding an optimal value. A cash-or-nothing option is taken to demonstrate the performance of the proposed algorithm for calculating the Greeks. The results show that the new treatment gives automatic and robust calculations for the Greeks.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"49 1","pages":"61-74"},"PeriodicalIF":0.6,"publicationDate":"2014-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76801503","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2014-03-25DOI: 10.12941/JKSIAM.2014.18.027
H. Lee, June-Yub Lee
In this paper, we present an efficient numerical method for multiphase image segmentation using a multiphase-field model. The method combines the vector-valued Allen- Cahn phase-field equation with initial data fitting terms containing prescribed interface width and fidelity constants. An efficient numerical solution is achieved using the recently developed hybrid operator splitting method for the vector-valued Allen-Cahn phase-field equation. We split the modified vector-valued Allen-Cahn equation into a nonlinear equation and a linear diffusion equation with a source term. The linear diffusion equation is discretized using an implicit scheme and the resulting implicit discrete system of equations is solved by a multigrid method. The nonlinear equation is solved semi-analytically using a closed-form solution. And by treating the source term of the linear diffusion equation explicitly, we solve the modified vector-valued Allen-Cahn equation in a decoupled way. By decoupling the governing equation, we can speed up the segmentation process with multiple phases. We perform some characteristic numerical experiments for multiphase image segmentation.
{"title":"A NUMERICAL METHOD FOR THE MODIFIED VECTOR-VALUED ALLEN–CAHN PHASE-FIELD MODEL AND ITS APPLICATION TO MULTIPHASE IMAGE SEGMENTATION","authors":"H. Lee, June-Yub Lee","doi":"10.12941/JKSIAM.2014.18.027","DOIUrl":"https://doi.org/10.12941/JKSIAM.2014.18.027","url":null,"abstract":"In this paper, we present an efficient numerical method for multiphase image segmentation using a multiphase-field model. The method combines the vector-valued Allen- Cahn phase-field equation with initial data fitting terms containing prescribed interface width and fidelity constants. An efficient numerical solution is achieved using the recently developed hybrid operator splitting method for the vector-valued Allen-Cahn phase-field equation. We split the modified vector-valued Allen-Cahn equation into a nonlinear equation and a linear diffusion equation with a source term. The linear diffusion equation is discretized using an implicit scheme and the resulting implicit discrete system of equations is solved by a multigrid method. The nonlinear equation is solved semi-analytically using a closed-form solution. And by treating the source term of the linear diffusion equation explicitly, we solve the modified vector-valued Allen-Cahn equation in a decoupled way. By decoupling the governing equation, we can speed up the segmentation process with multiple phases. We perform some characteristic numerical experiments for multiphase image segmentation.","PeriodicalId":41717,"journal":{"name":"Journal of the Korean Society for Industrial and Applied Mathematics","volume":"56 1","pages":"27-41"},"PeriodicalIF":0.6,"publicationDate":"2014-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81490466","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}