一种装配式建筑BIM模型构件生成方法及系统

未命名 10-25 阅读:71 评论:0

一种装配式建筑bim模型构件生成方法及系统
技术领域
1.本发明涉及装配式建筑领域,尤其是涉及装配式建筑bim模型构件生成方法及系统。


背景技术:

2.装配式建筑是指通过将预先生产好的建筑构件在现场进行组装而成的建筑。装配式建筑将大部分建筑工作转移到工厂内,将建筑构件在受控的环境中精确制造,然后在现场进行简单组装。这种方法大大提高了施工效率,节省了时间,相较于传统的逐项施工,装配式建筑可以将项目完成时间缩短数周甚至数月。而且工厂环境中的生产过程可监控和检测,减少了材料和工艺方面的缺陷,具有更高质量的建筑,减少了后期维修和修复的需求。
3.以往对建筑的建模大多是采用三维建模方法,但是实际施工可能会因材料、工艺或测量等因素而产生偏差,在复杂的建筑或工程项目中,三维建模可能无法完全捕捉复杂的几何形状和细节,在施工过程中需要不断的修改三维模型,尤其是实际生产的构件和模型中构件有较大差别时。点云是通过激光扫描或3d扫描采集的实际数据,因此提供了高度精确的现实世界反映,有助于准确捕捉建筑物构件的几何形状、尺寸和细节。但是如何快速、高效的生成构件是一个难点。


技术实现要素:

4.为了解决上述问题,本发明提供了一种装配式建筑bim模型构件生成方法,所述方法包括以下步骤:s1,获取构件的点云,将每个点分别作为基准点,统计以基准点为球心且半径r的球体内的点云数,根据所述点云数计算k值,计算距离所述基准点最近的k个点,并记录所述k个点中距离所述基准点最远距离d,对所述k个点进行平面拟合;计算距离所述基准点的距离小于2d的点,根据所述平面对距离小于2d的点进行过滤得到k个点,基于过滤得到的k个点计算基准点的法向量;s2,根据所述法向量对点云中的点进行聚类,并基于聚类中点云的点数对聚类过滤,从过滤后的每个聚类得到至少一个集合,对每个集合中的点进行平面拟合得到构件的平面,基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件。
5.优选地,所述根据所述平面对距离小于2d的点进行过滤得到k个点,具体为:将距离小于2d的点投影到所述平面上,计算以所述基准点为原点的坐标系中每个象限中投影点数,根据所述投影点数对小于2d的点进行过滤直到剩余k个点,使得剩余k个点在每个象限的投影数的方差最小。
6.优选地,所述从过滤后的每个聚类得到至少一个集合,具体为:s21,从聚类中选取一个关键点放入集合,迭代执行将聚类中距离集合中任一点的距离小于d的点从聚类中删除并放入集合中,直到聚类中没有点距离集合中任一点的距离小于d;
s22,对聚类中剩余点重复执行s21得到一个或多个集合,根据集合中的点数对集合过滤,得到聚类对应的至少一个集合。
7.优选地,所述每个聚类中选取一个关键点放入集合,具体为:对于聚类的每个点i,计算聚类中其他点和点i的距离小于d的点数,选取对应的点数最大的点作为关键点。
8.优选地,所述基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件,具体为:获取每个平面的边界点,根据边界点和所述点云剩余的点采用点云曲面重建方法进行三维重建;将重建后的三维文件导出为rfa文件。
9.另外,本发明还提供了一种装配式建筑bim模型构件生成系统,所述系统包括以下模块:法向量计算模块,用于获取构件的点云,将每个点分别作为基准点,统计以基准点为球心且半径r的球体内的点云数,根据所述点云数计算k值,计算距离所述基准点最近的k个点,并记录所述k个点中距离所述基准点最远距离d,对所述k个点进行平面拟合;计算距离所述基准点的距离小于2d的点,根据所述平面对距离小于2d的点进行过滤得到k个点,基于过滤得到的k个点计算基准点的法向量;三维重建模块,用于根据所述法向量对点云中的点进行聚类,并基于聚类中点云的点数对聚类过滤,从过滤后的每个聚类得到至少一个集合,对每个集合中的点进行平面拟合得到构件的平面,基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件。
10.优选地,所述根据所述平面对距离小于2d的点进行过滤得到k个点,具体为:将距离小于2d的点投影到所述平面上,计算以所述基准点为原点的坐标系中每个象限中投影点数,根据所述投影点数对小于2d的点进行过滤直到剩余k个点,使得剩余k个点在每个象限的投影数的方差最小。
11.优选地,所述从过滤后的每个聚类得到至少一个集合,具体为:s21,从聚类中选取一个关键点放入集合,迭代执行将聚类中距离集合中任一点的距离小于d的点从聚类中删除并放入集合中,直到聚类中没有点距离集合中任一点的距离小于d;s22,对聚类中剩余点重复执行s21得到一个或多个集合,根据集合中的点数对集合过滤,得到聚类对应的至少一个集合。
12.优选地,所述每个聚类中选取一个关键点放入集合,具体为:对于聚类的每个点i,计算聚类中其他点和点i的距离小于d的点数,选取对应的点数最大的点作为关键点。
13.优选地,所述基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件,具体为:获取每个平面的边界点,根据边界点和所述点云剩余的点采用点云曲面重建方法进行三维重建;将重建后的三维文件导出为rfa文件。
14.此外,本发明还提供了一种计算机可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序在被计算机执行时,实现如上所述的方法。
15.针对现有的利用点云生成构件的过程计算复杂、准确度低的问题,本发明在计算法向量时能够根据周围点的密度情况设置k值,通过k近邻的方式计算出法向量,同时能够尽可能均匀的得到基准点周围的k个点,防止在利用k近邻计算法线量中k个近邻集中到一起;然后根据法向量对点云的点进行聚类,得到构件的平面,这样减少了在利用构件的点云生成构件三维图过程中的计算量,尤其是减少了点云曲面重建中的点云数,提高了生成速度。
附图说明
16.为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
17.图1为实施例一的流程图;图2为点云的k近邻邻域示意图;图3为初始投影点分布示意图;图4为优化后投影点分布示意图;图5为一个具体实施例的流程图。
具体实施方式
18.在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个
……”
限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
19.下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
20.实施例一,本发明提供了一种装配式建筑bim模型构件生成方法,如图1所示,所述方法包括以下步骤:s1,获取构件的点云,将每个点分别作为基准点,统计以基准点为球心且半径r的球体内的点云数,根据所述点云数计算k值,计算距离所述基准点最近的k个点,并记录所述k个点中距离所述基准点最远距离d,对所述k个点进行平面拟合;计算距离所述基准点的距离小于2d的点,根据所述平面对距离小于2d的点进行过滤得到k个点,基于过滤得到的k个点计算基准点的法向量;点云(point cloud)是某个坐标系下的点的数据集,点包含了三维坐标x,y,z、颜色、分类值、强度值、时间等等,一般而言,点云至少包括三维坐标。法向量又称为法矢、法线
矢量,法向量是点云分析中的重要内容,法向量的计算准确度直接影响到点云的后续分析和处理,也影响这最后结果。k近邻法计算法向量是应用最为广泛的法向量计算方法,所谓k近邻是指距离某个点距离最近的k个点,这就会出现一种情况,如果一个点的k个近邻集中在这个点的一侧或者集中在这个点的某个小范围内,假设k=3,图2中1、2、3分别表示距离中心点最近的三个点,最后计算的法向量实际是这些点的法向量,没有利用到这个点其他方向上的点的信息,这进一步导致法向量的计算出现不准确。本发明在得到k值后,计算这k个点中和基准点距离最远的点,然后基于最远点与基准点的距离d,重新计算距离基准点的距离小于2d的所有点,根据距离和点的分布过滤得到k个点,然后通过过滤后的k个点计算基准点的法向量。
21.其中,所述根据所述点云数计算k值,具体为:k值和所述点云数成正比。在一个更为详细的实施例中,计算所述点云数与预设点云数的比值,根据所述比值与初始k值得到k值。
22.s2,根据所述法向量对点云中的点进行聚类,并基于聚类中点云的点数对聚类过滤,从过滤后的每个聚类得到至少一个集合,对每个集合中的点进行平面拟合得到构件的平面,基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件。
23.装配式建筑的构件很多都是对称或者较为规整的构件,相比于异形构件,装配式建筑的构件存在多个平面,而且平面的面积比较大,例如矩形砼梁实际是一个立方体,存在6个平面,在同一个平面的点的法向量相同或者接近,根据法向量对点云中的点进行聚类即可识别点云的点是否在一个平面中。矩形砼梁是一个理想的模型,完全有平面组成,装配式建筑中还会存在其他例如y型三通,这些构件也会存在非平面的地方。基于此,本发明根据点云的法向量对点云中的点进行聚类后,进一步对聚类进行过滤,当一个聚类中点云数小于预设点云数时,则将这个聚类过滤掉。其中,预设点云数和点云密度以及最小平面面积相关,点云密度通过点云采集设备获取,最小面积和构件相关,在一个更为具体的实施例中,最小面积是由用户输入。
24.每个聚类中的点的法向量相同或者相近,但是两个平行且朝向同一侧的两个平面在一个聚类中,对过滤后的每个聚类还需要再次分类,以区分上述情况中的不同平面,这样得到的每个集合对应一个平面。
25.当平面比较大时,例如一个墙体,则在一个平面中会存在很多点,在一个更为具体实施例中,所述对每个集合中的点进行平面拟合得到构件的平面,具体为:从集合中第一个点开始,按照预设间隔选取至少一个点,对所述至少一个点中每个点执行如下操作:计算得到点的k1个近邻,计算k1个近邻点x值的平均值、y值的平均值、z值的平均值和法向量的平均值,将上述平均值表示的点放入队列中;对队列中的点进行平面拟合得到构件的平面。
26.例如一个集合中含有100个点,则每隔十个点选取一个点,选取出1、11、21、31、41、51、61、71、81、91十个点,计算第1个点的k1个近邻,然后用k1个近邻的x值的平均值、y值的平均值、z值的平均值和法向量的平均值替代第1个点,以此类推,直到得到含有十个点的队列,然后对队列中的十个点进行平面拟合得到集合对应的平面。
27.为了能够使得一个点的k近邻在空间上能够分布均匀,所述根据所述平面对距离小于2d的点进行过滤得到k个点,具体为:
将距离小于2d的点投影到所述平面上,计算以所述基准点为原点的坐标系中每个象限中投影点数,根据所述投影点数对小于2d的点进行过滤直到剩余k个点,使得剩余k个点在每个象限的投影数的方差最小。
28.在原点确定后,在一个平面上建立一个坐标系有多种方式,例如x轴水平向右或者水平向左等,本发明中对于坐标系中x轴的方向和y轴的方向不作具体限定,在一个实施例中,x轴的方向为随机;当然,也可以根据距离小于2d的点在平面上的投影情况确定。
29.其中,所有象限中投影点数的方差最小,则每个象限中投影的点数分布最均匀。方差的具体计算过程是:计算每个象限中投影点数的期望,例如第一象限为5个,第二象限为10个,第三象限为2个,第四象限为2个,如图3所示,其中每个实线绘制的圆表示一个投影点,如果k=8,如果使得剩余k个点在每个象限的投影数的方差最小,则每个象限中最优的点数是接近2个,此时,从第一象限中删除距离原点最远的3个点,从第二象限中删除8个距离原点最远的点,得到的方差最小,最后结果如图4所示,其中每个虚线绘制的圆表示过滤掉的点。
30.所述根据所述投影点数对小于2d的点进行过滤直到剩余k个点,使得剩余k个点在每个象限的投影数的方差最小,更为具体的过程是:s11,计算得到k/4;s12,计算每个象限投影点数与k/4的差,按照差从大到小的顺序对象限排序,每次从排序第一的象限中删除距离原点最远的点;s13,执行s12 m次,所述m为距离小于2d的点投影到所述平面的个数与k的差值。
31.一个聚类是法向量相同或相近的所有点的集合,但是法向量相同或相近可能不在同一个平面上,所述从过滤后的每个聚类得到至少一个集合,如图5所示,具体为:s21,从聚类中选取一个关键点放入集合,迭代执行将聚类中距离集合中任一点的距离小于d的点从聚类中删除并放入集合中,直到聚类中没有点距离集合中任一点的距离小于d;假设共有4个聚类,对于第一个聚类,从第一个聚类中选取一个关键点放入集合c1中;将第一聚类中距离这个关键点的距离小于d的点从第一个聚类中删除并放入集合c1中;然后判断聚类中是否存在点到集合c1中的任意点的距离小于d,如果存在,则将这个点从第一个聚类中删除并放入集合c1中;不断重复这个过程,直到满足预设条件;预设条件也即聚类中不存在点到集合c1中的任意点的距离小于d。
32.s22,对聚类中剩余点重复执行s21得到一个或多个集合,根据集合中的点数对集合过滤,得到聚类对应的至少一个集合。
33.对于聚类中剩余的点不断重复执行s21可以得到一个或多个集合,然后过滤掉点数小于所述预设点云数的集合,如果一个聚类中的所有集合的点数都小于所述预设点云数,则将聚类删除,最后每个聚类对应至少一个集合。
34.为了加快集合的生成,在一个具体实施例中,所述每个聚类中选取一个关键点放入集合,具体为:对于聚类的每个点i,计算聚类中其他点和点i的距离小于d的点数,选取对应的点
数最大的点作为关键点。也即选择聚类中周围分布的点最密集的点作为关键点。
35.在构建平面后,点云中剩余的点数较少,通过曲面重建可以快速的对这些点进行三维重建,所述基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件,具体为:获取每个平面的边界点,根据边界点和所述点云剩余的点采用点云曲面重建方法进行三维重建;将重建后的三维文件导出为rfa文件。
36.在识别到点云的平面后,剩下的点云部分要么不是平面,要么是面积较小的平面,对于装配式建筑构件,剩下的点更多的是在平面的连接处,为了能够得到平滑的边界,本发明获取每个平面的边界点,将边界点和剩余的点进行点云曲面重建。点云曲面的方法包括但不限于基于voronoi图的点云曲面重建、基于点云插值的曲面重建。在重建后将文件导出为rfa文件即可导入到revit中使用。通过本发明生成构件的速度更快,这得益于先进行平面识别,然后对剩余的点进行三维重建,避免了对大量点云的三维重建。
37.实施例二,本发明还提供了一种装配式建筑bim模型构件生成系统,所述系统包括以下模块:法向量计算模块,用于获取构件的点云,将每个点分别作为基准点,统计以基准点为球心且半径r的球体内的点云数,根据所述点云数计算k值,计算距离所述基准点最近的k个点,并记录所述k个点中距离所述基准点最远距离d,对所述k个点进行平面拟合;计算距离所述基准点的距离小于2d的点,根据所述平面对距离小于2d的点进行过滤得到k个点,基于过滤得到的k个点计算基准点的法向量;三维重建模块,用于根据所述法向量对点云中的点进行聚类,并基于聚类中点云的点数对聚类过滤,从过滤后的每个聚类得到至少一个集合,对每个集合中的点进行平面拟合得到构件的平面,基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件。
38.优选地,所述根据所述平面对距离小于2d的点进行过滤得到k个点,具体为:将距离小于2d的点投影到所述平面上,计算以所述基准点为原点的坐标系中每个象限中投影点数,根据所述投影点数对小于2d的点进行过滤直到剩余k个点,使得剩余k个点在每个象限的投影数的方差最小。
39.优选地,所述从过滤后的每个聚类得到至少一个集合,具体为:s21,从聚类中选取一个关键点放入集合,迭代执行将聚类中距离集合中任一点的距离小于d的点从聚类中删除并放入集合中,直到聚类中没有点距离集合中任一点的距离小于d;s22,对聚类中剩余点重复执行s21得到一个或多个集合,根据集合中的点数对集合过滤,得到聚类对应的至少一个集合。
40.优选地,所述每个聚类中选取一个关键点放入集合,具体为:对于聚类的每个点i,计算聚类中其他点和点i的距离小于d的点数,选取对应的点数最大的点作为关键点。
41.优选地,所述基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件,具体为:获取每个平面的边界点,根据边界点和所述点云剩余的点采用点云曲面重建方法
进行三维重建;将重建后的三维文件导出为rfa文件。
42.实施例三,本发明还提供了一种计算机可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序在被计算机执行时,实现如实施例一所述的方法。
43.通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助加必需的通用硬件平台的方式来实现,当然也可以通过硬件和软件结合的方式来实现。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以计算机产品的形式体现出来,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
44.最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

