Reduced Graphene Oxide UWB Array Sensor: High Performance for Brain Tumor Imaging and Detection

A low cost, with high performance, reduced graphene oxide (RGO) Ultra-wide Band (UWB) array sensor is presented to be applied with a technique of confocal radar-based microwave imaging to recognize a tumor in a human brain. RGO is used to form its patches on a Taconic substrate. The sensor functioned in a range of 1.2 to 10.8 GHz under UWB frequency. The sensor demonstrates high gain of 5.2 to 14.5 dB, with the small size of 90 mm × 45 mm2, which can be easily integrated into microwave imaging systems and allow the best functionality. Moreover, the novel UWB RGO array sensor is established as a detector with a phantom of the human head. The layers’ structure represents liquid-imitating tissues that consist of skin, fat, skull, and brain. The sensor will scan nine different points to cover the whole one-sided head phantom to obtain equally distributed reflected signals under two different situations, namely the existence and absence of the tumor. In order to accurately detect the tumor by producing sharper and clearer microwave image, the Matrix Laboratory software is used to improve the microwave imaging algorithm (delay and sum) including summing the imaging algorithm and recording the scattering parameters. The existence of a tumor will produce images with an error that is lower than 2 cm.


Introduction
By 2030, it is estimated that almost 13 million individuals globally will be affected by cancer, which can be considered as one of the most complex illnesses that ever existed [1]. Early cancer detection might increase the survival ratio among the cancer patients due to the fact that less sever and more varied treatment options could be offered. Presently, Computed Tomography (CT) Scan, Magnetic Resonance Imaging (MRI), X-ray, microwave imaging, and the ultrasound machine have been used as current imaging equipment applied specifically to identify various types of cancers [2]. Among them, the microwave imaging technique for cancer detection offers numerous significant advantages as well as being firm, inexpensive, safe, and non-invasive [3,4]. The formation of tissue with an irregular cell development is well-defined as a brain tumor, which is considered to be one of the most complicated cancers ever to exist [5]. To visualize the internal human structure, microwave imaging is used in order to reveal electromagnetic fields that work in microwave frequencies, with a range of 300 MHz and 30 GHz [6]. There are three techniques used for microwave imaging, namely hybrid, active, and passive.
The novel development of a low-cost with high-performance reduced graphene oxide (RGO) UWB array sensor for brain tumor detection is described in this paper. The utilization of the combined RGO is pointed toward improving sensor gain and bandwidth. It enables the infiltration into the complex human head through interaction at various frequencies, as such cancers are normally situated inside the brain [21]. Moreover, the Improved Delay and Sum (IDAS) imaging algorithm is introduced by applying mslowess filter as for the smoothing technique by reducing the clutter to create more intensity of microwave imaging.   Figure 3a,b shows Field Emission Scanning Electron Microscope device and RGO SEM image respectively. The figure clearly illustrates the formed of scattered groups RGO sheets with a clear sign stacking. The RGO sheet looks like a layer sheet morphology in excess of outsized area due to large surface to volume ratio. Similar morphological nature of the RGO with sheets closely stacked due to elimination of oxygen group to form a closely associated stack arrangement has been reported in literature [22].    Figure 3a,b shows Field Emission Scanning Electron Microscope device and RGO SEM image respectively. The figure clearly illustrates the formed of scattered groups RGO sheets with a clear sign stacking. The RGO sheet looks like a layer sheet morphology in excess of outsized area due to large surface to volume ratio. Similar morphological nature of the RGO with sheets closely stacked due to elimination of oxygen group to form a closely associated stack arrangement has been reported in literature [22].  Figure 3a,b shows Field Emission Scanning Electron Microscope device and RGO SEM image respectively. The figure clearly illustrates the formed of scattered groups RGO sheets with a clear sign stacking. The RGO sheet looks like a layer sheet morphology in excess of outsized area due to large surface to volume ratio. Similar morphological nature of the RGO with sheets closely stacked due to elimination of oxygen group to form a closely associated stack arrangement has been reported in literature [22].     Figure 4a, the most intensive characteristic peak observed at 26.60 • corresponds to the (002) plane. It is clearly indicating the formation of pure graphite [23]. The reaction of oxidation in Figure 4b presents the broader diffraction peak at 12.70 • which tends to wider graphite unique peak and moves closely to 12.70 • . Moreover, it implies that the obtained diffraction shifts in GO due to the interface space higher than graphite [24]. In the case of RGO (Figure 4b), a new characteristic broad diffraction peak (broad hump) appears to be centered at 23.5 • . An increase of broad humps due to reduction of oxygen functionalities point out to decrease of intercalated layer sheets. The reductions of oxygen function groups induce an increase in interplanar distance. Similar work has been reported using reduction technique of functional group for oxygen removal by Wang et al. [25] where the stacking sheets are closely linked to each other due to the interplanar distance and d-spacing reduction. The formation of graphene single layer structures due to the removing of the multilayer component will increase the greater spreading peak. The decreasing of the size of sheets with a shapeless arrangement is formed from a broad hump at 2θ = 20-30 • [26]. Thus, the two-stepped method of this study in producing RGO is proven to have sufficient reduction properties as compared to previous techniques reported in the literature [27]. The removal or reduction of oxygen functionalities due to this chemical with the microwave assisted combination process reduction method leads to an increase in conductivity of RGO.

