Abstract

This study was aimed to explore the adoption value of magnetic resonance imaging (MRI) under convolutional neural networks (CNN) in the therapeutic effect of nasopharyngeal carcinoma (NPC) radiotherapy. A total of 54 NPC patients were recruited. CNN was employed to perform 3D visualization processing on magnetic resonance (MR) images of NPC patients. MRI changes were analyzed before and after the patient received radiotherapy. The image segmentation and radiotherapy effects of CNN were evaluated by the Recall, intersection over union (IOU), postoperative apparent diffusion coefficient (ADC), and diagnostic coincidence rate. Moreover, gradient vector flow (GVF) algorithm, fuzzy c-means (FCM), and SegNet were adopted for comparative evaluation. Recall of CNN was 94.89% and the IOU was 84.16%, which was remarkably different from other algorithms (P < 0.05). After analysis of the MRI images of patients receiving radiotherapy, ADC of local residual patients was 1.108 ± 0.097 measured by CNN, the ADC was 1.826 ± 0.115, and the missed diagnosis rate was only 7.14%. In summary, CNN had a good effect on the localization and segmentation of NPC patients, and can accurately evaluate the effect of patients receiving radiotherapy, which can assist clinical diagnosis and treatment of NPC.

1. Introduction

Nasopharyngeal carcinoma (NPC) is a malignant tumor in the mucosa of the nasopharynx. Its occurrence is obviously regional [1, 2]. Most NPC patients are often accompanied by neck mass. In addition, they are prone to hearing loss, tinnitus, earache, and recurrent ear infection [3]. The incidence of NPC is the first among head and neck malignant tumors. With the gradual development of radiotherapy technology, NPC has been controlled to a certain extent. Patients who receive radiation therapy typically have their lesions completely gone within three months. If the tumor does not disappear six months after surgery, it is residual tumor, and NPC residual or recurrence is the main reason that affects the prognosis and quality of life of patients. At present, the screening of NPC mainly includes imaging examination and detection of NPC virus infection in blood, including computerized tomography (CT) and magnetic resonance imaging (MRI) [4]. Its treatments are classified into local treatment (radiotherapy and surgery) and systemic treatment (chemotherapy, targeted therapy, immunotherapy, etc.) [5]. NPC is very sensitive to radiotherapy, and radiotherapy is the main treatment for NPC. In the early stage, only radiotherapy is utilized, and radiotherapy is mainly combined treatment in the local late stage [6]. With the popularization of intensity-modulated radiation therapy technology and the adoption of comprehensive treatment, the NPC local control rate is more than 90%, and the 5-year survival rate is more than 80%. In contrast to other tumors, NPC treatment has objective efficacy [7, 8].

At present, CT and MRI are commonly applied for the clinical diagnosis of NPC [9]. CT technology gives a lower resolution of soft tissue images, whereas MRI can clearly show the tumor site, invasion range, and lymph node metastasis [10]. In contrast to CT, MRI can remarkably increase the detection rate of soft tissue invasion around nasopharynx by about 30% and the detection rate of bone invasion at skull base by 20%. In addition, the detection rates of skull base invasion and parapharyngeal space invasion were remarkably improved [11, 12]. With breakthrough of medical imaging, the revolution of CT and MRI are enhanced, which provide good help for doctors in diagnosis and treatment [13]. CT has high density resolution, which can clearly show the patient’s anatomical relationship, and the clinical examination is quick. However, because of its large radiation amount and only cross-sectional scan, its use is limited. MRI is utilized to evaluate the efficacy of NPC radiotherapy due to its omni-directional and multi-parameter characteristics, which can accurately reflect the nerves and blood vessels at the lesion site. The quantification technology of image processing is actually very critical, which plays an important role in disease diagnosis and early prevention [14]. Now many quantitative analyses are a segmentation problem, and effective and automatic organ region segmentation is very challenging. The brightness difference between different parts is not obvious, connected with the surrounding parts, and the boundary is fuzzy. In addition, there are noise and other interference factors [15, 16]. Deep learning convolutional neural networks (CNN), cyclic neural networks, generative adversarial networks, etc. can focus on medical image imaging problems. In recent years, computer-aided medical image analysis has made major breakthroughs in technology and has obvious advantages in improving the efficiency and ability of medical services [17], which has also become an effective way to solve problems such as insufficient medical resources.

Therefore, it was attempted to build a CNN model, train the network to segment the MRI image of NPC, and conduct three-dimensional visualization of the MRI image, to make the image stereoscopic display. After training, CNN was utilized to segment NPC MRI images, to comprehensively study the adoption value of MRI based on CNN deep learning algorithm in the evaluation of therapeutic effect of diagnostic NPC radiotherapy. Moreover, further reference and basis were offered for the adoption of CNN in the diagnosis and curative effect evaluation of NPC.

