搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

用于多通道单分子定位的高精度图像配准方法

林丹樱 龚振权 黄黎琳 聂梦娇 于斌 屈军乐

引用本文:
Citation:

用于多通道单分子定位的高精度图像配准方法

林丹樱, 龚振权, 黄黎琳, 聂梦娇, 于斌, 屈军乐

A high-precision image registration method for multi-channel single molecule localization

Lin Dan-Ying, Gong Zhen-Quan, Huang Li-Lin, Nie Meng-Jiao, Yu Bin, Qu Jun-Le
PDF
HTML
导出引用
  • 单分子定位技术可以绕过光学系统的衍射限制, 在生物样品的单粒子追踪和超分辨显微成像中得到了广泛应用. 多通道单分子定位采用多个成像通道, 可以实现对不同目标的同时追踪或多色超分辨成像, 也可以提升单粒子追踪的轴向深度或实现更高的定位精度和密度. 但各通道图像间的差异会影响协同定位或定量分析, 因此图像配准是其图像数据预处理的关键环节; 且由于单分子定位精度高, 其对多通道图像配准精度的要求也很高. 现有技术一般采用基于控制点的配准方法, 且多采用复杂而精密的方式来获取基准物网格图像用于定位得到控制点对, 以实现高精度图像配准, 对样品或实验设备要求高, 难以直接推广. 为此, 本文基于局部非线性变换和误匹配点剔除, 发展了一种可以直接采用随机分布荧光珠样品作为基准物的高精度图像配准方法, 通过在特征匹配和变换模型参数估计的过程中对控制点进行监测和迭代筛选, 以剔除因单分子定位不准确或精度差而导致未精确匹配的控制点对, 从而消除以随机分布荧光珠样品作为基准物时对于控制点准确获取和精确匹配所带来的不良影响, 同时采用基于局部加权平均的二阶多项式拟合进行变换模型参数估计, 以更好地适用于不同通道间存在局部非线性形变的情形. 结果表明, 采用该方法只需要3次迭代, 就可以将未准确定位和精确匹配的控制点对找到并剔除, 从而实现更准确的变换模型参数估计, 将配准精度提高一个数量级, 在图像局部非线性形变情况严重的正交像散双通道单分子定位成像系统中实现了约6 nm的配准精度.
    Single-molecule localization technology has been widely used in single-particle tracking and super-resolution imaging of biological samples, as it can bypass the diffraction limit of optical systems. Multi-channel single-molecule localization uses multiple imaging channels to simultaneously track different targets or perform multi-color super-resolution imaging, and can also improve the axial depth of single-particle tracking or achieve higher localization precision and density for super-resolution imaging. However, the difference between images in each channel can affect collaborative localization or quantitative analysis, so image registration is a key step in its image data preprocessing. Moreover, due to the high precision of single-molecule localization, its requirements for multi-channel image registration accuracy are also high. Existing technologies generally use control point-based registration methods and often use complicated and precise methods to obtain fiducial images for locating control point pairs to achieve high-precision image registration, which involves high sample or experimental equipment requirements and is difficult to directly extend to other systems. Therefore, developed in this work, is a high-precision image registration method that can directly use randomly distributed fluorescent beads as fiducial samples based on local nonlinear transformation and elimination of mismatched points. By monitoring and iteratively filtering control points in the process of feature matching and transformation model parameter estimation to eliminate control point pairs that are not accurately matched due to inaccurate or poor precision of single-molecule localization, the adverse effects on accurate acquisition and precise matching of control points when using randomly distributed fluorescent beads as fiducial samples are eliminated. At the same time, a second-order polynomial fitting based on local weighted mean is used for estimating the transformation model parameter to better adapt to the existence of local nonlinear deformation between different channels. The results show that using this method only requires three iterations to find and eliminate control point pairs that are not accurately located and matched, thereby achieving more accurate transformation model parameter and improving the registration accuracy by an order of magnitude, achieving a registration accuracy of about 6 nm in a complex dual-channel single-molecule localization imaging system based on orthogonal astigmatism.
      通信作者: 林丹樱, dylin@szu.edu.cn ; 屈军乐, jlqu@szu.edu.cn
    • 基金项目: 国家重点研发计划(批准号: 2022YFF0706001)、国家自然科学基金(批准号: 62275165, 62235007, 61975131, 62175166, 61835009, 62127819)、深圳市光子学与生物光子学重点实验室(批准号: ZDSYS20210623092006020)和深圳市基础研究项目(批准号: JCYJ20200109105411133, JCYJ20220818100202005)资助的课题.
      Corresponding author: Lin Dan-Ying, dylin@szu.edu.cn ; Qu Jun-Le, jlqu@szu.edu.cn
    • Funds: Project supported by the National Key R&D Program of China (Grant No. 2022YFF0706001), the National Natural Science Foundation of China (Grant Nos. 62275165, 62235007, 61975131, 62175166, 61835009, 62127819), the Shenzhen Key Laboratory of Photonics and Biophotonics, China (Grant No. ZDSYS20210623092006020), and the Shenzhen Basic Research Project, China (Grant Nos. JCYJ20200109105411133, JCYJ20220818100202005).
    [1]

    Diezmann A V, Shechtman Y, Moerner W E 2017 Chem. Rev. 117 7244Google Scholar

    [2]

    Saxton M J, Jacobson K 1997 Annu. Rev. Biophys. Biomol. Struct. 26 373Google Scholar

    [3]

    Deich J, Judd E M, Mcadams H H, Moerner W E 2004 Proc. Natl. Acad. Sci. U. S. A. 101 15921Google Scholar

    [4]

    Rust M J, Bates M, Zhuang X 2006 Nat. Methods 3 793Google Scholar

    [5]

    Sigal Y M, Zhou R, Zhuang X 2018 Science 361 880Google Scholar

    [6]

    Driouchi A, Gray-Owen S D, Yip C M 2022 J. Biol. Chem. 298 102448Google Scholar

    [7]

    Louis B, Camacho R, Bresolí-Obach R, Abakumov S, Vandaele J, Kudo T, Masuhara H, Scheblykin I, Hofkens J, Rocha S 2020 Opt. Express 28 28656Google Scholar

    [8]

    Albrecht D, Winterflood C M, Ewers H 2015 Methods Appl. Fluoresc. 3 024001Google Scholar

    [9]

    Bates M, Dempsey G T, Chen K H, Zhuang X 2012 ChemPhysChem 13 99Google Scholar

    [10]

    Lehmann M, Gottschalk B, Puchkov D, Schmieder P, Schwagerus S, Hackenberger P, Haucke V, Schmoranzer J 2015 Angew. Chem. Int. Ed. 54 13230Google Scholar

    [11]

    Gu L, Sheng Y, Chen Y, Chang H, Zhang Y, Lv P, Ji W, Xu T 2014 Biophys. J. 106 2443Google Scholar

    [12]

    Min J, Holden S J, Carlini L, Unser M, Manley S, Ye J C 2014 Biomed. Opt. Express 5 3935Google Scholar

    [13]

    林丹樱, 武泽凯, 于斌, 黄黎琳, 张潇, 屈军乐 2022 物理学报 71 128701Google Scholar

    Lin D Y, Wu Z K, Yu B, Huang L L, Zhang X, Qu J L 2022 Acta Phys. Sin. 71 128701Google Scholar

    [14]

    Deschout H, Shivanandan A, Annibale P, Scarselli M, Radenovic A 2014 Histochem. Cell Biol. 142 5Google Scholar

    [15]

    Paul S, Pati U C 2021 Int. J. Remote Sens. 42 5396Google Scholar

    [16]

    Churchman L S, Okten Z, Rock R S, Dawson J F, Spudich J A 2005 Proc. Natl. Acad. Sci. U. S. A. 102 1419Google Scholar

    [17]

    Gahlmann A, Ptacin J L, Grover G, Quirin S, von Diezmann A R S, Lee M K, Backlund M P, Shapiro L, Piestun R, Moerner W E 2013 Nano Lett. 13 987Google Scholar

    [18]

    Huang B, Wang W Q, Bates M, Zhuang X 2008 Science 319 810Google Scholar

    [19]

    Goshtasby A 1988 Image Vison Comput. 6 255Google Scholar

    [20]

    Zagorchev L, Goshtasby A 2006 IEEE Trans. Image Process. 15 529Google Scholar

    [21]

    Huang F, Sirinakis G, Allgeyer E S, Toomre D, Booth M J, Bewersdorf J 2016 Cell 166 1028Google Scholar

    [22]

    Fischler M A, Bolles R C 1981 Commun. ACM 24 381Google Scholar

    [23]

    张岩, 孙世宇, 胡永江, 李建增, 范聪 2018 电子与信息学报 40 928Google Scholar

    Zhang Y, Sun S Y, Hu Y J, Li J Z, Fan C 2018 J. Electron. Inf. Technol. 40 928Google Scholar

    [24]

    赖焕杰, 孟祥印, 肖世德, 胡锴沣, 李召鑫 2023 传感器与微系统 42 135

    Lai H J, Meng X Y, Xiao S D, Hu K F, Li Z X 2023 Transducer Microsys. Technol. 42 135

  • 图 1  正交像散单分子定位及其图像配准需求示意图 (a)正交像散系统光路; (b)不同轴向位置物点对应的像散图像对; (c)正交像散导致的双通道图像形变及预期配准效果. ND, 衰减片; L, 透镜; M, 反射镜; DM, 二向色镜; OL, 物镜; F, 滤光片; A, 光阑; BS, 分束器; CL, 柱透镜; EMCCD, 电子倍增电荷耦合器件

    Fig. 1.  Illustration of the orthogonal astigmatism-based single molecule localization method and its requirements on image registration: (a) Optical path of the orthogonal astigmatism system; (b) astigmatic image pairs corresponding to object points at different axial positions; (c) dual-channel image distortion caused by orthogonal astigmatism and expected registration result. ND, attenuator; L, lens; M, mirror; DM, dichroic mirror; OL, objective lens; F, filter; A, aperture; BS, beam splitter; CL, cylindrical lens; EMCCD, electron-multiplying charge-coupled device.

    图 2  基于基准配准误差监测的误匹配点剔除(EMP-FRE)算法的基本原理和算法流程示意图 (a)基本原理; (b)算法流程图. I, 图像; P, 定位点; T, 配准参数矩阵; FRE, 基准配准误差; E, 误差阈值; 下标l/r, 左/右通道; 下标i, 控制点编号

    Fig. 2.  Schematic diagram (a) and algorithm flowchart (b) of the EMP-FRE (i.e., elimination of mismatched points based on fiducial registration error monitoring) algorithm. I, image; P, localizations; T, registration parameter matrix; FRE, fiducial registration error; E, error threshold; subscript l/r, left/right channel; subscript i, control point number.

    图 3  变换模型选择和误匹配点剔除对配准效果的影响 (a)左通道参考图像; (b)右通道待配准图像; (c) LWM法配准后的右通道图像; (d)采用不同模型配准后的定位点; (e)采用不同模型配准后的平均FRE; (f), (g)采用LWM模型和EMP-FRE算法的X, Y方向配准误差统计. 统计数据来自5帧不同视野荧光珠图像的140对定位点(控制点)

    Fig. 3.  Effects of transformation model selection and mismatched point elimination on registration performance: (a) Left channel, reference image; (b) right channel, image to be registered; (c) right channel, image after LWM registration; (d) localizations registered using different models; (e) average FRE after registration using different models; (f), (g) statistical data on registration errors in the X and Y directions using the LWM model and EMP-FRE algorithm, obtained from 140 pairs of localizations (control points) from 5 different field of view fluorescent bead images.

    图 4  实验荧光珠图像的配准结果对比 (a)传统LWM的配准结果; (b)结合LWM和EMP-FRE的配准结果; (c), (d)蓝色方框的放大图; (e), (f)红色方框的放大图 伪彩色: 紫色代表左通道图像, 绿色代表配准后的右通道图像, 两者准确重合时为灰白色

    Fig. 4.  Comparison of the registration results of experimental fluorescent bead images: (a) Registered image using traditional LWM; (b) registered image using LWM+EMP-FRE; (c), (d) enlarged views of the blue box; (e), (f) enlarged views of red box. Pseudo color: purple represents the left channel image, green represents the registered right channel image, and when they are accurately overlapped, the color should be grayish white.

    图 5  平行线模拟样品的超分辨重构结果比较 (a)不同情形数据的重构图像: 1-无误差数据; 2-有误差未配准数据; 3-AT配准数据; 4-2PF配准数据; 5-PWL配准数据; 6-LWM配准数据; 7-LWM+EMP-FRE配准数据; (b)—(d)图(a)中3个白色方框处对应的截面轮廓曲线

    Fig. 5.  Comparison of super-resolution reconstruction results for the simulated parallel line sample: (a) Reconstructed images from data in different cases: 1-error free data; 2-unregistered data with errors; 3-AT registered data; 4-2PF registered data; 5-PWL registered data; 6-LWM registered data; 7-LWM+EMP-FRE registered data; (b)–(d) the corresponding cross-section profiles of the three white boxes in panel (a).

    图 6  微管实验数据的超分辨重构结果比较 (a)宽场荧光图像; (b), (c)图(a)中白色方框区域的数据分别采用LWM法和LWM+EMP-FRE法配准后的重构结果; (d)图(b)和(c)中红色和蓝色方框处对应的截面轮廓曲线

    Fig. 6.  Comparison of super-resolution reconstruction results for microtubule experimental data: (a) Wide-field fluorescence image; (b), (c) reconstructed images from data of the white box region in panel (a), registered using LWM and LWM+EMP-FRE, respectively; (d) corresponding cross-section profiles of the red and blue boxes in panel (b) and (c), respectively.

    表 1  图3(a)图3(b)荧光珠图像的EMP-FRE迭代过程及利用剩余控制点对进行LWM配准所得的平均基准配准误差

    Table 1.  The EMP-FRE iteration process of the fluorescent bead images in Fig. 3(a), (b), and the average FRE obtained through LWM registration using residual control points.

    迭代次数剔除控制点对标号平均基准配准误差/ nm
    029.3
    13, 4, 6, 711.7
    21, 2, 3, 4, 53.4
    32, 3, 4, 52.9
    下载: 导出CSV

    表 2  误差阈值对迭代和平均配准误差的影响

    Table 2.  Effect of error threshold on iterative and average registration errors.

    误差阈值/nm 每次迭代剔除的控制点对数量 每次迭代后的平均配准误差/nm
    1 2 3 4 1 2 3 4
    10 4 16 10 9 11.7 5.8 5.6 5.4
    20 4 9 4 11.7 3.5 2.9
    30 4 5 4 11.7 3.4 2.9
    40 4 3 11.7 10.4
    50 4 2 11.7 11.6
    下载: 导出CSV

    表 3  误匹配点剔除算法的效果及效率对比

    Table 3.  Comparison of the effectiveness and efficiency of algorithms for eliminating mismatched points.

    误匹配点剔除算法迭代成功率迭代成功所需的
    平均迭代次数
    RANSAC0.84165
    RSCFDI0.49217
    RASCFDI+
    RANSAC2 match
    0.49260
    下载: 导出CSV

    表 4  平行线模拟样品重构图像的线间距和对比度定量比较

    Table 4.  Quantitative comparison of line spacing and contrast in reconstructed images of the simulated parallel line sample.

    线间距/nm 对比度/%
    b c d 平均 b c d 平均
    Control 58.6 59.4 57.2 58.4 82.1 87.2 76.2 81.8
    LWM 55.2 58.8 55.8 56.6 77.7 83.3 48.5 69.8
    LWM+EMP-FRE 57.2 58.4 59.0 58.2 74.7 79.4 82.1 78.7
    下载: 导出CSV
  • [1]

    Diezmann A V, Shechtman Y, Moerner W E 2017 Chem. Rev. 117 7244Google Scholar

    [2]

    Saxton M J, Jacobson K 1997 Annu. Rev. Biophys. Biomol. Struct. 26 373Google Scholar

    [3]

    Deich J, Judd E M, Mcadams H H, Moerner W E 2004 Proc. Natl. Acad. Sci. U. S. A. 101 15921Google Scholar

    [4]

    Rust M J, Bates M, Zhuang X 2006 Nat. Methods 3 793Google Scholar

    [5]

    Sigal Y M, Zhou R, Zhuang X 2018 Science 361 880Google Scholar

    [6]

    Driouchi A, Gray-Owen S D, Yip C M 2022 J. Biol. Chem. 298 102448Google Scholar

    [7]

    Louis B, Camacho R, Bresolí-Obach R, Abakumov S, Vandaele J, Kudo T, Masuhara H, Scheblykin I, Hofkens J, Rocha S 2020 Opt. Express 28 28656Google Scholar

    [8]

    Albrecht D, Winterflood C M, Ewers H 2015 Methods Appl. Fluoresc. 3 024001Google Scholar

    [9]

    Bates M, Dempsey G T, Chen K H, Zhuang X 2012 ChemPhysChem 13 99Google Scholar

    [10]

    Lehmann M, Gottschalk B, Puchkov D, Schmieder P, Schwagerus S, Hackenberger P, Haucke V, Schmoranzer J 2015 Angew. Chem. Int. Ed. 54 13230Google Scholar

    [11]

    Gu L, Sheng Y, Chen Y, Chang H, Zhang Y, Lv P, Ji W, Xu T 2014 Biophys. J. 106 2443Google Scholar

    [12]

    Min J, Holden S J, Carlini L, Unser M, Manley S, Ye J C 2014 Biomed. Opt. Express 5 3935Google Scholar

    [13]

    林丹樱, 武泽凯, 于斌, 黄黎琳, 张潇, 屈军乐 2022 物理学报 71 128701Google Scholar

    Lin D Y, Wu Z K, Yu B, Huang L L, Zhang X, Qu J L 2022 Acta Phys. Sin. 71 128701Google Scholar

    [14]

    Deschout H, Shivanandan A, Annibale P, Scarselli M, Radenovic A 2014 Histochem. Cell Biol. 142 5Google Scholar

    [15]

    Paul S, Pati U C 2021 Int. J. Remote Sens. 42 5396Google Scholar

    [16]

    Churchman L S, Okten Z, Rock R S, Dawson J F, Spudich J A 2005 Proc. Natl. Acad. Sci. U. S. A. 102 1419Google Scholar

    [17]

    Gahlmann A, Ptacin J L, Grover G, Quirin S, von Diezmann A R S, Lee M K, Backlund M P, Shapiro L, Piestun R, Moerner W E 2013 Nano Lett. 13 987Google Scholar

    [18]

    Huang B, Wang W Q, Bates M, Zhuang X 2008 Science 319 810Google Scholar

    [19]

    Goshtasby A 1988 Image Vison Comput. 6 255Google Scholar

    [20]

    Zagorchev L, Goshtasby A 2006 IEEE Trans. Image Process. 15 529Google Scholar

    [21]

    Huang F, Sirinakis G, Allgeyer E S, Toomre D, Booth M J, Bewersdorf J 2016 Cell 166 1028Google Scholar

    [22]

    Fischler M A, Bolles R C 1981 Commun. ACM 24 381Google Scholar

    [23]

    张岩, 孙世宇, 胡永江, 李建增, 范聪 2018 电子与信息学报 40 928Google Scholar

    Zhang Y, Sun S Y, Hu Y J, Li J Z, Fan C 2018 J. Electron. Inf. Technol. 40 928Google Scholar

    [24]

    赖焕杰, 孟祥印, 肖世德, 胡锴沣, 李召鑫 2023 传感器与微系统 42 135

    Lai H J, Meng X Y, Xiao S D, Hu K F, Li Z X 2023 Transducer Microsys. Technol. 42 135

  • [1] 方靖岳, 文之豪, 朱海碧涛, 李欣幸, 邓联文. 基于集成法布里-珀罗微腔阵列的16通道快照式多光谱成像. 物理学报, 2024, 73(7): 074205. doi: 10.7498/aps.73.20231775
    [2] 王辉林, 廖艳林, 赵艳, 章文, 谌正艮. 基于多激光束驱动准单能高能质子束模拟研究. 物理学报, 2023, 72(18): 184102. doi: 10.7498/aps.72.20230313
    [3] 田晓俊, 孔繁芳, 经士浩, 郁云杰, 张尧, 张杨, 董振超. 单分子瞬时带电态中电子-振动耦合特性的亚纳米荧光成像研究. 物理学报, 2022, 71(6): 063301. doi: 10.7498/aps.71.20212003
    [4] 林丹樱, 武泽凯, 于斌, 黄黎琳, 张潇, 屈军乐. 正交像散高密度三维单分子定位显微的数值模拟. 物理学报, 2022, 71(12): 128701. doi: 10.7498/aps.71.20212091
    [5] 霍勇刚, 严江余, 张全虎. 缪子多模态成像图像质量分析. 物理学报, 2022, 71(2): 021401. doi: 10.7498/aps.71.20211083
    [6] 李四维, 林丹樱, 邹小慧, 张炜, 陈丹妮, 于斌, 屈军乐. 基于双螺旋点扩散函数工程的多焦点图像扫描显微. 物理学报, 2021, 70(3): 038701. doi: 10.7498/aps.70.20200640
    [7] 霍勇刚, 严江余, 张全虎. 缪子多模态成像图像质量分析. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211083
    [8] 苏博, 陶芬, 李可, 杜国浩, 张玲, 李中亮, 邓彪, 谢红兰, 肖体乔. 同步辐射纳米CT图像配准方法研究. 物理学报, 2021, 70(16): 160704. doi: 10.7498/aps.70.20210156
    [9] 王琼, 王凯歌, 孟康康, 孙聃, 韩仝雨, 高爱华. 基于单分子成像技术研究λ-DNA分子穿越微米通道端口的电动力学特性. 物理学报, 2020, 69(16): 168202. doi: 10.7498/aps.69.20200074
    [10] 王殿伟, 韩鹏飞, 范九伦, 刘颖, 许志杰, 王晶. 基于光照-反射成像模型和形态学操作的多谱段图像增强算法. 物理学报, 2018, 67(21): 210701. doi: 10.7498/aps.67.20181288
    [11] 李四维, 吴晶晶, 张赛文, 李恒, 陈丹妮, 于斌, 屈军乐. 用于大景深单分子定位显微的多功能全息相位片的设计及数值模拟. 物理学报, 2018, 67(17): 174202. doi: 10.7498/aps.67.20180569
    [12] 王岩, 王飞, 王挺峰, 谢京江. 基于自适应阈值的阵列激光三维点云配准. 物理学报, 2016, 65(24): 249501. doi: 10.7498/aps.65.249501
    [13] 赵辽英, 吕步云, 厉小润, 陈淑涵. 基于尺度不变特征变换和区域互信息优化的多源遥感图像配准. 物理学报, 2015, 64(12): 124204. doi: 10.7498/aps.64.124204
    [14] 周先春, 汪美玲, 周林锋, 吴琴. 基于Demons算法改进的图像去噪模型研究. 物理学报, 2015, 64(2): 024205. doi: 10.7498/aps.64.024205
    [15] 陈鹤, 于斌, 陈丹妮, 李恒, 牛憨笨. 超衍射成像中双螺旋点扩展函数的三维定位精度. 物理学报, 2013, 62(14): 144201. doi: 10.7498/aps.62.144201
    [16] 刘国忠, 周哲海, 邱钧, 王晓飞, 刘桂礼, 王瑞康. 幅值和相位配准技术及其在光学相干层析血流成像中的应用. 物理学报, 2013, 62(15): 158702. doi: 10.7498/aps.62.158702
    [17] 陈丹妮, 刘磊, 于斌, 牛憨笨. HeLa细胞突起中微丝束的纳米分辨荧光成像. 物理学报, 2010, 59(10): 6948-6954. doi: 10.7498/aps.59.6948
    [18] 胡昕, 张继彦, 杨国洪, 刘慎业, 丁永坤. 基于布拉格反射镜的X射线多色单能成像谱仪. 物理学报, 2009, 58(9): 6397-6402. doi: 10.7498/aps.58.6397
    [19] 王琛, 王桂英, 徐至展. 全内反射荧光显微术应用于单分子荧光的纵向成像. 物理学报, 2004, 53(5): 1325-1330. doi: 10.7498/aps.53.1325
    [20] 张连水, 李云静, 李晓苇, 傅广生. 碘分子单频及双频多光子光谱. 物理学报, 1997, 46(6): 1088-1095. doi: 10.7498/aps.46.1088
计量
  • 文章访问数:  1505
  • PDF下载量:  125
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-10-25
  • 修回日期:  2023-12-01
  • 上网日期:  2023-12-15
  • 刊出日期:  2024-03-20

/

返回文章
返回