Structure Study from XRD
X-Ray Diffraction (XRD) device and XRD patterns of pure graphite, GO an RGO are presented in Figure 4a, Figure 4b, and Figure 4c respectively. From Figure 4a, the most intensive characteristic peak observed at 26.60° corresponds to the (002) plane. It is clearly indicating the formation of pure graphite [23]. The reaction of oxidation in Figure 4b presents the broader diffraction peak at 12.70° which tends to wider graphite unique peak and moves closely to 12.70°. Moreover, it implies that the obtained diffraction shifts in GO due to the interface space higher than graphite [24]. In the case of RGO (Figure 4b), a new characteristic broad diffraction peak (broad hump) appears to be centered at 23.5°. An increase of broad humps due to reduction of oxygen functionalities point out to de-crease of intercalated layer sheets. The reductions of oxygen function groups induce an increase in interplanar distance. Similar work has been reported using reduction technique of functional group for oxygen removal by Wang et al. [25] where the stacking sheets are closely linked to each other due to the interplanar distance and d-spacing reduction. The formation of graphene single layer structures due to the removing of the multilayer component will increase the greater spreading peak. The decreasing of the size of sheets with a shapeless arrangement is formed from a broad hump at 2θ = 20-30° [26]. Thus, the twostepped method of this study in producing RGO is proven to have sufficient reduction properties as compared to previous techniques reported in the literature [27]. The removal or reduction of oxygen functionalities due to this chemical with the microwave assisted combination process reduction method leads to an increase in conductivity of RGO.    Figure 5a,b shows RGO in the liquid form and Four Point Probe device respectively. RGO recorded 3.38 × 10 8 s/m (conductivity), 2.01 × 10 −3 m (resistivity), 0.2 V (power), and 9.48 × 10 −1 A (current) according to the experiments. The measured high conductivity and low resistivity values are similar as compared to the early findings on conductivity and resistivity of RGO by several researchers [25][26][27]. The observed electrical properties are sufficient to be utilized as conductive material to replace copper and other metals for various applications. The chemical with microwave assisted combination process-based reduction method is playing a prime role to enhance the electrical properties.  Figure 5a,b shows RGO in the liquid form and Four Point Probe device respectively. RGO recorded 3.38 × 10 8 s/m (conductivity), 2.01 × 10 −3 m (resistivity), 0.2 V (power), and 9.48 × 10 −1 A (current) according to the experiments. The measured high conductivity and low resistivity values are similar as compared to the early findings on conductivity and resistivity of RGO by several researchers [25][26][27]. The observed electrical properties are sufficient to be utilized as conductive material to replace copper and other metals for various applications. The chemical with microwave assisted combination process-based reduction method is playing a prime role to enhance the electrical properties.

Sensor Design and Development
It is computed how much RGO is needed to make circular patches with a radius of 15 mm and thicknesses of 0.035. As shown in Figure 6, the patch positions are subsequently taken on top of the substrate by a bespoke mould. Once all four moulds are filled with liquid RGO, are heated for 5 h in hot oven with 50 °C. Excessive heating with longer time and higher temperature could compromise the structure of the RGO as the process needs to be performed with care as it may lead to cracked RGO layers.   Figure 7a. The patches are equal, having diameter of 15 mm and the square parasitic element dimensions of 8.0 × 32.0 mm 2 , which is located adjacent to the transmission lines. It has Ultra-Wide Band characteristic through focusing the wave within desired direction by terminating the wave propagated in unwanted directions [28]. This contributes towards the novel design of the UWB RGO