2. Methods

2.1. Samples

In this study, 42 patients with NPC who visited hospital from March 2019 to March 2021 were recruited. There were 23 male patients and 19 female patients in total, with the average age of 55.2 ± 11.3 years old. All patients were diagnosed pathologically and graded into T1N0M0 of 7 cases, T1-2N1M0 of 7 cases, T2N0M0 of 5 cases, T1-3N2M0 of 7 cases, T3N0-2M0 of 5 cases, T3N1M0 of 4 cases, T4N2M0 of 4 cases, and T2N2M0 of 3 cases according to the Tumor Node Metastasis (TNM) staging established by American Joint Committee on Cancer. The details of TNM standard are shown in Table 1. All patients underwent radiotherapy, followed by MRI examination before and after radiotherapy. The target area and normal organs were delineated based on MRI images, and radiotherapy was implemented based on the delineation results. All patients had signed the informed consent forms and this study had been approved by the ethics committee of hospital.

Inclusion criteria: (i) those confirmed as NPC by pathological biopsy before radiotherapy; (ii) patients aged between 18 and 70 years; (iii) the patient had no previous tumor history; (iv) patients who were able to participate in the study and follow-up; (v) no serious metabolic diseases such as heart, brain, and kidney.

Exclusion criteria: (i) biopsy pathology was not clear; (ii) the patient had contraindications to radiotherapy and chemotherapy and MRI.

2.2. MRI Examination

Magnetic Resonance Imaging System was employed to perform MRI scans on 42 NPC patients. The patient was in a supine position, and the cross-sectional, sagittal, and coronal planes of the patient were scanned using standard head coils and spin echo sequences. The T1 weighted imaging (T1WI) scan sequence was utilized for the transverse and coronal planes, and the T1WI and T2 weighted imaging (T2WI) scan sequences were utilized for the sagittal plane. Time of inversion (TI): 110 ms, time of echo (TE): 48 ms, time of repetition (TR): 3,800∼5,100 ms, layer spacing: 6 mm, layer thickness: 5 mm, and flip angle: 90°.

2.3. Radiotherapy

Radiation therapy was performed on all NPC patients participating in this study. Before treatment, the clinical target volume (CTV) of the patients was delineated according to MRI. High-risk CTV1 included the parapharyngeal space, slope, rupture hole, and lymphatic drainage area where the positive lymph nodes in the neck were located. Low-risk CTV2 included low-risk lymphatic drainage areas where subclinical metastasis may occur. After the delineation was over, each target area was expanded back by 2–3 mm, and expanded 3–4 mm outside of other directions to be the planned target volume. During radiotherapy, 6MV-X linear accelerator was utilized for external irradiation, the irradiated skin area was 100 mm × 125 mm∼130 mm × 160 mm, the dose was 6,000∼7,000 Gy, and the course of treatment was 45∼65 days.

2.4. Construction of NPC Radiotherapy Curative Effect Model under CNN

NPC radiotherapy efficacy evaluation including data preprocessing and CNN model was completed. Data need to be normalized as the dimension between the features in the actual prediction is different. Normalization of maximum and minimum values is adopted to standardize the data.

In equation (1), x is data after normalization of maximum and minimum, is pre-processed data, is the maximum of dimensional data, and is the minimum. It is necessary to reduce the dimension of data to speed up the training. The data redundancy information is presented by principal component analysis. If matrix H contains a samples, and b is row vector of different attributes of each sample, then the matrix H is expressed as follows:

Column vector of H is P, then . The covariance is obtained regarding different P, and a matrix G is formed.

and are expectations of and , respectively. Matrix element values reflect linear correlation degree between two random variables. Two random variables are linearly independent with covariance of 0. Large covariance reflects strong correlation. H is transformed to harvest transformation matrix of H to calculate eigenvalues and eigenvectors of H. Large feature value means large numerical range and variance of it in data compression. The transformation matrix is constructed by using the eigenvectors corresponding to large eigenvalues to obtain new features of high importance.

After data treatment, CNN model is utilized to evaluate the curative effect of NPC radiotherapy. The convolution layer includes convolution, activation, and pooling. After the original data is input, feature mapping is performed on the data obtained after the original data is convolved with the convolution kernel (CK). The calculation of neural network on 2D tensors is shown as below.

is output 2D tensor, is input 2D tensor, and is CK’s weight. Each layer in a complete CNN contains several feature images. If each layer contains a CK, and each CK is composed of massive parameters, the convolution is implemented for each layer as follows.

