Next Article in Journal
How Are Sports-Trackers Used by Runners? Running-Related Data, Personal Goals, and Self-Tracking in Running
Previous Article in Journal
Temperature Sensing Optimization for Home Thermostat Retrofit
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Method of Rice Moisture Content Determination Using Voxel Weighting-Based from Radio Tomography Images

by
Nurul Amira Mohd Ramli
1,2,
Mohd Hafiz Fazalul Rahiman
1,2,*,
Latifah Munirah Kamarudin
2,3,
Latifah Mohamed
1,2,
Ammar Zakaria
1,2,
Anita Ahmad
4 and
Ruzairi Abdul Rahim
4
1
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis, Pauh Putra Campus, Perlis 02600, Malaysia
2
Centre of Excellence for Advanced Sensor Technology (CEASTech), Universiti Malaysia Perlis, Pauh Putra Campus, Perlis 02600, Malaysia
3
Faculty of Electronic Engineering Technology, Universiti Malaysia Perlis, Pauh Putra Campus, Perlis 02600, Malaysia
4
School of Electrical Engineering, Faculty of Engineering, Universiti Teknologi Malaysia, Johor 81310, Malaysia
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(11), 3686; https://doi.org/10.3390/s21113686
Submission received: 6 April 2021 / Revised: 17 April 2021 / Accepted: 21 April 2021 / Published: 26 May 2021
(This article belongs to the Section Sensor Networks)

Abstract

:
This manuscript presents a new method to monitor and localize the moisture distribution in a rice silo based on tomography images. Because the rice grain is naturally hygroscopic, the stored grains’ quality depends on their level of moisture content. Higher moisture content leads to fibre degradation, making the grains too frail and possibly milled. If the moisture is too low, the grains become brittle and are susceptible to higher breakage. At present, the single-point measurement method is unreliable because the moisture build-up inside the silo might be distributed unevenly. In addition, this method mostly applies gravimetric analysis, which is destructive. Thus, we proposed a radio tomographic imaging (RTI) system to address these problems. Four simulated phantom profiles at different percentages of moisture content were reconstructed using Newton’s One-Step Error Reconstruction and Tikhonov Regularization algorithms. This simulation study utilized the relationship between the maximum voxel weighting of the reconstructed RTI image and the percentage of moisture content. The outcomes demonstrated promising results, in which the weighting voxel linearly increased with the percentage of moisture content, with a correlation coefficient higher than 0.95 was obtained. Therefore, the results support the possibility of using the RTI approach for monitoring and localizing the moisture distribution inside the rice silo.

1. Introduction

In Malaysia, rice is a symbol of the local culture, and is consumed daily either as cooked rice or indirectly in the form of rice flour. Despite various initiatives that have been implemented to help expand the national rice production [1], a scenario of increasing demand and declining supply has been forecasted for the next ten years [2,3]. In analysis of the challenges of this agricultural industry, researchers have found that appropriate silo facilities and effective farm management practices should be emphasized to guarantee an adequate supply of domestic production [4,5].
In general, the primary function of a silo is to protect grain crops from any elements, particularly moisture build-up, fungal load, and pest infestation [6,7,8]. However, statistics record that 10–30% of total paddy values (quality and quantity) are degraded due to poor on- and off-farm storage operations [9]. Because the paddy and rice grains are naturally hygroscopic, their value post-harvesting primarily depends on their level of moisture content. The literature has also acknowledged that moisture content is the prominent factor that determines both the proper time for harvesting and the requirement for safe storage conditions [10,11,12].
Following the mechanization of agricultural production, heated forced air dryers were developed to maintain the moisture level inside the silo as a means to preserve the quality and longevity of the grains. Nonetheless, in Malaysia, the storage of rice grains is challenging due to the hot, humid climate, which has a relative humidity of approximately 80%. According to [9,13,14], to preserve the rice quality and allow long-term storage, grains must be dried to 11–14% moisture content. Grains with higher moisture content tend to be frail and may be pulverized, and encourage the development of fungal invasion. In contrast, if the moisture content is too low, the grains become brittle and are exposed to a higher rate of breakage [15,16].
At present, moisture content is usually measured using a straightforward, fast, and convenient moisture meter [17,18]. The sensing application depends on the relationship of the material’s dielectric properties and the frequency of interest [19]. Although this type of instrument provides rapid measurement, it is not practical for moisture sensing of bulk grain in a silo. This failure is due to the single-point measurement, meaning that the measured moisture content does not represent the location and distribution of moisture at the bulk grain or silo scale [20]. Moreover, this method usually applies gravimetric analysis. This destructive process eventually reduces the grains’ quantity throughout the storage period.
Therefore, a comprehensive sensing method is essential to improve the efficiency of moisture distribution assessment in rice silos. Our systematic study presents a simulation of an RTI technique, which can help to visualize the location and percentage of the moisture content. The proposed method focuses on image analysis by employing simulation procedures for data collection, a tomographic technique based on algorithms, and linear regression to evaluate the effectiveness of the proposed method. Consequently, the findings allowed us to explore the feasibility of a tested RTI system as a tool for sensing changes of moisture content, with the aim of establishing a new approach for online, off-site, reliable, and non-invasive sensing methods.

2. Methodology

Tomography is a unique approach that is used to reconstruct a cross-sectional image of an object’s internal structure through sensory data obtained by the system [21]. The tomography system does not require invasion of the object of interest (OI), and was designed to analyze the monitored area’s internal composition using penetrating waves to calculate virtual cross-sections. Traditionally, the system’s structure can be visualized as three main components, which are the sensory system, data acquisition, and image reconstruction and display system [22]. The deployed sensor structure acts as the core input to the system, which is built based on the suitability of surroundings, nature of the OI, and the resolution quality of the information required by the user [23,24].
Various sensory systems in tomography applications have been used to localize moisture distribution, for example, electrical capacitance tomography (ECT), electrical resistance tomography (ERT), and microwave tomography (MWT). ECT and ERT are known as electrical tomography, which was developed to visualize the permittivity and conductivity distribution inside the OI by measuring a set of capacitance and resistance between inter-electrodes mounted around the periphery of the OI [25]. In contrast, the MWT technique reconstructs the tomogram based on the measurement of the scattered electromagnetic field from the arrayed sensor nodes [26]. Despite having different sensing elements, researchers have pursued the optimal tomography method to measure the moisture distribution in numerous frameworks.
In the construction industry, previous works as in [27,28] provided a non-destructive tool for approximately imaging the flow, shape, and position of the moisture distribution inside a cement-based material using ERT and ECT techniques. Refrence [29] implemented the ERT system in their geophysics study, based on the detection of moisture availability in the Mediterranean area. In addition, study by [30] found that the ECT system is feasible for imaging the moisture content information in wood. Overall, it can be concluded that, although electrical tomography is laborious, it can successfully deliver the required information about the OI.
However, in the agricultural sector, and particularly for sensing of grain moisture, the MWT system is usually preferred. The MWT system is based on the relationship between electromagnetic behaviour and the material’s relative permittivity [26]. Regardless of the successful application of the MWT technique [31], it is well known to be complicated and costly. This is because the system must be calibrated for dissimilar materials or if more than one variable is encountered in the framework [26].