技术特征:
1.一种装配式建筑bim模型构件生成方法,其特征在于,所述方法包括以下步骤:s1,获取构件的点云,将每个点分别作为基准点,统计以基准点为球心且半径r的球体内的点云数,根据所述点云数计算k值,计算距离所述基准点最近的k个点,并记录所述k个点中距离所述基准点最远距离d,对所述k个点进行平面拟合;计算距离所述基准点的距离小于2d的点,根据所述平面对距离小于2d的点进行过滤得到k个点,基于过滤得到的k个点计算基准点的法向量;s2,根据所述法向量对点云中的点进行聚类,并基于聚类中点云的点数对聚类过滤,从过滤后的每个聚类得到至少一个集合,对每个集合中的点进行平面拟合得到构件的平面,基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件。2.如权利要求1所述的方法,其特征在于,所述根据所述平面对距离小于2d的点进行过滤得到k个点,具体为:将距离小于2d的点投影到所述平面上,计算以所述基准点为原点的坐标系中每个象限中投影点数,根据所述投影点数对小于2d的点进行过滤直到剩余k个点,使得剩余k个点在每个象限的投影数的方差最小。3.如权利要求1所述的方法,其特征在于,所述从过滤后的每个聚类得到至少一个集合,具体为:s21,从聚类中选取一个关键点放入集合,迭代执行将聚类中距离集合中任一点的距离小于d的点从聚类中删除并放入集合中,直到聚类中没有点距离集合中任一点的距离小于d;s22,对聚类中剩余点重复执行s21得到一个或多个集合,根据集合中的点数对集合过滤,得到聚类对应的至少一个集合。4.如权利要求3所述的方法,其特征在于,所述每个聚类中选取一个关键点放入集合,具体为:对于聚类的每个点i,计算聚类中其他点和点i的距离小于d的点数,选取对应的点数最大的点作为关键点。5.如权利要求1所述的方法,其特征在于,所述基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件,具体为:获取每个平面的边界点,根据边界点和所述点云剩余的点采用点云曲面重建方法进行三维重建;将重建后的三维文件导出为rfa文件。6.一种装配式建筑bim模型构件生成系统,其特征在于,所述系统包括以下模块:法向量计算模块,用于获取构件的点云,将每个点分别作为基准点,统计以基准点为球心且半径r的球体内的点云数,根据所述点云数计算k值,计算距离所述基准点最近的k个点,并记录所述k个点中距离所述基准点最远距离d,对所述k个点进行平面拟合;计算距离所述基准点的距离小于2d的点,根据所述平面对距离小于2d的点进行过滤得到k个点,基于过滤得到的k个点计算基准点的法向量;三维重建模块,用于根据所述法向量对点云中的点进行聚类,并基于聚类中点云的点数对聚类过滤,从过滤后的每个聚类得到至少一个集合,对每个集合中的点进行平面拟合得到构件的平面,基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件。7.如权利要求6所述的系统,其特征在于,所述根据所述平面对距离小于2d的点进行过
滤得到k个点,具体为:将距离小于2d的点投影到所述平面上,计算以所述基准点为原点的坐标系中每个象限中投影点数,根据所述投影点数对小于2d的点进行过滤直到剩余k个点,使得剩余k个点在每个象限的投影数的方差最小。8.如权利要求6所述的系统,其特征在于,所述从过滤后的每个聚类得到至少一个集合,具体为:s21,从聚类中选取一个关键点放入集合,迭代执行将聚类中距离集合中任一点的距离小于d的点从聚类中删除并放入集合中,直到聚类中没有点距离集合中任一点的距离小于d;s22,对聚类中剩余点重复执行s21得到一个或多个集合,根据集合中的点数对集合过滤,得到聚类对应的至少一个集合。9.如权利要求8所述的系统,其特征在于,所述每个聚类中选取一个关键点放入集合,具体为:对于聚类的每个点i,计算聚类中其他点和点i的距离小于d的点数,选取对应的点数最大的点作为关键点。10.如权利要求6所述的系统,其特征在于,所述基于所述平面对所述构件进行三维重建得到适用于bim中的格式文件,具体为:获取每个平面的边界点,根据边界点和所述点云剩余的点采用点云曲面重建方法进行三维重建;将重建后的三维文件导出为rfa文件。