Sensor Design and Development
It is computed how much RGO is needed to make circular patches with a radius of 15 mm and thicknesses of 0.035. As shown in Figure 6, the patch positions are subsequently taken on top of the substrate by a bespoke mould. Once all four moulds are filled with liquid RGO, are heated for 5 h in hot oven with 50 • C. Excessive heating with longer time and higher temperature could compromise the structure of the RGO as the process needs to be performed with care as it may lead to cracked RGO layers.  Figure 5a,b shows RGO in the liquid form and Four Point Probe device respectively. RGO recorded 3.38 × 10 8 s/m (conductivity), 2.01 × 10 −3 m (resistivity), 0.2 V (power), and 9.48 × 10 −1 A (current) according to the experiments. The measured high conductivity and low resistivity values are similar as compared to the early findings on conductivity and resistivity of RGO by several researchers [25][26][27]. The observed electrical properties are sufficient to be utilized as conductive material to replace copper and other metals for various applications. The chemical with microwave assisted combination process-based reduction method is playing a prime role to enhance the electrical properties.

Sensor Design and Development
It is computed how much RGO is needed to make circular patches with a radius of 15 mm and thicknesses of 0.035. As shown in Figure 6, the patch positions are subsequently taken on top of the substrate by a bespoke mould. Once all four moulds are filled with liquid RGO, are heated for 5 h in hot oven with 50 °C. Excessive heating with longer time and higher temperature could compromise the structure of the RGO as the process needs to be performed with care as it may lead to cracked RGO layers.   Figure 7a. The patches are equal, having diameter of 15 mm and the square parasitic element dimensions of 8.0 × 32.0 mm 2 , which is located adjacent to the transmission lines. It has Ultra-Wide Band characteristic through focusing the wave within desired direction by terminating the wave propagated in unwanted directions [28]. This contributes towards the novel design of the UWB RGO  On the other hand, RGO construct the patches while parasitic element and transmission lines are constructed by common copper as shown in Figure 7a. The patches are equal, having diameter of 15 mm and the square parasitic element dimensions of 8.0 × 32.0 mm 2 , which is located adjacent to the transmission lines. It has Ultra-Wide Band characteristic through focusing the wave within desired direction by terminating the wave propagated in unwanted directions [28]. This contributes towards the novel design of the UWB RGO sensor. Meanwhile, Figure 7b shows the partial ground made of copper, which measures 18 × 90 mm 2 , as well as a Sub Miniature connection placement at the back of the substrate plane. This structure's total dimensions are 90 mm × 45 mm 2 . In order to improve gain and directivity, a reflector made of copper has been placed 20 mm away behind the sensor to eliminate back lobes [29], which is proven by the high gains obtained and shown in Figure 7c. The reflector dimensions are similar to the sensor dimensions. sensor. Meanwhile, Figure 7b shows the partial ground made of copper, which measures 18 × 90 mm 2 , as well as a Sub Miniature connection placement at the back of the substrate plane. This structure's total dimensions are 90 mm × 45 mm 2 . In order to improve gain and directivity, a reflector made of copper has been placed 20 mm away behind the sensor to eliminate back lobes [29], which is proven by the high gains obtained and shown in Figure  7c. The reflector dimensions are similar to the sensor dimensions. The physical dimensions of the sensor patch (r), a radius of 15 mm and thickness of 0.035, could be determined as follows [30]: where r = patch radius, h = substrate thickness, fr = resonant frequency, Ɛr = substrate dielectric constant, Ɛeff = effective dielectric constant. Through side feeding, the proposed sensor was supplied by a 50 Ω SMA connection. The feed line is 2.9 mm in diameter to exactly match the sensor's source feed of 50 Ω. The following is the formula for determining microstrip feed lines dimension to obtain impedance of 50 Ω [30]: The parasitic element is electromagnetically attached to the radiating element to form passive resonator to increase bandwidth and efficiency by suppressing the surface current at the transmission line of the sensor which has a length of 32 mm and width of 8 mm for two closely situated radiating elements [31]. To establish the sensor's radiation and ability to cut off the waves in an unwanted direction, the waves from the parasitic elements should not interfere. The quarter-wave converter technique for current matching purpose ensures the same amount of current is received by the entire RGO patches. Then, 70.71 Ω quarter-wave converters were utilized to perfectly match the impedance lines of 100 and 50 Ω. The quarter-wave converter's formula is presented in the following equation [32]: The physical dimensions of the sensor patch (r), a radius of 15 mm and thickness of 0.035, could be determined as follows [30]: where r = patch radius, h = substrate thickness, fr = resonant frequency, ε r = substrate dielectric constant, ε eff = effective dielectric constant. Through side feeding, the proposed sensor was supplied by a 50 Ω SMA connection. The feed line is 2.9 mm in diameter to exactly match the sensor's source feed of 50 Ω. The following is the formula for determining microstrip feed lines dimension to obtain impedance of 50 Ω [30]: where Wm = microstrip feed line width, h = substrate thickness, ε r = substrate dielectric constant, Z 0 = input impedance. The parasitic element is electromagnetically attached to the radiating element to form passive resonator to increase bandwidth and efficiency by suppressing the surface current at the transmission line of the sensor which has a length of 32 mm and width of 8 mm for two closely situated radiating elements [31]. To establish the sensor's radiation and ability to cut off the waves in an unwanted direction, the waves from the parasitic elements should not interfere. The quarter-wave converter technique for current matching purpose ensures the same amount of current is received by the entire RGO patches. Then, 70.71 Ω quarter-wave converters were utilized to perfectly match the impedance lines of 100 and 50 Ω. The quarter-wave converter's formula is presented in the following equation [32]: where Wm = microstrip feed line width, h = substrate thickness, ε r = substrate dielectric constant, Z 0 = input impedance, Z 1 = transformer characteristic impedance, R in = edge resistance at resonance. Equations (5) and (6) are used to calculate the width of the impedance feedlines (50 and 100) as well as quarterwave transformer (70.71) [32]: where h = substrate height while A and B could be evaluated as follows:

