磁共振酰胺质子转移成像磁化转移率检测方法及相关设备与流程

专利检索2022-05-10  3



1.本技术涉及医疗影像技术领域,特别涉及一种磁共振酰胺质子转移成像磁化转移率检测方法及相关设备。


背景技术:

2.化学交换饱和转移是一种新型的磁共振分子成像技术,其具有检测组织、细胞中游离蛋白质和多肽等优点,在生物医院领域迅速引起广泛的探索研究。其中,酰胺质子转移(amide proton transfer,apt)技术的新成像技术引入了新颖图像对比度,apt技术基于化学交换饱和转移(cest)效应,通过选择性预饱和游离蛋白质及多肽中的酰胺质子的信号,探测经过与酰胺质子交换之后的周围自由水信号,通过采集自由水饱和前后信号的变化,间接反应游离蛋白质及多肽浓度的改变,使得磁共振酰胺质子转移成像成为识别肿瘤组织的重要指标。因此,希望可以全面的获取磁共振酰胺质子转移成像的转移率数据,以提高磁共振酰胺质子转移成像的评估价值。


技术实现要素:

3.本技术要解决的技术问题在于,针对现有技术的不足,提供一种磁共振酰胺质子转移成像磁化转移率可视量化的检测方法及相关设备。
4.为了解决上述技术问题,本技术实施例第一方面提供了一种磁共振酰胺质子转移成像磁化转移率检测方法,所述方法包括:
5.获取磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型;
6.确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域;
7.基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据,其中,所述转移率数据包括所述感兴趣区域中所有像素点的酰胺质子转移成像磁化转移率数值,酰胺质子转移成像磁化转移率数值的统计数据、酰胺质子转移成像磁场均匀度的z谱图、酰胺质子转移成像磁化转移率数值的直方图、酰胺质子转移成像的映射图像、主磁场的映射图像以及酰胺质子转移成像与原始图像叠加的映射图像中的一种或者多种。
8.所述磁共振酰胺质子转移成像磁化转移率检测方法,其中,所述确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域具体包括:
9.将所述磁共振酰胺质子转移成像的图像数据集中的各图像数据进行归一化,以得到归一化后的各图像数据;
10.在归一化后的各图像数据中选取基础图像数据,并且确定所述基础图像数据对应的感兴趣区域;
11.将所述感兴趣区域作为所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域。
12.所述磁共振酰胺质子转移成像磁化转移率检测方法,其中,所述基础图像数据为归一化后的各图像数据中按照成像顺序位于最前的归一化后的图像数据;所述确定所述基础图像数据对应的感兴趣区域具体为:
13.将所述基础图像数据中的各像素点的像素值乘以255,以将各像素点的像素值映射至0

255之间,以得到映射后的基础图像数据;
14.基于映射后的基础图像数据,确定所述基础图像数据对应的感兴趣区域。
15.所述磁共振酰胺质子转移成像磁化转移率检测方法,其中,所述基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据具体包括:
16.基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域的转移率图像,其中,所述在z频谱中以水的频率作为中心频率,将水的频率设置为0ppm,从z频谱中频率为3.3