技术总结
本发明涉及装配式建筑领域,具体地,提供了一种装配式建筑BIM模型构件生成方法及系统,将每个点分别作为基准点,统计以基准点为球心且半径R的球体内的点云数,根据点云数计算K值,计算距离基准点最近的K个点,并记录K个点中距离基准点最远距离D,对K个点进行平面拟合;计算距离基准点的距离小于2D的点,根据平面对距离小于2D的点进行过滤得到K个点,基于过滤得到的K个点计算基准点的法向量;根据法向量对点云中的点进行聚类,从过滤后的每个聚类得到至少一个集合,对每个集合中的点进行平面拟合得到构件的平面,基于平面对所述构件进行三维重建得到适用于BIM中的格式文件。本发明能够识别构件中的平面,提高构件模型生成速度。度。度。


技术研发人员:尹国安 宋秉权 于本福 刘梦莹 毛小婷
受保护的技术使用者:临沂大学
技术研发日:2023.09.08
技术公布日:2023/10/20
版权声明

本文仅代表作者观点,不代表航家之家立场。
本文系作者授权航家号发表,未经原创作者书面授权,任何单位或个人不得引用、复制、转载、摘编、链接或以其他任何方式复制发表。任何单位或个人在获得书面授权使用航空之家内容时,须注明作者及来源 “航空之家”。如非法使用航空之家的部分或全部内容的,航空之家将依法追究其法律责任。(航空之家官方QQ:2926969996)

航空之家 https://www.aerohome.com.cn/

航空商城 https://mall.aerohome.com.cn/

航空资讯 https://news.aerohome.com.cn/

分享:

扫一扫在手机阅读、分享本文

评论

相关推荐