Convolutional Neural Network-Processed MRI Images in the Diagnosis of Plastic Bronchitis in Children

Objective The study focused on the features of the convolutional neural networks- (CNN-) processed magnetic resonance imaging (MRI) images for plastic bronchitis (PB) in children. Methods 30 PB children were selected as subjects, including 19 boys and 11 girls. They all received the MRI examination for the chest. Then, a CNN-based algorithm was constructed and compared with Active Appearance Model (AAM) algorithm for segmentation effects of MRI images in 30 PB children, factoring into occurring simultaneously than (OST), Dice, and Jaccard coefficient. Results The maximum Dice coefficient of CNN algorithm reached 0.946, while that of active AAM was 0.843, and the Jaccard coefficient of CNN algorithm was also higher (0.894 vs. 0.758, P < 0.05). The MRI images showed pulmonary inflammation in all subjects. Of 30 patients, 14 (46.66%) had complicated pulmonary atelectasis, 9 (30%) had the complicated pleural effusion, 3 (10%) had pneumothorax, 2 (6.67%) had complicated mediastinal emphysema, and 2 (6.67%) had complicated pneumopericardium. Also, of 30 patients, 19 (63.33%) had lung consolidation and atelectasis in a single lung lobe and 11 (36.67%) in both two lung lobes. Conclusion The algorithm based on CNN can significantly improve the segmentation accuracy of MRI images for plastic bronchitis in children. The pleural effusion was a dangerous factor for the occurrence and development of PB.


Introduction
Plastic bronchitis (PB) is a rare respiratory disease arising from pathogenic infections. e endogenous foreign matter produced then leads to bronchial congestion, ventilatory disorder, and air exchange dysfunction [1]. Clinically, it manifests as dyspnea, wheezing, chest pain, fever, and even life-threatening respiratory circulatory failure in severe cases. Children aged 4-12 are predominantly affected [2]. According to the pathological characteristics, it mainly falls into type I and type II plastic bronchitis. Type I refers to the inflammatory type, where fibrous protein and a great number of inflammatory cells are visible in the bronchus, thanks to infectious diseases of the respiratory tract. Type II is cell-free, where the bronchus is mainly blocked with mucoprotein, without or with just a few infiltration cells, and a great amount of lymph fluid is visible in the alveolus [3]. e incidence of PB is increasing in China, and it is associated with the Mycoplasma pneumoniae infection or virus infection of the respiratory tract [4].
Electronic computed tomography (CT) is a routine examination to diagnose plastic bronchitis in children. With the continuous maturity of CT technology, CT has gradually become the first choice for lung examination and is a common examination method to determine various diseases. However, the effective radiation dose involved in CT examination is relatively harmful to children, and CT examination can only observe the macroscopic signs of the lesion site, unable to accurately determine the shape, size, and density, and the conclusion is relatively simple, which has certain limitations in the diagnosis of plastic bronchitis in children. Magnetic resonance imaging (MRI) belongs to positron emission tomography (PET), during which the radio-frequency pulse at a specific frequency is applied to the human body in a magnetostatic field, to motivate the hydrogen proton in the human body, triggering magnetic resonance. Next, the electromagnetic signal of the human body is used to reconstruct the image [5]. MRI mainly consists of three systems, namely, the magnet system, the spectrometer system, and the computer image reconstruction system [6,7]. During the examination, the patient is in a magnetic field and the orientation of magnetic moment is consistent with the direction of magnetic induction line. e radio-frequency pulses enable the lower-energy-level nuclei to jump to the higher-energy level by absorbing the radiofrequency energy and disrupting the proton movement in the tissue. When the pulses are stopped, the energy level and phase position of protons return to the premotivated state.
is process is called the relaxation [8]. T1-weighted image (T1WI) can display the anatomical details, while T2weighted image (T2WI) can display lesions. e black-white contrast of MRI stems from the signal differences in different tissues. Short T1 presents white, and long T1 presents black; long T2 presents white, and short T2 presents black. e air and bone cortex presents black no matter in T1 or T2. e infarction, inflammation, tumor, and fluid present lowsignal T1WI but high-signal T2WI [9]. MRI can display the anatomical structure of any section of the human body, with high resolution. It can clearly display the lesions in the spinal cord, brainstem, fossa cranii posterior, and bronchi. Additionally, MRI is free of ionizing radiation and carcinogenic risk. It is a very safe examination. However, MRI examination also has certain limitations, and image artifacts are often produced by the influence of respiration and blood flow, which hinder the application of MRI in lung imaging. High-quality images are beneficial to the diagnosis of lung diseases, which can fully reflect the details of the image and significantly improve the diagnostic coincidence rate. e convolutional neural network (CNN) is a feedforward neural network. It has a deep structure and contains convolutional calculation and is representative of deep learning [10]. In the medical field, the CNN is usually used to learn original MRI images and then applied in image segmentation and image classification. Sathish et al. (2019) [11] put forward in their study that the pixel information of different dimensions can be combined to extract the optimal size information. Wang et al. (2019) [12] believed that reducing the size of convolution kernels can improve the operating speed of the neural networks.
In the study, the CNN-based algorithm and Active Appearance Model-(AAM-) based algorithm were compared for the segmentation results of MRI images in 30 PB children.
e objective of the study was to explore the features of MRI images for PB patients, expected to provide reliable reference for its clinical diagnosis.

