2.1. Description of HP Based Imaging Procedure
It has been discussed in a study on breast cancer detection [
14] that microwave technology relies on detecting the contrast at the interface between two tissues with differing dielectric properties (between healthy and malignant tissues). UWB technology has recently been applied to medical applications for detection and monitoring purposes. The use of UWB offers numerous benefits, for instance high resolution (owing to the high bandwidth spectrum) in order to indicate the inclusions in resulting images, high capacity and reliability, and low power transmission. The same principle has been applied to determine the contrast between blood and brain matter to identify stroke in head mimicking phantoms [
15]. In this paper, an algorithm based on HP is applied to forward propagate the waves [
16]. The use of HP removes the necessity for having to solve complex inverse problems. The scattered electric field
E is achieved through Equation (
1) by summing the signals
collected at the points
displaced along a circular surface having radius
, where
is the number of receiving points (from 1 to
):
In the above equation,
≡
is the observation point;
represents the wave number for the media in the imaging zone;
is the spatial sampling (representing the distance between two adjacent receiving points that can be calculated by Equation (2), and
indicates the number of transmitting antennas operating at frequency
f. The “reconstructed” internal field has been indicated by the string “rcstr”, whilst the string HP indicates that Huygens based procedure will be employed in Equation (1). S21 is the parameter representing the complex transfer function from the transmitting antenna to the receiving antenna.
In Equation (1), the Green’s function
G is applied to propagate the field using Equation (3);
The proposed HP-based procedure differs from the Kirchhoff migration method whose algorithms usually perform time reversal and back-propagation to find the phase, that is, time, traces. Conversely, here we are not interested in finding the phase traces: in fact, we use HP with the aim of reconstructing the field. More details can be found in [
17]. It should be pointed out that the reconstructed electric field depends on both illuminating source and frequency. If we use
frequencies
, then the intensity of the final image
I can be obtained mathematically using Equation (4).
To remove artefacts, transmitter image and the reflections of the layers, signal pre-processing methods are employed. An accurate detection will not be feasible to obtain when the artefact is not totally eliminated and the inclusions (e.g., the strokes) are masked. In this paper, we have applied different artefact removal methods with the data from the measurements in order to find the best method of artefact removal, specifically for the detection of brain haemorrhage.
2.2. Artefact Removal Methods
In this approach of artefact removal, two sets of experimental measurements are performed for the imaging scenario; First, the “no target” scenario, for a phantom with no target. Secondly, the “target” scenario a repeat of the procedure using a phantom with a millimetric cylindrically shaped inclusion to emulate a haemorrhagic brain stroke. The inclusion will be indicated in the image obtained through Equation (4), by subtracting the recorded S21 of the “no target” scenario from that of the “target” scenario “”. This expression can be used in Equation (1). An “ideal” case is then generated to prove the concept of the technology. This would mean that the “ideal” image method could be used as a reference for comparison with the resulting images from other artefact removal methods. This will be presented in the next section.
It is important to point out that, when dealing with a real scenario, there is no possibility of applying this artefact removal method to medical imaging. We apply this method to show the strong feasibility of detecting the brain haemorrhage in an ideal way. Moreover, one of the efforts of the current paper is to find an algorithm which generates as close of a match as possible to the ideal result. Applying several different artefact removal algorithms has the potential to vary the resulting image. Therefore, it is not feasible to test the algorithm in situations where the ideal response is not calculated or known (i.e., measured data from a real human head with a brain stroke). Hence, in clinical trials this artefact removal method cannot be helpful and effective.
The rotation subtraction strategy has been performed for artefact removal by employing a setup to replicate a signal from two transmitters, which have been positioned 4.5
apart on the perimeter of the cylinder. In the rotation subtraction method [
18,
19], artefacts are eliminated by subtracting two measurements collected using two slightly displaced transmitting positions. The RS strategy is described mathematically as “
”. This expression can be used in Equation (1). The current procedure has been implemented by subtracting the transmitting position
from the transmitting position
m, where
m and
correspond to the same set of recorded measurements.
In this approach of artefact removal [
19], we endeavour to successfully remove the artefact by using the average values of the complex S21, which have been received from many transmitting antennas. These transmitting antennas are positioned 4.5
apart from each other. We apply the subtraction between received signals
and the mean of received signals
mathematically through “
”. This expression can be employed in Equation (1).
The medium we eventually will be using (the brain) has a natural symmetry. By exploiting the (eventual) object symmetry, it may be possible to apply another artefact removal method using the difference between the receivers placed symmetrically opposite. This method is initially derived from the literature of [
20,
21,
22] used by Mustafa et al. In the DSR Type method, artefacts are removed by performing the subtraction between each receiver value from its symmetrically opposite receiver. The artefact removal can be obtained mathematically through DSR method using the following equation:
where
is the maximum number of receiving antennas,
is an index to indicate the maximum transmitting positions and
is the original recorded complex S21. This creates a Differential (Symmetric Receiver Type) matrix
S. The AS or RS method is then applied to this matrix.
The DSR method mentioned above is based on the natural symmetry of some media (e.g., the brain), across the left and right halves. It is worthwhile pointing out that the images which have used the symmetric method may contain mirrored artefacts. To eliminate the mirrored section, SSD could be successful. The ellipsoidal shape of the human head has a distinct left-right line of symmetry. The front-back sections of the brain also contain similar densities of tissue. Whilst not absolutely symmetrical, the resemblance in shape and density could be utilised to offer an artefact removal method by summing a differential matrix formed from the left-right differential and a second matrix formed from a front-back differential. Hence, this can provide a more intense peak at the area of inclusion, and subsequently the mirrored artefacts will have a reduced intensity. A differential matrix
S is constructed the same way as before by applying Equation (5). A second matrix
R is constructed across the front-back receivers and defined as follows:
Then the combined matrix C is composed of summing matrices S and R and used as a Differential (Summed Symmetric Receiver Type) matrix. The AS or RS method is then applied to this matrix.
2.3. Image Quantification Metrics
Imaging performance has been investigated using image quantification. As a portion of this research, the calculation of certain parameters will be required to quantify the imaging’s detection capabilities, compare the proposed artefact removal methods, and provide a quantifiable measurement system for comparing images. There are several metrics used for quantifying stroke detection capability.
Based on two scenarios discussed in this paper (the “no target” scenario and the “target” scenario), we introduced the calculation of six metrics which allow us to perform a quantification of detection accuracy. These metrics are categorised into those that rely on a reference image called “ideal” image and those independent of that, and are comprised of Area Difference (ArD) index, Polyshape Construction, Centroid Difference (CD), Signal-to-Noise Ratio (SNR), Signal-to-Clutter ratio (S/C) and Structural Similarity Index Metric (SSIM). Amongst these, SSIM, ArD, and CD are dependent on an “ideal” image.
For the purposes of this experiment, an “ideal” image has been considered as shown in the results section. Further details with reference to these metrics are explained below.
When evaluating an image, the inclusion in the “ideal” image is approximately made up of the normalised values which are over 75% of the maximum value. To evaluate the shape of the inclusion, the image is adjusted through expanding the values greater than 0.75 to 1 and enforcing to 0 the values less than 0.75. By applying MATLAB’s polyboundary and polyshape functions, the resulting shape can be achieved [
23].
Having set a threshold and defined two areas as the “target” area and the “background” area, a comparison can then be drawn between the size of the target area for an “ideal” image and that of the test image. Mathematically speaking, it can be obtained through the following equation:
In the above, and are the number of values over the target threshold in the experiment image and “ideal” image, respectively. It is important to highlight that this is a useful measure of the precision of the target area, but not the accuracy. For this, another metric is required.
Evaluation of detection accuracy is carried out by assessing the Euclidean difference between the centroid of an “ideal” image polyshape and that of the test image. Assuming the shape has a constant density; the centroid function will take a polyshape and find the exact centre of mass. Using Cartesian coordinates, if the centroid of the “ideal” image is positioned at
and the centroid of the test image is positioned at
, then the distance,
, is calculated in metres using the Pythagorean formula. This is done using the MATLAB pdist function [
23].
This method of image quantification has basically been derived from the literature of [
19]. The SNR is a useful metric in determining how clear any detected inclusion is by providing an assessment of the ratio between the background noise and the desired signal. It is calculated in dBs by using the above mentioned threshold for computing the Polyshape, with the aim of determining the target and background area. The calculations of the SNR are performed through Equation (9):
In the above equation, and are the mean values of the detected target and background regions, respectively, and is the standard deviation of the background.
By applying the SSIM method, a value is calculated between 0 and 1, which indicates the structural similarity between two images (with 1 meaning the images are identical) [
24]. We proceed by writing the following equation:
where
x represents the reference image,
y represents the test image,
and
are the corresponding mean,
and
represent the corresponding variance,
is the covariance of the reference and test images. In above equation,
, and
are small constants. SSIM can be calculated on the basis of two input images using the SSIM function in MATLAB [
23]. This will output both a value and a monochrome mapping, which is a useful visual assessment of the structural similarity between the images.
To evaluate the performance of the imaging algorithms in this research, a quantitative metric is used, which is S/C ratio. Clutter comes from undesired scatter from objects within the radar beam that are not targets and has been characterised by several distributions. As the resulting images might possibly contain some clutter even after applying artefact removal procedures, it is applicable to present a parameter with the intention of quantifying and comparing the performance in detection when applying different artefact removal algorithms. Typically, the S/C ratio has been described as the ratio of the maximum brain haemorrhage response to the maximum clutter response. Here, S/C is described as the ratio between the maximum intensity evaluated in the region of the lesion divided by the maximum intensity outside the region of the lesion [
25]. This metric evaluates the accuracy and strength of the results and would be highly effective to employ in comparisons.
2.4. Description of Phantom
With the aim of validating the MWI system and finding the best method of artefact removal for detection of brain haemorrhagic strokes, a simple phantom is presented. A simple phantom is constructed mimicking the dielectric properties of the human brain, to which a millimetric cylindrically shaped inclusion is applied to emulate a brain haemorrhagic stroke. The two layers mimic: (I) average brain tissues with realistic dielectric properties of white matter and grey matter; and (II) blood inclusion (the inclusion has been placed inside the first layer). There is a dominating region of emulated white and grey matter with a mean conductivity and permittivity of 1.01 S/m and 44, respectively at 1.5 GHz. The emulated blood provides a contrast through its conductivity and permittivity of 1.79 S/m and 60, respectively at 1.5 GHz. The dielectric constant values of both the first layer and the inclusion have been derived from [
26]. In order to simulate the dielectric properties of the human brain and haemorrhage tissues, different materials were applied with the purpose of approximating the values described. The mentioned phantom contains different combinations of liquids emulating grey/white matter and blood (brain haemorrhage). Each of the combinations comprises a mixture of deionised water and glycerine at different ratios. Furthermore, choosing the materials was driven by many beneficial aspects such as the stability over long periods of storage time, the low cost of materials, and the ease of access (off-the-shelf). The first layer, which imitates the brain layer, has been fabricated by creating a mixture of deionised water and glycerine with ratios of 60% and 40%, respectively as given in [
26]. The inclusion has been mimicked using a combination of 15% glycerine and 85% water [
26]. In the mixture proposed, deionised water is employed as the principal substance or source of permittivity, as it shows high dielectric values over a wide bandwidth. Also, a pinch of salt is applied to control the relative permittivity of the compounds. The Epsilon dielectric measurement device (Biox Epsilon Model E100, fabricated at Biox System Company Ltd.) has been used for measuring the electrical properties of the combinations of blood and brain emulating tissue, at room temperature. It should be highlighted that such dielectric properties can be considered representative of brain and an inclusion constituted of blood [
26]. A round-bottom cylindrical ABS-Plastic mould with diameter of 110 mm and height of 115 mm has been designed and fabricated to maintain the brain equivalent material. A cylindrically shaped tube with a diameter of 4 mm has been employed to contain an inclusion.
2.5. Description of Microwave Imaging (MWI) Device
This research has provided a new application direction for a portable MWI device which has been previously used for breast imaging through phantom measurements [
27], and verified in preliminary clinical trials [
28]. We have adapted the methods used in this research to detect brain haemorrhage and have produced promising research using the technology for stroke detection.
The MWI device consists of an aluminium cylindrical hub (radius equal to 50 cm) containing two antennas, one transmitting (tx) and one receiving (rx). The hub is internally covered by microwave absorbers, and is equipped with a hole with a cup, allowing the insertion of the object to be imaged. The antennas are installed at the same height, in free space and can rotate around the azimuth in order to collect microwave signals in a multi-bistatic fashion from different angular positions [
13].
The tx and rx are connected to a 2-port VNA (S5065, Copper Mountain, Indianapolis, IN, USA). Concerning the transmitting positions, all the experiments have been performed by employing 15 transmitting positions, displaced in 5 triplets centred at 0
, 72
, 144
, 216
, and 288
. In each triplet, the transmitting positions are displaced by 4.5
. For each transmitting position, we recorded the complex S21 at
= 80 receiving positions, uniformly displaced along a circular surface having radius
= 7 cm (the receiving positions are denoted with the index
np = 1, 2, ..., 80, and the transmitting positions with the index
m = 1, 2, ..., 15). Specifically, we measure complex S21 at the points
≡
≡
, displaced along a circular surface having radius
. For each transmitting and receiving position, the complex S21 is collected from 1 to 1.5 GHz, as this band has been demonstrated as being ideal and optimal for brain imaging [
1], with 5 MHz sampling.
Figure 1a shows a sketch of the configuration of the MWI device, where the light green dashed circle indicates the perimeter where the receiving antenna can be moved circularly, in order to receive the signals from different positions. The transmitting antennas positions are displaced along a circular surface having radius of 35 cm. The brain phantom and inclusion (imitating blood) are presented in light blue and red respectively.
Figure 1b shows the MWI device.