Radio Tomographic Imaging (RTI) Approach

As a promising device-free localization technology, the concept of the RTI technique for monitoring and locating targets was initially proposed by [32]. This technology enables the location of targeted objects using image reconstruction based on changes in received signal strength (RSS) of the radio frequency (RF) signals between each stationary sensor node link in a wireless network area [33,34]. It is understood that when an object obstructs the transmission links, the RSS quality of the associated links experiences significant loss, whereas unblocked links are unaffected. Thus, within a monitored RF sensor network, the RTI system determines the targeted object’s location by reconstructing the RSS attenuation map across the sensor network. Figure 1 below illustrates the conceptual view of the RTI system.
In this work, the RTI system was studied in a finite element modelling (FEM) simulation platform. The structure of the geometrical dimension, number of RF nodes, signal frequency, and material dielectric properties were considered thoroughly. The sensor network consists of 20 RF nodes deployed on a conventional 2D silo to test the capability of the RTI technique for moisture distribution sensing in a rice silo. The moisture profiles were reconstructed using appropriate image reconstruction algorithms applied to the entire attenuation maps. The mapping process was relatively simple because it depends on the linearization and weighting model for RTI. In solving the forward problem, the measuring strategy was based on the linearization of its normalized weighted back-projected (Jacobian) matrix from each measurement [25,35].

3. Finite Element Modelling of RTI System

It is essential to study the behaviour of RF energy in a rice medium to control for the confounding variables and ensure the repeatability of obtaining RF attenuation for each scenario in the study. Figure 2 demonstrates the configuration setup of the FEM simulation. The RF sensor network comprising 20 nodes was configured on a 2D square silo. Each node is equally arrayed 40 mm apart, along the perimeter of a 0.5 m × 0.5 m modelled silo with an operating frequency of 2.4 GHz. The 2.4 GHz frequency was chosen to utilize the low-cost and standard ISM band of the wireless sensor network (WSN). The nodes were arranged in such a way to receive the scattered electric fields in multiple directions. The number of pixels that reflect the image’s resolution was set to 300 × 300 pixels. The following subsections comprehensively explain the design and selection of the simulation parameters.

3.1. Modelling Electromagnetic Behaviour

The interaction between an electromagnetic wave and the material’s dielectric properties has led to the development of a wide range of practical solutions in remote sensing. The remarkable features of electromagnetic waves when propagating in a dielectric medium include the carrying by the travelling energy of measurable electric ( E ) and magnetic ( H ) fields [36,37]. The description of electromagnetism is expressed by Maxwell’s equation below.
Faraday s   Law :   × E ( r ) = j ω B ( r )
Ampere s   Law :   × H ( r ) = j ω D ( r ) + J ( r )
Gauss s   Law :   · B ( r ) = 0
Gauss s   Law :   · D ( r ) = ρ
where E is the electric field intensity ( V m ) , H is the magnetic field intensity ( A m ) , B is the magnetic flux density ( T ) , D is the electric flux density ( C m 2 ) , J is the electric current density ( A m 2 ) , ρ is the electric charge density ( C m 3 ) , and r indicates the position vector.
These equations depend indirectly on e j w t , where j = 1 denotes the imaginary units. The constitutive relations have been applied in the existence of materials in which electromagnetic phenomena take place.
D ( r ) = ε o ε r E ( r )
J ( r ) = σ E ( r )
B ( r ) = μ o μ r H ( r )
where ε o is the permittivity of free space ( F m ) , ε r is the material’s relative permittivity, μ o is the permeability of free space ( H m ) , and μ r is the material’s relative permeability.
Below is the derivation of the equation used in the FEM study. Equation (7) is substituted into Equation (1) and curl on both sides, as shown in Equation (8).
× ( × E ( r ) ) = j ω μ o μ r ( × H ( r ) )
Substituting Equations (2), (5), and (6) into Equation (8):
× ( × E ( r ) ) = j ω μ o μ r ( j ω ε o ε r E ( r ) + σ E ( r ) )
× ( × E ( r ) ) = ω 2 μ o μ r ε o ε r E ( r ) j ω μ o μ r σ E ( r )
Equation (10) is simplified into Equations (11) and (12) as below.
× ( × E ( r ) ) = k o 2 μ r ( ε r j ( σ ε o ω ) )
× μ r 1 ( × E ( r ) ) = k o 2 ( ε r j ( σ ε o ω ) )
where k o = ω ε o μ o is the wavenumber in free space.

3.2. Modelling and Simulation: Rice Grains

As discussed above, the complex dielectric properties of rice grains, ε r , are the parameters that vary significantly with exposed electromagnetic waves. Dielectric properties, or permittivity, are the electrical characteristics of poorly conducting materials. Their dependence on the material’s moisture content (MC) has been of interest due to their practicality for rapid moisture sensing [38]. When interacting with the RF, these properties are studied as the parameters that can be polarized by the electric field. Thus, this RTI system is based on the relationship between electromagnetic energy and the medium’s dielectric properties. The relative complex permittivity of the rice is expressed in Equation (13) [19].
ε r = ε j ε
where ε is the dielectric constant that describes the capability of a material to store energy in the electric field, and ε is the dielectric loss factor that indicates the capability of a material to dissipate energy from the electric field, which is then converted into heat energy.
Furthermore, the conductivity of the material’s dielectric σ cannot be neglected, and is defined by Equation (14).
σ = ω ε o ε
where ω = 2 π f is the angular frequency with f in Hz, and ε o is the permittivity of free space, which is equivalent to 8.854 × 10 12 F/m.
The nature of the correlation between dielectric properties and changes in MC of rice grains is well known [19,39]. Apparently, at any applied frequency, the dielectric constant, ε , is almost linearly correlated with rice grain MC. The loss factor, ε , has been found to be less predictable, and may increase or decrease, depending on the particular range of MC.
The modelling of the rice grains was separated into two parts, which are the dried rice as a reference medium and the phantoms that indicate the rice with ‘X’ MC. The parameters for both conditions, such as dielectric properties, size, and shape, were defined with certanity to understand the ordinary deterministic electromagnetic scattering. Because the nature of moisture distribution is unpredictable [40], this initial study presents a simple geometric shape for the rice phantoms using a built-in primitive object. This commonly used shape can be modified based on the desired dimension and position.
The dried rice medium was assumed to be at 14% MC based on the standard level of MC for rice storage following the drying process. For the rice phantoms, the MC was set incrementally by 2% MC, starting from 16% up to 28%. The upper value of 28% MC was chosen to demonstrate the worst condition for rice grains, because the maximum MC of paddy grains post-harvesting process is about 25% [9]. Equations (15) and (16) below represent the correlation between complex dielectric properties of rice grains, ε and ε , with MC, based on the regression method modelled by [41,42,43], using the average of ten common types of rice grain in Malaysia’s marketplace.
ε = a 0 + a 1 M C + a 2 M C 2
ε = a 0 + a 1 M C + a 2 M C 2
where the value of coefficients a 0 , a 1 , a 2 and a 0 , a 1 , a 2 are for the polynomial expression of the dielectric constant, ε , and dielectric loss factor, ε , respectively.
The coefficients in Equations (15) and (16) are expressed as a function of frequency, f , as tabulated in Table 1 below [41,42,43].
The particular electrical conductivity, σ , for the dry medium and rice phantom ‘X’ MC, was calculated using Equation (14). In this simulation study, the temperature was presumed and fixed at 26 °C, because a change in temperature causes the grains’ dielectric properties to vary slightly [44]. In addition, in Malaysia, to preserve the quality of paddy and rice grains, the domestic crops are temporarily stored in warehouses under an ambient temperature range between 26 and 28 °C [3].