3.7ppm处提取酰胺质子的磁化转移率图像;
17.对所述转移率图像进行去噪,并基于去噪后的转移率图像确定所述感兴趣区域对应的转移率数据。
18.所述磁共振酰胺质子转移成像磁化转移率检测方法,其中,所述对所述转移率图像进行去噪具体包括:
19.分别所述转移率图像中的各像素点的像素值的像素绝对值;
20.将像素绝对值小于预设阈值的像素点过滤掉,以得到去噪后的转移率图像。
21.所述磁共振酰胺质子转移成像磁化转移率检测方法,其中,所述统计数据包括数值的平均值、中位数、最大值、最小值、标准偏差、峰度系数、偏度系数以及第90%位数中的一种或者多种。
22.本技术实施例第二方面提供了一种磁共振酰胺质子转移成像磁化转移率检测系统,所述的检测系统包括:
23.获取模块,用于获取磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型;
24.第一确定模块,用于确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域;
25.第二确定模块,用于基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据。
26.所述磁共振酰胺质子转移成像磁化转移率检测系统所述的检测系统,其中,所述件系统为基于python语言建立的。
27.本技术实施例第三方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如上任一所述的磁共振酰胺质子转移成像磁化转移率检测方法中的步骤。
28.本技术实施例第四方面提供了一种终端设备,其包括:处理器、存储器及通信总线;所述存储器上存储有可被所述处理器执行的计算机可读程序;
29.所述通信总线实现处理器和存储器之间的连接通信;
30.所述处理器执行所述计算机可读程序时实现如上任一所述的磁共振酰胺质子转
移成像磁化转移率检测方法中的步骤。
31.有益效果:与现有技术相比,本技术提供了一种磁共振酰胺质子转移成像磁化转移率检测方法及相关设备,所述方法包括获取磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型;确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域;基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据。本技术通过对头颈部肿瘤的磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型进行分析,得到确定所述感兴趣区域对应的转移率数据,并且提高了转移率数据的丰富性,从而提高了磁共振酰胺质子转移成像的应用价值。
附图说明
32.为了更清楚地说明本技术实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本技术的一些实施例,对于本领域普通技术人员而言,在不符创造性劳动的前提下,还可以根据这些附图获得其他的附图。
33.图1为本技术提供的磁共振酰胺质子转移成像磁化转移率检测方法的流程图。
34.图2为z谱图。
35.图3为主磁场的映射图像。
36.图4为酰胺质子转移成像与原始图像叠加的映射图像。
37.图5为酰胺质子转移成像的映射图像。
38.图6为酰胺质子转移成像磁化转移率数值的直方图。
39.图7为本技术提供的磁共振酰胺质子转移成像磁化转移率检测装置的结构原理图。
40.图8为本技术提供的终端设备的结构原理图。
具体实施方式
41.本技术提供一种磁共振酰胺质子转移成像磁化转移率检测方法及相关设备,为使本技术的目的、技术方案及效果更加清楚、明确,以下参照附图并举实施例对本技术进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本技术,并不用于限定本技术。
42.本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本技术的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或无线耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的全部或任一单元和全部组合。
43.本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语),具有与本技术所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语,应该被理解为具有与现有技术的上下文中的
意义一致的意义,并且除非像这里一样被特定定义,否则不会用理想化或过于正式的含义来解释。
44.应理解,本实施例中各步骤的序号和大小并不意味着执行顺序的先后,各过程的执行顺序以其功能和内在逻辑确定,而不应对本技术实施例的实施过程构成任何限定。
45.发明人经过研究发现,化学交换饱和转移成像(cest)是一种新型的磁共振分子成像技术,具有无创伤检测组织、细胞中游离蛋白质和多肽等优点,在生物医学领域迅速引起广泛的探索研究。其中,酰胺质子转移(apt)基于化学交换饱和转移效应,通过选择性预饱和游离蛋白质及多肽中的酰胺质子的信号,探测经过与酰胺质子交换之后的周围自由水信号,通过采集自由水饱和前后信号的变化,间接反应游离蛋白质及多肽浓度的改变,使得磁共振酰胺质子转移成像成为识别肿瘤组织的重要指标。因此,希望可以全面的获取磁共振酰胺质子转移成像的转移率数据,以提高磁共振酰胺质子转移成像的评估价值。
46.为了解决上述问题,在本技术实施例中,获取磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型;确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域;基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据。本技术通过对头颈部肿瘤的磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型进行分析,得到确定所述感兴趣区域对应的转移率数据,并且提高了转移率数据的丰富性,从而提高了磁共振酰胺质子转移成像的应用价值。
47.下面结合附图,通过对实施例的描述,对申请内容作进一步说明。
48.本实施例提供了一种磁共振酰胺质子转移成像磁化转移率检测方法,所述磁共振酰胺质子转移成像磁化转移率检测方法可以作为一个功能模块,通过该功能模块来实现执行本实施例提供的磁共振酰胺质子转移成像磁化转移率检测方法的步骤。在一个具体实现方式中,该功能模块可以通过python程序设计语言实现的,并且在通过python程序设计语言实现该功能,调用了simpleitk,opencv,scipy等开源库。
49.如图1所示,本实施例提供的磁共振酰胺质子转移成像磁化转移率检测方法具体包括:
50.s10、获取磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型。
51.具体地,磁共振酰胺质子转移成像的图像数据集包括酰胺质子转移成像序列,例如,酰胺质子转移成像序列包括通过磁共振扫描仪扫描,饱和度偏移量从