Subjects and Grouping.
In the study, 30 PB patients diagnosed in xx Hospital from Jan. 2019 to Dec. 2020 were selected as the subjects, including 19 boys and 11 girls, and they all underwent the MRI examination. ey were between 6 and 12 years old, with an average age of 8.21 ± 1.25 years old. e study was approved by the Ethics Committee of xx Hospital, and all the subjects had signed an informed consent form.
Inclusion criteria were as follows: (1) children diagnosed with plastic bronchitis by hospital pathological examination; (2) children with good coordination during the examination; (3) children with complete general data and imaging data.
Exclusion criteria were as follows: (1) children with suspected clinical diagnosis of plastic bronchitis and no plastic material found in imaging examination; (2) children who had contraindications; (3) children with mental diseases; (4) children who could not normally cooperate with the examination. Figure 1 shows the structure of the CNN model. LeNet-5 is used in the study.

Structure of the CNN Model.
e input data include a matrix formed by 32 × 32 pixels. e first feature layer includes 6 feature maps. After a 5 × 5 window is used to convolve the input image, a 28 × 28 feature map can be obtained. en, downsampling is performed on the first feature layer to obtain 6 feature maps with a size of 14 × 14. C stands for the convolution layer, and S stands for the pooling layer. e C3 layer is a convolutional layer, and the size of the convolution kernel is 5 × 5, as with C1. e downsampling is performed again in S4. e C5 layer performs a convolution operation on the S4 layer, and each convolution kernel in C5 layer is convolved on the basis of the S4 feature map. en, on the basis of the C5 layer, through a fully connected network, a result of 1 × 10 is finally obtained as output. In the vector of the output result of 1 × 10, the image of the final output model is classified by ReLu function and the classification result of the network output is the corresponding position of the maximum component.

CNN-Based PB Segmentation Model.
e convolutional neural network is mainly composed of the convolutional layer, the pooling layer, the fully connected layer, and the deconvolution layer. e convolutional layer detects the local features of MRI images, which can be expressed as follows: where k l ij represents the convolution kernel, i indicates the feature map in the l − 1 layer, j represents the feature map in the l layer, B l j indicates the feature map in the convolutional layer, and B l−1 i represents the feature map adjacent to B l j . M(l − 1) indicates the number of feature maps in the l − 1 layer, • indicates the convolution, c l j is the bias, and f represents the nonlinear activation function. e activation function is to process linear nonseparable data.
e study uses ReLu as the activation function, expressed as Equation (2), and the relevant derivative function can be expressed in Equation (3).
As for an MRI image A, X is the set of all pixels in image A: X � (x 1 , x 2 , . . . , x n ) and B indicates segmentation cate- , where m � 4; the probability to output b j of pixel x i in the j channel is expressed as follows: where , and the relevant loss function is expressed as follows: Adam is a commonly used method to optimize neural networks, and it is employed in the study to optimize the CNN model. e update parameter W and the offsets S dw and V dw can be expressed as follows: where P is the circulation number of Adam, α 1 and α 2 indicate the hyperparameters to control average weighted value of the two coefficients, β is the learning efficiency, and ε is a small constant, which is set to avoid that the denominator is zero. e parameters of the convolutional neural network are shown in Table 1.

