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].
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 , 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).
where
and
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,
, higher than 0.95 regardless of the location and size of the rice phantoms used in this study.
Although the correlation coefficient,
, 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, . 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, , 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.