In equation (5), feature image of s-th layer convolution is , and CK parameter of i-th feature image of convolution layer mapped to the feature image of s-th layer convolution is . i-th feature image of CL is , and offset of feature image of s-th CL is . Convolution refers to linear transformation. Usually, multi-layer nonlinear feature is formed with a nonlinear activation function after convolution of all layers. Taking the sigmoid function as an example, the expression is as follows:

Without performing weighting or requiring activation functions, pooling performs maximum or average operations to express statistical characteristics of the pooling results. It is only involved in network institution design to reduce data dimensions, speed up calculations, and avoid over-fitting (Figure 1).

Through the full connection layer, image is transmitted to the output layer after original features extraction through the CL and main features from the sub-sampling layer. The full connection layer maps the information processed by CL and sub-sampling layer to the label space of the sample. 2D feature graph is converted into 1D feature vector, and the information is summarized for full-neural network training (Figure 2). Different from CL, the fully connected layer captures the nonlinear relationship between CL’s comprehensive feature information and sub-sampling layer to learn feature information, to accurately classify features.

In the process of calculating CNN, in order to prevent the model from over-fitting due to the excessive feature information contained in the full connection layer, it is usually necessary to introduce the normalized function to carry out the final classification. The final classification layer utilized is Softmax function in this work. The classification function is as follows:

In equation (7), probability of a sample in a category is , output of the neuron is u, and number of categories is K.

2.5. Evaluation Indexes

When NPC lesions were detected, the results obtained by subjective judgment alone were unreliable. Therefore, segmentation results of CNN were objectively evaluated regarding recall and intersection over union (IOU). GVF, FCM, and SegNet were adopted for comparison. The function expressions of Recall and IOU are as follows:

In equations (8) and (9), TP represents the pixels of the positive lesions in the area delineated by the algorithm, FN represents the pixels of the positive lesions outside the area delineated by the algorithm, S is the area delineated by the algorithm, and S0 is the area delineated by the expert [1820].

2.6. Statistical Treatment

SPSS 22.0 was employed. Count data were expressed in the form of %, quantitative data are represented by (), tested by χ2 test. t-test was utilized between independent samples. P < 0.05 meant considerable differences.

3. Results

3.1. MRI Image Preprocessing of NPC Patients

Due to the limitation of MRI imaging, as well as the inhomogeneity of the magnetic field or the movement of the patient during the imaging process, it is easy to cause the image to produce motion artifacts, noise, and blurred edges. Therefore, the image was pre-processed for image enhancement to highlight image features and prepare for the next step of lesion segmentation. The preprocessing result is shown in Figure 3. The patient’s nasopharyngeal cavity was narrowed, and the nasopharyngeal mucosa on the right side was swollen and distributed asymmetrically. Due to the noise artifacts in the original image, the image was smoothed and gray-scale adjusted.

3.2. NPC Lesion Detection Model Based on CNN

The CNN was utilized to locate the NPC lesion area, then this area was delineated and segmented, and the segmented image area was reduced at the same time. When the lesion area was delineated, manual delineation required the doctor to observe the scope of the lesion carefully and comprehensively on the MRI image and delineate the tumor target area layer by layer. Therefore, the results obtained only by the subjective judgment of clinicians were often unstable. When CNN was adopted for detection, the detection results were adjusted to make the results stable. In Figure 4, the blue area was the NPC lesion area outlined by the radiologist, and the red area was the area automatically segmented by CNN. Notably, the lesion area marked by CNN was relatively more accurate, which was similar to the expert’s outline of the lesion.

In addition to subjective evaluation, segmentation effect of CNN was objectively evaluated regarding Recall and IOU. The results are shown in Figure 5. In contrast to other algorithms, it was found that the Recall and IOU of CNN were 94.89% and 84.16%, respectively. The Recall and IOU of the GVF were 67.48% and 52.15%, respectively. Those of the FCM were 84.12% and 68.71%, respectively. The Recall and IOU of the SegNet were 78.44% and 65.01%, respectively. After the image was pre-processed in the initial stage of model training, the Recall and IOU obtained by CNN segmentation were remarkably higher than other algorithms (P < 0.05).

3.3. Evaluation of the Curative Effect of NPC Radiotherapy Based on CNN