3.3. Modelling and Simulation: Rice Silo

Based on Figure 2, the rice silo was designed as a square shape made of an acrylic plate with a wall thickness of 10 mm. Because the scattering mechanism may include wave diffraction, the silo size was considered to be sufficiently substantial relative to the wavelength of the propagating frequency to avoid diffraction of the electromagnetic waves between the sensor nodes [45]. However, the attenuation of the electromagnetic waves increases with the increasing size of the medium [46].
As mentioned above, the dimensions of the silo used in this simulation study were 0.5 × 0.5 m. These dimensions were specifically chosen to investigate the capability of the RTI technique to image the moisture distribution inside the silo. The dielectric properties of the silo material were chosen from the predefined FEM input model. Polymethyl methacrylate (PMMA) with relative permittivity and conductivity of 2.5 and 0.0068 S/m, respectively, was selected due to its popularity in various applications [47]. In addition, the same material will be used for our future experimental study.

3.4. Modelling and Simulation: RF Nodes

It is requisite to have a mathematical model for the RF sensor nodes so that the RTI system can accommodate simulation of silos of any size with a flexible number of sensor nodes. In this study, we considered a rectangular waveguide design for the RF nodes operating at 2.4 GHz and TE10 mode. To apply a scalar model for the RTI problem, it was assumed that under the TE10 mode, an incident wave is polarized along the z-axis where E ( x , y , z ) = E ( x , y ) e i k z z , which signifies that E x and E y are equal to zero [48]. The width of the sensor node, w , was modelled and calibrated through Equation (17) [48].
f c = c o 2 w ε r
where f c is the cut-off frequency depending on the width of the waveguide, w . c o is the speed of light in a vacuum, which is equivalent to 3 ×   10 8 ms−1. If the cut-off frequency is lower than the operating frequency, RF signals can propagate in the dried rice medium.
According to [49], in TE10 mode, the cut-off frequency must be higher than 1.22 GHz and less than 2.45 GHz. Herein, the width of the sensor node, w , was 60 mm at a cut-off frequency of 1.31 GHz. In addition, the node density also plays a key role in the accuracy of the RTI results. In the same area of geometry, the imaging result is expected to be more accurate where nodes are placed more closely together than in areas where the nodes are spaced at a larger distance [35,50]. Our FEM simulation used 20 RF nodes in the square network within an imaging area of 0.25 m2. Although a greater number of nodes could attain better images, we considered the size of the silo.
Furthermore, the transmission power was restricted to 1 mW. When dealing with electromagnetic energy, heating of dielectric materials is subject to the material density and the power source [37]. Therefore, the designated transmission power should be carefully considered, because it may lead to changes in the physicochemical properties of the rice in practice [51].

4. Solving Forward Problem

4.1. Measuring Strategy

The basic linear formulation model was applied to determine the scattered electric field measurement values for known attenuation coefficients on each discretized pixel weight. The discretized electric field was solved based on the distribution of dielectric properties of the dried rice, rice phantoms, and rice silo. This technique is accessible because the dielectric properties are capable of manipulating the behaviour of the electric field.
In the RTI system, the electric field was excited within the network region through the rice phantom with ‘X’ MC. An array of 20 RF nodes were utilized to determine the signal strength from each transmitting and receiving node. The transmission links between node 1 and 2, and between node 2 and 1, in the network area are the same because any pair of nodes is regarded as a link, regardless of whether communication exists between them. Therefore, if K is the number of sensor nodes deployed in the RTI system, then the total number of unique two-ways links, M , is expressed as Equation (18) [32].
M = K 2 + K 2
If node 1 is assigned to function as transmitter, then the remaining nodes, 2, 3, 4, …, K = 20 , are allocated to serve as receivers. The procedure was sequentially repeated until node K = 20 was transmitting, and nodes 1 to 19 were receiving. Hence, with a K number of 20 nodes, the maximum number of measured radio links, M , would be M = 190 .

4.2. Weighting Model

From the local linearization, 190 sensitivity maps, also known as the Jacobian matrix, J , were computed from the electric field strength distribution, E ( x , y ) . One element, J j i , of this matrix defines the electric field, E ( x , y ) , and changes at a measurement area due to a slight perturbation in the dielectric properties of one element, ε r ( x , y ) , of the model [48]. This is stated as in Equation (19).
J j i = E ( x , y ) ε r ( x , y )
Each Jacobian matrix was derived from the FEM mesh of the forward model without the elements established in the dried rice (reference medium) and RF nodes. In this study, the electric field distribution was computed on the 11,000 finer tetrahedral mesh. The grid size for each Jacobian matrix was 0.002 m, and the number of grid pixels in the computational domain was 300 × 300. Figure 3 shows the normalized sensitivity maps of excitation at RF nodes, j = 3 , with the corresponding adjacent and opposite measurement at RF node i = 4 and i = 13 , respectively. It can be seen that, as the distance between two nodes increases, the sensitivity decreases.
In the RTI technique, the signals’ attenuation due to the phantom’s presence can be integrated by uniformly divided the imaging area into many small square pixels, N [52]. Because the attenuation of a link contributes differently to each pixel value, a weighting model was applied. Therefore, in this simulation study, the weight balanced map could be approximated as the sum of multiplication of the Jacobian matrix, J j i , by its corresponding signal loss, L j i , as shown in Equation (20).
W N = j = 1 T i = 1 R L j i × J j i
where L j i = S i j t o t S i j i n c , T denotes the number of excitation nodes, R denotes the number of measuring nodes, S i j t o t is the electric field measured when there is a rice phantom with ‘X’ MC, and S i j i n c is the electric field measured before the presence of the rice phantom.

5. Solving Inverse Problem

Image reconstruction aims to solve the inverse problem by generating the real image, which determines the unknown electric field distribution vector in the FEM domain. Although this technique does not provide precise information about the phantom’s location (moisture distribution), the electric field does represent some changes in the distribution of the different dielectric properties, which later enable reconstruction of rice phantoms [22]. This study applied two qualitative imaging algorithms, namely Newton’s One-Step Error Reconstruction (NOSER) and Tikhonov Regularization.

5.1. NOSER Algorithm

NOSER is a fast and practical image reconstruction algorithm due to its stability and direct linear method [53]. The formulation of NOSER in this RTI system is as shown in Equation (21).
G N O S E R ( x , y ) = [ H ] T × j = 1 T i = 1 R L S j i × J j i ( x , y )
where   L S j i = [ S i j t o t S i j i n c ] 2
and   H = J j i ( x , y ) [ J j i ( x , y ) ] T
where H is the Hessian matrix and J j i ( x , y ) is the computed Jacobian matrix corresponding to its signal loss, L j i as described above.
L j i is calculated using the least square method. The least square method aims to minimize the variation between the S i j i n c and S i j t o t for a given electric field distribution in the computational domain. In addition, the transposed Hessian matrix is utilized as a rough approximation instead of using inverse because it is impossible to compute the direct inverse [54].

