We describe a method of processing an EEG and/or MEG signal to generate image data representing a 3D current distribution, J, within the brain, the method comprising: capturing a plurality of electric and/or magnetic measurements from the exterior of the head； solving an integral equation for a part of said current distribution to generate said image data representing said 3D current distribution, wherein said integral equation comprises an integral of a first function representing said part of said current distribution and of a second function (∇τvs (r, τ)) representing the geometry and conductivity of the head independent of said current distribution； wherein said solving comprises: modelling the head as at least two regions separated by at least one internal boundary, and solving a set of partial differential equations, one for each said internal region, each partial differential equation comprising a geometry-conductivity function (w(r, τ)) representing the geometry and conductivity of the respective region, wherein said solving is subject to a boundary condition that either i) the gradients of the functions across the or each said internal boundary are smooth when conductivity is taken into account, or ii) a normal component of the electric field of said part of said current distribution is continuous across the or each said internal boundary, and wherein said geometry-conductivity function for an outermost said region of said head defines said second function (∇Tvs (r, τ))