After the radiotherapy, the radiotherapy effect was evaluated by brain MRI images of the patients. Among the 54 NPC patients, 47 cases had no nasopharyngeal local residual, and 7 cases had nasopharyngeal local residual. After manual outlining and judgment by radiologists, there were 38 cases of nasopharyngeal cavity deformation, 44 cases of nasopharyngeal wall thickening (3 cases with parapharyngeal soft tissue thickening), 4 cases of sphenoid wing plate destruction, and 2 cases of slope bone destruction. The thickening of nasopharyngeal wall was shown as irregular thickening of nasopharyngeal soft tissue on MRI, and the thickening of parapharyngeal tissue showed slightly high signal (Figure 6). Sphenoid wing plate destruction and slope bone destruction showed partial bone destruction with uneven high intensity on MRI.

After radiotherapy, the ADC of NPC patients was evaluated, and the ADC measured based on different algorithm outlines was compared with the ADC measured by radiologists. The results showed that the ADC measured by CNN for the local residual patients after radiotherapy was 1.108 ± 0.097. For locally residual patients, the ADC measured by CNN was 1.826 ± 0.115 (Figure 7). The measurement result calculated by CNN was relatively closer to the expert manual measurement result, which was statistically different compared with other algorithms (P < 0.05).

3.4. Contrast of Coincidence Rate of Postoperative Diagnosis Based on CNN

After 54 NPC patients were treated with radiotherapy, the MRI images of the patients were analyzed for diagnosis. The results showed that 50 patients (92.59%) were accurately diagnosed by CNN for the first time, and 4 patients (7.41%) were missed for the first time (Figure 8). In contrast to other algorithms, the diagnostic accuracy rate was remarkably different (P < 0.05). The four cases of missed diagnosis were all patients with skull base invasion, and MRI showed abnormal skull base signals, leading to missed diagnosis.

4. Discussion

The incidence of NPC is the first among head and neck malignant tumors. With the gradual development of radiotherapy technology, the local control rate of NPC has been improved to a certain extent. Patients who receive radiation therapy typically have their lesions completely gone within three months. If the tumor does not disappear six months after surgery, it is residual tumor, and NPC residual or recurrence is the main reason that affects the prognosis and quality of life of patients [21]. Clinical imaging examination is still needed for local residual or recurrence of MPC. Currently, CT and MRI are two commonly utilized examination methods. Compared with CT examination, MRI has the omni-directional and multi-parameter characteristics, which can accurately reflect the nerves and blood vessels of the lesion. Thus, MRI can be applied in the efficacy evaluation of NPC radiotherapy [22]. Many quantitative analyses are a segmentation problem. The brightness difference between different parts is not obvious, connected with the surrounding parts, and the boundary is fuzzy. CNN, cyclic neural networks, and generative adversarial networks focus on medical image imaging problems [23]. In recent years, computer-aided medical image analysis has made major breakthroughs in technology and has obvious advantages in improving the efficiency and ability of medical services, which has also become an effective way to solve insufficient medical resources [24].

NPC MRI images were pre-processed, and then CNN was utilized to accurately delineate and segment the lesions in the MRI images of NPC patients. Then, the proposed CNN was compared with the GVF, FCM, and SegNets, and it was found that the Recall and IOU of CNN were 94.89% and 84.16%, respectively. The similarity with the expert’s delineation of the lesions was high, suggesting that CNN had ideal segmentation effect on NPC MRI images. Moreover, the MRI images of NPC patients after chemotherapy were analyzed to evaluate the ADC of the patient. ADC is an important way to predict the sensitivity of the tumor to radiotherapy and chemotherapy and whether the tumor remains in recurrence. The higher the ADC, the better the anti-tumor treatment effect. It was found that for patients with local residual tumors and local residual tumors after radiotherapy, the ADCs measured by CNN were 1.108 ± 0.097 and 1.826 ± 0.115, respectively. These data were in line with the manual measurement results of experts. The results were in line with the study of Peng et al. [25]. Diagnosis of patients after surgery suggested that the diagnostic accuracy of CNN was 92.59%, which was remarkably different from the other three algorithms (P < 0.05), suggesting ideal diagnostic effect of CNN on radiotherapy efficacy.

5. Conclusion

CNN was utilized to outline and segment the lesion sites on NPC MRI images, and the NPC radiotherapy efficacy was judged by MRI images. It was proved that CNN could identify the lesion sites on MRI images of NPC patients and accurately evaluate the radiotherapy effect of patients. The algorithm presented in this study can assist clinicians, which is suitable to be applied in imaging diagnosis of NPC patients. It can be concluded that MRI based on CNN yielded high clinical adoption value for evaluation of the NPC radiotherapy effect on patients. However, there are still some deficiencies in this study, which are manifested as a small number of included cases and that the truth all patients within six months of radiotherapy were just confirmed by biopsy pathology without long-term follow-up, which needs to be further studied by increasing the sample size in future studies.

Data Availability

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

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Haifeng Qian and Yancheng Fang contributed equally to this work.