A Conjugate Direction Method for Geophysical Inversion Problems

Mark S. Frank, Constantine Balanis

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

In geophysical tomography, algebraic methods are often used to linearize the nonlinear problem of determining the characteristics of an underground region given measurements of the earth's attenuation to electromagnetic or seismic waves. In this way,a of linear equations is developed such that the unknowns are the picture elements (pixels) of the region being scanned. Classically, these linear equations have been solved using the algebraic reconstructiontechnique (ART) algorithm. In this paper, a new algorithm that is a member of the set of conjugate direction (CD) methods is developed and comparisons are made between this algorithm and the ART algorithm for data arising from simulated electromagnetic probing. This new method, which we call the constrained conjugate gradient (CCG) algorithm, is shown to have a much faster convergence to a final solution than the ART algorithm. In addition, for applications involving high-contrast anomalies (for example, tunnel detection) the CCG is shown to have superior performance in locating the anomalous region for almost all test cases considered.

Original languageEnglish (US)
Pages (from-to)691-701
Number of pages11
JournalIEEE Transactions on Geoscience and Remote Sensing
VolumeGE-25
Issue number6
DOIs
StatePublished - 1987

Fingerprint

Linear equations
Seismic waves
electromagnetic wave
seismic wave
Electromagnetic waves
tomography
Tomography
inversion
method
Tunnels
pixel
tunnel
Pixels
Earth (planet)
anomaly
detection
comparison
test

ASJC Scopus subject areas

  • Earth and Planetary Sciences(all)
  • Electrical and Electronic Engineering

Cite this

A Conjugate Direction Method for Geophysical Inversion Problems. / Frank, Mark S.; Balanis, Constantine.

In: IEEE Transactions on Geoscience and Remote Sensing, Vol. GE-25, No. 6, 1987, p. 691-701.

Research output: Contribution to journalArticle

@article{d8b24114fb654c8c829f02ab723d7127,
title = "A Conjugate Direction Method for Geophysical Inversion Problems",
abstract = "In geophysical tomography, algebraic methods are often used to linearize the nonlinear problem of determining the characteristics of an underground region given measurements of the earth's attenuation to electromagnetic or seismic waves. In this way,a of linear equations is developed such that the unknowns are the picture elements (pixels) of the region being scanned. Classically, these linear equations have been solved using the algebraic reconstructiontechnique (ART) algorithm. In this paper, a new algorithm that is a member of the set of conjugate direction (CD) methods is developed and comparisons are made between this algorithm and the ART algorithm for data arising from simulated electromagnetic probing. This new method, which we call the constrained conjugate gradient (CCG) algorithm, is shown to have a much faster convergence to a final solution than the ART algorithm. In addition, for applications involving high-contrast anomalies (for example, tunnel detection) the CCG is shown to have superior performance in locating the anomalous region for almost all test cases considered.",
author = "Frank, {Mark S.} and Constantine Balanis",
year = "1987",
doi = "10.1109/TGRS.1987.289738",
language = "English (US)",
volume = "GE-25",
pages = "691--701",
journal = "IEEE Transactions on Geoscience and Remote Sensing",
issn = "0196-2892",
publisher = "Institute of Electrical and Electronics Engineers Inc.",
number = "6",

}

TY - JOUR

T1 - A Conjugate Direction Method for Geophysical Inversion Problems

AU - Frank, Mark S.

AU - Balanis, Constantine

PY - 1987

Y1 - 1987

N2 - In geophysical tomography, algebraic methods are often used to linearize the nonlinear problem of determining the characteristics of an underground region given measurements of the earth's attenuation to electromagnetic or seismic waves. In this way,a of linear equations is developed such that the unknowns are the picture elements (pixels) of the region being scanned. Classically, these linear equations have been solved using the algebraic reconstructiontechnique (ART) algorithm. In this paper, a new algorithm that is a member of the set of conjugate direction (CD) methods is developed and comparisons are made between this algorithm and the ART algorithm for data arising from simulated electromagnetic probing. This new method, which we call the constrained conjugate gradient (CCG) algorithm, is shown to have a much faster convergence to a final solution than the ART algorithm. In addition, for applications involving high-contrast anomalies (for example, tunnel detection) the CCG is shown to have superior performance in locating the anomalous region for almost all test cases considered.

AB - In geophysical tomography, algebraic methods are often used to linearize the nonlinear problem of determining the characteristics of an underground region given measurements of the earth's attenuation to electromagnetic or seismic waves. In this way,a of linear equations is developed such that the unknowns are the picture elements (pixels) of the region being scanned. Classically, these linear equations have been solved using the algebraic reconstructiontechnique (ART) algorithm. In this paper, a new algorithm that is a member of the set of conjugate direction (CD) methods is developed and comparisons are made between this algorithm and the ART algorithm for data arising from simulated electromagnetic probing. This new method, which we call the constrained conjugate gradient (CCG) algorithm, is shown to have a much faster convergence to a final solution than the ART algorithm. In addition, for applications involving high-contrast anomalies (for example, tunnel detection) the CCG is shown to have superior performance in locating the anomalous region for almost all test cases considered.

UR - http://www.scopus.com/inward/record.url?scp=0011412431&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0011412431&partnerID=8YFLogxK

U2 - 10.1109/TGRS.1987.289738

DO - 10.1109/TGRS.1987.289738

M3 - Article

AN - SCOPUS:0011412431

VL - GE-25

SP - 691

EP - 701

JO - IEEE Transactions on Geoscience and Remote Sensing

JF - IEEE Transactions on Geoscience and Remote Sensing

SN - 0196-2892

IS - 6

ER -