Sensor Principle of Operation
The UWB RGO sensor exhibits operated frequency of less than −10 dB ranging from 1.2 GHz-10.8 GHz for measured while 2.1-10.2 GHz for simulated with enormous energy generated ranging from 5.2-14.5 dB for measured gain and 5.0-14.0 dB for simulated gain, as demonstrated by Figure 8. The different results of measured and simulated reflection coefficient (s11) are due to deficiency of RF cable which connected between sensor and vector network analyzer (VNA) during the measurement. Yet, the operated frequencies are still within the desired frequencies for both simulated and measured. The simulations were performed in Microwave Studio by CST, a full wave electromagnetic solver, specifically using time solver. For the radiation safety part, the sensor radiated microwave energy is set to be less than 1.6 W/Kg for Specific Absorption Rate standardized by Federal Communications Commission (FCC). The value recorded by the sensor is 1.2 (W/Kg).

Sensor Principle of Operation
The UWB RGO sensor exhibits operated frequency of less than −10 dB ranging from 1.2 GHz-10.8 GHz for measured while 2.1-10.2 GHz for simulated with enormous energy generated ranging from 5.2-14.5 dB for measured gain and 5.0-14.0 dB for simulated gain, as demonstrated by Figure 8. The different results of measured and simulated reflection coefficient (s11) are due to deficiency of RF cable which connected between sensor and vector network analyzer (VNA) during the measurement. Yet, the operated frequencies are still within the desired frequencies for both simulated and measured. The simulations were performed in Microwave Studio by CST, a full wave electromagnetic solver, specifically using time solver. For the radiation safety part, the sensor radiated microwave energy is set to be less than 1.6 W/Kg for Specific Absorption Rate standardized by Federal Communications Commission (FCC). The value recorded by the sensor is 1.2 (W/Kg).   Figure 9 exhibits the sensor's polar energy form in the azimuthal direction for microwave brain imaging frequency resonance at 2 GHz and 3 GHz. Basically, both simulated and measured radiation patterns do radiate in a wide frequency band and could be considered as a unidirectional sensor due to the radiated wave continuously radiating toward the main lobe as compared to the side and rear lobes [33,34]. Moreover, the results of simulation and measurement indicate a satisfactory agreement. All the characteristics such as sensor wave shape, wide operating frequency, and high gain have made the sensor an ideal sensor for brain tumor detection. Figure 9 exhibits the sensor's polar energy form in the azimuthal direction for microwave brain imaging frequency resonance at 2 GHz and 3 GHz. Basically, both simulated and measured radiation patterns do radiate in a wide frequency band and could be considered as a unidirectional sensor due to the radiated wave continuously radiating toward the main lobe as compared to the side and rear lobes [33,34]. Moreover, the results of simulation and measurement indicate a satisfactory agreement. All the characteristics such as sensor wave shape, wide operating frequency, and high gain have made the sensor an ideal sensor for brain tumor detection.

