• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.329-336
• /
• 2007

For scaling of the gradient of misfit function, we develop a new pseudo-Hessian matrix constructed by combining amplitude field and pseudo-Hessian matrix. Since pseudo- Hessian matrix neglects the calculation of the zero-lag auto-correlation of impulse responses in the approximate Hessian matrix, the pseudo-Hessian matrix has a limitation to scale the gradient of misfit function compared to the approximate Hessian matrix. To validate the new pseudo- Hessian matrix, we perform frequency-domain elastic full waveform inversion using this Hessian matrix. By synthetic experiments, we show that the new pseudo-Hessian matrix can give better convergence to the true model than the old one does. Furthermore, since the amplitude fields are intrinsically obtained in forward modeling procedure, we do not have to pay any extra cost to compute the new pseudo-Hessian. We think that the new pseudo-Hessian matrix can be used as an alternative of the approximate Hessian matrix of the Gauss-Newton method.

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.191-196
• /
• 2004

Seismic waveform inversion can be solved by using the classical Gauss-Newton method, which needs to construct the huge Hessian by the directly computed Jacobian. The property of Hessian mainly depends upon a source and receiver aperture, a velocity model, an illumination Bone and a frequency content of source wavelet. In this paper, we try to invert the Marmousi seismic data by controlling the huge Hessian appearing in the Gauss-Newton method. Wemake the two kinds of he approximate Hessian. One is the banded Hessian and the other is the approximate Hessian with automatic gain function. One is that the 1st updated velocity model from the banded Hessian is nearly the same of the result from the full approximate Hessian. The other is that the stability using the automatic gain function is more improved than that without automatic gain control.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.14 no.4
• /
• pp.22-31
• /
• 1977

A general equation by which the Hessian matrix of an error function can be determined directly, has been derived. It was verified to be useful in optimization processes that include the Hessian matrix. A few design examples had shown that this method had accelerated the processes of finding the minimums. The advantage of this technique is the possibility of optimizing functions that composed of both the phases and magnitudes.

• Journal of the Korean Data and Information Science Society
• /
• v.23 no.5
• /
• pp.1027-1035
• /
• 2012

This paper proposes a new Levenberg-Marquardt algorithm that is accelerated by adjusting a Jacobian matrix and a quasi-Hessian matrix. The proposed method partitions the Jacobian matrix into block matrices and employs the inverse of a partitioned matrix to find the inverse of the quasi-Hessian matrix. Our method can avoid expensive operations and save memory in calculating the inverse of the quasi-Hessian matrix. It can shorten the training time for fast convergence. In our results tested in a large application, we were able to save about 20% of the training time than other algorithms.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.9 no.12
• /
• pp.5073-5086
• /
• 2015

The rock fracture detection by image analysis is significant for fracture measurement and assessment engineering. The paper proposes a novel image segmentation algorithm for the centerline tracing of a rock fracture based on Hessian Matrix at Multi-scales and Steger algorithm. A traditional fracture detection method, which does edge detection first, then makes image binarization, and finally performs noise removal and fracture gap linking, is difficult for images of rough rock surfaces. To overcome the problem, the new algorithm extracts the centerlines directly from a gray level image. It includes three steps: (1) Hessian Matrix and Frangi filter are adopted to enhance the curvilinear structures, then after image binarization, the spurious-fractures and noise are removed by synthesizing the area, circularity and rectangularity; (2) On the binary image, Steger algorithm is used to detect fracture centerline points, then the centerline points or segments are linked according to the gap distance and the angle differences; and (3) Based on the above centerline detection roughly, the centerline points are searched in the original image in a local window along the direction perpendicular to the normal of the centerline, then these points are linked. A number of rock fracture images have been tested, and the testing results show that compared to other traditional algorithms, the proposed algorithm can extract rock fracture centerlines accurately.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.159-162
• /
• 2002

We tried to obtain an initial velocity model for prestack depth migration via waveform inversion. For application of any field data we chose a smooth background layered velocity model (v=v0 + k x z) as an initial velocity model. Newton type waveform inversion needs to invert huge Hessian matrix. In order to compute full Hessian matrix arising from full aperture data and full illumination zone, we meet insurmountable difficulties of paying astronomical computing cost. For the layered media, approximate Hessian emerging from single shot aperture data can be used repeatedly for split spread source configuration. In our work of using this Hessian characteristic of layered media we attempted to obtain the approximate velocity model as close as possible to the true velocity model in first iteration.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.44-48
• /
• 2009

In this study, a vortex visualization method based on the vorticity magnitude is developed. One of the simplest models for a vortex is a vortex filament with the maximum vorticity on its center. The proposed method is based on the observation of this ideal distribution of vorticity magnitude. Laplacian and Hessian matrix of vorticity magnitude are tested for detecting the local maximum of vorticity magnitude. These ideas were applied to wake flow past a sphere. It was found that the Laplacian method is not able to distinguish vortices from the underlying shear layer clearly, while the Hessian matrix method does not suffer from this problem.

• Journal of Korea Multimedia Society
• /
• v.12 no.8
• /
• pp.1082-1088
• /
• 2009

Fundus images are highly useful in evaluating patients' retinal conditions in diagnosing eye diseases. In particular, vessel regions are essential in diagnosing diabetes and hypertension. In this paper, we used top-hat filter to compensate for non-uniform background. Image contrast was enhanced by using contrast limited adaptive histogram equalization (CLAHE) method. Hessian matrix was next applied to detect vessel regions. Results indicate that our method is 1.3% more accurate than matched filter method. Our proposed method is expected to contribute to diagnosing eye diseases.

• The Pure and Applied Mathematics
• /
• v.2 no.2
• /
• pp.133-140
• /
• 1995

The D-optimal design criterion for precise parameter estimation in nonlinear regression analysis is called the determinant criterion because the determinant of a matrix is to be maximized. In this thesis, we derive the gradient and the Hessian of the determinant criterion, and apply a QR decomposition for their efficient computations. We also propose an approximate form of the Hessian matrix which can be calculated from the first derivative of a model function with respect to the design variables. These equations can be used in a Gauss-Newton type iteration procedure.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.1
• /
• pp.47-57
• /
• 2010

Subgrid scale stabilization method is applied to Woo and Liu(2004)'s extended Boussinesq FEM numerical model to eliminate the 2dx wiggles. In order to optimize the computational efficiency, Hessian operator is introduced and the matrix of velocity vector is combined to one matrix for solving matrix equations. The mass lumping technique is also applied to the matrix equations of auxiliary variables. The newly developed code is applied to simulate Vincent and Briggs(1989)' wave transformation experiments and the results show that the numerical solution is almost wiggle-free and it matches very well with experimental data. Due to improvement of computational efficiency and wiggle reduction, it is plausible to apply this model to a realistic problem such as harbor oscillation problems.