5.2. Tikhonov Regularization Algorithm

Tikhonov Regularization is a commonly used method that provides a simple framework for integrating desired characteristics into the RTI reconstruction. In this study, the summation of the Hessian matrix was decomposed into three constituent matrices, which were U , V , and Σ . U and V are the unitary matrices, and Σ is a diagonal matrix. The diagonal elements of Σ are singular values, ρ , as expressed in Equation (24).
ρ = d i a g ( Σ )
Using Equation (23), regularization involves introducing additional information into the mathematical model to handle these small singular values and stabilize the inverse problem [55].
T t i k h , κ = ρ ρ 2 + κ
where κ is the regularization parameter and κ > 0 .
The formula to reconstruct the image using the Tikhonov Regularization algorithm is shown in Equation (26) for this RTI system.
G t i k h ( x , y ) = U × V × d i a g ( T t i k h , κ ) × j = 1 T i = 1 R L S j i × J j i ( x , y )

6. Results and Discussion

The image reconstruction and analysis were based on the pixel values. These pixel values represent the density related to the scattered electric field distribution inside the imaging area. In this section, we discuss the imaging results of the static moisture phantoms, which involved different MC percentages, investigate the approximation of MC percentage based on the regression plot of maximum pixel values, and examine the feasibility of the developed RTI system for rice MC imaging. We include an assessment of image quality based on different image reconstruction algorithms, i.e., NOSER and Tikhonov Regularization.
The results of all reconstructed images were quantified using a Mean Structural Similarity Index (MSSIM) image quality assessment by comparing the reconstructed image with the reference image. The evaluation aimed to analyze the similarity between the two images in terms of structure, luminance, and contrast, which gives an output index in the range of 0 to 1 [56]. An output index approaching one indicates the reconstructed image is close to the reference image.
The FEM simulation studies were conducted using the selected frequency, phantom properties, and background values. Four different scenarios of moisture distribution with seven different sets of MC percentage each were considered to evaluate the performance of the RTI concept, as shown in Figure 4. The reconstructed images using NOSER and Tikhonov Regularization algorithms for Phantoms A, B, C, and D are presented in Table 2 and Table 3, where the highest and lowest constrasts are indicated in red and blue, respectively.
Using the Tikhonov Regularization algorithm, different regularization parameters yield different tomogram results. Therefore, the chosen parameter should provide good results based on the ability to distinguish the changes in electric field attenuation and the limitation of the iteration number. In this study, the regularization parameter, κ of κ = 10 16 , was selected as the optimum parameter to reconstruct the Phantoms A, B, C, and D at MC percentages of 16%, 18%, 20%, 22%, 24%, 26%, and 28%.
Table 2 and Table 3 illustrate that the simulation results exhibit successful image reconstruction, which supports the possibility of using RTI to monitor and localize the moisture distribution in a rice silo. The rice phantoms reconstructed using both NOSER and Tikhonov Regularization algorithms possessed a position corresponding to the original simulated phantoms’ profiles for seven different MC percentages. Despite the presence of smearing-effect artefacts, the phantom shape and position can be rocognized and distinguished from the dried rice context.
Although both NOSER and Tikhonov Regularization approaches applied the least square method, NOSER does not perform an iterative procedure, which leads to shorter reconstruction time. Tikhonov Regularization iteratively involves additional regularization parameters into the mathematical model, thus required extra calculation time. In general, we noticed that the image reconstruction results using Tikhonov Regularization were preferable to those of NOSER, for which the regularization leads to smoothing of the reconstructed profiles. Hence, the selection of the algorithm should be task-specific.

6.1. Analysis of Maximum Pixel Value

Before the measurement and image reconstruction procedures, the reference medium was represented by dried rice at its respective dielectric properties, and the reference measurement was conducted. Then, a different permittivity of rice phantom at ‘X’ MC was introduced into the reference medium. The reconstruction algorithm was used to analyze the difference between the pixels’ density data for the scenario with and without phantoms, and to calculate a qualitative approximation of the permittivity distribution. Therefore, the maximum (max) weighting pixel of the reconstructed images should indicate its corresponding ‘X’ MC.
Based on this understanding, we utilized the max pixel value from the reconstructed image’s grid to form a linear regression plot. Because NOSER and Tikhonov Regularization algorithms offered different baselines of the max pixel value, the data was quantified using the percentage of the max pixel value from the tomogram of the reference medium in which the rice phantom was absent. As in the dried rice medium, the max pixel value was equivalent to 16.89. After the calculation, the max pixel value in terms of percentage for both algorithms can be appropriately compared.
Table 4 and Table 5 below show the percentage of max pixel value as a function of MC using NOSER and Tikhonov Regularization for MC at 16% to 28%, respectively. The tabulated data are plotted, as depicted in Figure 5 and Figure 6.
All of the reconstructed images in Table 2 and Table 3 show the existence of rice phantoms with ‘X’ MC by indicating changes in the color scale. However, these tomograms only exhibit positive results on the location and size of the moisture distribution. In addition, due to the smearing effects, it was difficult to determine the changes in terms of MC intensity. Therefore, the analyses, as illustrated in Figure 5 and Figure 6 are at least measurable.
The relationship between max pixel value and rice MC percentage can be quantitatively reconstructed using a linear regression approach. Both reconstruction algorithms, NOSER and Tikhonov Regularization, indicated that the max pixel value increases steadily as the rice MC increases. Because the max pixel value increases linearly with the percentage of MC, a first-order polynomial regression model was preliminary predicted, as shown in Equation (27).
P V m a x = a M C b
where a and b are the coefficients specific to each regression.
The statistical regression analyses are summarized in Table 6 and Table 7 for NOSER and Tikhonov Regularization algorithms, respectively. Based on the corresponding tables, the results indicate that the percentage of rice MC is the most significant parameter affecting the max pixel value, with a correlation coefficient, R 2 , higher than 0.95 regardless of the location and size of the rice phantoms used in this study.
Although the correlation coefficient, R 2 , results using the NOSER algorithm were slightly higher than those of the Tikhonov Regularization algorithm, we concluded that the Tikhonov Regularization is the preferred approach because it yielded the most stable linear regression model, as seen in Figure 5 and Figure 6. It is clear that the ill-posed inverse problem of RTI is stabilized by introducing a regularization parameter into the mathematical framework.

6.2. MSSIM Analysis