MRI Images Processed by the CNN-Based
Model. e experimental equipment is a workstation equipped with Quadro K40c graphics card, 128 G memory, Intel Xeon(R) CPU E5-2660 V3, and Ubuntu 14.04 LTS (64 bit) system. e entire training is completed under the Tensorflow framework. e testing phase can be completed either in the MATLAB 2015B platform or in the Tensorflow framework.
Occurring simultaneously than (OST) factor is used to evaluate the quality of MRI image. It is the area of the intersection to the area of the union. A higher OST indicates a more accurate target detection bounding box.
It is calculated as follows: OST � area(S) ∩ area(P) area(S) ∪ area(P) , where S represents the standard area segmented by the doctor and Prepresents the result segmented by the CNN model. e shape of the bronchus is mostly variable and irregular, so the Dice coefficient, precision, and Jaccard coefficient were used to evaluate the segmentation effects. e Dice coefficient is a measure of the overlap between the segmentation result and the gold standard area. It is calculated as follows: where S represents the standard area segmented by the doctor and e Jaccard coefficient reflects the difference between the data, and it is calculated as follows: where E is the left bronchus predicted by the CNN and F represents the right bronchus area delineated by the doctor.

MRI Scanning
Method. MR Prisma 3.0 (SIEMENS, Germany) is used to conduct the scanning of 30 patients. During the scanning, the child lied on his/her back and was instructed to breathe steadily. For the child who did not cooperate, taking chloral hydrate orally was required and the doctor should inject chlorpromazine and Phenergan intramuscularly to sedate him/her. During examination, it is a must to ensure that the child maintains one posture all the time. First, cross-section and sagittal-section scans were performed to find the tracheal plane. en, coronal-section scan was performed with the point lightly above the crotch of the trachea and bronchus as the central point. Scanning parameters were as follows: matrix is 251 × 251; layer thickness is 3.5 mm; vision is 25 × 25 cm; the flip angle is 15°; and the slice gap is 6.1 mm. After the scan, two experienced professional physicians in the imaging department performed manual segmentation and analysis of the image.

Statistical
Method. e data were processed by SPSS20.0. e measurement data were expressed as mean-± deviation (x ± s), and t-test was adopted. e count data were expressed by percentage, and χ 2 test was used. P < 0.05 was the threshold for significance. Figure 2(a) shows the OST values before and after correction. Before correction, the mean OST value of target box and real box was 0.658. After correction, it was 0.921, increasing by 0.263. Figure 2(b) shows the Euclidean distance between the central point predicted by the CNN-based model and the real central point. It was between 0 and 20 mm and mainly ranged from 4 to 10 mm. Further statistics showed that the Euclidean distance less than 10 mm accounted for 72% among those ranging from 0 to 20 mm.

Evaluation of MRI Image Quality of Children Based on
Convolutional Neural Network Processing. AAM is an image feature extraction method, which can be used to extract features from target images. Figure 3(a) shows the Dice coefficients of CNN-processed MRI images and AAMprocessed MRI images under different training cycles. As the number of training cycles increased, the Dice coefficients of both the two algorithms rose, but the Dice coefficient of CNN-based algorithm was always higher than that of AAM model. e maximum Dice coefficient of CNN-based algorithm reached 0.946, while that of AAM was 0.843. It indicated that the results predicted by the CNN model had smaller differences from the real results versus the AAM. Figure 3(b) shows the precision and Jaccard coefficient under the two algorithms compared. e precision of CNN algorithm was 0.952, while that of the AAM model was 0.934, showing no obvious difference (P > 0.05). e Jaccard coefficient of CNN algorithm was 0.894, while that of the AAM model was 0.758, and the difference was statistically significant (P < 0.05). Figures 3(c) and 3(d) show the segmentation results of the CNN model and AAM model, respectively.

