pith. sign in

arxiv: 1907.01504 · v1 · pith:K5L7GC42new · submitted 2019-06-21 · 🧬 q-bio.NC · physics.comp-ph

FEM-based Scalp-to-Cortex data mapping via the solution of the Cauchy problem

Nikolay Koshev , Nikolay Yavich , Mikhail Malovichko , Ekaterina Skidchenko , Maxim Fedorov This is my paper

Pith reviewed 2026-05-25 18:33 UTC · model grok-4.3

classification 🧬 q-bio.NC physics.comp-ph
keywords EEGCauchy problemLaplace equationfinite element methodscalp-to-cortex mappingspatial resolutionelectrocorticography
0
0 comments X

The pith

Finite element solution of the Cauchy problem maps EEG potentials from scalp to cortex surface.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents a numerical method to solve the ill-posed Cauchy problem for Laplace's equation in order to map measured potentials from the scalp to the brain surface. This pre-processing step for EEG data uses tetrahedral finite elements with linear approximation. A sympathetic reader would care because it promises to boost the spatial resolution of non-invasive EEG to levels approaching invasive ECoG recordings. The approach treats the head as a conductor where potential satisfies Laplace's equation and solves for the potential on the cortex.

Core claim

The proposed FEM-based algorithm for solving the Cauchy problem generates an accurate mapping of the electric potential from the scalp to the brain surface, sufficiently increasing the spatial resolution of EEG to make it comparable with ECoG.

What carries the argument

The numerical solution of the ill-posed Cauchy problem for Laplace's equation using tetrahedral finite element linear approximation, which extrapolates the potential from scalp measurements inward to the cortex.

If this is right

  • The method enables higher resolution EEG without invasive procedures.
  • EEG data can be pre-processed to generate cortex surface potentials directly.
  • Resolution gain makes EEG competitive with ECoG in spatial detail.
  • Applicable to standard EEG setups for improved source localization.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the mapping is accurate, it could reduce the need for invasive recordings in some clinical settings.
  • Combining this with existing source localization techniques might further improve brain activity mapping.
  • Testing on real patient data with known cortical activity would validate the resolution claims.
  • The stability under noise suggests potential for real-time applications if computation is optimized.

Load-bearing premise

The head can be modeled as a homogeneous or piecewise-homogeneous conductor where the potential obeys Laplace's equation, and the numerical solution stays stable and accurate enough with realistic noise to achieve the resolution improvement.

What would settle it

If simulations or experiments show that the mapped cortical potentials do not match direct ECoG measurements within the claimed accuracy under typical noise levels, the resolution gain would not hold.

Figures

Figures reproduced from arXiv: 1907.01504 by Ekaterina Skidchenko, Maxim Fedorov, Mikhail Malovichko, Nikolay Koshev, Nikolay Yavich.