To quantitatively compare different algorithms, the MSSIM Index was calculated for the reconstructed images concerning the real distribution. The data is tabulated in Table 8 and presented in Figure 7.
As presented in Figure 7, we can observe that the Tikhonov Regularization algorithm recorded the highest MSSIM Index for all four moisture distribution profiles. The highest and lowest MSSIM Indexes scored by the Tikhonov Regularization algorithm were 0.7211 for Phantom A at MC of 26%, and 0.3092 for Phantom D at MC of 18%, respectively. For the MSSIM Index scored by the NOSER algorithm, the highest and lowest indexes were 0.5389 for Phantom A at MC of 28%, and 0.3040 for Phantom D at MC of 20%, respectively. Therefore, it can be deduced that the Tikhonov Regularization technique is better in reconstructing the moisture distribution in the rice silo.
Overall, moderate MSSIM indexes were obtained by the RTI system. This might be attributed to the difference in dielectric properties of dried rice and rice with ‘X’ MC, which causes a relatively small re-polarization (to the initial field) of an incident field, E i n c . However, the results demonstrate the possibility of implementing the RTI system for rice moisture monitoring and localization.

7. Conclusions

In conclusion, monitoring and localizing the moisture distribution in a rice silo is essential in the agricultural industry to ensure the sustainability of domestic production. In this preliminary study, the feasibility of an RTI system for rice moisture detection was successfully simulated, and the analysis of its performance was appropriately carried out. The comprehensive studies on the proposed FEM simulation data confirmed the findings and ensure that the results were not artefacts due to simulation errors.
The results indicate that the percentage of rice MC is a parameter that has a significant effect on the max pixel value, with a correlation coefficient, R 2 , higher than 0.95. Therefore, this study provides evidence to support the possibility of using the RTI technique for monitoring the rice moisture distribution, based on the study of maximum pixel values from the reconstructed tomogram. Furthermore, we also evaluated MSSIM Indexes to quantitatively compare the two algorithms. The findings showed that both NOSER and Tikhonov Regularization algorithms performed good image reconstruction with moderate MSSIM Indexes.
Our assessments presented in this paper demonstrate that the Tikhonov Regularization algorithm possessed a more stable technique in monitoring and localizing the rice moisture distribution in the RTI system. This stability is due to the incorporation of the regularization parameter into the mathematical model. Nevertheless, the results of this initial study could perhaps serve as a non-destructive tool to monitor and localize the moisture distribution in rice silos.
Finally, regarding an actual case study, further research work should be conducted to explore the reliability, sensitivity, and effectiveness of the proposed image reconstruction technique. It is recommended that future studies focus on the development of pre-processing of the sensory data [57,58] because the reconstructed images may be affected by the intrinsic behavior of the images due to factors related to noise and the environment.

Author Contributions

Conceptualization, N.A.M.R. and M.H.F.R.; methodology, N.A.M.R. and M.H.F.R.; software, N.A.M.R. and M.H.F.R.; writing—original draft preparation, N.A.M.R. and M.H.F.R.; writing—review and editing, N.A.M.R., M.H.F.R., L.M.K., L.M., and A.Z.; supervision, M.H.F.R.; funding acquisition, N.A.M.R., M.H.F.R., L.M.K., L.M., A.Z., A.A., and R.A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was funded by the Ministry of Higher Education (MOHE) Malaysia under the grant Transdisciplinary Research Grant Scheme (Ref: TRGS/1/2018/UNIMAP/02/4/2) titled “A Localization of Moisture Distribution in Agricultural Silos via Sensor Network and Artificial Intelligence Assisted”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated in this manuscript are available upon request, with provided approval by the corresponding author. In addition, the models and code used during the research work cannot be shared at present because the data forms part of an ongoing study.

Acknowledgments