7ppm

7ppm,其中,饱和度偏移量的偏移步长为0.2ppm。所述磁共振酰胺质子转移成像序列的零阶对称模型为磁共振酰胺质子转移成像序列需添加的匀场模型,通过该零阶对称模型可以提高主磁场的均匀性,并把水校正为0ppm,使水左右两边的信号呈对称分辨。在本实施例的一个实现方式中,磁共振酰胺质子转移成像的图像数据集为头颈部肿瘤磁共振酰胺质子转移成像的图像数据集,当然,在实际应用中,磁共振酰胺质子转移成像的图像数据集可以是其他部分的磁共振酰胺质子转移成像的图像数据集。
52.s20、确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域。
53.具体地,所述感兴趣区域为基于磁共振酰胺质子转移成像的图像数据集中的一张磁共振酰胺质子基础图像数据确定,并且感兴趣区域包括于该磁共振酰胺质子所有图像数
据中。在选取到感应器区域后,可以拟合感兴趣区域的z谱,通过分析z谱图观察最低点是否出现在0ppm处,两边的信号是否对称性分布,只有最低点在0ppm才说明主磁场的非均匀性得到了较好的解决,这是质量控制的需要,也是保证apt值准确性的前提。另外,使用了较宽的偏移频率扫描范围,避免或减少z谱偏移后出现apt信号的丢失。其中,感兴趣区域为磁共振酰胺质子图像数据中的肿块或病变图像区域,避开血管、空气以及骨头的图像区域,不包含正常组织。
54.在本实施例的一个实现方式中,所述确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域具体包括:
55.将所述磁共振酰胺质子转移成像的图像数据集中的各图像数据进行归一化,以得到归一化后的各图像数据;
56.在归一化后的各图像数据中选取基础图像数据,并且确定所述基础图像数据对应的感兴趣区域;
57.将所述感兴趣区域作为所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域。
58.具体地,所述归一化指的是将图像数据中的像素点的像素取值范围映射至0

1范围之内。基础图像数据为磁共振酰胺质子转移成像的图像数据集中的一图像数据,其中,基础图像数据可以随机选取,也可以是根据使用需求选取。在本实施例的一个实现方式中,基础图像数据可以为磁共振酰胺质子转移成像的图像数据集中初始饱和度对应的图像数据,也就是说,基础图像数据可以为该磁共振酰胺质子转移成像的图像数据集中按照图像采集顺序位于最前的图像数据。在获取到基础图像数据后,可以将基础图像数据显示给用户,以使得用户可以在所述基础图像数据上选取感兴趣区域,其中,感兴趣区域根据可以为肿块或病变区域,其可以任意形状,例如,为圆形,正方形或者其他不规则形状等。
59.在本实施例的一个实现方式中,所述确定所述基础图像数据对应的感兴趣区域具体为:
60.将所述基础图像数据中的各像素点的像素值乘以255,以将各像素点的像素值映射至0

255之间,以得到映射后的基础图像数据;
61.基于映射后的基础图像数据,确定所述基础图像数据对应的感兴趣区域。
62.具体地,将所述基础图像数据中的各像素点的像素值乘以255指的是将归一化后的基础图像像素值乘以255,以将基础图像数据中的像素点的像素取值映射至0