Figure 1
Figure 1. Figure 1: Computational domains and boundaries. In the EEG technique, the measurements are being performed in a bounded area of the surface located commonly at the top of a head. In further consider￾ation, this part of the boundary, which we call accessible part of the boundary, is represented by the subdomain of the outer surface Γ1 ⊂ ∂Ω1 (see [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Spherical model (sources located under the surface of the active domain): a) sim [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Spherical model (sources located on the surface of the inner sphere): a) The current [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: MRI-based model of the head: case with piecewise-constant current distribution on [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: MRI-based model of the head: case with random current, distributed within big [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
read the original abstract

We propose an approach and the numerical algorithm for pre-processing of the electroencephalography (EEG) data, enabling to generate an accurate mapping of the potential from the measurement area - scalp - to the brain surface. The algorithm based on the solution of ill-posed Cauchy problem for the Laplace's equation using tetrahedral finite elements linear approximation. Application of the proposed algorithm sufficiently increases the spatial resolution of the EEG technique, making it comparable with much more complicated electrocorticography (ECoG) method.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The paper proposes a numerical algorithm based on linear tetrahedral finite-element discretization to solve the ill-posed Cauchy problem for Laplace's equation, thereby mapping measured scalp potentials to the cortical surface. It asserts that this pre-processing step raises the effective spatial resolution of non-invasive EEG to a level comparable with invasive electrocorticography.

Significance. A validated, stable implementation would supply a practical, non-invasive route to higher-resolution cortical potential estimates from standard EEG arrays, which could be useful for source-localization pipelines and clinical monitoring. The approach rests on well-known FEM machinery applied to a standard forward model, but its claimed resolution gain hinges on unshown numerical properties.

major comments (2)
  1. [Abstract] Abstract: the central claim that the method 'sufficiently increases the spatial resolution ... making it comparable with ... ECoG' is unsupported by any reconstruction error, noise-sensitivity study, or comparison against ground-truth cortical data; the exponential ill-posedness of the Cauchy problem makes such quantification load-bearing for the resolution assertion.
  2. [Abstract] Abstract / method description: no regularization strategy, stability bound, or convergence analysis is indicated for the linear FEM solution of the Cauchy problem under the 10–20 dB SNR levels typical of EEG; without this, the mapping cannot be guaranteed to preserve high-frequency cortical detail.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript. We address each major comment below, indicating where revisions will be made to the manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that the method 'sufficiently increases the spatial resolution ... making it comparable with ... ECoG' is unsupported by any reconstruction error, noise-sensitivity study, or comparison against ground-truth cortical data; the exponential ill-posedness of the Cauchy problem makes such quantification load-bearing for the resolution assertion.

    Authors: We agree that the abstract's claim regarding spatial resolution becoming comparable to ECoG lacks supporting quantitative evidence such as reconstruction errors, noise-sensitivity studies, or ground-truth comparisons. The manuscript presents the FEM-based mapping algorithm and illustrative numerical results but does not contain the requested analyses. We will revise the abstract to remove this specific claim about comparability to ECoG, limiting the description to the method's purpose of mapping scalp potentials to the cortical surface. revision: yes

  2. Referee: [Abstract] Abstract / method description: no regularization strategy, stability bound, or convergence analysis is indicated for the linear FEM solution of the Cauchy problem under the 10–20 dB SNR levels typical of EEG; without this, the mapping cannot be guaranteed to preserve high-frequency cortical detail.

    Authors: The described algorithm applies a direct linear tetrahedral FEM discretization to the Cauchy problem without an explicit regularization strategy, stability bounds, or convergence analysis. We acknowledge that the ill-posedness of the problem implies potential instability at typical EEG noise levels, and the current text does not address preservation of high-frequency details. In the revised manuscript we will add a dedicated discussion of these aspects, including the role of the discretization in providing implicit regularization and suggestions for explicit regularization approaches. revision: yes

Circularity Check

0 steps flagged

No circularity: standard numerical PDE solver

full rationale

The paper presents a direct numerical algorithm (linear tetrahedral FEM) for solving the Cauchy problem for Laplace's equation to map scalp to cortical potentials. No equations, parameters, or self-citations are shown that reduce the output to the input by construction, nor any fitted-input-called-prediction or ansatz-smuggled steps. The derivation is self-contained as a forward solver of a standard elliptic PDE; the resolution-gain claim is an empirical assertion about the solver's output rather than a definitional equivalence.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the standard bioelectromagnetic modeling assumption that the head volume satisfies Laplace's equation; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Electric potential inside the head volume obeys Laplace's equation in source-free regions.
    Invoked implicitly when the Cauchy problem for Laplace's equation is posed on the head domain.

pith-pipeline@v0.9.0 · 5622 in / 1228 out tokens · 33912 ms · 2026-05-25T18:33:59.270482+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    Pascual-Marqui, C.M

    R.D. Pascual-Marqui, C.M. Michel, and D. Lehman. Low resolution elec- tromagnetic tomography: a new method for localizing electrical activity in the brain. Int.J. Psychophysiol., 18:49–65, 1994

    work page 1994

  2. [2]

    Pascual-Marqui

    R.D. Pascual-Marqui. Standardized low-resolution brain electromagnetic tomography (sloreta): technical details. Methods Find. Exp. Clin. Phar- macol. (Suppl. D) , 24:5–12, 2002

    work page 2002

  3. [3]

    H+m+l+inen and R

    M. H+m+l+inen and R. J. Ilmoniemi. Interpreting magnetic fields of the brain: minimum norm estimates. Medical and biological engineering and computing, 32(1):35–42, 1994

    work page 1994

  4. [4]

    F.H. Lin, T. Witzel, S.P. Ahlfors, S.M. Stufflebeam, J. W. Belliveau, and M.S. H+m+l+inen. Assessing and improving the spatial accuracy in meg source localization by depth-weighted minimum-norm estimates. NeuroIm- age, 31(1):160 – 171, 2006

    work page 2006

  5. [5]

    Baillet, J

    S. Baillet, J. C. Mosher, and R. M. Leahy. Electromagnetic brain mapping. IEEE Signal Processing Magazine , 18(6):14–30, Nov 2001

    work page 2001

  6. [6]

    Camilleri, Si- mon G

    Roberta Grech, Tracey Cassar, Joseph Muscat, Kenneth P. Camilleri, Si- mon G. Fabri, Michalis Zervakis, Petros Xanthopoulos, Vangelis Sakkalis, and Bart Vanrumste. Review on solving the inverse problem in eeg source 14 analysis. Journal of NeuroEngineering and Rehabilitation , 5(1):25, Nov 2008

    work page 2008

  7. [7]

    B. D. Van Veen, W. van Drongelen, M. Yuchtman, and A. Suzuki. Localiza- tion of brain electrical activity via linearity constrained minimum variance spatial filtering. IEEE Trans. Biomed. Eng., 44(9), 1997

    work page 1997

  8. [8]

    Mosher and R.M

    J.C. Mosher and R.M. Leahy. Recursive music: a framework for eeg and meg source localization. IEEE Trans. Biomed. Eng., 45(11), 1998

    work page 1998

  9. [9]

    Arjan Hillebrand and Gareth R. Barnes. Beamformer analysis of meg data. In Magnetoencephalography, volume 68 of International Review of Neuro- biology, pages 149 – 171. Academic Press, 2005

    work page 2005

  10. [10]

    Sekihara and S.S

    K. Sekihara and S.S. Nagarajan. Adaptive Spatial Filters for Electromag- netic Brain Imaging . Springer: Berlin, 2008

    work page 2008

  11. [11]

    Cortical mapping by laplace–cauchy trans- mission using a boundary element method

    Maureen Clerc and Jan Kybic. Cortical mapping by laplace–cauchy trans- mission using a boundary element method. Inverse Problems, 23(6):2589, 2007

    work page 2007

  12. [12]

    De Munck, C.H

    J.C. De Munck, C.H. Wolters, and M. Clerc. EEG and MEG: forward modeling, chapter Chapter 6, pages 203–272. CUP/BRRR, 2012

    work page 2012

  13. [13]

    Evaluation of numerical techniques for solving the cur- rent injection problem in biological tissues

    Fredrik Andersson. Evaluation of numerical techniques for solving the cur- rent injection problem in biological tissues. Proceedings - International Symposium on Biomedical Imaging , 2016-June:876–880, 2016

    work page 2016

  14. [14]

    Convergence rates for the quasi-reversibility method to solve the cauchy problem for laplace’s equation

    Bourgeois L. Convergence rates for the quasi-reversibility method to solve the cauchy problem for laplace’s equation. Inverse Problems, 22(2):413– 430, 2006

    work page 2006

  15. [15]

    A mixed formulation of quasi-reversibility to solve the cauchy problem for laplace’s equation

    Bourgeois L. A mixed formulation of quasi-reversibility to solve the cauchy problem for laplace’s equation. Inverse Problems, 21(3):1087–1104, 2005

    work page 2005

  16. [16]

    Conductivity of living intracranial tissues

    Juha Latikka, Timo Kuurne, and Hannu Eskola. Conductivity of living intracranial tissues. Physics in Medicine and Biology , 46:1611–1616, 2001

    work page 2001

  17. [17]

    J. Sarvas. Basic mathematical and electromagnetic concepts of the bio- magnetic inverse problem. Physics in Medicine and Biology , 32(1):11–22, 1987

    work page 1987

  18. [18]

    Biomagnetic Phenomena

    Plonsey R. Biomagnetic Phenomena. New York: McGraw-Hill, 1969

    work page 1969

  19. [19]

    H¨ am¨ al¨ ainen, R

    M. H¨ am¨ al¨ ainen, R. Hari, R. J. Ilmoniemi, J. Knuutila, and O. V. Lounas- maa. Magnetoencephalography: theory, instrumentation, and applications to noninvasive studies of the working human brain. Rev. Modern Phys. , 1993. 15

    work page 1993

  20. [20]

    Zakharova, P.I

    T.V. Zakharova, P.I. Karpov, and V.M. Bugaevskii. Localization of the activity source in the inverse problem of magnetoencefalography. Compu- tational Mathematics and Modeling , 28(2):148–157, 2017

    work page 2017

  21. [21]

    Oostendorp

    Thom F. Oostendorp. The conductivity of the human skull: Results of in vivo and in vitro measurements. IEEE TRANSACTIONS ON BIOMEDI- CAL ENGINEERING, 47:1487–1492, 2000

    work page 2000

  22. [22]

    Engl and A

    H.W. Engl and A. Leit˜ ao. A mann iterative regularization method for elliptic cauchy problems. Numerical Functional Analysis and Optimization, 22(7-8):861–884, 2001

    work page 2001

  23. [23]

    Numerical solution of a cauchy prob- lem for the laplace equation

    Frederik Berntsson and Lars Eld ’en. Numerical solution of a cauchy prob- lem for the laplace equation. Inverse Problems, 17:839–853, 2001

    work page 2001

  24. [24]

    H.J. Li, X.S. Feng, J. Xiang, and P.B. Zuo. New approach for solving the inverse boundary value problem of laplace’s equation on a circle: Technique renovation of the grad-shafranov (gs) reconstruction.Journal of Geophisical Research: space physics, 118:1–7, 2013

    work page 2013

  25. [25]

    M.Klibanov, N

    V. M.Klibanov, N. A. Koshev, L.I.Jingzhi, and A. G. Yagola. Numerical solution of an ill-posed cauchy problem for a quasilinear parabolic equa- tion using a carleman weight function. Journal of Inverse and Ill-posed Problems, 2016

    work page 2016

  26. [26]

    Anatoly. B. Bakushinsky, Michael V. Klibanov, and Nikolaj A. Koshev. Carleman weight functions for a globally convergent numerical method for ill-posed cauchy problems for some quasilinear pdes. Nonlinear Analysis: Real World Applications, 34:201–224, 2017

    work page 2017

  27. [27]

    Pereverzev

    Hui Cao and Sergei V. Pereverzev. The balancing principle for the regular- ization of elliptic cauchy problems. Inverse Problems, 23:19431961, 2007

    work page 2007

  28. [28]

    Klibanov, and Sergei V

    Hui Cao, Michael V. Klibanov, and Sergei V. Pereverzev. A carleman estimate and the balancing principle in the quasi-reversibility method for solving the cauchy problem for the laplace equation. Inverse Problems , 25:25pp, 2009

    work page 2009

  29. [29]

    Anhdiv-based mixed quasi-reversibility method for solving elliptic cauchy problems.SIAM Journal on Numerical Analysis , 51:2123–2148, 2013

    J´ er´ emi Dard´ e, Antti Hannukainen, and Nuutti Hyv¨ onen. Anhdiv-based mixed quasi-reversibility method for solving elliptic cauchy problems.SIAM Journal on Numerical Analysis , 51:2123–2148, 2013

    work page 2013

  30. [30]

    Lavrentyev, V.G

    M.M. Lavrentyev, V.G. Romanov, and S.P. Shishatsky. Ill-posed Problems of Mathematical Physics and Analysis . M. Nauka, 1980

    work page 1980

  31. [31]

    Bourgeois

    L. Bourgeois. Convergence rates for the quasi-reversibility method to solve the cauchy problem for laplaces equation. Inverse Problems, 22:413–430, 2006. 16

    work page 2006

  32. [32]

    Kuznetsov, O.V

    Yu.A. Kuznetsov, O.V. Boiarkine, I.V. Kapyrin, and N.B. Yavich. Numer- ical analysis of a two-level preconditioner for the diffusion equation with an anisotropic diffusion tensor. Russ. J. Numer. Anal. Math. Modelling , 22(4):377–391, 2007

    work page 2007

  33. [33]

    http://iso2mesh.sourceforge.net/cgi-bin/index.cgi?Home. 17

    work page