Imaging Algorithm-Based Real-Time Shear Wave Elastography Combined with Thyroglobulin Antibodies in Diagnosis of Thyroid Cancer

Department of Ultrasonography, e First Affiliate Hospital of Qiqihar Medical College, Qiqihar 161041, Heilongjiang, China Department of Price Division, e First Affiliate Hospital of Qiqihar Medical College, Qiqihar 161041, Heilongjiang, China Department of Nuclear Medicine, e First Affiliate Hospital of Qiqihar Medical College, Qiqihar 161041, Heilongjiang, China Department of Nuclear Medicine Department, e First Affiliate Hospital of Qiqihar Medical College, Qiqihar 161041, Heilongjiang, China


Introduction
Differentiated thyroid carcinoma (DTC) is one of the common endocrine malignancies. In recent years, its incidence has shown an obvious upward trend [1,2]. Clinically, the diagnosis of DTC has a high rate of misdiagnosis due to the largest diameter of the tumor, the distribution of special subtypes, atypical clinical symptoms, and small lesions. Although the progress of papillary thyroid carcinoma (PTC) is slow and the biological behavior tends to be benign, it can spread in the glands and metastasize to regional lymph nodes, leading to further lesions, so preoperative diagnosis is very important. At present, advanced and effective imaging techniques are used for inspection to improve the diagnosis rate of DTC [3]. In recent years, shear wave elastography (SWE), which is implemented based on conventional ultrasound examination, has gradually become the preferred method for thyroid disease. It can detect tiny nodules about 1 to 2 mm and clearly shows the gravel-like calcification in the tiny nodules and the changes in the nature of the surrounding lymph nodes [4]. Due to the limitation of image clarity in conventional ultrasound images, the imaging effect of thyroid nodules is not good, which is not easy to distinguish from other normal tissues, resulting in missed diagnosis. erefore, improving the corresponding quality indicators of ultrasound images is of great significance for the early diagnosis of patients with thyroid cancer.
As a new technology in recent years, SWE uses ultrahighspeed imaging technology to generate various forms of sound waves through a multiwave imaging platform and a probe chip. It adopts acoustic radiation pulse control technology, continuous focusing at various depths of the tissue, and forming a "Mach cone" phenomenon to promote the generation of shear waves [5]. Since SWE has a quality mode to control the quality of the elastic image during the inspection process, it has good objectivity and strong repeatability. It can display the diseased tissue in time and does not require artificial pressure on the patient's tissue, which can reduce the impact on the patient. erefore, it has become an important auxiliary diagnosis for the differential diagnosis of benign and malignant thyroid nodules [6]. erefore, it has become an important auxiliary diagnosis for the differential diagnosis of benign and malignant thyroid nodules. Real-time SWE can visually display the color coding of elastic images under the real-time guidance of twodimensional images and directly measures the absolute value of Young's modulus of thyroid nodules, which is expressed quantitatively in the form of Young's value (kPa). SWE uses color images to intuitively reflect the elastic distribution of the organization in real time. SWE can measure many Young's moduli, including the maximum value (E max ), average value (E mean ), minimum value (E min ), and standard deviation (E sd ) of the measured structure. It can comprehensively and accurately quantitatively analyze the hardness of the lesion tissue and the changes with the surrounding tissues. Moreover, it can well describe the morphological characteristics of the lesion and can clarify the state of the lesion, reflecting the characteristics of the blood supply of the tumor and the relationship with the surrounding tissue structure. It can also objectively display the size and structure of the tumor. Unfortunately, due to the noise and artifacts in SWE imaging, there is a certain degree of missed diagnosis rate in clinical diagnosis, and other methods are needed for auxiliary diagnosis [2]. yroglobulin antibodies (TGAb) are an emerging diagnostic marker for thyroid cancer tumors which have emerged in recent years. TGAb is the corresponding antibody of thyroglobulin (Tg). As one of the main serum tumor markers in DTC patients, Tg can monitor the tumor recurrence and metastasis after treatment. High titer TGAb will affect the accuracy of Tg determination and reduces the sensitivity of Tg as a tumor marker [7][8][9]. Some studies suggested that TGAb alone can also be used as a potential predictor of DTC, which has become a DTC substitute tumor marker. TGAb is used to assist physicians in the effective diagnosis of thyroid cancer to improve the accuracy of DTC diagnosis [10,11]. erefore, combination of realtime SWE and TGAb under the imaging algorithm was used to explore the rapid clinical diagnosis of thyroid cancer.
In summary, a hybrid displacement estimation algorithm with weighted phase separation and two-dimensional cross-correlation was proposed in the traditional imaging algorithm design, which was applied to the image analysis of 102 patients with DTC in real-time SWE. e optimization effect of the algorithm was evaluated through indicators such as sensitivity, specificity, and AUC. Receiver operating characteristic (ROC) curve, Young's value E max , etc. were used to evaluate the diagnosis performance for thyroid cancer, so as to comprehensively evaluate the application value of real-time SWE under the imaging algorithm combined with TGAb in the diagnosis of DTC.