General Statistics on Pediatric Patients.
Of 30 PB patients, the boys (19) outnumbered the girls (11) by a wider margin and the difference was statistically significant (P < 0.05) (Figure 4). e subjects aged from 6 to 12 years, with an average age of 8.21 ± 1.25 years. e onset time of PB was concentrated in spring and summer, and there were 13 cases in spring (43.33%), 11 cases in summer (36.67%), 4 cases in autumn (13.33%), and 2 cases in winter (6.67%). Obviously, the incidence in autumn and winter was lower than in spring (P < 0.05) ( Figure 5). e average hospital stay was 13.2 ± 7.6 days. Of the 30 patients, 17 had a hospital stay of 7-14 days (56.67%), 8 had a hospital stay of 0-6 days (26.67%), and 5 had a hospital stay of 15-21 days (16.66%) ( Figure 6). Figure 7 shows MRI images processed by CNN-based algorithm. Figure 7(a) presents chronic emphysema and pleural thickening and adhesions; Figure 7(b) presents aortic dissection; Figure 7(c) presents calcification of the pleura beside the right upper mediastinum; and Figure 7(d) presents the double cavities. e MRI images clearly display the trend of the trachea and bronchus, location, morphology, size, and density of foreign matters. It manifested that the solidified sticky secretion was tubular in shape and its density was close to or slightly higher than the density of liquid. Furthermore, there were high-density dot or strip shadows.

Discussion
PB arises from the block of bronchus, leading to abnormal ventilation of the lung tissue and even acute dyspnea in severe cases. Generally, foreign matters are in a tree-like shape [13]. PB may occur at any stages, but children are predominantly affected, especially those aged 4-12 years [14]. e PB subjects in the study were from 6 to 12 years old, with an average age of 8.21 ± 1.25 years. Of 30 PB patients, boys outnumbered girls and the difference was statistically significant (P < 0.05), which was consistent with the research results of Xiong et al. (2019) [15]. e outset time was concentrated in spring and summer, suggesting that the onset of PB is seasonal. Before correction, the average OST value of target box and real box was 0.658. After correction, it was 0.921. e Euclidean distance from the predicted central point to the real central point ranged from 0 to 20 mm and the Euclidean distance less than 10 mm accounted for 72.5%. It indicated that the CNN-based algorithm can accurately segment MRI images for PB. e Dice coefficient has been widely applied in the verification of the segmentation effects [16]. Under different number of training cycles, the maximum Dice coefficient of CNN-based algorithm reached 0.946, while that of AAM algorithm was 0.843. It indicated that the algorithm of the study had better segmentation effects, which was in line with the research results of Fujioka et al. (2021) [17].
Clinically, PB mainly manifests as fever and cough [18,19]. A majority of subjects included had fever. e average hospital stay was 13.2 ± 7.6 days. Of the 30 patients, 17 had a hospital stay of 7-14 days (56.67%), 8 had a hospital stay of 0-6 days (26.67%), and 5 had a hospital stay of 15-21 days (16.66%). e MRI images for PB presented chronic emphysema, pleural thickening and adhesions, aortic dissection, calcification of the pleura beside the right upper mediastinum, and double-cavity signs. MRI can clearly display the trend of the trachea and bronchus, location, morphology, size, and density of foreign matters. It manifested that the solidified sticky secretion was tubular in shape and its density was close to or slightly higher than the density of liquid. Furthermore, there were high-density dot or strip shadows.

Conclusion
In this study, a CNN-based algorithm was constructed and applied to segment MRI images for 30 PB patients, to explore the imaging features of PB patients. e algorithm based on CNN can significantly improve the segmentation accuracy of MRI images for plastic bronchitis in children. MRI images prompted that pleural effusion was a dangerous factor in the occurrence and development of PB. However, the shortcomings of this study are that the sample size is small, the MRI examination time is long, and the cost is high, which is limited in the actual operation. In short, the CNN-based MRI image segmentation model greatly raises the segmentation accuracy of the bronchus, providing a theoretical basis for the diagnosis of PB.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest.