A method of image reconstruction in EIT (Electrical Impedance Tomography) to determine inhomogeneities of a sample, the method comprising: - positioning contacts at an edge of an investigated domain of the sample; - inducing flow of current through the contacts ; and - reading potentials on the contacts , by using a system of connections connected between the contacts and a data acquisition system ; the method characterized by : - determining a conformal map ( f ) that represents a transformation of a shape of a disk to a shape of the investigated domain; - wherein the contacts are positioned at the edge of the investigated domain at positions which are calculated from equally spaced contacts located on an edge of the disk using the conformal map ( f ); - providing a plurality of hole positions within the disk; - by means of a processor connected to the data acquisition system , processing the currents and potentials read at the contacts by: - selecting four ( α , β ,γ, δ ) of the contacts and measuring a four-probe resistance between the selected contacts ( R αβγδ = V γδ / j αβ ); - determining what radiuses of holes at the provided hole positions provide a best fit for the measured four-probe resistance; and - imaging the determined radiuses of holes that are the best fit by means of a visualisation system , by transforming the positions and radiuses of these holes to the investigated domain using the conformal map ( f ). - positioning contacts at an edge of an investigated domain of the sample; - inducing flow of current through the contacts ; and - reading potentials on the contacts , by using a system of connections connected between the contacts and a data acquisition system ; - determining a conformal map ( f ) that represents a transformation of a shape of a disk to a shape of the investigated domain; - wherein the contacts are positioned at the edge of the investigated domain at positions which are calculated from equally spaced contacts located on an edge of the disk using the conformal map ( f ); - providing a plurality of hole positions within the disk; - by means of a processor connected to the data acquisition system , processing the currents and potentials read at the contacts by: - selecting four ( α , β ,γ, δ ) of the contacts and measuring a four-probe resistance between the selected contacts ( R αβγδ = V γδ / j αβ ); - determining what radiuses of holes at the provided hole positions provide a best fit for the measured four-probe resistance; and - imaging the determined radiuses of holes that are the best fit by means of a visualisation system , by transforming the positions and radiuses of these holes to the investigated domain using the conformal map ( f ). - selecting four ( α , β ,γ, δ ) of the contacts and measuring a four-probe resistance between the selected contacts ( R αβγδ = V γδ / j αβ ); - determining what radiuses of holes at the provided hole positions provide a best fit for the measured four-probe resistance; and - imaging the determined radiuses of holes that are the best fit by means of a visualisation system , by transforming the positions and radiuses of these holes to the investigated domain using the conformal map ( f ). <EOS>