临床研究
成人股骨CT建模及解剖学参数三维自动测量
中华医学杂志, 2019,99(39) : 3093-3099. DOI: 10.3760/cma.j.issn.0376-2491.2019.39.011
摘要
目的

通过成人股骨CT数据精确重建三维模型,并实现股骨解剖学参数自动化测量,研究国人股骨解剖学参数特点。

方法

在148例成年患者中应用Mimics17.0软件根据CT数据实行建模,把建成的模型导入Geomagic Studio软件中,提取股骨解剖学标志建立三维坐标系,统一坐标系;利用Matlab软件自编程序,采用最近点迭代ICP算法和解剖学标志点的三维自动提取算法进行精确定位,自动测量股骨解剖学参数,最后对数据进行分析比较。研究数据采用独立样本t检验进行比较。

结果

本组股骨髁左右径男性为(84.1±3.6)mm,女性为(74.8±3.3)mm;股骨髁前后径男性为(66.5±3.7)mm,女性为(61.2±3.5)mm;球头直径男性为(48.8±2.1)mm,女性为(43.4±2.2)mm,差异均有统计学意义(t=16.21、8.84、15.20,均P<0.05);而股骨矩弧度半径男性为(112.5±24.5)mm,女性为(124.7±19.2)mm(t=3.30,P=0.002)。颈干角男性为124.9°±4.0°,女性为126.1°±5.5°(t=1.40,P>0.05)。随着身高增长,大部分解剖参数均随着增加。

结论

股骨解剖参数三维自动测量数据较手动测量重复性高、精度高。建立和丰富股骨解剖数据库,设计研发符合国人的内固定产品十分必要。

引用本文: 薛磊, 李乾, 胡刚峰, 等.  成人股骨CT建模及解剖学参数三维自动测量 [J] . 中华医学杂志, 2019, 99(39) : 3093-3099. DOI: 10.3760/cma.j.issn.0376-2491.2019.39.011.
参考文献导出:   Endnote    NoteExpress    RefWorks    NoteFirst    医学文献王
扫  描  看  全  文

正文
作者信息
基金 0  关键词  0
English Abstract
评论
阅读 0  评论  0
相关资源
引用 | 论文 | 视频

版权归中华医学会所有。

未经授权,不得转载、摘编本刊文章,不得使用本刊的版式设计。

除非特别声明,本刊刊出的所有文章不代表中华医学会和本刊编委会的观点。

近年来,关节置换手术的开展不断增多。在髋关节置换手术中,为了获得满意的效果,术者期望股骨柄与股骨近端达到完美的初始匹配,这样可减少术中股骨劈裂及术后假体松动等并发症[1]。而在全膝关节置换术中一些并发症,例如假体松动和位置不良、髌股关节不稳及术后疼痛、屈伸活动受限[2],也可能与初始匹配不佳有关。有文献报道,股骨相关参数及形状受人种、环境、年龄、性别等多种因素的影响[3],并且国人股骨解剖参数与西方人存在差异[4,5]。因此,对股骨相关参数的精确测量显得尤为重要。本研究采用成人股骨CT数据进行精确建模,同时开发自动测量的程序,提升测量精度和速度,为建立国人的骨骼数据库提供方法学的参考,并在完善股骨解剖数据库的同时为假体的设计提供数据支持。

对象与方法
一、研究对象

11例股骨尸体骨标本;148例随机入选患者来源于苏州大学附属第二医院放射科2015年1月至2018年8月的行股骨全长CT扫描的成年患者。其中男83例,女65例,年龄(66±16)岁。纳入标准:健康的成年人,对试验方案知情同意。排除标准:股骨有任何骨折及骨缺损者、CT数据不全者、股骨头坏死者、骨关节退变严重者、既往有股骨手术病史者、股骨合并有肿瘤者。

二、研究方法
1.CT数据要求:

CT为薄层扫描(包含11例尸体骨标本),CT数据层厚为0.63 mm。处理软件:Mimics 17.0软件系统、Geomagic Studio软件系统和SPSS 18.0统计软件系统。

2.测量方法:

(1)建模:原始数据为包含股骨全长的CT图像,导入Mimics 17.0软件系统。建模采用阈值调整等技术,对未成像松质髓腔进行实体填充,骨骼边界精细化分割,重建股骨模型,对松质骨和髓腔的部分,通过实体填充与皮质骨求差,精确建立髓腔和松质骨的模型。并输出STL模型(图1[6]。(2)后处理:利用Geomagic Studio软件,导入STL模型,对模型进行松弛、删除钉状物、去噪、快速平滑处理,去除建模中的伪影和人为误差。同时对导入的所有模型进行预处理,初步统一模型坐标系,为下一步自动测量建立良好的初始条件[7]。并输出点云VTX格式(图2)。(3)编程自动测量参数:利用Matlab软件自编程序,导入VTX点云格式,采用最近点迭代ICP算法和解剖学标志点的三维自动提取算法进行精确定位,最后采用解剖学参数的三维自动测量(包括空间距离和角度)[8],精确获得每个股骨样本的解剖学参数。(4)用游标卡尺对11例尸体骨股骨头数据进行两次测量。手动第一次测量值记为A,第二次记为B;自动测量记为C;手动测量均值记为D。(5)把本组测量的生理前后径(AP值)、假体左右径(ML值)和市面上常见的股骨髁假体数据在坐标系中标出,并根据数据分布拟合直线,比较分析不同。

点击查看大图
图1
精确建立股骨和髓腔(含松质骨)三维模型:A为导入dicom图像;B为阈值分割填充等建立股骨表面模型;C为求差法建立髓腔和松质骨模型
点击查看大图
图1
精确建立股骨和髓腔(含松质骨)三维模型:A为导入dicom图像;B为阈值分割填充等建立股骨表面模型;C为求差法建立髓腔和松质骨模型
点击查看大图
图2
股骨三维模型的后处理过程:A为原始模型;B、C分别为后处理前后图像;C、D为坐标系统一
点击查看大图
图2
股骨三维模型的后处理过程:A为原始模型;B、C分别为后处理前后图像;C、D为坐标系统一
3.参数定义:

(1)股骨总长:为股骨头上缘到股骨髁下缘的距离。(2)球头直径:应用Geomagic软件中的最佳拟合球命令,将股骨头拟合成球体,找到球心后,测量该球体的直径记为股骨头直径。(3)颈干角:即为股骨颈的长轴与股骨干纵轴之间形成的夹角。(4)前倾角:为人体股骨颈的中轴线与股骨内外髁中点间的连线形成的夹角。(5)近端髓腔直径:截面0是股骨颈方向(股骨头中心和股骨颈中心连线)和近端髓腔的方向(股骨柄安装方向)的交点,以此为基准做水平面的辅助线,改线即为截面0,以后每隔20 mm为一个截面测量位置。(6)股骨髁AP及ML:为股骨髁内外侧缘和前后缘的大小。(7)生理内外侧后髁半径:把内外侧后髁拟合为球体,测量球体的半径即为后髁半径。(8)冠状面角度、髓腔冠状面角度:在股骨冠状面观测,股骨近端髓腔的方向和股骨干髓腔的方向所形成的夹角。(9)矢状面弧度半径、髓腔矢状面弧度半径:把骨干作为圆的一部分,找出其圆心后测量半径。(10)股骨距弧度半径:在股骨小粗隆侧测量,方法同股骨矢状面弧度。(11)颈长:即为股骨头中心至髓腔和股骨颈中线交点的距离。(12)假体ML:模拟按照股骨远端9 mm截骨后测量股骨髁左右侧距离即为假体ML(图3)。

点击查看大图
图3
股骨表面及其髓腔的参数定义示意图:A为股骨各长度参数;B、C、D为各角度参数;E为假体左右径(ML)
点击查看大图
图3
股骨表面及其髓腔的参数定义示意图:A为股骨各长度参数;B、C、D为各角度参数;E为假体左右径(ML)
4.结果比较:

对尸体骨股骨头自动测量和手动测量结果进行比较,分析本组结果与文献结果的差异,同时分析不同性别和体质指数(BMI)股骨测量参数的差异。

三、统计学方法

实验结果均采用SPSS 18.0软件进行统计分析,分析单个变量与解剖参数的相关性。测量数据资料用±s表示(长度参数保留1位小数),采用独立样本t检验进行比较。以P<0.05为差异有统计学意义。

结果
一、尸体骨股骨头自动测量和手动测量比较分析

测量值A-B平均绝对值大小为(1.7±2.1)mm,A-C为(0.7±0.5)mm,B-C为(1.8±2.2)mm,D-C为(1.1±1.1)mm。理论上,3个测量值应一致,但统计分析得出自动测量与手动测量存在差异,且手动测量重复性差,误差较大,精度低。

二、股骨解剖参数测量结果与国内外学者测量结果比较

在股骨头直径及颈干角方面,本组数据和国外数据无明显差异[9,10],但在总长及颈长这两项参数上,可看到国外学者[11,12]测得的数据与本组数据相差较大。与国内其他学者[13,14]测得的数据比较,在颈干角、颈长、球头直径方面,本组测量的数据无明显大小差异。而在生理ML、生理AP方面,本组与国内平均数值差异较大。

点击查看表格
表1

本研究与国内外学者数据比较(±s,mm)

表1

本研究与国内外学者数据比较(±s,mm)

项目本研究国外学者[9, 10, 11, 12]国内学者[13, 14]
总长428.5±24.9436.8±35.3-
球头直径46.5±3.546.1±4.847.1±3.5
颈长45.7±4.035.6±2.747.0±3.0
颈干角(°)125.5±4.8124.7±7.4128.1±5.7
前倾角(°)18.2±7.812.4±3.817.0±4.2
生理左右径80.0±5.881.2±6.075.6±3.8
生理前后经64.2±4.563.6±5.860.7±4.0

注:-:无数据

三、股骨解剖参数测量结果男女比较

在股骨近端方面,球头直径、颈长以及髓腔近端截面都是男性大于女性(均P<0.05)。而在颈干角、前倾角及股骨距弧度方面,男性则比女性小,其中颈干角、前倾角比较均P>0.05,而股骨距弧度半径比较P<0.05。股骨远端参数中,男性生理ML、假体ML、生理AP、生理内侧后髁半径、生理外侧后髁半径均大于女性(均P<0.05)。在股骨总长、冠状面角度、髓腔冠状面角度、矢状面弧度、髓腔矢状面弧度均为男性高于女性(均P<0.05),而冠状面角度、矢状面弧度半径方面两性虽有差异,但均P>0.05。

四、股骨解剖参数测量结果与身高的关系

本组受试者身高(163±8)cm(150~180 cm)。股骨干总长、球头直径、颈长、生理ML、生理AP、生理内外后侧髁半径、矢状面弧度以髓腔近端截面大小均随身高的增加而增大。但颈干角、前倾角、冠状面角度及髓腔冠状面角度未显示出与身高存在线性关系。其中身高在180~189 cm的病例仅2例,对于前倾角过大无明显统计学意义。

五、股骨解剖参数测量结果与体质指数的关系

在BMI从低于正常到肥胖的过程中,可看到球头直径、矢状面的弧度大致呈上升趋势,与之相反,颈干角、股骨距弧度则呈下降趋势。由于过度肥胖人数较少,统计意义不大,因此在统计过程中予以剔除,在球头直径、矢状面弧度方面,BMI大者高于BMI小者(t=2.15、2.29,均P<0.05);而颈干角、股骨矩弧度则相反(t=0.30、1.09,均P>0.05)。

六、股骨解剖参数中生理内外后侧髁大小比较

本组148例数据中股骨生理内侧后髁的大小为(19.1±1.5)mm,而生理外侧后髁为(18.3±1.6)mm,总体上内侧稍大(t=4.03,P=0.001)。并且生理内、外侧后髁并不是完全呈正相关(r=0.80,P<0.05),也就意味着内外侧后髁的伴随增长幅度不相同(图4)。

点击查看大图
图4
生理内外侧髁线性相关性
点击查看大图
图4
生理内外侧髁线性相关性
七、本组股骨生理AP、假体ML和常见假体数据比较分析
点击查看表格
表2

男女股骨参数数据比较(±s,mm)

表2

男女股骨参数数据比较(±s,mm)

性别球头直径颈长颈干角(°)前倾角(°)截面0截面20 mm截面40 mm截面60 mm截面80 mm假体ML
48.9±2.146.7±4.1125.0±4.117.9±7.257.5±6.435.0±4.322.8±2.817.6±2.415.3±2.175.2±2.8
43.4±2.244.5±3.5126.1±5.618.6±8.847.0±5.827.6±4.119.0±2.815.4±2.113.8±2.067.0±2.9
t15.23.41.45.410.4----16.21
P0.0010.010.170.590.001----0.001
性别总长生理ML生理AP生理内侧 后髁半径生理外侧 后髁半径冠状面 角度(°)矢状面弧 度半径髓腔冠状 面角度(°)髓腔矢状面 弧度半径股骨距 弧度半径
443.2±18.184.1±3.666.5±3.719.9±1.219.4±1.2174.5±4.4988.3±199.8176.5±2.4943.2±193.5112.5±24.5
409.7±18.874.8±3.361.2±3.518.1±1.217.1±1.1173.3±3.2927.9±193.3174.0±6.0878.5±183.0124.7±19.2
t10.8516.218.848.8411.961.841.853.152.073.30
P0.0010.0010.0010.0010.0010.070.070.0020.040.001

注:男83例,女65例;-:未统计;ML:左右径;AP:前后径

本组148例数据中股骨髁生理AP为(64.2±4.5)mm,假体ML值为(71.6±5.0)mm。本实验测得的数据与市面上常见的膝关节假体进行比较,可以看到每组数据之间存在一些差异,本组数据与施乐辉、史赛克、强生三组产品数据的差异均无统计学意义(均P>0.05,图5)。

点击查看大图
图5
本组数据与常见膝关节假体比较
点击查看大图
图5
本组数据与常见膝关节假体比较
点击查看表格
表3

不同身高受试者股骨解剖参数(±s,mm)

表3

不同身高受试者股骨解剖参数(±s,mm)

项目150~159 cm160~169 cm170~179 cm180~189 cm
总长405.4±17.8428.9±18.1450.9±14.6484.9±1.7
球头直径43.4±2.546.7±2.949.3±2.351.8±0.7
颈长43.6±3.745.4±3.248.1±3.655.0±0.6
颈干角(°)125.3±4.8126.0±4.9124.8±4.8125.4±1.1
前倾角(°)19.5±8.717.4±7.317.5±7.426.8±2.4
生理左右径74.9±3.880.3±4.884.7±3.989.8±3.0
生理前后径61.1±3.564.1±3.967.1±3.476.7±1.0
生理内侧后髁半径18.0±1.219.1±1.220.2±1.322.0±1.5
生理外侧后髁半径17.1±1.418.4±1.319.4±1.320.9±0.6
冠状面角度(°)172.8±3.1174.9±2.6173.8±5.9170.6±6.0
矢状面弧度半径887.1±164.8970.4±203.91 025.3±230.81 093.8±179.1
髓腔冠状面角度(°)173.4±6.8176.7±2.1175.6±3.4174.5±1.4
髓腔矢状面弧度半径845.3±159.6925.0±189.8973.3±206.0980.8±253.8
股骨距弧度半径123.2±21.2114.7±22.0117.1±26.4117.2±22.3
截面0 mm48.0±7.353.1±7.257.6±7.560.0±4.1
截面20 mm27.8±4.332.1±5.135.2±4.641.1±5.2
截面40 mm19.2±2.521.1±3.223.0±2.728.0±7.6
截面60 mm15.4±1.816.5±2.317.9±2.321.8±6.7
截面80 mm13.8±1.714.5±2.115.4±2.019.5±5.7
截面100 mm12.9±1.713.4±2.014.1±1.817.7±4.8
点击查看表格
表4

不同体质指数受试者股骨解剖参数(±s,mm)

表4

不同体质指数受试者股骨解剖参数(±s,mm)

项目<18.5≥18.5~24.0≥24.0~28.0≥28.0~32.0≥32.0
总长421.5±19.5427.9±28.9430.4±20.4429.2±21.8415.6±7.1
球头直径44.9±3.046.0±3.647.1±3.147.4±4.147.1±1.8
颈长44.0±3.346.0±4.445.8±3.544.9±3.643.2±2.2
颈干角(°)127.0±3.7125.4±5.2125.5±3.8123.8±5.6131.4±7.6
前倾角(°)26.0±5.517.3±7.919.0±7.418.2±7.27.1±9.4
生理左右径75.8±5.379.2±5.781.3±5.481.3±6.278.4±9.6
生理前后径63.0±3.463.6±4.764.8±4.165.5±4.963.3±5.9
生理内侧后髁半径18.4±1.719.0±1.519.2±1.419.4±1.819.1±1.9
生理外侧后髁半径17.8±1.718.2±1.718.5±1.418.7±1.918.3±1.7
冠状面角度(°)175.0±2.0173.6±5.0174.5±2.8173.3±2.0173.6±4.1
矢状面弧度半径834.7±169.7935.3±199.11 004.5±202.01 018.0±163.5819.5±52.3
髓腔冠状面角度(°)176.0±1.6174.8±5.8176.1±2.7175.2±2.6178.0±2.4
髓腔矢状面弧度半径785.5±169.4894.0±193.7952.4±192.8959.1±156.5781.8±19.1
股骨距弧度半径120.9±16.9119.7±24.0115.7±23.6116.6±20.4107.0±20.0
截面049.8±9.051.7±7.954.9±7.652.4±9.758.2±7.8
截面60 mm15.7±2.816.6±2.716.8±2.016.7±2.717.0±2.1
截面100 mm12.9±2.213.7±2.213.4±1.713.3±2.013.8±3.5

注:体质指数单位kg/m2

讨论

此前国外有学者已对股骨解剖学进行了相当多的研究,国内也有相当一部分的学者对国人的股骨解剖学参数进行了报道[15],并且对国内外数据进行了对比[16]。但是国内很多测量数据是通过X线直接测量获得的[17],因患者的拍摄姿势和测量医师的水平差异,和国外的三维精确测量难以精确对比分析,测量精度有待提高,而且测量方法也多为人工测量,效率较低,可重复新差,因此测量的误差也比较大。也有学者进行三维数据的手动测量,这种测量可比X线在一定程度上提高精度,但相比于三维自动测量仍有一些不足。既往有学者已提出自动测量法比手动测量方法更加精确[8],从统计数据中可看到手动测量和自动测量的平均误差在1.1 mm,但两次手动测量的误差较大,平均误差在1.7 mm,这也验证了自动测量的精度高于手动测量。本研究采用对股骨CT数据建模,在Geomagic Studio软件中统一坐标系,对解剖参数进行自动化测量,对于解剖标志点的提取是固定的,并且可对数据进行大样本及重复性测量检验,这是手动测量不能实现的。因此测量的精度得到了很大的提高、测量的效率也得到了很大提升,并且此种测量方法可减少人为误差。

本研究测量所得结果与国外测量数据存在一定差异。在股骨总长及股骨颈长的数据上,国内外存在较大差别,尤其在前倾角上,本组实验数据为18.2°±7.8°,而国外测量的数据为12.4°±3.8°,两者数据差异较大,在总长和颈长上两者也存在差异,这可能与欧美人与国人在体格以及种族上存在差异,并且本研究发现这两者的数据随着身高的增加而增加,因此国内外在这两者数据上存在差异也符合本研究结果。同时在本研究中,我们发现男性和女性在多项参数上存在差异,例如生理ML男性平均为84.1 mm,女性平均为74.8 mm;生理AP男性平均为66.5 mm,女性平均为61.2 mm,差异均有统计学意义,但目前对于骨科内植物产品暂未区分男女,这可能导致一部分患者的假体与其本身不相匹配,进而导致后续内固定失败或假体松动。在本研究中,我们加入了BMI这个分类,但并未发现明显的测量数据变化规律;同时,我们在自动三维测量方法中测量了股骨矩弧度半径及股骨干矢状面弧度半径,其中男女在股骨矩弧度半径上存在显著差异,这两项数据对于髓内钉及股骨柄的设计有着重要的参考意义。我们还发现,股骨内外侧后髁大小并不完全一致,生理内侧后髁的大小为(19.1±1.5)mm,而生理外侧后髁的大小为(18.3±1.6)mm,我们根据数据做了拟合直线图,因此在假体的设计中若把内外侧髁大小设计成一样大小就可能带来假体的不匹配。并且我们分析了生理AP及假体ML值上与关节外科中常见股骨髁假体数据进行对比,发现假体数据与本组统计数据存在差异,没有完全符合的数据,这也可以为设计符合国人的股骨髁假体提供参考。同时研究发现股骨近端髓腔值越往远端其个体差异越来越小,这也意味着我们在设计股骨柄假体时应重点考虑靠近股骨头侧数据变化。

本研究的也存在一些不足之处,比如样本数量还不够大,缺少我国其他地区的样本,并且缺少少数民族人群统计,而这些因素可能会对股骨相关参数产生一定的影响,因此还需进一步研究。

目前是信息化高度发展的社会,精准医疗也是越来越受到大家的关注,设计出更加符合国人股骨解剖结构的内固定器材是未来骨科发展的一个趋势。并且此项测量技术同样也可用在术前准备上,在术前可以精确测量出术中所需的解剖参数,为假体的选择及其放置角度提供一定参考,以减少术中医源性损伤。

利益冲突
利益冲突

所有作者均声明不存在利益冲突

参考文献
[1]
IssaK, PivecR, WuestemannT, et al. Radiographic fit and fill analysis of a new second-generation proximally coated cementless stem compared to its predicate design[J]. J Arthroplasty, 2014, 29(1):192-198. DOI: 10.1016/j.arth.2013.04.029.
[2]
EbraheimNA, RayJR, WandtkeME, et al. Systematic review of periprosthetic tibia fracture after total knee arthroplasties[J]. World J Orthop, 2015,6(8):649-654. DOI: 10.5312/wjo.v6.i8.649.
[3]
陈成旺温宏潘孝云,.年龄及性别对股骨近端解剖参数的影响及其临床意义[J].浙江医学,2014,36(1):5-8.
[4]
张晟余斌罗吉伟,.股骨近端髓腔三维解剖测量[J].中国临床解剖学杂志,2010,28(3):258-261.
[5]
YangZ, JianW, LiZH, et al. The geometry of the bone structure associated with total hip arthroplasty[J]. PLoS One, 2014,9(3):e91058. DOI: 10.1371/journal.pone.0091058.
[6]
ZhuZ, LiG. Construction of 3D human distal femoral surface models using a 3D statistical deformable model[J]. J Biomech, 2011, 44(13):2362-2368. DOI: 10.1016/j.jbiomech.2011.07.006.
[7]
朱忠林丁辉王广志.利用可变形统计模型进行膝关节建模与运动测量[J].清华大学学报(自然科学版),2013,53(1): 139-144.
[8]
GottschlingH, SchröderM, ReimersN, et al. A system for performing automated measurements on large bone databases[G]. Munich: World Congress on Medical Physics and Biomedical Engineering, 2009: 910-913.
[9]
MassinP, GeaisL, AstoinE, et al. The anatomic basis for the concept of lateralized femoral stems: a frontal plane radiographic study of the proximal femur[J]. J Arthroplasty, 2000,15(1):93-101.
[10]
NoblePC, AlexanderJW, LindahlLJ, et al. The anatomic basis of femoral component design[J]. Clin Orthop Relat Res, 1988, (235):148-165.
[11]
GriffinFM, MathK, ScuderiGR, et al. Anatomy of the epicondyles of the distal femur: MRI analysis of normal knees[J]. J Arthroplasty, 2000,15(3):354-359.
[12]
KuoTY, SkedrosJG, BloebaumRD. Measurement of femoral anteversion by biplane radiography and computed tomography imaging: comparison with an anatomic reference[J]. Invest Radiol, 2003,38(4):221-229. DOI: 10.1097/01.RLI.0000059542.90854.EF.
[13]
范吉星李宁龚晓峰,.老年髋部骨折患者股骨近端形态学参数的评估[J].中华医学杂志,2017,97(31):2443-2446. DOI: 10.3760/cma.j.issn.0376-2491.2017.31.012.
[14]
FanL, XuT, LiX, et al. Morphologic features of the distal femur and tibia plateau in Southeastern Chinese population: A cross-sectional study[J]. Medicine (Baltimore), 2017, 96(46): e8524. DOI: 10.1097/MD.0000000000008524.
[15]
孙月芳莫展豪程凯亮,.股骨远端及后髁曲率半径的性别差异及其意义[J].吉林大学学报(医学版),2010,36(1): 159-162.
[16]
汪轶平张恒辉王燎,.股骨近端解剖参数的自动化三维测量[J].医用生物力学,2016, 31(1): 1-7. DOI: 10.3871/j.1004-7220.2016.01.001.
[17]
皮益刚赵耀超王万春,. 500例湖南地区股骨近端形态样本测量分析[J].中南大学学报(医学版),2013,(9):925-930. DOI: 10.3969/j.issn.1672-7347.2013.09.010.
 
 
展开/关闭提纲
查看图表详情
回到顶部
放大字体
缩小字体
标签
关键词