The authors would like to thank Ministry of Higher Education (MOHE) Malaysia for funding this research under the Transdisciplinary Research Grant Scheme (Ref: TRGS/1/2018/UNIMAP/02/4/2). Special acknowledgement indeed goes to Universiti Malaysia Perlis for providing the facilities and technical support.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Siwar, C.; Mohd Idris, N.D.; Yasar, M.; Morshed, G. Issues and challenges facing rice production and food security in the granary areas in the East Coast Economic Region (ECER), Malaysia. Res. J. Appl. Sci. Eng. Technol. 2014, 7, 711–722. [Google Scholar] [CrossRef]
  2. Rajamoorthy, Y.; Rahim, K.b.A.; Munusamy, S. Rice industry in Malaysia: Challenges, policies and implications. Procedia Econ. Financ. 2015, 31, 861–867. [Google Scholar] [CrossRef] [Green Version]
  3. Che Omar, S.; Shaharudin, A.; Tumin, S.A. The Status of the Paddy and Rice Industry in Malaysia. Available online: www.krinstitute.org (accessed on 2 December 2019).
  4. Espino, L.A.; Greer, C.A.; Mutters, R.G.; Thompson, J.F. Survey of rice storage facilities identifies research and education needs. Calif. Agric. 2014, 68, 38–46. [Google Scholar] [CrossRef] [Green Version]
  5. Pekmez, H. Cereal storage techniques: A review. J. Agric. Sci. Technol. B 2016, 6, 67–71. [Google Scholar]
  6. Bhandari, G.; Achhami, B.B.; Karki, T.B.; Bhandari, B.; Bhandari, G. Survey on maize post-harvest losses and its management practices in the Western Hills of Nepal. J. od Maize Res. Dev. 2015, 1, 98–105. [Google Scholar] [CrossRef] [Green Version]
  7. Jyoti. Effect of treatment, packing material and storage on viability in paddy (Oryza sativa L.). J. Pharmacogn. Phytochem. 2017, 6, 962–964. [Google Scholar]
  8. Poojari, P.; Krishna, A. Influence of seed storage condition on seed moisture content and germination in impatiens talbotii hook. J. Plant Dev. Sci. 2018, 10, 435–443. [Google Scholar]
  9. Postharvest Unit, C. Grain Storage and Pest Management. Available online: www.knowledgebank.irri.org (accessed on 3 December 2019).
  10. Putri, R.E.; Yahya, A.; Adam, N.M.; Aziz, S.A. Related fracture resistance with moisture content in different grain orientation of paddy grain. J. Biol. Agric. Healthc. 2015, 5, 64–70. [Google Scholar]
  11. Furahisha, K.; Chove, L.M.; Chaula, D. Effect of final moisture content, cooling time and paddy variety on milling quality of rice (Oryza sativa L.). J. Agric. Sci. Food Technol. 2016, 2, 169–179. [Google Scholar]
  12. Naik, D.S.; Chetti, M.B. Influence of packaging and storage conditions on the moisture content and its effect on fungal load of paddy. Res. J. Agric. Sci. 2017, 8, 370–374. [Google Scholar]
  13. Fu, H.; Cao, D.D.; Hu, W.M.; Guan, Y.J.; Fu, Y.Y.; Fang, Y.F.; Hu, J. Studies on Optimum Harvest Time for Hybrid Rice Seed. J. Sci. Food Agric. 2017, 97, 1124–1133. [Google Scholar] [CrossRef]
  14. Alhendi, A.S.; Al-Rawi, S.H.; Jasim, A.M. Effect of moisture content of two paddy varieties on the physical and cooked properties of produced rice. Braz. J. Food Technol. 2019, 22. [Google Scholar] [CrossRef] [Green Version]
  15. Alizadeh, M.R.; Dabbaghi, A.; Rahimi-Ajdadi, F.R. Effect of final paddy moisture content on breaking force and milling properties of rice varieties. Elixir Int. J. Agric. 2011, 36, 3186–3189. [Google Scholar]
  16. Shafiekhani, S.; Wilson, S.A.; Atungulu, G.G. Impacts of storage temperature and rice moisture content on color characteristics of rice from fields with different disease management practices. J. Stored Prod. Res. 2018, 78, 89–97. [Google Scholar] [CrossRef]
  17. Mbie, N.Z. APLMF Guide Document on Rice Moisture Measurement. Available online: http://www.aplmf.org/uploads/5/7/4/7/57472539/aplmf_guide_-_rice_moisture_measurement.pdf (accessed on 22 January 2020).
  18. Vera, M.; Dutta, B.; Mercer, D.G.; Maclean, H.L.; Touchie, M.F. Assessment of moisture content measurement methods of dried food products in small-scale operations in developing countries: A review. Trends Food Sci. Technol. 2019, 88, 484–496. [Google Scholar] [CrossRef]
  19. Nelson, S.O. Chapter 13—Dielectric properties models for grain and seed. In Dielectric Properties of Agricultural Materials and Their Application; Academic Press: Cambridge, MA, USA, 2015; pp. 175–193. [Google Scholar]
  20. Abdullah, M.S.M.; Rahiman, M.H.F.; Zakaria, A.; Kamarudin, L.M.; Mohamed, L. A review on moisture measurement technique in agricultural silo. IOP Conf. Ser. Mater. Sci. Eng. 2019, 705, 47–51. [Google Scholar] [CrossRef]
  21. Jamaludin, J.; Zawawir, M.Z.; Abdul Rahim, R.; Mohd Yunos, F.R.; Nor Ayob, N.M.; Ridzuan Aw, S.; Fadzil, N.S.; Zakaria, Z.; Fazalul Rahiman, M.H. A review of tomography system. J. Teknol. Sci. Eng. 2013, 64, 47–51. [Google Scholar] [CrossRef] [Green Version]
  22. Beck, M.S.; Williams, R.A. Process tomography: A European innovation and its applications. Meas. Sci. Technol. 1996, 7, 215–224. [Google Scholar] [CrossRef]
  23. Abd Rahman, N.A.; Hong, L.E.; Abdul Rahim, R.; Abdul Rahim, H.; Ahmad, N.; Bunyamin, S.; Abas, K.H.; Nor Ayob, N.M.; Mohd Yunos, F.R.; Badri Mansor, M.S. A review: Tomography systems in medical and industrial processes. J. Teknol. 2015, 73, 1–11. [Google Scholar] [CrossRef] [Green Version]
  24. Yao, J.; Takei, M. Application of process tomography to multiphase flow measurement in industrial and biomedical fields: A review. IEEE Sens. J. 2017, 17, 8196–8205. [Google Scholar] [CrossRef]
  25. Cui, Z.; Wang, Q.; Xue, Q.; Fan, W.; Zhang, L.; Cao, Z.; Sun, B.; Wang, H.; Yang, W. A review on image reconstruction algorithms for electrical capacitance/resistance tomography. Sens. Rev. 2016, 36, 429–445. [Google Scholar] [CrossRef]
  26. Rahiman, M.H.F.; Kiat, T.T.W.; Jack, S.P.; Rahim, R.A. Microwave tomography application and approaches—A review. J. Teknol. 2015, 73, 133–138. [Google Scholar] [CrossRef] [Green Version]
  27. Hallaji, M.; Seppänen, A.; Pour-Ghaz, M. Electrical resistance tomography to monitor unsaturated moisture flow in cementitious materials. Cem. Concr. Res. 2015, 69, 10–18. [Google Scholar] [CrossRef]
  28. Voss, A.; Pour-Ghaz, M.; Vauhkonen, M.; Seppänen, A. Electrical capacitance tomography to monitor unsaturated moisture ingress in cement-based materials. Cem. Concr. Res. 2016, 89, 158–167. [Google Scholar] [CrossRef]
  29. Nijland, W.; Van der Meijde, M.; Addink, E.A.; De Jong, S.M. Detection of soil moisture and vegetation water abstraction in a mediterranean natural area using electrical resistivity tomography. Catena. 2010, 81, 209–216. [Google Scholar] [CrossRef]
  30. Voss, A.; Seppänen, A.; Siltanen, S.; Salokangas, L.; Baroudi, D. Imaging of moisture content in wood using electrical capacitance tomography. In Proceedings of the World Confrence on Timber Engineering, Vienna, Austria, 22–25 August 2016; pp. 22–25. [Google Scholar]
  31. Semenov, S.N. Imaging the dielectric objects by microwave tomography method. St. Petersbg. Polytech. Univ. J. Phys. Math. 2015, 1, 315–318. [Google Scholar] [CrossRef] [Green Version]
  32. Wilson, J.; Patwari, N. Radio tomographic imaging with wireless networks. IEEE Trans. Mob. Comput. 2010, 9, 621–632. [Google Scholar] [CrossRef] [Green Version]
  33. Yigitler, H.; Jantti, R.; Kaltiokallio, O.; Patwari, N. Detector based radio tomographic imaging. IEEE Trans. Mob. Comput. 2018, 17, 58–71. [Google Scholar] [CrossRef]
  34. Dharmadasa, M.A.I.M.; Gamage, C.D.; Keppitiyagama, C.I. Radio tomographic imaging (RTI) and privacy implications. In Proceedings of the 2018 18th International Conference on Advances in ICT for Emerging Regions, Colombo, Sri Lanka, 27–28 September 2018; pp. 413–419. [Google Scholar]
  35. Chiu, C.; Dujovne, D. Experimental characterization of radio tomographic imaging using tikhonov’s regularization. In Proceedings of the 2014 IEEE Biennial Congress of Argentina, Bariloche, Argentina, 11–13 June 2014; pp. 468–472. [Google Scholar]
  36. Baker-jarvis, J. The interaction of radio-frequency fields with dielectric materials at macroscopic to mesoscopic scales. J. Res. Natl. Inst. Stand. Technol. 2012, 117. [Google Scholar] [CrossRef]
  37. Mazima, J.K.; Johnson, A.; Manasseh, E. An overview of electromagnetic radiation in grain crops. Int. J. Food Sci. Technol. 2018, 1, 21–32. [Google Scholar]
  38. Nelson, S.O. Fundamentals of dielectric properties measurements and agricultural applications. J. Microw. Power Electromagn. Energy 2010, 44, 98–113. [Google Scholar] [CrossRef] [PubMed]
  39. Hou, K.M.; Kok, Y.Y.; Dimon, M.N. Rice grain moisture determination using microstrip wide-ring and microstrip coupled-line sensors. Am. J. Appl. Sci. 2015, 12, 112–120. [Google Scholar]
  40. Bartosik, R. An inside look at the silo-bag system. In Proceedings of the 9th International Conference on Controlled Atmosphere and Fumigation in Stored Products CAF 2012, Antalya, Turkey, 15–19 October 2012; pp. 117–128. [Google Scholar]
  41. You, K.Y.; Mun, H.K.; You, L.L.; Salleh, J.; Abbas, Z. A small and slim coaxial probe for single rice grain moisture sensing. Sensors 2013, 13, 3652–3663. [Google Scholar] [CrossRef] [Green Version]
  42. You, K.Y.; You, L.L.; Yue, C.S.; Mun, K.; Yue, C.S.; Mun, H.K.; Lee, C.Y. Physical and chemical characterization of rice using microwave and laboratory methods. In Rice—Technology and Production; Fahad, A.S., Ed.; InTech: Houston, TX, USA, 2017; pp. 81–99. [Google Scholar]
  43. Salleh, J. Development of Monopole Sensors for Rice Quality Characterization in Malaysia. Master’s Thesis, Universiti Teknologi Malaysia, Johor, Malaysia, July 2012. [Google Scholar]
  44. Park, C.E.; Kim, Y.S.; Park, K.J.; Kim, B.K. Changes in physicochemical characteristics of rice during storage at different temperatures. J. Stored Prod. Res. 2012, 48, 25–29. [Google Scholar] [CrossRef]
  45. Hu, G. Study on the signal transmission characteristics of 2.4GHz wireless network in dorms. Int. J. Online Eng. 2016, 12, 58–63. [Google Scholar] [CrossRef] [Green Version]
  46. Kiat, T.T.W.; Rahiman, M.H.F.; Jack, S.P.; Rahim, R.A. Initial study on microwave technique to find the relationship between attenuation of microwave and thickness of a subject. J. Teknol. 2015, 77, 99–103. [Google Scholar]
  47. Pawar, E.A. Review article on acrylic PMMA. ISOR J. Mech. Civ. Eng. 2016, 13, 1–4. [Google Scholar]
  48. Mallach, M.; Gebhardt, P.; Musch, T. 2D microwave tomography system for imaging of multiphase flows in metal pipes. Flow Meas. Instrum. 2016, 53, 80–88. [Google Scholar] [CrossRef]
  49. Kittiamornkul, N.; Jirasereeamornkul, K.; Kiattisin, S.; Chamnongthai, K.; Higuchi, K. The electric field distribution of a hybrid rectangular and circular waveguide resonator for use in granular material dielectric measurements. IEEJ Trans. Fundam. Mater. 2015, 135, 439–449. [Google Scholar] [CrossRef]
  50. Kiat, T.T.W. Simulation Study of Tomography for Agarwood Evaluation. Master’s Thesis, Universiti Malaysia Perlis, Perlis, Malaysia, September 2017. [Google Scholar]
  51. Ding, C.; Khir, R.; Pan, Z.; Wood, D.F.; Venkitasamy, C.; Tu, K.; El-Mashad, H.; Berrios, J. Influence of infrared drying on storage characteristics of brown rice. Food Chem. 2018, 264, 149–156. [Google Scholar] [CrossRef]
  52. Wang, Q.; Yiǧitler, H.; Jäntti, R.; Huang, X. Localizing Multiple Objects Using Radio Tomographic Imaging Technology. IEEE Trans. Veh. Technol. 2016, 65, 3641–3656. [Google Scholar] [CrossRef]
  53. Mallach, M.; Gevers, M.; Gebhardt, P.; Musch, T. Fast and precise soft-field electromagnetic tomography systems for multiphase flow imaging. Energies 2018, 11, 1199. [Google Scholar] [CrossRef] [Green Version]
  54. Yunos, Y.M.; Ruzairi Abd Rahim, R.G.; Ruzairi Abd Rahim, R.G.; Hafiz, M.; Rahiman, F. Image reconstruction using iterative transpose algorithm for optical tomography. J. Teknol. 2007, 47, 91–102. [Google Scholar]
  55. Beck, B.; Ma, X.; Baxley, R. Ultrawideband tomographic imaging in uncalibrated networks. IEEE Trans. Wirel. Commun. 2016, 15, 6474–6486. [Google Scholar] [CrossRef]
  56. Zhang, L.; Zhang, L.; Mou, X.; Zhang, D. FSIM: A feature similarity index for image quality assessment. IEEE Trans. Image Process. 2011, 20, 2378–2386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Kaur, P.; Chaira, T.A. Novel fuzzy approach for segmenting medical images. Soft Comput. 2021, 25, 3565–3575. [Google Scholar] [CrossRef]
  58. Versaci, M.; Morabito, F.C. Image edge detection: A new approach based on fuzzy entropy and fuzzy divergence. Int. J. Fuzzy Syst. 2021, 1–19. [Google Scholar] [CrossRef]