Sensor Integration and Experimental Test
Interaction of electromagnetic waves in the frequency band of 0.3 to 30 GHz with biological tissues could lead to signal absorption by the tissues or reflection back to the source [35]. There are many types of interactions that can possibly take place for different applications. In microwave imaging, reflection or scattering signals have been considered for cancer detection purposes [36]. As microwaves propagate through free space into a particular tissue, there is an interface layer between free space and the tissue medium. Different impedance of both mediums resulted in different microwave interaction for both tissues. These two mediums have different impedance values, and the microwaves will interact with both tissues differently. Some of the microwave signals will transmit through into the tissue and some will be reflected to the source or scattered away from it. The transmitted signal (Et) and reflected signal (Er) can be expressed as [37]: where Ei = incident wave, τ = transmission coefficient and Γ = reflection coefficient. Two different mediums of reflection coefficient could be determined from their impedance values. Transmission line theory stated that wave reflection could occur when incoherence incident wave undergoes impedance characteristics changes. Thus, reflection coefficient is given as [38]: μo = permeability of free space, εo = permittivity of free space. For a medium 'm', permittivity εm is given as: where εm = relative permittivity of the material.

Sensor Integration and Experimental Test
Interaction of electromagnetic waves in the frequency band of 0.3 to 30 GHz with bio-logical tissues could lead to signal absorption by the tissues or reflection back to the source [35]. There are many types of interactions that can possibly take place for different applications. In microwave imaging, reflection or scattering signals have been considered for cancer detection purposes [36]. As microwaves propagate through free space into a particular tissue, there is an interface layer between free space and the tissue medium. Different impedance of both mediums resulted in different microwave interaction for both tissues. These two mediums have different impedance values, and the microwaves will interact with both tissues differently. Some of the microwave signals will transmit through into the tissue and some will be reflected to the source or scattered away from it. The transmitted signal (Et) and reflected signal (Er) can be expressed as [37]: where Ei = incident wave, τ = transmission coefficient and Γ = reflection coefficient. Two different mediums of reflection coefficient could be determined from their impedance values. Transmission line theory stated that wave reflection could occur when incoherence incident wave undergoes impedance characteristics changes. Thus, reflection coefficient is given as [38]: µo = permeability of free space, εo = permittivity of free space. For a medium 'm', permittivity εm is given as: where εm = relative permittivity of the material. Material impedance could be expressed as: Similarly, for medium 'n', Replacing Equations (13) and (14) for Equation (15), the reflection coefficient is given as: On the other hand, transmission coefficient is given as: These Equations and analyses could be applied to all types of microwave imaging techniques including UWB radar imaging. As for the initial step in microwave imaging detection, UWB sensor reflections can determine the existence of abnormalities at tissue boundaries [37,38]. Transmitted and reflected signal ratio to the incident signal are known as scattering parameters. These parameters may be measured or calculated. They are measured in microwave imaging when a source sensor or an array sensor sequentially transmits microwave signals into a tissue of interest receive the reflected or transmitted signals [33]. Scattering parameters are better understood from the perspective of multiport networks or the array system of sensors. The simulated and experimental arrangement of the microwave imaging method to identify the brain tumor is illustrated Figure 10. The tumor is detected by the arrangements comprised of human head phantom, RGO sensor, and VNA.
Material impedance could be expressed as: Similarly, for medium 'n', Replacing Equations (13) and (14) for Equation (15), the reflection coefficient is given as: On the other hand, transmission coefficient is given as: These Equations and analyses could be applied to all types of microwave imaging techniques including UWB radar imaging. As for the initial step in microwave imaging detection, UWB sensor reflections can determine the existence of abnormalities at tissue boundaries [37,38]. Transmitted and reflected signal ratio to the incident signal are known as scattering parameters. These parameters may be measured or calculated. They are measured in microwave imaging when a source sensor or an array sensor sequentially transmits microwave signals into a tissue of interest receive the reflected or transmitted signals [33]. Scattering parameters are better understood from the perspective of multiport networks or the array system of sensors. The simulated and experimental arrangement of the microwave imaging method to identify the brain tumor is illustrated Figure 10. The tumor is detected by the arrangements comprised of human head phantom, RGO sensor, and VNA. To achieve an accurate result, the measurements are performed in an anechoic chamber in order to remove the unwanted signal that might interfere with the result. The RGO sensor is placed parallelly 10 mm away from the head phantom for every reading in order to have optimum signal penetration [39]. The human head phantom used in this work consists of a four-layered rectangular structure sized at 194 × 164 × 260 mm 3 . Its shell is made of Plexiglas (εr = 2.3) and is partitioned into sections of brain, skull, fat, and skin. The brain tissue is placed in the innermost (center) of this phantom with an area of 170 × 140 mm 2 , followed by the skull (10 mm thick), fat (5 mm thick), and skin (5 mm thick). To achieve an accurate result, the measurements are performed in an anechoic chamber in order to remove the unwanted signal that might interfere with the result. The RGO sensor is placed parallelly 10 mm away from the head phantom for every reading in order to have optimum signal penetration [39]. The human head phantom used in this work consists of a four-layered rectangular structure sized at 194 × 164 × 260 mm 3 . Its shell is made of Plexiglas (εr = 2.3) and is partitioned into sections of brain, skull, fat, and skin. The brain tissue is placed in the innermost (center) of this phantom with an area of 170 × 140 mm 2 , followed by the skull (10 mm thick), fat (5 mm thick), and skin (5 mm thick). Each tissue is separated by a 1-mm-thick Plexiglas layer. A rectangular shape of the head phantom has been used due to the bending limitation of the Plexiglass structure. There are two different conditions in emitting the signal towards head phantom, namely with and without tumor, to acquire signal different values to be applied toward the confocal microwave radar-based algorithm to generate an image. Figure 11 shows the sensor-emitted microwave signal for nine specific regions in order to have full scanning of the entire one-sided region of the phantom where the tumor is located at region 5. The brain tumor is made of mixture of deionized water, polyethylene powder, aluminium powder and hydroxy-ethylcellulose to become a semi-solid structure to resemble the actual tumor, which has the permittivity value of 63.2 (similar to an actual brain tumor). The received backscattered signal (S11) from nine different regions of presence and absence of tumor are calculated, where the difference values acquired are then being converted using image processing for generating the image.
phantom has been used due to the bending limitation of the Plexiglass structure. There are two different conditions in emitting the signal towards head phantom, namely with and without tumor, to acquire signal different values to be applied toward the confocal microwave radar-based algorithm to generate an image. Figure 11 shows the sensor-emitted microwave signal for nine specific regions in order to have full scanning of the entire one-sided region of the phantom where the tumor is located at region 5. The brain tumor is made of mixture of deionized water, polyethylene powder, aluminium powder and hydroxy-ethylcellulose to become a semi-solid structure to resemble the actual tumor, which has the permittivity value of 63.2 (similar to an actual brain tumor). The received backscattered signal (S11) from nine different regions of presence and absence of tumor are calculated, where the difference values acquired are then being converted using image processing for generating the image. To acquire a clear image, averaging and interpolation steps have been added to the available original delay-and-sum algorithm for the enhancement. The purpose of averaging is to reduce the effects of messes originating from dominating reflections of environment, sensor, and equipment. At each frequency sample f and each sensor location (v_k) in (f, k ∈ N), the complex reflection coefficient matrix Γ (j, k) is calculated. A theoretical interpolation approach, on the other hand, is used to produce new data points in a series of known data points. In the meantime, interpolation helps to fill up the unfilled elements among the adjacent pixels by predicting based on the value of certain adjacent pixels [40].
Essentially, the energy at the brain's focus point may be computed by [41]: where = total amount of received signal energy, = simulated received energy signal, = time-delay of each focal point, and = integration range. Assuming that the point of convergence is at ( , , ) and distance ( ) to sensor ( ), the time-delay can be communicated as [41]: where v is the brain tissue spread velocity that could be additionally determined with [41]: Since the distance between the port and the aperture centre at the sensor surface is limited, a time-delay must account for both values determined by the equation above and the time-delay taken by the port to aperture sensor surface. On the other hand, because of the tumor response constantly coming at low amplitude and always being associated with To acquire a clear image, averaging and interpolation steps have been added to the available original delay-and-sum algorithm for the enhancement. The purpose of averaging is to reduce the effects of messes originating from dominating reflections of environment, sensor, and equipment. At each frequency sample f and each sensor location (v_k) in (f, k ∈ N), the complex reflection coefficient matrix Γ (j, k) is calculated. A theoretical interpolation approach, on the other hand, is used to produce new data points in a series of known data points. In the meantime, interpolation helps to fill up the unfilled elements among the adjacent pixels by predicting based on the value of certain adjacent pixels [40].
Essentially, the energy at the brain's focus point may be computed by [41]: where M = total amount of received signal energy, y = simulated received energy signal, T = time-delay of each focal point, and τ = integration range. Assuming that the point of convergence is at (x, y, z) and distance (d) to sensor (i), the time-delay can be communicated as [41]: where v is the brain tissue spread velocity that could be additionally determined with [41]: Since the distance between the port and the aperture centre at the sensor surface is limited, a time-delay must account for both values determined by the equation above and the time-delay taken by the port to aperture sensor surface. On the other hand, because of the tumor response constantly coming at low amplitude and always being associated with responses from other tissue, there is a need to extract the tumor response from other tissues' response in order to have tumor response only for signal processing purposes. Priori method is used for the signal extracting purpose where signals are collected under two different conditions; absence of tumor within the human head phantom (reference signal) and when the tumor is placed within the head phantom. Both conditions of signals are subtracted from each other to have only the backscattered signal originating from the tumor located inside the human head phantom.
An averaging approach is used to eliminate artefacts originating from environmental reflections and sensor reverberation. Averaging is performed by taking an average of all the tumor-only backscattered responses and then subtracting this average from all the tumor responses, which could be stated mathematically by [42]: A sensor with number of positions (M) was used in collecting backscattered signals and b'[n] is the backscattered signal at each sensor position. On the other hand, the Hamming Window was applied in removing lobes at both side of the tumour response, and the Inverse Fast Fourier Transform was applied in converting the domain from frequency to time to obtain a temporal equivalent of the tumor response's frequency-dependent scattering profiles. The introduction of a filter for smoothing technique is the modification proposed in IDAS imaging algorithm. The smoother lowess filter was applied for the purpose of improving visual information in a scatter plot with minimal additional cost on the plotting function [41]. Lowess, also known as "locally weighted scatter plot smooth", applied a linear regression method using local weights to smooth the data. The smoothing process uses neighboring data points and a regression weighting function to determine smoothed value. This smoother was implemented in the IDAS using the MATLAB built-in mslowess function.
Synthetic focusing computes the intensity value of the focal points which are directly related to the pixel values in microwave images. This is achieved by first computing the time delays of the smoothed tumor response signals for time shifting purposes. The time delays are calculated for sensor distance toward each pixel points in image projection area and evaluated by [42]: where h is the distance between sensor and tumor, x i and y j represents pixel points coordinate and X and Y define the sensor positions. Time delay at each pixel point is then: D(x,y) is the sensor distance to each pixel point, c is the speed of light, εr is the average relative permittivity of the medium, given as 40 for the propagation of UWB pulses in a head model [43].
The intensity value at each pixel point is derived through calculating smoothed da-ta values for each time, t(x,y). This is calculated for each pixel point as: where S[n] is the tumor response on which smoothing has been performed, r is the total number of sensor positions in a row and c is the total number of sensor positions in a column.