255区间上,然后通过调用opencv窗口的鼠标事件,展示基础图像数据,并通过人机交互式的方式进行感兴趣区域的选取。在通过人机交互式的方式进行感兴趣区域后,将选取到的感兴趣区域提取出来以形成感兴趣区域图。
63.s30、基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据。
64.具体地,所述转移率数据包括所述感兴趣区域中所有像素点的酰胺质子转移成像磁化转移率数值,酰胺质子转移成像磁化转移率数值的统计数据、酰胺质子转移成像磁场均匀度的z谱图、酰胺质子转移成像磁化转移率数值的直方图、酰胺质子转移成像的映射图像、主磁场的映射图像以及酰胺质子转移成像与原始图像叠加的映射图像中的一种或者多种。例如,如图2所示的z谱图,通过z谱图观察最低点是否出现在0ppm处,两边的信号是否对
称性分布,以体现主磁场的均匀性;如图3所示的主磁场的映射图像,通过主磁场的映射图像说明感兴趣区主磁场的均匀性;如图4所示的头颈部肿瘤的酰胺质子转移成像与原始图像叠加的映射图像,用于直观观察感兴趣区是否存在肿瘤蛋白,根据颜色深浅,分析肿瘤的增殖程度,颜色越红或越蓝,apt绝对值越高,说明肿瘤细胞增殖旺盛;如图5所示的酰胺质子转移成像的映射图像,用于直观观察感兴趣区是否存在肿瘤蛋白,及其肿瘤细胞增殖旺盛程度;如图6所示的酰胺质子转移成像磁化转移率数值的直方图,用于感兴趣区apt值的可视量化,尤其肿瘤治疗前后的对比,可直观的观察肿瘤细胞的存活率。
65.所述统计数据包括数值的平均值、中位数、最大值、最小值、标准偏差、峰度系数、偏度系数以及第90%位数中的一种或者多种。在本实施例的一个实现方式中,所述转移率数据包括第90%位数,由于酰胺质子磁化传递率数值正负值,且数值多数呈非正态分布,并且通过60例头颈部酰胺质子转移成像磁化传递率分析,感兴趣区域的第90%位数相对于平均值可以更好的提供生物学信息。
66.在实施例的一个实现方式中,所述基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据具体包括:
67.基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域的转移率图像;
68.对所述转移率图像进行去噪,并基于去噪后的转移率图像确定所述感兴趣区域对应的转移率数据。
69.具体地,所述感兴趣区域的转移率图像为以水的频率作为中心频率,在z谱图中将水的频率设置为0ppm,从z频谱中频率为3.3

3.7ppm(以3.5ppm为中心,外扩
±
0.2ppm)处获取到的酰胺质子的磁化转移率的图像。此外,转移率图像中的每个像素点为感兴趣区域中的各像素点对应的转移率。也就是说,在获取到感兴趣区域后,基于磁共振酰胺质子转移成像的图像数据集中的各图像数据确定感兴趣区域的z谱,然后从z频谱中频率为3.3

3.7ppm处获取酰胺质子的磁化转移率的图像,以得到感兴趣区域对应的转移率图像。
70.在本实施的一个实现方式中,不对称磁化转移百分比的计算如下面公式:mtrasym(%)=(s

δω

s δω)
÷
s0其中,s

δω代表负向频率偏移时的信号强度,s δω代表同样对称位置的正向频率偏移时的信号强度;s0代表没有射频饱和时的信号强度。而在aptw加权像中,不对称磁化转移百分比mtrasym(%)为转化率apt的权重,aptweight,并且酰胺质子是在 3.5ppm,从而以3.5ppm计算的mtrasym为apt的权重。因此,转移率的计算公式可以为:aptw%=mtrasym(%)=(s

δω

s δω)
÷
s0,令δω= 3.5ppm。也就是说,s

δω代表