Research
Objects. 102 patients with DTC who were pathologically confirmed by pathology in the hospital were selected as the research objects and were divided into a test group (TGAb-positive) and a control group (TGAb-negative). Among them, 59 were males and 45 were females, aged 29-65 years. e patients were randomly divided into test group (TGAb positive) and control group (TGAb negative), with 51 people in each group. In the experimental group (TGAb-positive), TGAb ≥ 115 IU/mL, there were a total of 56 patients, TGAb level was 498.0 (356.5-998.6) IU/mL, and average age was 41.74 ± 11.87 years, of which 15 were males and 36 were females. In the control group (TGAb-negative), TGAb < 115 IU/mL, TGAb level was 9.78 (8.79-12.96) IU/mL, and average age was 39.67 ± 9.88 years, including 21 males and 30 females. e study had been approved by the ethics committee of the hospital, and the patients and their families included in the study knew the study process and signed the informed consent.
Inclusion criteria were as follows: (i) patients without history of thyroid surgery before the experiment, (ii) patients with complete basic clinical data, and (iii) patients with preoperative serological examination results. Exclusion criteria were as follows: (i) patients whose surgery was performed in a different hospital or whose surgical method was partial thyroidectomy or lymph node dissection, (ii) patients with incomplete pathological data, and (iii) patients who had not been tested for serological indicators or had incomplete inspection indicators.

SWE Image Acquisition.
e Mindray Resona 7 ultrasonic color Doppler diagnostic apparatus equipped with SWE software was used, the probe model was L11-3U, and the frequency was 5.6, 10.0 MHz. e patient took a supine position and breathed calmly. e thyroid imaging mode was set to perform routine ultrasound examination of thyroid nodules, such as location, size, number, blood flow, cystic change, and calcification. e largest longitudinal section of the thyroid nodule was taken. After waiting for 3-5 seconds, the image was frozen after the image was stable. e quality mode of SWE was observed, and it was switched to the measurement mode after the credibility index was as close as possible to 100% and the lesion area was uniformly green. If the reliability index was low and there was a lot of yellow in the lesion area, the measurement was not accurate and needed to be measured again. For the area of interest, the measurement method that comes with the ultrasonic instrument was adopted for further measurement. During the measurement, the liquid anechoic area and the large calcification should be avoided as much as possible, and the lesion should be covered as much as possible. ree repeated operations were performed for each lesion. e measurement data of the measurement lesion should include the maximum value (E max ), average value (E mean ), and minimum value (E min ) of the elastic modulus value. Statistics of routine indicators such as tumor size and number of lesions of patients were conducted, and all operations were performed by physicians who had received SWE examination technology training.

Real-Time Shear Wave Model Construction Based on
Imaging Algorithm. Imaging algorithm was employed to optimize real-time SWE, which was of great application value in clinical tumor diagnosis, and can greatly reduce the degree of dependence on the subjective judgment of physicians. e amplitude modulation correction (AMC) can reduce the amplitude modulation noise in the traditional calculation algorithm structure, which also estimated the displacement through the fluctuation of the signal amplitude. Its computing efficiency was very high, and its performance was superior to stretching technology. e strain estimation based on the amplitude modulation theory was expressed by the phase difference as follows: In equation (1), a 2 and a 1 represented the displacement of the continuous window, and N 2 and N 1 represented the relative position of the displacement.
First, the signal model was set as follows: In equations (2) and (3), b 1 (q) and b 2 (q) were the analytical signals of presignal and postsignal, respectively, which were obtained through Hilbert transform or IQ modulation.
en, the AMC technique was used in the phase estimation, which was expressed as follows: . (4) In equation (4), q was the axial sampling point, and q � 0 was the surface of the transducer. Δq was the sampling frequency, Q was the window length, and d(q) was the real displacement at q. e position of h m was related to the weight, such as amplitude modulation. e weight can not only be the amplitude, but all signal attributes can be used as the weight. d(q) was expressed as In equation (5), z was the tissue strain, and the position estimate τ m at h m was defined as follows: Displacement estimation based on phase information was as follows.
e phase method processed the analytic signal and obtained the analytic signal through the Hilbert transform. e complex cross correlation was expressed as follows: In equation (7), k 1 and k 2 were the analytic signals before and after compression, k * 1 was the conjugate complex number of k 1 , mΔq was the start of the signal of the mth window, Q was the window length, and X was the estimated displacement value. e phase expression was as follows: Weight phase separation (WPS) was adopted for precise estimation of axial displacement. WPS was set as the matching criterion to obtain the subsampling point-level axial displacement accuracy. e following equation was used for iteration: In equation (9), d x and d x+1 were the displacements, φ pre and φ post were the phase of the RF signal before and after compression, and u 0 was the center frequency. U(i, j) was the weight, and its expression was as follows: In equation (10), q 1 and q 2 were the amplitude signals corresponding to presignal and postsignal. In this way, the algorithm was divided into two levels. e first-level displacement estimation was rough estimation, which was accurate to the sampling point, and the second displacement estimation was the net estimation, which improved the axial resolution. Moreover, the displacement data was smoothed. In the strain estimation, to avoid the gradient calculation to further amplify the displacement estimation noise, the leastsquares strain estimation method was adopted.

Real-Time Shear
Wave Elastic Imaging Feature Acquisition and Label Establishment. After processing of the realtime SWE image with a hybrid displacement estimation algorithm based on weighted phase separation and two-Scientific Programming dimensional cross-correlation, the specific process and effect changes are shown in Figures 1 and 2. e collected SWE images of 102 patients were optimized by algorithm step by step. Before and after the algorithm processing during the collection optimization period, the processing algorithms included one-dimensional elastic algorithm, two-dimensional elastic algorithm, and two-dimensional cross-correlation hybrid algorithm.
e SWE image quality processed by different algorithms was compared and analyzed, and the SWE image optimization effects of each algorithm were summarized.

e Tumor Pathological Characteristics and Preoperative Serological Index Detection of the Two Groups of Patients.
Preoperative serological indexes (thyroid-stimulating hormone (TSH) and thyroid peroxidase antibody (TPOAb)) and pathological characteristics of tumors were compared between the test group and the control group. e following pathological features of the surgically resected specimen were evaluated by at least two experienced pathologists: the size of the tumor (the maximum diameter of the largest tumor was measured) and whether there was peripheral invasion, bilateral, multifocal (≥2 tumor lesions in one thyroid lobe), lymph node metastasis rate (number of lymph node metastases/number of lymph nodes surgically dissected), lymph node staging (N0: no regional lymph node metastasis; N1A: lymph node metastasis in VI group, namely, pretracheal, paratracheal, and prelaryngeal lymph nodes; N1B: unilateral, bilateral, or contralateral cervical or superior mediastinal lymph node metastases).

Statistical Methods.
e data processing of this study was performed by SPSS 19.0. Counting data was described by percentage (%), and measurement data was represented by mean ± standard deviation (x ± s). e Mann-Whitney U test method was used to compare the differences in the macroscopic characteristics and imaging omics characteristics of the two groups of patients. e ROC was used to calculate the area under the curve (AUC) and sensitivity, specificity, and accuracy of each feature. P < 0.05 indicated that the difference was statistically considerable. e t-test was used to compare the difference between the two groups for continuous variables that met the normal distribution and the homogeneity of variance. Continuous variables that did not meet the normal distribution (preoperative TSH, TGAb, TPOAb, tumor size, and lymph node metastasis rate) were tested by the rank-sum test (Kruskal-Wallis). Categorical variables (tumor size and multifocality) were tested by chi-square test to compare their constituent ratios between the two groups. P < 0.05 indicated that the difference was considerable. Figure 3 shows a comparison diagram of SWE effects under different algorithm processing. Figures 3(a)-3(d) show SWE original image, one-dimensional elastic algorithm processed image, two-dimensional elastic algorithm processed image, and two-dimensional cross-correlation hybrid algorithm processed image, respectively. e image quality of the SWE image processed by the algorithm was improved, and the image improvement effect of the SWE image processed by the two-dimensional cross-correlation hybrid algorithm was the most obvious, followed by the two-dimensional elastic algorithm and the one-dimensional elastic algorithm. Figure 4 shows the comparison of signal-to-noise ratios of SWE strain diagrams processed by different algorithms. Compared with the one-dimensional elastic algorithm and the two-dimensional elastic algorithm, the two-dimensional cross-correlation hybrid algorithm had better signal-tonoise optimization effect. e signal-to-noise ratio curves of the two-dimensional elastic algorithm and the two-dimensional cross-correlation hybrid algorithm crossed, but the overall SNR optimization effect of the two-dimensional cross-correlation hybrid algorithm was better. Figure 5 showed the comparison of operation time of different algorithms.

Quality Evaluation of Algorithm Optimization.
e calculation time of the two-dimensional cross-correlation hybrid algorithm was longer than that of the two-dimensional elastic algorithm in the early stage, and the calculation time was shorter than that of the two-dimensional elastic algorithm in the later stage. e computing times of the two-dimensional cross-correlation hybrid algorithm and the two-dimensional elastic algorithm were significantly shorter than that of the one-dimensional elastic algorithm.

Analysis of Preoperative Serum Indexes of Two Groups of Patients.
rough the comparisons of the levels of TSH, TGAb, and TPOAb in the preoperative serum indexes of the two groups of patients, it was found that the preoperative TSH, TGAb, and TPOAb of TGAb-positive patients were all greater than those of the TGAb-negative group (P < 0.05), as illustrated in Figure 6.

Comparison of Tumor Pathological
Characteristics between the Two Groups. It was found that there were remarkable differences in the tumor size and the number of lesions between the two groups after the tumor pathological characteristics of the two groups of patients were compared. e preoperative tumors of experimental group were multifocal, and the tumor size was larger than that of control group (P < 0.05), with considerable differences (Figures 7  and 8).

Comparison of E max Value between Two Groups of Patients and Analysis of the Correlation between E max Value and TGAb
Level. Young's modulus E max , E mean value, and E min value of the two groups of patients were measured and compared, and the results showed that maximum Young's modulus E max of experimental group was superior to that of control group (P < 0.05) (Figure 9). Moreover, the patient's TGAb level and the corresponding data of maximum Young's modulus E max were analyzed, and it was found that there was a very significant positive correlation between the two (P < 0.001), which indicated that the level of TGAb in the patient's body greatly affected the value of E max , as presented in Figure 10.

SWE Image Quality Analysis Based on Imaging Algorithm.
Patient's SWE image optimized by the elastography algorithm was analyzed, and the optimization effect of the algorithm was evaluated regarding the detection sensitivity, specificity, and AUC. Figures 11 and 12 show that the sensitivity, specificity, and AUC of the joint detection of patient real-time SWE E max and TGAb for the malignant degree of thyroid cancer were dramatically greater than those of a single real-time SWE and TGAb, with statistical significance (P < 0.05).

Discussion
yroid cancer is a malignant tumor with a very high incidence. Due to the complex anatomical structure of the human thyroid gland, which is closely connected with multiple organs, and the fact that thyroid cancer is highly erosive, it is easy to invade the surrounding tissue structure   and lymph node metastasis, making the implementation of surgical resection very difficult [12,13]. As an emerging technology, real-time SWE is widely used in assisted medical imaging diagnosis and analysis of various cancers [14,15]. At present, the preoperative diagnosis of thyroid cancer using SWE is relatively common [16,17]. However, the traditional method of using the physician's subjective judgment will easily lead to misdiagnosed cases. erefore, how to accurately locate the tumor lesions and lymph nodes still needs further research.
In this work, the real-time SWE based on the hybrid displacement estimation algorithm was adopted to automate the image segmentation to improve the image processing effect, by which the macroscopic and microscopic imaging omics characteristics were obtained. TGAb detection was used as an auxiliary diagnosis method for DTC [18]. ROC curve was used to analyze the sensitivity, specificity, and AUC of joint detection of real-time SWE combined with TGAb to detect the malignant degree of thyroid cancer. It was found that, in 102 patients with thyroid cancer, the   Scientific Programming preoperative TSH, TGAb, and TPOAb of experimental group were greater than those of the control group (P < 0.05). e preoperative tumors of experimental group were multifocal, and the tumor size was larger than that of TGAb-positive patients (P < 0.05). Such results were similar to the results of Patel et al. [19], indicating that TGAb levels directly affected the tumor size. E max of Young's modulus of experimental group was greater than that of the control group (P < 0.05). ere was a very significant positive correlation between the TGAb level and E max of Young's modulus (P < 0.001). Moreover, the sensitivity, specificity, and AUC of the joint detection of real-time SWE E max and TGAb for the malignant degree of thyroid cancer were significantly greater than those of the single real-time SWE and TGAb, and the difference was substantial (P < 0.05). e joint detection of SWE E max and TGAb had high prediction accuracy and AUC and can be used as an improved model for predicting thyroid cancer. In short, joint detection of real-time SWE based on hybrid displacement estimation algorithm combined with TGAb had high sensitivity, specificity, and AUC for the diagnosis of DTC, which was suitable for clinical application.

Conclusion
In this study, real-time SWE was optimized based on mathematical imaging algorithms, combined with TGAb, and applied to the clinical imaging diagnosis and analysis of 102 cases of DTC. It was found that the joint detection of real-time SWE based on the optimization of the imaging algorithm combined with TGAb showed high sensitivity, specificity, and AUC for the diagnosis and prediction of DTC. However, the selection of patient samples in this study is small and the source is single, and the different manifestations are not discussed in detail, which makes it impossible to verify the influence of these characteristics on the accuracy of diagnosis. In the future, it should be considered to increase the sample size of patients with thyroid cancer and further adopt a multicenter collaborative analysis method for research. In conclusion, the results provide a good theoretical basis for the optimized real-time SWE algorithm combined with TGAb in the clinical diagnosis of DTC.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.