Result and Discussion
Averaging and interpolation steps were applied in the original delay-and-sum algorithm to achieve better image clarity. In this manuscript, we are comparing the performance of the DAS algorithm technique between copper-based and Reduced Graphene Oxide (RGO) sensors for brain tumor detection where Refs. [13,40] used only the conventional copper sensor. Figure 12 shows the effects of interpolation and averaging algorithm during the process of producing the image. Averaging eliminates the messes effect originating from the direct reflections originated from the equipment and surrounding. Data points are made from the construction of interpolation of different data points from a discrete set. It assists the filling of the empty elements in between the adjacent pixels [34]. Figure 13 shows the images processed applied by the delay and sum algorithm in presence and absence of tumor.
conventional copper sensor. Figure 12 shows the effects of interpolation and averaging algorithm during the process of producing the image. Averaging eliminates the messes effect originating from the direct reflections originated from the equipment and surrounding. Data points are made from the construction of interpolation of different data points from a discrete set. It assists the filling of the empty elements in between the adjacent pixels [34]. Figure 13 shows the images processed applied by the delay and sum algorithm in presence and absence of tumor.  Figure 13 clearly demonstrates that a microwave imaging system equipped with RGO sensor manages to produce a sharper and clearer image for both condition with tumor and without tumor as compared with another microwave imaging system equipped with a common sensor. This is due to the higher conductivity of graphene that leads to a higher signal produced by the RGO sensor to be radiated towards the tumor. The images demonstrated the efficiency of RGO sensor by perfect and clear images with better pixel intensity. The color with the most intensity (reddish-orange area) represents the presence of a tumor in the image because of the strongest scattering area that matches up to tumor presence in the human head phantom. Table 1 demonstrates the localization error of the detected tumor in the image as compared with the tumor's actual position with the value of 0.2 cm or 2% of error which could be considered a tolerable error to prove the system is efficiently and accurately detecting the desired tumor. conventional copper sensor. Figure 12 shows the effects of interpolation and averaging algorithm during the process of producing the image. Averaging eliminates the messes effect originating from the direct reflections originated from the equipment and surrounding. Data points are made from the construction of interpolation of different data points from a discrete set. It assists the filling of the empty elements in between the adjacent pixels [34]. Figure 13 shows the images processed applied by the delay and sum algorithm in presence and absence of tumor.  Figure 13 clearly demonstrates that a microwave imaging system equipped with RGO sensor manages to produce a sharper and clearer image for both condition with tumor and without tumor as compared with another microwave imaging system equipped with a common sensor. This is due to the higher conductivity of graphene that leads to a higher signal produced by the RGO sensor to be radiated towards the tumor. The images demonstrated the efficiency of RGO sensor by perfect and clear images with better pixel intensity. The color with the most intensity (reddish-orange area) represents the presence of a tumor in the image because of the strongest scattering area that matches up to tumor presence in the human head phantom. Table 1 demonstrates the localization error of the detected tumor in the image as compared with the tumor's actual position with the value of 0.2 cm or 2% of error which could be considered a tolerable error to prove the system is efficiently and accurately detecting the desired tumor.  Figure 13 clearly demonstrates that a microwave imaging system equipped with RGO sensor manages to produce a sharper and clearer image for both condition with tumor and without tumor as compared with another microwave imaging system equipped with a common sensor. This is due to the higher conductivity of graphene that leads to a higher signal produced by the RGO sensor to be radiated towards the tumor. The images demonstrated the efficiency of RGO sensor by perfect and clear images with better pixel intensity. The color with the most intensity (reddish-orange area) represents the presence of a tumor in the image because of the strongest scattering area that matches up to tumor presence in the human head phantom. Table 1 demonstrates the localization error of the detected tumor in the image as compared with the tumor's actual position with the value of 0.2 cm or 2% of error which could be considered a tolerable error to prove the system is efficiently and accurately detecting the desired tumor.

Conclusions
This paper presents a reduced graphene oxide (RGO) array sensor to detect a tumor in a human brain phantom via a confocal radar-based microwave imaging technique. The RGO patches are fed using transmission lines and coupled to a parasitic element etched using conventional copper on a Taconic substrate. The sensor operated at frequencies of 1.2 GHz to 10.8 GHz within the ultra-wide band (UWB) frequency. The sensor demonstrates a high gain at 14.5 dB and has a size of 90 mm × 45 mm 2 . This system is operating experimentally by using an in-house-developed head phantom. The quality of the resulting image can be achieved successfully by improved delay and sum algorithm (IDAS). For future work, the presence of multiple tumors should be clearly investigated in a larger experimental region.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.