<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="zh-CN">
 <head>
  <meta http-equiv="Content-Type" content="text/html; charset=utf-8" pages_book="12-15" pages_ebook="19-22" column_name="专论" author_name="梁永刚，付丽媛，陈坚，黄志峰，熊晖，李良杰，钟群，肖慧，许尚文，陈自谦">
  <link rel="stylesheet" type="text/css" href="css/style.css">
  <title>基于SMR170体模磁共振质量控制线性度自动测算程序设计</title>
 </head>
 <body>
  <h1 class="title" data-meta-name="篇名">基于SMR170体模磁共振质量控制线性度自动测算程序设计</h1>
  <p class="author_info" data-meta-name="作者信息"><span title="pagenumber_ebook=19,pagenumber_book=12" class="emphasis_bold">梁永刚<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，付丽媛<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，陈坚<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，黄志峰<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，熊晖<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，李良杰<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，钟群<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，肖慧<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，许尚文<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>，陈自谦<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1,2</span></span></p>
  <p class="author_info" data-meta-name="作者信息">1.联勤保障部队第900医院（原南京军区福州总医院）放射诊断科，福建福州350025；2.厦门大学附属东方医院放射诊断科，福建福州350025</p>
  <p class="meta" data-meta-name="摘要"><span data-meta-name="摘要.标题" class="key">[摘要]</span>目的设计基于SMR170体模磁共振成像系统质量控制检测的线性度自动测量程序，实现线性度参数的自动测算。方法在Matlab平台设计编写了图像处理算法，对线性度参数进行自动测量，并采用Bland-Altman统计学方法对手工测量和自动测量两种方法的一致性进行了分析。结果利用Matlab平台设计编写的自动检测方法可实现线性度的自动检测，并且与手工测量方法有较好的一致性。结论线性度自动检测程序与手工测量可相互替代。</p>
  <p class="meta" data-meta-name="关键词"><span data-meta-name="关键词.标题" class="key">[关键词]</span>磁共振成像系统；SMR170体模；质量控制检测；线性度；自动测量</p>            
  <div id="article_content">
   <h2 class="sect2">引言</h2>
   <p class="content">为了确保磁共振成像（Magnetic Resonance Imaging，MRI）设备参数准确、可靠及良好的运行状态，就需要定期对其进行质量控制检测<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript" tag="1-4">[1-4]</span>。MRI系统构成复杂，影响MRI影像质量的因素很多，对整个MRI系统进行全面测试十分困难，所以用户仅进行常规参数进行检测<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript" tag="5-6">[5-6]</span>。线性度是MRI系统质控检测必检的参数之一，又称为几何畸变，是描述MRI系统所产生的图像变形程度的参数，体现了MRI重现物体几何尺寸的能力。图像的线性度不好，即所得图像有几何扭曲，就不能真实反映成像物体的几何结构<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript" tag="7-8">[7-8]</span>。</p>
   <p class="content">线性度的检测需要经历体模扫描和参数评估两个阶段，采用体模进行质量控制检测时，真正费时费力的不是体模扫描本身，而是后续的线性度参数计算评估。线性度参数计算与评估所需时间会受检测操作人员经验程度的影响<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript" tag="9-11">[9-11]</span>，结果也会随检测者不同因人而异，部分受到主观影响。为了最大程度提高效率，同时去除不同检测人员带来的主观影响，本研究基于SMR170体模设计了自动测量线性度系统，以期实现对MRI质量控制的线性度参数进行自动、客观的评价。</p>
   <h2 class="sect2">1 材料与方法</h2>
   <h3 title="pagenumber_ebook=19,pagenumber_book=12" class="sect3">1.1 检测体模及受检设备</h3>
   <p class="content">采用Magphan SMR 170性能测试体模，该体模由美国体模实验室研制，内部加注硫酸铜溶液。受检设备为Siemens Trio 3.0 T MRI（德国西门子公司制造），2008年10月投入使用。</p>
   <h3 title="pagenumber_ebook=20,pagenumber_book=13" class="sect3">1.2 体模扫描方法及参数</h3>
   <p class="content">将体模稳定置于安装好的头线圈内，并使用水平仪检查左右、前后水平，将激光灯对准体模中心后送入磁体等中心位置，静置5 min后进行扫描。首先进行三平面定位像扫描，然后在所得的三平面定位像上确定对体模轴位图像的扫描位置层面，从上到下依次为1～5层（图1）。扫描采用自旋回波成像脉冲序列 （Spin-Echo，SE），TR=500 ms，TE=30 ms，FOV=24 cm×24 cm，Matrix=256×256，激励次数为 1次，单层扫描层厚为5 mm，接收带宽为156 Hz/pixel。</p>
   <div class="picture_figure">
    <img title="pagenumber_ebook=20,pagenumber_book=13" class="picture_figure" id="image_20_305_960_1139_1254" width="305" height="107" src="images/18ed8400cbd9b17de3dda2d17b13c7cc.jpg">
    <p class="picture_figure_title" id="image_20_305_960_1139_1254">图1 在三平面定位像上确定对体模轴位图像的扫描位置和层数示意图</p>
   </div>
   <h3 title="pagenumber_ebook=20,pagenumber_book=13" class="sect3">1.3 线性度计算方法</h3>
   <p class="content">MRI图像的线性度要从两个方向进行测算，即频率编码方向和相位编码方向，要求测量多组数据进行平均得到图像线性度的测量结果。一般用畸变百分率表示图像线性的好坏，这是因为畸变越大图像线性越差，畸变越小图像线性越好。畸变百分率能较好地反映图像线性度的好坏。畸变百分率（线性度）的定义为：<img title="pagenumber_ebook=20,pagenumber_book=13" class="picture_line_figure" width="121" height="31" src="images/ae82309c20098966d8354a6249675d77.jpg">，式中，LR为实际距离，单位mm，LM为测量距离，单位mm。轴位图像第三层和第四层，采用测距功能测量图像中的小孔间的间距（主要是交叉方向、横纵方向）（图2），并与实际的距离比较，用上述公式进行计算。</p>
   <div class="picture_figure">
    <img title="pagenumber_ebook=20,pagenumber_book=13" class="picture_figure" id="image_20_226_2086_1184_2532" width="351" height="163" src="images/3dc3339e24c9062a4324cdda627afa2c.jpg">
    <p class="picture_figure_title" id="image_20_226_2086_1184_2532">图2 SMR170体模空间线性度测量示意图</p>
    <p class="picture_figure_note" id="image_20_226_2086_1184_2532">注：a.交叉方向；b.横纵方向。</p>
   </div>
   <h3 title="pagenumber_ebook=20,pagenumber_book=13" class="sect3">1.4 线性度自动测量软件实现方法</h3>
   <p class="content">本自动测量系统基于Matlab 2012b平台进行编写，可对DICOM格式的文件进行解析与处理，能够详细的获取扫描参数。通过Matlab图像处理，自动测算出线性度参数，并输出结果。系统框图如图3所示。</p>
   <p class="content">图像预处理：结合SMR170常规检测方法，发现要计算线性度等参数，必须先对图像中特征区域进行提取。通过对图像进行分析，发现体模中央区域除体模四个支撑柱处为低信号外，整体信号强度远远高于周围背景区域，可利用自动获取阈值的方法进行图像二值化。选取二值化图像中心点作为种子点进行区域生长，并用合适的参数进行膨胀腐蚀，使得图形边缘平滑规则，可以得到图像中心的高信号方形区域，然后采用Hough直线检测，通过Matlab直线检测相关函数Houghlines，可以得到检测出的四条直线的极坐标表示以及每条直线经过的两点，根据这些参数，求得直线在图像笛卡尔坐标系中的斜率K，作为后续部分相关参数的提取的基础。同时，由于SMR170体模为圆柱体，在摆放时可能会有旋转角度偏差，通过K值可以对允许偏差范围内的图像进行校准，而不必放弃图像重新扫描。</p>
   <div class="picture_figure">
    <img title="pagenumber_ebook=20,pagenumber_book=13" class="picture_figure" id="image_20_1548_1012_1987_1604" width="160" height="217" src="images/d301d8420e8457a3e38352b9d7bc0769.jpg">
    <p class="picture_figure_title" id="image_20_1548_1012_1987_1604">图3 线性度自动测算系统框图</p>
   </div>
   <p class="content">SMR170体模轴位图像第三层和第四层图像均可用于线性度的检测，但由于第四层中还存在空间分辨率的测量部件，图像成分较为复杂，不易对线性度检测部件单独提取，因此选取第三层图像进行检测。将图像视为一个大的矩阵，只需确定线性度测量部件的各小圆点在矩阵中的坐标位置，即可求出各点间的两两间距，与间距的相应标准值作比较并代入公式，即可求得线性度，实现流程如图4所示。</p>
   <div class="picture_figure">
    <img title="pagenumber_ebook=20,pagenumber_book=13" class="picture_figure" id="image_20_1588_2122_1961_2781" width="136" height="241" src="images/9305612f26739489d884449f5c23b69b.jpg">
    <p class="picture_figure_title" id="image_20_1588_2122_1961_2781">图4 线性度测量实现流程</p>
   </div>
   <p class="content">根据图像中央方形特征区域相关参数，可以得到其中心坐标，从而得到如图5所示的横纵两条直线，沿着两条直线提取二值图中对应位置的行与列，并分别对其求一阶导，如图6所示的位置特征图。</p>
   <div class="picture_figure">
    <img title="pagenumber_ebook=21,pagenumber_book=14" class="picture_figure" id="image_21_459_279_1002_823" width="199" height="199" src="images/7b3a09fb6ddeba2d6a892bf53fb8cb85.jpg">
    <p class="picture_figure_title" id="image_21_459_279_1002_823">图5 第三层图像二值图</p>
   </div>
   <div class="picture_figure">
    <img title="pagenumber_ebook=21,pagenumber_book=14" class="picture_figure" id="image_21_223_895_1247_1299" width="375" height="148" src="images/ea189a757b502693d168f1c7f42b64e4.jpg">
    <p class="picture_figure_title" id="image_21_223_895_1247_1299">图6 位置特征图</p>
   </div>
   <p class="content">原二值图上高低信号变化时在特征图中就会出现1或-1，则通过一组1和-1，求其位置中点，即可确定原图中一个小圆点的位置。</p>
   <p class="content">原图中纵向为x方向，横向为y方向。以横向为例，特征图中从左至右、从右至左2～4组1和-1分别可以确定横向特征线穿过的6各点的y坐标，其x坐标由特征线本身确定，然后分别以左侧和右侧第二个点的y坐标作纵向特征线，类似地求出上下两点的x坐标，其他小圆点坐标也可以用类似方法求出，得到小圆点坐标如图7所示。</p>
   <div class="picture_figure">
    <img title="pagenumber_ebook=21,pagenumber_book=14" class="picture_figure" id="image_21_483_1923_981_2420" width="182" height="182" src="images/cb06fa71b7a6c314de04e2d029d483ab.jpg">
    <p class="picture_figure_title" id="image_21_483_1923_981_2420">图7 小圆点提取位置示意图</p>
   </div>
   <p class="content">在得到所有圆点坐标后，根据体模使用说明，求出标准距离为2、8、10、12 cm的圆点间的图上测量距离，代入公式即可求得线性度。</p>
   <h3 title="pagenumber_ebook=21,pagenumber_book=14" class="sect3">1.5 统计学处理</h3>
   <p class="content">采用Bland-Altman统计学方法对手工测量和自动测量两种方法的一致性进行统计分析，以测量结果差值为纵轴，以两种方法计算结果的均值为横轴，绘制出散点坐标图，并标出95%一致性界限，通常当图中位于区间内的点占到95%以上且不超出专业上可接受的临界值范围时，则认为两种检测方法的一致性满足相互替代的要求<span title="pagenumber_ebook=21,pagenumber_book=14" class="superscript" tag="12-13">[12-13]</span>。</p>
   <h2 class="sect2">2 结果</h2>
   <p class="content">以某次基于SMR170体模的MR设备质量控制得到的第4层图像作为检测对象，根据图2a所示对体模2、8、10、12 cm的圆点间的距离进行测量，采用手工测量和自动测量两种方法对线性度进行测算，每个距离测量四次，并绘制Bland-Altman图，见图8。手工测量和自动测量两种方法对线性度测量得到的数据均位于一致性界限范围内（-0.42%，0.68%），两种方法的一致性较好，可相互替代。</p>
   <div class="picture_figure">
    <img title="pagenumber_ebook=21,pagenumber_book=14" class="picture_figure" id="image_21_1334_779_2196_1386" width="316" height="222" src="images/f1eb1f6ff3584d516d94db4da908b23d.jpg">
    <p class="picture_figure_title" id="image_21_1334_779_2196_1386">图8 Bland-Altman图</p>
   </div>
   <h2 class="sect2">3 讨论</h2>
   <p class="content">国内关于MRI设备质量控制检测规范中所采用的体模多数为SMR170体模，采用该体模可实现信噪比、均匀度、层厚、线性度、空间分辨率和低对比度分辨率等重要参数的检测与评估，线性度是不可或缺的MRI系统质控检测必检参数之一<span title="pagenumber_ebook=21,pagenumber_book=14" class="superscript" tag="14-16">[14-16]</span>。</p>
   <p class="content">目前线性度的测算多采用手工测算，这种方法的缺点是效率低下，且易受操作人员测量习惯及主观影响，如果一直是同一人进行评估，其连续评估结果具有较好的连续性与一致性，有利于多次检测结果对比。如果每次测量均不是同一人，则测量结果有可能出入较大，不利于前后对比。为了有效解决此问题，本研究在Matlab平台上，设计了基于SMR170体模MRI质量控制检测的线性度自动测量程序，实现了对MRI质量控制的线性度参数进行自动测量，大幅缩短检测结果分析时间，提高了工作效率，同时对结果的评判更加客观准确，减少了不同检测人员的主观差异。</p>
   <p class="content">在系统编写完成后，对某次基于SMR170体模的MR设备质量控制检测得到的图像进行了测试，分别采用手工测量和自动测量对线性度进行测量，并采用Bland-Altman统计学方法对手工测量和自动测量两种方法的一致性进行统计分析，发现采用两种方法对线性度测量得到的数据均位于一致性界限范围内，两种方法的一致性较好，可相互替代。</p>
   <p class="content">在程序运行中发现，如果线性度测量层面的图像有伪影，往往会造成图像识别出现错误，同时导致了错误的测量结果，因此要想确保结果准确，必须在扫描体模的过程中做好质量保证，确保测量线性度层面图像没有伪影。同时，本研究也存在一些需要改进和提高的方面，在确保质量控制图像本身的准确性方面，本研究采用轴向定位像中的局部特征元素进行比对来判定图像定位是否准确，然后在检测图像通过中央特征区域参数提取，来判断体模是否有旋转偏差。这种判定方法存在一定局限性，应当对体模的多个层面扫描后分析图像，得到更多的评判条件以确保图像定位的准确性。</p>
   <p class="content">另外，无论是手工测算还是自动测量，由于缺乏“金标准”，因此关于测量的精度仍存在一些争议，为此，一方面需要在测量过程中做好质控，另一方面需要通过长时间定期检测建立设备运行基线，以基线及检测规范为参考标准，进一步验证测量的精确度。在大多数情况下，通过图像预处理，可以获得图像中期望的图像点或图像像素。但是由于图像数据质量较差或者检测程序存在缺陷等原因，可能存在点或像素缺失等情况，此时可以考虑采用其他算法，如基于标准Hough变换的圆检测算法或者随机Hough变换，从而保证参数检测准确性<span title="pagenumber_ebook=22,pagenumber_book=15" class="superscript" tag="17">[17]</span>。最后，系统还未进行大批量多次的性能测试，未来应该扩大测试范围，进一步确定系统的可靠性。</p>
   <p title="pagenumber_ebook=22,pagenumber_book=15" class="meta" data-meta-name="参考文献"><span data-meta-name="参考文献.标题" class="key">[参考文献]</span></p>
   <p class="content" name="ref1" data-meta-name="参考文献.条目">[1]杨正汉,王振常.重视MRI设备临床质量控制——发挥设备最优性能[J].中国医学影像技术,2017,33(11):1605-1606.</p>
   <p class="content" name="ref2" data-meta-name="参考文献.条目">[2]尹红霞,杨娉娉,刘雅文,等.MRI设备中心频率和发射增益稳定性检测及处置界限建立[J].中国医学影像技术,2017,33(8):1260-126.</p>
   <p class="content" name="ref3" data-meta-name="参考文献.条目">[3]尹红霞,刘雅文,杨娉娉,等.正交鸟笼头线圈质量控制检测与建立处置界限的初步实践[J].中国医学影像技术,2018,34(1):123-127.</p>
   <p class="content" name="ref4" data-meta-name="参考文献.条目">[4]刘雅文,尹红霞,杨娉娉,等.MRI设备临床质量控制流程和制度的现状与发展趋势[J].中国医学影像技术,2017,33(11):1615-1619.</p>
   <p class="content" name="ref5" data-meta-name="参考文献.条目">[5]刘隺是,李彦,严福华,等.垂体磁共振临床图像质量评价研究[J].中国医疗设备,2019,34(7):34-40.</p>
   <p class="content" name="ref6" data-meta-name="参考文献.条目">[6]邓洁.核磁共振成像系统的质量保证与质量控制[J].中国医疗设备,2011,26(7):13-20.</p>
   <p class="content" name="ref7" data-meta-name="参考文献.条目">[7]倪萍,孙钢.医用磁共振成像设备质量控制检测技术[M].北京:中国质检出版社,中国标准出版社,2016:9.</p>
   <p class="content" name="ref8" data-meta-name="参考文献.条目">[8]储呈晨,王龙辰,毕帆,等.磁共振图像质量控制中的若干评价指标探讨[J].中国医疗设备,2016,31(7):124-127.</p>
   <p class="content" name="ref9" data-meta-name="参考文献.条目">[9]熊晖,陈自谦,钟群,等.磁共振质量控制参数中均匀度检测方法的对比研究[J].中国医疗设备,2018,33(10):15-18.</p>
   <p class="content" name="ref10" data-meta-name="参考文献.条目">[10]林迪逵.磁共振质量控制的智能评价系统研究[D].福州:福州大学,2017.</p>
   <p class="content" name="ref11" data-meta-name="参考文献.条目">[11]熊晖.基于ACR体模的医用磁共振质量控制检测智能评价系统研究[D].福州:福州大学,2018.</p>
   <p class="content" name="ref12" data-meta-name="参考文献.条目">[12]闻浩,陆梦洁,刘玉秀,等.定量测量Bland-Altman一致性评价方法研究及临床应用[J].医学研究生学报,2015,28(10):1107-1111.</p>
   <p class="content" name="ref13" data-meta-name="参考文献.条目">[13]萨建,刘桂芬.定量测量结果的一致性评价及Bland-Altman法的应用[J].中国卫生统计,2011,28(4):409-411.</p>
   <p class="content" name="ref14" data-meta-name="参考文献.条目">[14]汪立文,刘子龙,代鲲鹏,等.医用磁共振模体空间分辨力测试模块校准方法研究[J].计量与测试技术,2018,(8):59-61.</p>
   <p class="content" name="ref15" data-meta-name="参考文献.条目">[15]玉碧坚,苏毅.医用磁共振成像系统(MRI)成像原理及其计量参数检测的实验研究[J].中国科技纵横,2016,(4):194-196.</p>
   <p class="content" name="ref16" data-meta-name="参考文献.条目">[16]万宇.医用磁共振成像(MRI)系统的计量检测[J].计量与测试技术,2018,45(11):52-54.</p>
   <p class="content" name="ref17" data-meta-name="参考文献.条目">[17]林迪逵,付丽媛,陈自谦,等.基于随机Hough变换圆检测算法的MRI QC空间线性度自动检测程序设计[J].中国医疗设备,2018,33(10):19-22.</p>
  </div> 
  <form>
   <input type="hidden" id="html_version" value="3">
  </form>  
  <div id="en">
   <p class="title_translate" data-meta-name="译文篇名">Automatic Measurement of Linearity in MRI Quality Control Test Based on SMR170 Phantom</p>
   <p class="author_info_translate" data-meta-name="译文作者信息">LIANG Yonggang<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,FU Liyuan<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,CHEN Jian<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,HUANG Zhifeng<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,XIONG Hui<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,LI Liangjie<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,ZHONG Qun<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,XIAO Hui<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,XU Shangwen<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1</span>,CHEN Ziqian<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">1,2</span><br>1.Department of Radiology,900<span title="pagenumber_ebook=19,pagenumber_book=12" class="superscript">th</span>Hospital of Joint Logistics Support Force(Formerly Fuzhou General Hospital of Nanjing Military Region),Fuzhou Fujian 350025,China;2.Department of Radiology,The Affiliated Dongfang Hospital of Xiamen University,Fuzhou Fujian 350025,China</p>
   <p class="meta" data-meta-name="译文摘要"><span data-meta-name="译文摘要.标题" class="key">Abstract:</span>Objective To design an automatic linearity measurement program based on the quality control test of the SMR170 phantom MRI system in order to realize the automatic measurement of linearity parameters.Methods The image processing algorithm was designed and written on the Matlab platform,and the linearity parameters were automatically measured.The consistency between manual measurement method and automatic measurement method was analyzed by Bland-Altman statistical method.Results The automatic measurement program based on Matlab platform could measure linearity automatically,which had good consistency with manual measurement method.Conclusion The automatic linearity detection program and manual measurement can be substituted for each other.</p>
   <p class="meta" data-meta-name="译文关键词"><span data-meta-name="译文关键词.标题" class="key">Key words:</span>magnetic resonance imaging system;SMR170 phantom;quality control test;linearity;automatic measurement</p>
  </div>
  <div id="other">
   <p class="meta" data-meta-name="中图分类号"><span data-meta-name="中图分类号.标题" class="key">[中图分类号]</span>R445.2</p>
   <p class="meta" data-meta-name="doi"><span data-meta-name="doi.标题" class="key">doi：</span>11.3969/j.issn.1674-1633.2020.11.004</p>
   <p class="meta" data-meta-name="文献标识码"><span data-meta-name="文献标识码.标题" class="key">[文献标识码]</span>A</p>
   <p class="meta" data-meta-name="文章编号"><span data-meta-name="文章编号.标题" class="key">[文章编号]</span>1674-1633(2020)11-0012-04</p>
   <p class="meta" data-meta-name="收稿日期"><span data-meta-name="收稿日期.标题" class="key">收稿日期：</span>2020-08-27</p>
   <p class="meta" data-meta-name="基金项目"><span data-meta-name="基金项目.标题" class="key">基金项目：</span>国家重点研发计划数字诊疗装备研发重点专项子课题（2016 YFC0103103）；军区重大课题（14ZX23）；医院科研计划（2019Q02）。</p>
   <p class="meta" data-meta-name="通信作者"><span data-meta-name="通信作者.标题" class="key">通信作者：</span>陈自谦，主任医师，教授，博士生导师，主要研究方向为影像医学与核医学。</p>
   <p class="meta" data-meta-name="通信作者"><span data-meta-name="通信作者.标题" class="key">通信作者</span>邮箱：chenziqianfz@sina.com</p>
   <p class="author_info_right" data-meta-name="作者信息">本文编辑 皮志超</p>
  </div>
      <script type="text/javascript">
      var css={
       cssBox:'position: fixed; left: 0; top: 0; width: 98%; height: 98%; padding: 1%; background: rgba(0,0,0,0.6);',
       cssBody:'display: -webkit-box; -webkit-box-orient:vertical; -webkit-box-pack: center; -webkit-box-align:center; width: 100%; height: 100%;',
       cssImg:'max-width: 100%; max-height: 100%;'
      }
      var body=document.querySelector('body');
   function ImgShow(body1,css1){
    this.show1=function(){
     var imgBox=document.createElement("div");
        var imgBody=document.createElement("div");
        var newImg=document.createElement("img");
        imgBox.style.cssText=arguments[2].cssBox;
        imgBody.style.cssText=arguments[2].cssBody;
        newImg.style.cssText=arguments[2].cssImg;
        newImg.setAttribute('src',arguments[1].getAttribute('src'));
        imgBody.appendChild(newImg);
        imgBox.appendChild(imgBody);
        document.querySelector('body').appendChild(imgBox);
        imgBox.onclick=function(){
         this.style.display='none';
         this.remove();
        }
    }
    this.imgList=body1.getElementsByTagName('img');
    for(var i=0;i<this.imgList.length;i++){
     var _self=this;
     this.imgList[i].addEventListener('click',function(event){
      _self.show1(event,this,css1);
     },false);
    }
   }
   var zImg=new ImgShow(body,css);
     </script>
</body>
</html>