3.5ppm频率偏移时的信号强度;s δω代表 3.5ppm频率偏移时的信号强度。相应的,aptw%的计算公式可以改写为:aptw%=(s

3.5ppm

s 3.5ppm)
÷
s0。
71.在本实施例的一个实现方式中,所述对所述转移率图像进行去噪具体包括:
72.分别所述转移率图像中的各像素点的像素值的像素绝对值;
73.将像素绝对值小于预设阈值的像素点过滤掉,以得到去噪后的转移率图像。
74.具体地,预设阈值为预先设置的,用于过滤转移率图像中的像素点,其中,转移率图像中的各像素点与感兴趣区域中的各像素点相对应,并且各像素点的像素值用于反映其对应的像素点的转移率。像素绝对值小于预设阈值的像素点过滤掉可以将转移率图像中具
有不合理的酰胺质子转移成像值,拟合优度差、非均质性、非常大或部分(不完全)饱和的像素排出。在本实施例的一个实现方式中,所述预设阈值可以基于转移率图像的平均数确定,例如,预设阈值为平均数的一半,或者预设阈值与平均数的差值绝对值为预设的。此外,在获取到去噪后的转移率图像后,可以基于去噪后的转移率图像计算平均值、中位数、最大值、最小值、标准偏差、峰度系数、偏度系数以及第90%位数,同时可以将转移率图像叠加到原始图上,其中,原始图为用于确定感兴趣区域的基础图像数据。此外,在获取到转移率数据后,可以将转移率数据通过显示窗口显示给用户,以使得转移率数据可视化。当然,在基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域的转移率图像时,可以采用现有方法确定得到转移率图像,这里不具体说明。
75.本实施例提供了一种磁共振酰胺质子转移成像磁化转移率检测方法及相关设备,所述方法包括获取磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型;确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域;基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据。本技术通过对头颈部肿瘤的磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型进行分析,得到确定所述感兴趣区域对应的转移率数据,并且提高了转移率数据的丰富性,从而提高了磁共振酰胺质子转移成像的应用价值。
76.基于上述磁共振酰胺质子转移成像磁化转移率检测方法,本实施例提供了一种磁共振酰胺质子转移成像磁化转移率检测系统,如图7所示,所述的检测系统包括:
77.获取模块100,用于获取磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型;
78.第一确定模块200,用于确定所述磁共振酰胺质子转移成像的图像数据集对应的感兴趣区域;
79.第二确定模块300,用于基于所述磁共振酰胺质子转移成像的图像数据集以及磁共振酰胺质子转移成像序列的零阶对称模型,确定所述感兴趣区域对应的转移率数据。
80.所述磁共振酰胺质子转移成像磁化转移率检测系统,其中,所述的检测系统基于python语言建立的。
81.基于上述磁共振酰胺质子转移成像磁化转移率检测方法,本实施例提供了一种计算机可读存储介质,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如上述实施例所述的磁共振酰胺质子转移成像磁化转移率检测方法中的步骤。
82.基于上述磁共振酰胺质子转移成像磁化转移率检测方法,本技术还提供了一种终端设备,如图8所示,其包括至少一个处理器(processor)20;显示屏21;以及存储器(memory)22,还可以包括通信接口(communications interface)23和总线24。其中,处理器20、显示屏21、存储器22和通信接口23可以通过总线24完成相互间的通信。显示屏21设置为显示初始设置模式中预设的用户引导界面。通信接口23可以传输信息。处理器20可以调用存储器22中的逻辑指令,以执行上述实施例中的方法。
83.此外,上述的存储器22中的逻辑指令可以通过软件功能单元的形式实现并作为独
立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。
84.存储器22作为一种计算机可读存储介质,可设置为存储软件程序、计算机可执行程序,如本公开实施例中的方法对应的程序指令或模块。处理器20通过运行存储在存储器22中的软件程序、指令或模块,从而执行功能应用以及数据处理,即实现上述实施例中的方法。
85.存储器22可包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据终端设备的使用所创建的数据等。此外,存储器22可以包括高速随机存取存储器,还可以包括非易失性存储器。例如,u盘、移动硬盘、只读存储器(read

only memory,rom)、随机存取存储器(random access memory,ram)、磁碟或者光盘等多种可以存储程序代码的介质,也可以是暂态存储介质。
86.此外,上述存储介质以及移动终端中的多条指令处理器加载并执行的具体过程在上述方法中已经详细说明,在这里就不再一一陈述。
87.最后应说明的是:以上实施例仅用以说明本技术的技术方案,而非对其限制;尽管参照前述实施例对本技术进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本技术各实施例技术方案的精神和范围。
转载请注明原文地址:https://win.8miu.com/read-50176.html

最新回复(0)