Figure 1. The conceptual view of the RTI system consisting an arrangement of 20 RF nodes, data acquisition process, and image reconstruction system.
Figure 1. The conceptual view of the RTI system consisting an arrangement of 20 RF nodes, data acquisition process, and image reconstruction system.
Sensors 21 03686 g001
Figure 2. Configuration of the RTI system.
Figure 2. Configuration of the RTI system.
Sensors 21 03686 g002
Figure 3. Normalized sensitivity maps of excitation at node, j = 3, with corresponding measurement at node i = 4 and i = 13.
Figure 3. Normalized sensitivity maps of excitation at node, j = 3, with corresponding measurement at node i = 4 and i = 13.
Sensors 21 03686 g003
Figure 4. Simulation of rice phantom (moisture distribution).
Figure 4. Simulation of rice phantom (moisture distribution).
Sensors 21 03686 g004
Figure 5. Percentage of maximum pixel value using NOSER algorithm.
Figure 5. Percentage of maximum pixel value using NOSER algorithm.
Sensors 21 03686 g005
Figure 6. Percentage of maximum pixel value using Tikhonov Regularization algorithm.
Figure 6. Percentage of maximum pixel value using Tikhonov Regularization algorithm.
Sensors 21 03686 g006
Figure 7. MSSIM Index graphs computed on the reconstructed images of Phantom A, B, C, and D.
Figure 7. MSSIM Index graphs computed on the reconstructed images of Phantom A, B, C, and D.
Sensors 21 03686 g007
Table 1. Quadratic coefficients of Equations (15) and (16).
Table 1. Quadratic coefficients of Equations (15) and (16).
Coefficient Functions
Dielectric constant, ε a 0 = 0.4605 f 2 + 8.5289 f 5.1708
a 1 = 0.0544 f 2 0.9928 f 0.2982
a 2 = 0.0016 f 2 + 0.0281 f + 0.0224
Dielectric loss factor, ε a 0 = 0.1347 f 2 + 0.6010 8.4663
a 1 = 0.0166 0.0693 f + 0.8215
a 2 = 0.0005 f 2 + 0.0023 f 0.0149
Table 2. Reconstructed images for Phantoms A and B at different percentages of MC.
Table 2. Reconstructed images for Phantoms A and B at different percentages of MC.
MC (%)Phantom APhantom B
NOSERTikhonov RegularizationNOSERTikhonov Regularization
16 Sensors 21 03686 i001 Sensors 21 03686 i002 Sensors 21 03686 i003 Sensors 21 03686 i004High
18 Sensors 21 03686 i005 Sensors 21 03686 i006 Sensors 21 03686 i007 Sensors 21 03686 i008 Sensors 21 03686 i009
20 Sensors 21 03686 i010 Sensors 21 03686 i011 Sensors 21 03686 i012 Sensors 21 03686 i013
22 Sensors 21 03686 i014 Sensors 21 03686 i015 Sensors 21 03686 i016 Sensors 21 03686 i017
24 Sensors 21 03686 i018 Sensors 21 03686 i019 Sensors 21 03686 i020 Sensors 21 03686 i021
26 Sensors 21 03686 i022 Sensors 21 03686 i023 Sensors 21 03686 i024 Sensors 21 03686 i025
28 Sensors 21 03686 i026 Sensors 21 03686 i027 Sensors 21 03686 i028 Sensors 21 03686 i029Low
Table 3. Reconstructed images for Phantoms C and D at different percentages of MC.
Table 3. Reconstructed images for Phantoms C and D at different percentages of MC.
MC (%)Phantom CPhantom D
NOSERTikhonov RegularizationNOSERTikhonov Regularization
16 Sensors 21 03686 i030 Sensors 21 03686 i031 Sensors 21 03686 i032 Sensors 21 03686 i033High
18 Sensors 21 03686 i034 Sensors 21 03686 i035 Sensors 21 03686 i036 Sensors 21 03686 i037 Sensors 21 03686 i038
20 Sensors 21 03686 i039 Sensors 21 03686 i040 Sensors 21 03686 i041 Sensors 21 03686 i042
22 Sensors 21 03686 i043 Sensors 21 03686 i044 Sensors 21 03686 i045 Sensors 21 03686 i046
24 Sensors 21 03686 i047 Sensors 21 03686 i048 Sensors 21 03686 i049 Sensors 21 03686 i050
26 Sensors 21 03686 i051 Sensors 21 03686 i052 Sensors 21 03686 i053 Sensors 21 03686 i054
28 Sensors 21 03686 i055 Sensors 21 03686 i056 Sensors 21 03686 i057 Sensors 21 03686 i058Low
Table 4. Percentage of maximum pixel value using NOSER algorithm.
Table 4. Percentage of maximum pixel value using NOSER algorithm.
MC (%)Phantom APhantom BPhantom CPhantom D
16 1.092 × 10 4 1.107 × 10 4 1.248 × 10 4 1.283 × 10 4
18 1.579 × 10 4 1.607 × 10 4 1.886 × 10 4 1.954 × 10 4
20 2.167 × 10 4 2.199 × 10 4 2.233 × 10 4 2.266 × 10 4
22 2.840 × 10 4 2.873 × 10 4 2.893 × 10 4 2.927 × 10 4
24 3.053 × 10 4 3.074 × 10 4 3.096 × 10 4 3.110 × 10 4
26 3.429 × 10 4 3.454 × 10 4 3.471 × 10 4 3.482 × 10 4
28 3.509 × 10 4 3.533 × 10 4 3.552 × 10 4 3.558 × 10 4
Table 5. Percentage of maximum pixel value using Tikhonov Regularization algorithm.
Table 5. Percentage of maximum pixel value using Tikhonov Regularization algorithm.
MC (%)Phantom APhantom BPhantom CPhantom D
16 6.334 × 10 14 6.393 × 10 14 6.542 × 10 14 6.630 × 10 14
18 1.290 × 10 13 1.302 × 10 13 1.415 × 10 13 1.540 × 10 13
20 3.114 × 10 13 3.114 × 10 13 3.208 × 10 13 3.108 × 10 13
22 4.682 × 10 13 4.724 × 10 13 4.730 × 10 13 4.753 × 10 13
24 5.452 × 10 13 5.475 × 10 13 5.505 × 10 13 5.529 × 10 13
26 6.570 × 10 13 6.630 × 10 13 6.630 × 10 13 6.689 × 10 13
28 6.511 × 10 13 6.570 × 10 13 6.570 × 10 13 6.570 × 10 13
Table 6. Correlation coefficients of the regression model using NOSER algorithm.
Table 6. Correlation coefficients of the regression model using NOSER algorithm.
PhantomRegression Model
P V m a x = a M C b
Correlation   Coefficient ,   R 2
A P V m a x = 2.113 E 3 M C 2.125 E 4 0.9550
B P V m a x = 2.115 E 3 M C 2.103 E 4 0.9534
C P V m a x = 1.954 E 3 M C 1.674 E 4 0.9583
D P V m a x = 1.915 E 3 M C 1.559 E 4 0.9541
Table 7. Correlation coefficients of the regression model using Tikhonov Regularization algorithm.
Table 7. Correlation coefficients of the regression model using Tikhonov Regularization algorithm.
PhantomRegression Model
P V m a x = a M C b
Correlation   Coefficient ,   R 2
A P V m a x = 5.452 E 14 M C 7.959 E 13 0.9529
B P V m a x = 5.502 E 14 M C 8.039 E 13 0.9536
C P V m a x = 5.447 E 14 M C 7.883 E 13 0.9545
D P V m a x = 5.449 E 14 M C 7.871 E 13 0.9553
Table 8. MSSIM Indexes computed on the reconstructed images.
Table 8. MSSIM Indexes computed on the reconstructed images.
MC (%)Phantom APhantom BPhantom CPhantom D
NOSERTikhonov RegularizationNOSERTikhonov RegularizationNOSERTikhonov RegularizationNOSERTikhonov Regularization
160.52210.66440.35730.41610.33160.36650.31180.3294
180.49120.59290.34190.35770.32460.33470.30910.3092
200.45960.61940.33590.38730.31710.36260.30400.3266
220.44770.60140.33870.38350.33400.38600.31890.3490
240.46730.67830.34950.39330.34640.40020.33080.3628
260.53490.72110.39070.47760.36400.43190.34970.3869
280.53890.70490.39470.45890.37030.42370.35580.3754
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mohd Ramli, N.A.; Fazalul Rahiman, M.H.; Kamarudin, L.M.; Mohamed, L.; Zakaria, A.; Ahmad, A.; Rahim, R.A. A New Method of Rice Moisture Content Determination Using Voxel Weighting-Based from Radio Tomography Images. Sensors 2021, 21, 3686. https://doi.org/10.3390/s21113686

AMA Style

Mohd Ramli NA, Fazalul Rahiman MH, Kamarudin LM, Mohamed L, Zakaria A, Ahmad A, Rahim RA. A New Method of Rice Moisture Content Determination Using Voxel Weighting-Based from Radio Tomography Images. Sensors. 2021; 21(11):3686. https://doi.org/10.3390/s21113686

Chicago/Turabian Style

Mohd Ramli, Nurul Amira, Mohd Hafiz Fazalul Rahiman, Latifah Munirah Kamarudin, Latifah Mohamed, Ammar Zakaria, Anita Ahmad, and Ruzairi Abdul Rahim. 2021. "A New Method of Rice Moisture Content Determination Using Voxel Weighting-Based from Radio Tomography Images" Sensors 21, no. 11: 3686. https://doi.org/10.3390/s21113686

APA Style

Mohd Ramli, N. A., Fazalul Rahiman, M. H., Kamarudin, L. M., Mohamed, L., Zakaria, A., Ahmad, A., & Rahim, R. A. (2021). A New Method of Rice Moisture Content Determination Using Voxel Weighting-Based from Radio Tomography Images. Sensors, 21(11), 3686. https://doi.org/10.3390/s21113686

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop