# Dosiomics: Extracting 3D Spatial Features From Dose Distribution to Predict Incidence of Radiation Pneumonitis

- Department of Radiation Oncology, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China

Radiation pneumonitis (RP) is one of the major toxicities of thoracic radiation therapy. RP incidence has been proven to be closely associated with the dosimetric factors and normal tissue control possibility (NTCP) factors. However, because these factors only utilize limited information of the dose distribution, the prediction abilities of these factors are modest. We adopted the dosiomics method for RP prediction. The dosiomics method first extracts spatial features of the dose distribution within ipsilateral, contralateral, and total lungs, and then uses these extracted features to construct prediction model via univariate and multivariate logistic regression (LR). The dosiomics method is validated using 70 non-small cell lung cancer (NSCLC) patients treated with volumetric modulated arc therapy (VMAT) radiotherapy. Dosimetric and NTCP factors based prediction models are also constructed to compare with the dosiomics features based prediction model. For the dosimetric, NTCP and dosiomics factors/features, the most significant single factors/features are the mean dose, parallel/serial (PS) NTCP and gray level co-occurrence matrix (GLCM) contrast of ipsilateral lung, respectively. And the area under curve (AUC) of univariate LR is 0.665, 0.710 and 0.709, respectively. The second significant factors are V_{5} of contralateral lung, equivalent uniform dose (EUD) derived from PS NTCP of contralateral lung and the low gray level run emphasis of gray level run length matrix (GLRLM) of total lungs. The AUC of multivariate LR is improved to 0.676, 0.744, and 0.782, respectively. The results demonstrate that the univariate LR of dosiomics features has approximate predictive ability with NTCP factors, and the multivariate LR outperforms both the dosimetric and NTCP factors. In conclusion, the spatial features of dose distribution extracted by the dosiomics method effectively improves the prediction ability.

## Introduction

Radiation pneumonitis (RP) is one of the major toxicities of thoracic radiation therapy. The clinical symptoms range from fever, cough to pulmonary function failure, which may occur during the first 6 months after irradiation. Reducing the prescription dose could lower the risk of RP incidence, but also impairs tumor control. An accurate RP predictor (or prediction model) is desired to “safely” irradiate the tumor target without increasing the risk of RP incidence.

RP incidence is directly associated with the dose distribution within lung volume. Dosimetric factors, such as mean lung dose (MLD) and the lung volume within which the dose is greater than xGy (V_{x}), are widely used for RP prediction. Boonyawan et al. reported that RP incidence increases with V_{10} and V_{20} (1). Ramella et al. found that RP incidence is associated with V_{20} and V_{30} (2). Briere et al. reported that RP incidence significantly increases if the sparing lung volume is <1852cc (receiving dose≤40Gy) (3). Pinnix et al. reported that V_{5} has better prediction capability than V_{10}, V_{15}, and V_{20} (4). Palma et al. analyzed 836 patient cases from international institutions and concluded that symptomatic RP is associated with V_{20}, and fatal RP associated with the mean dose per day during treatment (5). Those studies demonstrate that the dosimetric factors are associated with RP incidence. Although enlightening for understanding the causes of RP incidence, the conclusion of those studies differs from individual institution or dataset.

The dosimetric factors only utilize partial information contained in the dose distribution. For instance, V_{x} is only a discrete point on the dose volume histogram (DVH) curve. Compared with dosimetric factors, the normal tissue complication probability (NTCP) model utilizes all information of the DVH curve by compressing the entire curve to a single factor with dose response functions (DRFs). Better prediction performance could be achieved by fitting the DRF parameters (6–8). And the standard deviation of the fitted parameters is 16.5% between different institutions (9). The improved prediction ability of NTCP factors are partly contributed by the introduction of DRF, but more importantly because the NTCP factors utilize all information of DVH curve. However, since the DVH curve does not take the spatial information into consideration, completely different dose distributions may result in identical DVH curve, and then identical NTCP factors.

By investigating the prediction ability of dosimetric and NTCP factors, it is reasonable to hypothesize that if the spatial information of dose distribution were utilized properly, the prediction ability should be further improved. The recently emerged radiomics method extracts numerous spatial features from medical images and uses these features to predict therapeutic responses (10–12). Enlightened by those works, the “**dosiomics**” method has been proposed, which attempts to extract the spatial features from dose distribution for radiotherapy response prediction (13–15). In this work, we adopted the dosiomics method for RP incidence prediction. The prediction ability of the dosiomics features is validated using 70 non-small cell lung cancer (NSCLC) patients treated with volumetric modulated arc therapy (VMAT) radiotherapy, and further evaluated by comparing with the dosimetric and NTCP factors.

## Methods and Materials

### Patient Data

Seventy NSCLC patients treated in our institution from 2013 to 2016 are used in this study. All patients were treated with 6MV VMAT without surgical operation. Treatment plans were designed using Pinnacle treatment planning system (v 9.0). The slice spacing of planning CT image was 5 mm, and the grid spacing of dose calculation was 4 × 4 × 4 mm. The dose was prescribed to 95% of the planning target volume (PTV). RP was graded from 0 to 5 according to Common Terminology Criteria for Adverse Events (CTCAE v3.0). Bootstrap method is adopted to address the issue of limited dataset. Bootstrap of original dataset is performed 1,000 times, and the resulting 1,000 bootstrap samples are used as training datasets for both the univariate and multivariate LR.

### Feature Extraction

The dosiomics features are extracted from the dose distribution within ipsilateral, contralateral and total lungs, separately. The extracted features are a set of indices, such as autocorrelation, sum of squares and cluster prominence etc, derived from the gray level co-occurrence matrix (GLCM) and gray level run length matrix (GLRLM). The calculation formulas can be referred in (16). All the extracted features are normalized to zero mean and unit variance with z-score normalization before further processing.

### Univariate Analysis

The endpoint of this study is the occurrence of grade ≥2 RP. Prediction model is built using LR, of which the coefficient is derived using maximum likelihood estimation method. For each extracted feature, univariate LR is performed 1,000 times based on the bootstrap samples of original dataset. The mean AUC of training datasets (bootstrap samples) is calculated. The most significant feature is determined as the feature with maximal mean training AUC. The final coefficient is determined as the median of the 1,000 resulting coefficients. The predictive ability of each single feature is evaluated on the original entire dataset by the area under curve (AUC) of receiver operating characteristic (ROC) (17).

### Multivariate Analysis

According to the one tenth rule, the number of predictors of multivariate LR should be one tenth of the outcome events. Recent studies show that the one tenth rule is generally too conservative, and can be relaxed for the logistic and Cox regression (18). The outcome event (grade ≥2 RP occurrence) of this study is 15, and the number of predictors of multivariate LR is relaxed to 2.

In this study, all possible two-feature combinations are traversed to search for the optimal combination. For each combination, multivariate LR and Spearman test between the two features are performed 1,000 times using bootstrap samples. The optimal combination is determined as the combination with maximal mean training AUC while the mean Spearman correlation within [−0.8, 0.8]. The Spearman correlation threshold excludes the combination of strongly correlated features to prevent overfitting. The final coefficient is also determined as the median of the 1,000 resulting coefficients, and the multivariate LR model is validated on the entire dataset.

### Dosimetric and NTCP Factors

The predictive ability of dosiomics features is further evaluated by comparing with the dosimetric and NTCP factors. Both univariate and multivariate analyses are performed using: 1. dosimetric factors and 2. NTCP factors. As listed in Table 1, the dosimetric factors include V_{5}, V_{10}, V_{15}, V_{20}, and MLD. The NTCP factors are two sets of equivalent uniform dose (EUD) and NTCP factors of Lyman (19) and parallel/serial (PS) models (20). All factors are calculated for the ipsilateral, contralateral and total lungs, separately.

With DVH reduction technique, the heterogeneous dose distribution is reduced to a single EUD_{L} parameter, which by definition has the same complication probability with the original heterogeneous dose distribution. Lyman model assumes that all sub-volumes are of the same contribution to side effects, and the NTCP factor (NTCP_{L}) is calculated as:

where *TD*_{50} is the tolerance dose, under which the complication probability of normal tissue is 50%. *m* is the curve slope at *TD*_{50}.

The PS model considers the normal tissues are composed of a large number of microscopic functional sub-volumes. The parallel model supposes the sub-volumes are independent. Side effects occur only if a significant number of sub-volumes are destroyed, and small portion of the damaged sub-volumes will not cause side effects. In contrast, the serial model assumes the sub-volumes are dependent. Side effects occur even only a small portion of sub-volumes are damaged. The entire volume is an arbitrary mixture of both serial and parallel sub-volumes. The NTCP (NTCP_{PS}) factor is calculated as:

where *k* is the ratio of serial to parallel sub-volume. Like Lyman model, *m* is the slope of NTCP curve at *TD*_{50}. EUD_{PS} is calculated according to formula proposed in (21):

The study presented in (22) optimized the parameters of NTCP models for grade≥2 RP prediction based on 382 thoracic patient cases. As the endpoint is the same and patient population is similar with our study, we directly used the same parameters, which are listed in Table 2.

The dose distribution within lung volumes is derived from the treatment plan data in DicomRT format. This procedure is implemented with Matlab software (MathWorks, Natick, MA). Feature extraction is implemented with the python pyradiomics package (v2.0.0) (16). LR is implemented in R language, with the stats (v3.4.1) package (23).

## Results

Clinical factors of analyzed patients are listed in Table 3. Figure 1A shows the dose distribution within total lungs, and Figures 1B,C shows the corresponding GLCM and GLRLM in logarithmic scale for the sake of clearance. Twenty and sixteen features, which are commonly used in radiomics studies, are derived from GLCM and GLRLM, respectively. All the features are calculated for the dose distribution within ipsilateral, contralateral and total lungs, separately. In total, 129 [(27 + 16) × 3] features are extracted for each patient.

**Figure 1**. Intermediate results of dosiomics method. **(A)** 3D dose distribution, **(B)** GLCM, and **(C)** GLRLM of total lung.

### Univariate Analysis

Table 4 lists the results of univariate analysis. The most significant predictors of the three sets of factors/features are the MLD (MLD^{I}), NTCP_{PS} (NTCP${}_{\text{PS}}^{\text{I}}$) and GLCM (GLCM^{I}) contrast of ipsilateral lung, respectively. The GLCM contrast measures local variation of dose distribution. It is interesting to notice that the most significant predictors are all derived from the dose distribution within ipsilateral lung. The median and 10–90th% odds ratio (OR) of the most significant factor/feature and constant term are listed in Table 4. The OR can be interpreted as: increasing the corresponding factor by one unit the probability of RP incidence increases by *OR* times. Increasing MLD^{I}, NTCP${}_{\text{PS}}^{\text{I}}$ and GLCM^{I} contrast by one unit will increase the probability of RP incidence by 1.667, 2.041 and 2.010 times. The range of 10–90th% OR measures the repeatability of the derived predictive model. The 10–90th% OR range of GLCM^{I} contrast is greater than MLD^{I} but lower than NTCP${}_{\text{PS}}^{\text{I}}$. The OR of constant term is the ratio of RP probability over non-RP probability of the dataset without using any predictor. For the three LR models, the median and 10–90th% ORs of constant term are almost identical. The AUC values of NTCP and dosiomics factor/feature are approximate, but higher than dosimetric factors.

### Multivariate Analysis

Table 5 lists the results of multivariate analysis. The optimal combinations all contain the most significant single predictors. The second significant predictors are V_{5} of contralateral lung (V${}_{5}^{C}$), EUD_{PS} (EUD${}_{\text{PS}}^{C}$) of contralateral lung and GLRLM (GLRLM^{T}) low gray level run emphasis of total lungs, respectively. V${}_{5}^{C}$ represents the lower dose within contralateral lung. EUD${}_{\text{PS}}^{C}$ by definition is the uniform dose that has the same complication probability with the original heterogeneous dose distribution. The GLRLM^{I} low gray level run emphasis measures the region of low dose, with a higher value indicating a greater concentration of low dose distribution. The most significant single factors/features are all extracted from ipsilateral lung, while the second from either contralateral or total lungs. As shown in Figure 2, the factors/features extracted from same dose distributions are more correlated, while the factors/features extracted from different dose distributions are less correlated, especially the factors/features of ipsilateral and contralateral lungs. In order to prevent overfitting, the strongly correlated factors/features are excluded. This explains why the second predictors are derived from either contralateral or total lungs.

**Figure 2**. Mean Spearman correlation of 1,000 bootstrap samples. **(A–C)** mean Spearman correlation of dosimetric factors, NTCP factors, and dosiomics features. **(D)** zoomed correlation of the features extracted from ipsilateral and total lungs. All the 4 figures are diagonal symmetric. For **(A–C)**, the features are sorted in the order of ipsilateral, contralateral, and total lungs from left to right and from top to down. In **(D)**, the features are sorted the order of GLCM and GLRLM. The correlation of the features of the optimal combination are denoted with red circle. Cor, correlation; Emp, emphasis; GL, gray level; Homo, homogeneity; LRun, long run; SRun, short run; N-uni, non-uniformity; RL, run length.

The mean Spearman correlation between the two selected factors/features of 1,000 bootstrap samples and its standard deviation are also listed in Table 5. MLD^{I} and V${}_{5}^{C}$ is positive correlated. This is because increasing the MLD of ipsilateral lung would increase the scatter dose delivered to contralateral lung thus increase the value of V_{5}. For NTCP and dosiomics factors/features, the Spearman correlation is negative and of lower magnitude, indicating that the selected predictors are weakly negative correlated.

For dosiomics features, the increase of AUC is obvious when switching from univariate LR to multivariate LR. On the other hand, the increase of AUC for dosimetric and NTCP factors is limited. This is because either the dosimetric or NTCP factors describe the dose distribution from the similar perspective. Adding another predictor will not significantly improve the predictive ability. On the contrary, the dosiomics features display a rich diversity, which is benefit for revealing the hidden correlation with RP incidence.

## Discussion

We investigated the published studies on the correlation between dosimetric factors and RP incidence, and found the conclusions differ from individual institution or dataset. The quantitative analysis of normal tissue effects in the clinic (QUANTEC) summarized available published data and performed a logistic regression between MLD and RP (9). Despite the differences in patient selection and RP grade of published data, an obvious trend could be observed: the probability of RP incidence increases with MLD. This conclusion supports our finding: MLD^{I} is the most significant dosimetric predictor. Most published studies on the correlation of NTCP factors and RP incidence focus on fitting the parameters of NTCP models to better predict RP incidence. In this study, we directly used the optimized parameters presented in (22), and found that NTCP${}_{PS}^{\text{I}}$ is the most significant predictor. Tsougos et al. (7) also reported that PS model outperforms the rest NTCP models for RP (grade 2) prediction of breast cancer radiotherapy. Both studies demonstrate that RP occurs if significant sub-volumes are damaged. This conclusion is further validated by the study reported in (3), which found RP incidence significantly increases if the sparing lung volume (dose ≤ 40Gy) is less than 1852cc.

The results of multivariate LR demonstrate that the prediction ability of dosiomics features outperform dosimetric and NTCP factors. Meanwhile the NTCP factors has better performance than the dosimetric factors. The results validate the hypothesis that the predictive ability improves with more information of the dose distribution are used by the prediction model. The application of dosiomics method is not limited to RP prediction. It is suitable for any radiotherapy outcome, either positive (like survival, control rate) or negative (like normal tissue damage, complication).

We have to admit that the patient dataset of this study is limited. In order to address this issue, the bootstrap approach is adopted. The coefficients of univariate and multivariate LR is determined as the median of the fitting results using bootstrap samples. The predictor number of multivariate LR is set to 2 to avoid overfitting. For dosiomics features, the improvement on predictive ability is significant when switching from univariate LR to multivariate LR. With the diverse dosiomics features, it is reasonable to assume that the predictive ability could be further improved when the patient dataset is enlarged with more positive cases, in which case more predictors could be included in multivariate LR.

Except for enlarging the patient dataset, another promising direction to further improve the predictive ability is to dig deeper of the patient dataset. Dose distribution, even strongly correlated with RP incidence, is not the unique factor. Other clinical factors, such as the age, smoking history, chemotherapy, are also found to be associated with RP incidence. In addition, other “omics” features, such as the features extracted from radiomics, genome and proteomics, also provide insightful information. All these factors/features could be integrated into the model to improve its prediction ability and robustness.

Although the dosiomics features demonstrate good prediction ability, the understanding of these features is still qualitative. The main reason is the process of transform dose distribution into GLCM and GLRLM cannot be accurately described with analytic function. Therefore the features based on GLCM and GLRLM are not as simple and straightforward as dosimetric factors. Based on the finding of this study, we qualitatively conclude that the higher local dose variation within ipsilateral lung and the greater low dose region of total lungs, the greater probability of RP incidence. Furthermore, how to utilize the features for treatment plan design is not quite clear. In other words, currently the dosiomics method is limited to RP incidence prediction. How to use the dosiomics method to prevent RP incidence during treatment planning is one goal of our future study.

## Ethics Statement

This study was carried out in accordance with the declaration of Helsinki and approved with exemption from informed consent by the independent ethics committee of cancer hospital, Chinese academy of medical sciences (No. NCC2015 G-15).

## Author Contributions

BL, HY, and JD conceived the project and wrote the paper. XC, YT, and LY collected and analyzed the data. TZ, ZZ, and LW provided expert clinical knowledge. All authors edited the manuscript.

## Funding

This work was supported by the National Natural Science Foundation of China (11475261, 81801799 and 81502649) and the National Key R&D Program of China (2016YFC0904600).

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

1. Boonyawan K, Gomez DR, Komaki R, Xu Y, Nantavithya C, Allen PK, et al. Clinical and dosimetric factors predicting grade≥ 2 radiation pneumonitis after postoperative radiotherapy for patients with non-small cell lung carcinoma. *Int J Radiat Oncol Biol Phys.* (2018) 101:919–26. doi: 10.1016/j.ijrobp.2018.04.012

2. Ramella S, Trodella L, Mineo TC, Pompeo E, Stimato G, Gaudino D, et al. Adding ipsilateral V20 and V30 to conventional dosimetric constraints predicts radiation pneumonitis in stage IIIA-B NSCLC treated with combined-modality therapy. *Int J Radiat Oncol Biol Phys*. (2010) 76:110–5. doi: 10.1016/j.ijrobp.2009.01.036

3. Briere TM, Krafft S, Liao Z, Martel MK. Lung size and the risk of radiation pneumonitis. *Int J Radiat Oncol Biol Phys*. (2016) 94:377–84. doi: 10.1016/j.ijrobp.2015.10.002

4. Pinnix CC, Smith GL, Milgrom S, Osborne EM, Reddy JP, Akhtari M, et al. Predictors of radiation pneumonitis in patients receiving intensity modulated radiation therapy for Hodgkin and non-Hodgkin lymphoma. *Int J Radiat Oncol Biol Phys*. (2015) 92:175–82. doi: 10.1016/j.ijrobp.2015.02.010

5. Palma DA, Senan S, Tsujino K, Barriger RB, Rengan R, Moreno M, et al. Predicting radiation pneumonitis after chemoradiation therapy for lung cancer: an international individual patient data meta-analysis. *Int J Radiat Oncol Biol Phys*. (2013) 85:444–50. doi: 10.1016/j.ijrobp.2012.04.043

6. Avanzo M, Trovo M, Furlan C, Barresi L, Linda A, Stancanello J, et al. Normal tissue complication probability models for severe acute radiological lung injury after radiotherapy for lung cancer. *Phys Med*. (2015) 31:1–8. doi: 10.1016/j.ejmp.2014.10.006

7. Tsougos I, Mavroidis P, Rajala J, Theodorou K, Jarvenpaa R, Pitkanen MA, et al. Evaluation of dose-response models and parameters predicting radiation induced pneumonitis using clinical data from breast cancer radiotherapy. *Phys Med Biol.* (2005) 50:3535–54. doi: 10.1088/0031-9155/50/15/004

8. Yorke ED, Jackson A, Rosenzweig KE, Merrick SA, Gabrys D, Venkatraman ES, et al. Dose-volume factors contributing to the incidence of radiation pneumonitis in non-small-cell lung cancer patients treated with three-dimensional conformal radiation therapy. *Int J Radiat Oncol Biol Phys*. (2002) 54:329–39. doi: 10.1016/S0360-3016(02)02929-2

9. Marks LB, Bentzen SM, Deasy JO, Kong F-M, Bradley JD, Vogelius IS, et al. Radiation dose-volume effects in the lung. *Int J Rad Oncol Biol Phys.* (2010) 76:S70–S6. doi: 10.1016/j.ijrobp.2009.06.091

10. Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More than Pictures, They Are Data. *Radiology*. (2016) 278:563–77. doi: 10.1148/radiol.2015151169

11. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. *Nat Commun*. (2014) 5:4006. doi: 10.1038/ncomms5006

12. Lambin P, Riosvelazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, et al. Radiomics: extracting more information from medical images using advanced feature analysis. *Eu J Cancer.* (2012) 48:441–6. doi: 10.1016/j.ejca.2011.11.036

13. Rossi L, Bijman R, Schillemans W, Aluwini S, Cavedon C, Witte M, et al. Texture analysis of 3D dose distributions for predictive modelling of toxicity rates in radiotherapy. *Rad Oncol.* (2018) 129:548–53. doi: 10.1016/j.radonc.2018.07.027

14. Gabryś HS, Buettner F, Sterzing F, Hauswald H, Bangert M. Design and selection of machine learning methods using radiomics and dosiomics for normal tissue complication probability modeling of xerostomia. *Front Oncol.* (2018) 8:35. doi: 10.3389/fonc.2018.00035

15. Zhen X, Chen J, Zhong Z, Hrycushko B, Zhou L, Jiang S, et al. Deep convolutional neural network with transfer learning for rectum toxicity prediction in cervical cancer radiotherapy: a feasibility study. *Phys Med Biol.* (2017) 62:8246–63. doi: 10.1088/1361-6560/aa8d09

16. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational radiomics system to decode the radiographic phenotype. *Cancer Res.* (2017) 77:e104–e7. doi: 10.1158/0008-5472.CAN-17-0339

17. Fawcett T. An introduction to ROC analysis. *Patt Recogn Lett.* (2006) 27:861–74. doi: 10.1016/j.patrec.2005.10.010

18. Vittinghoff E, McCulloch CE. Relaxing the Rule of Ten Events per Variable in Logistic and Cox Regression. *Am J Epidemiol.* (2007) 165:710–8. doi: 10.1093/aje/kwk052

19. Lyman JT. Complication probability as assessed from dose-volume histograms. *Radiat Res Suppl*. (1985) 8:S13–9.

20. Kallman P, Agren A, Brahme A. Tumour and normal tissue responses to fractionated non-uniform dose delivery. *Int J Radiat Biol*. (1992) 62:249–62.

21. Mavroidis P, Lind BK, Brahme A. Biologically effective uniform dose (D) for specification, report and comparison of dose response relations and treatment plans. *Phys Med Biol*. (2001) 46:2607–30. doi: 10.1088/0031-9155/46/10/307

Keywords: dosiomics, radiomics, dose distribution, pneumonitis prediction, logistic regression

Citation: Liang B, Yan H, Tian Y, Chen X, Yan L, Zhang T, Zhou Z, Wang L and Dai J (2019) Dosiomics: Extracting 3D Spatial Features From Dose Distribution to Predict Incidence of Radiation Pneumonitis. *Front. Oncol.* 9:269. doi: 10.3389/fonc.2019.00269

Received: 23 December 2018; Accepted: 25 March 2019;

Published: 12 April 2019.

Edited by:

John Varlotto, University of Massachusetts Medical School, United StatesReviewed by:

Hongqing Zhuang, Peking University Third Hospital, ChinaWenzheng Sun, Shandong University, China

Copyright © 2019 Liang, Yan, Tian, Chen, Yan, Zhang, Zhou, Wang and Dai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Lvhua Wang, wlhwq@yahoo.com

Jianrong Dai, dai_jianrong@cicams.ac.cn