期刊文献+
共找到30篇文章
< 1 2 >
每页显示 20 50 100
地球空间参考网格系统建设初探 被引量:23
1
作者 程承旗 吴飞龙 +3 位作者 王嵘 秦永刚 童晓冲 陈波 《北京大学学报(自然科学版)》 EI CAS CSCD 北大核心 2016年第6期1041-1049,共9页
为了弥补用经纬度作为地理位置标识存在的描述复杂、无区域特征、计算复杂度高等不足,以北京大学提出的GeoSOT网格为基础,构建一套全球空间参考网格系统。该系统以等度、等分、等秒的完美四叉树剖分网格为基础,建立适合协调空地联合行... 为了弥补用经纬度作为地理位置标识存在的描述复杂、无区域特征、计算复杂度高等不足,以北京大学提出的GeoSOT网格为基础,构建一套全球空间参考网格系统。该系统以等度、等分、等秒的完美四叉树剖分网格为基础,建立适合协调空地联合行动的网格体系,设计了简单实用、可距离计算的定位编码方法,可实现多部门、多源异构地理空间数据的统一检索,并发展出高效编码代数运算、地理空间管理与计算框架、三维地球空间剖分框架等新方法,将在未来的地理空间大数据应用中发挥重要作用。 展开更多
关键词 地球空间参考网格系统 GeoSOT 经纬度 四叉树
在线阅读 下载PDF
多尺度空间填充曲线空间连续性研究 被引量:5
2
作者 翟卫欣 陈波 +1 位作者 童晓冲 程承旗 《北京大学学报(自然科学版)》 EI CAS CSCD 北大核心 2018年第2期331-335,共5页
将二维Hilbert编码和Z编码拓展到以尺度维作为第三维的三维填充曲线:多尺度Hilbert曲线和Z曲线。在多尺度数据条件下,这两种曲线能够提高空间填充曲线的空间连续性,适应多尺度的需求。依托四叉树模型,将多尺度的Hilbert曲线与按照相同... 将二维Hilbert编码和Z编码拓展到以尺度维作为第三维的三维填充曲线:多尺度Hilbert曲线和Z曲线。在多尺度数据条件下,这两种曲线能够提高空间填充曲线的空间连续性,适应多尺度的需求。依托四叉树模型,将多尺度的Hilbert曲线与按照相同思路设计的多尺度Z曲线进行两类对比试验,验证了多尺度Hilbert曲线相对于Z曲线在空间连续性方面的优势,提高的比例在15%~30%之间。 展开更多
关键词 多尺度 HILBERT曲线 空间连续性
在线阅读 下载PDF
基于GPT2w模型化加权平均温度反演可降水量 被引量:6
3
作者 翟树峰 吕志平 +3 位作者 李林阳 吕浩 邝英才 王方超 《大地测量与地球动力学》 CSCD 北大核心 2019年第7期733-737,共5页
提出一种基于GPT2w模型化加权平均温度反演大气可降水量的方法,并分析附加系统偏差改正的模型化加权平均温度对可降水量的影响。结果表明,基于GPT2w模型化加权平均温度反演的大气可降水量的精度与基于Bevis公式计算的加权平均温度反演... 提出一种基于GPT2w模型化加权平均温度反演大气可降水量的方法,并分析附加系统偏差改正的模型化加权平均温度对可降水量的影响。结果表明,基于GPT2w模型化加权平均温度反演的大气可降水量的精度与基于Bevis公式计算的加权平均温度反演的大气可降水量的精度相当;对GPT2w模型化加权平均温度进行系统偏差改正后,大气可降水量的精度有一定改善,但改善率不到1%。 展开更多
关键词 天顶对流层湿延迟 GPT2w 加权平均温度 可降水量
在线阅读 下载PDF
结合向量化和FFT技术的模型扰动引力快速计算 被引量:3
4
作者 邢志斌 李姗姗 +2 位作者 范雕 张驰 马越原 《中国惯性技术学报》 EI CSCD 北大核心 2019年第3期321-326,共6页
用地球重力场模型计算扰动引力矢量时,随着模型阶次的提高,计算效率会迅速降低。针对此问题,研究了扰动引力矢量计算的矩阵表示形式、计算过程的向量化技术以及顾及起始计算点坐标的快速傅里叶变换(FFT)计算表达式。利用EGM2008地球重... 用地球重力场模型计算扰动引力矢量时,随着模型阶次的提高,计算效率会迅速降低。针对此问题,研究了扰动引力矢量计算的矩阵表示形式、计算过程的向量化技术以及顾及起始计算点坐标的快速傅里叶变换(FFT)计算表达式。利用EGM2008地球重力场模型进行数值实验,结果表明:结合向量化和FFT技术可实现快速计算区域和全球格网模型扰动引力矢量,计算全球扰动引力矢量时,当重力场模型展开至2160阶次对应5′×5′的分辨率,计算效率可提高5倍,展开至360 阶次对应30′×30′的分辨率,计算效率可提高近10倍。所提出的方法对于快速计算模型重力场元具有一定的参考作用。 展开更多
关键词 扰动引力矢量 矩阵运算 向量化技术 快速傅里叶变换
在线阅读 下载PDF
面向越野通行分析的土壤数据分类方法研究 被引量:1
5
作者 李坤伟 游雄 +1 位作者 张欣 汤奋 《系统仿真学报》 CAS CSCD 北大核心 2019年第1期158-165,共8页
土壤是影响部队越野机动的重要因素,土壤与气象因素相结合使得越野通行性分析变得异常复杂。基于统一土壤分类系统(Unified Soil Classification System,USCS)的土壤数据是定量分析土壤通行性的基础,采用随机森林方法,利用土壤的多种属... 土壤是影响部队越野机动的重要因素,土壤与气象因素相结合使得越野通行性分析变得异常复杂。基于统一土壤分类系统(Unified Soil Classification System,USCS)的土壤数据是定量分析土壤通行性的基础,采用随机森林方法,利用土壤的多种属性信息来预测土壤的USCS类型。从已有的USCS土壤数据中提取样本数据构建8种随机森林模型,对不同随机森林模型的精度和特征变量的重要性进行分析,根据土壤数据集的特点采用第三种随机森林模型对其进行处理,构建适合越野通行分析的土壤数据。相比以往的方法,这种方法精度更高,更能满足越野通行分析的需要。 展开更多
关键词 越野机动 土壤数据 随机森林 圆锥指数
在线阅读 下载PDF
地球重力场模型计算点向径取值对模型重力场元影响 被引量:1
6
作者 邢志斌 李姗姗 +1 位作者 范雕 张驰 《中国惯性技术学报》 EI CSCD 北大核心 2019年第2期167-175,共9页
利用地球重力场模型计算重力场元时,为提高计算速度,有时将计算点的高程近似为0,为提高精度,在球面对计算点高程进行泰勒级数展开。但是,考虑到地球表面真实情况,这种近似将会引入巨大的距离误差,同时也会影响计算点的地心纬度值。针对... 利用地球重力场模型计算重力场元时,为提高计算速度,有时将计算点的高程近似为0,为提高精度,在球面对计算点高程进行泰勒级数展开。但是,考虑到地球表面真实情况,这种近似将会引入巨大的距离误差,同时也会影响计算点的地心纬度值。针对此,分析了计算点向径不同取值之间的差异及对计算点地心纬度、模型垂线偏差计算结果的影响,引入等效高程讨论了按照泰勒级数展开快速计算模型垂线偏差方法的适用性,采用实测数据比较了计算点向径不同取值计算模型垂线偏差的精度。实验结果表明:将计算点向径近似为地球长半轴与计算点高程之和,造成的距离误差可达20km,而近似为参考椭球面向径与地面点高程之和造成的距离误差最大不超过6cm,对模型垂线偏差的影响则不足0.05″,这为精确快速计算模型重力场元提供了参考。 展开更多
关键词 地球重力场模型 向径 模型垂线偏差 等效高程
在线阅读 下载PDF
基于三维Taylor多项式建立局域海洋地磁场模型 被引量:2
7
作者 种洋 柴洪洲 《大地测量与地球动力学》 CSCD 北大核心 2016年第10期854-858,共5页
在海洋磁力测量中,波浪和风等外界因素的变化会对磁力测量数据在垂直空间上产生影响,而二维Taylor多项式在建立局域海洋地磁场模型时没有考虑高度变化对模型的影响。针对这一问题,提出使用三维Taylor多项式完全展开法来构建局域海洋地... 在海洋磁力测量中,波浪和风等外界因素的变化会对磁力测量数据在垂直空间上产生影响,而二维Taylor多项式在建立局域海洋地磁场模型时没有考虑高度变化对模型的影响。针对这一问题,提出使用三维Taylor多项式完全展开法来构建局域海洋地磁场模型,通过模拟某海域水下球体磁场,仿真船载磁力测量数据,分析两种模型截断阶数的选择、完全与不完全展开的区别和模型边界效应等问题。实验表明,在顾及高度变化的情况下,三维Taylor多项式模型能更详细地反映海洋磁力异常的分布情况,且较好地改善了边界效应。 展开更多
关键词 局域海洋地磁场模型 三维Taylor多项式 截断阶数 边界效应
在线阅读 下载PDF
福建省CORS系统在沿海大地控制网测量中的应用 被引量:1
8
作者 欧阳明达 李军 张海东 《海洋测绘》 2013年第5期27-30,共4页
福建省CORS系统是地球空间信息服务的重要基础设施,是当代GNSS发展的热点之一。结合工程实际,并综合利用陆态网络观测数据,设计出三种不同方案,对沿海大地控制网进行解算。结果表明,三种方案均具有较好的可靠性,坐标中误差数据显示,由... 福建省CORS系统是地球空间信息服务的重要基础设施,是当代GNSS发展的热点之一。结合工程实际,并综合利用陆态网络观测数据,设计出三种不同方案,对沿海大地控制网进行解算。结果表明,三种方案均具有较好的可靠性,坐标中误差数据显示,由于利用1σ约束CORS站点作为控制的解算方案,大大缩短了基线长度,并使图形结构更加合理,能够获得较强约束陆态网络基准站作为控制方案更好的测量精度。 展开更多
关键词 CORS GAMIT GLOBK 大地控制网
在线阅读 下载PDF
基于Filtersim的多源降水数据融合方法研究
9
作者 李坤伟 游雄 +1 位作者 张欣 汤奋 《系统仿真学报》 CAS CSCD 北大核心 2019年第6期1232-1238,共7页
降雨是影响车辆越野机动的重要因素,高精度的降雨数据是定量评估战场越野通行能力,合理制定行动方案的前提。针对传统插值方法的不足,论文采用flitersim多点统计法来提高降水数据的融合精度。将降雨量分解为局部均值与局部残差之和,以... 降雨是影响车辆越野机动的重要因素,高精度的降雨数据是定量评估战场越野通行能力,合理制定行动方案的前提。针对传统插值方法的不足,论文采用flitersim多点统计法来提高降水数据的融合精度。将降雨量分解为局部均值与局部残差之和,以卫星降水数据的局部残差作为"训练数据",以气象站点的局部残差数据作为"硬数据",以卫星降水的局部残差数据和气象站点的局部残差数据的插值数据作为"软数据"进行flitersim多点地统计模拟,最终得到1km分辨率的降水数据。通过比较数据的平均绝对值误差、均方根误差和相关系数,表明flitersim明显优于协同克里金插值,能有效提高降水数据的精度。 展开更多
关键词 降水数据 融合方法 均方根误差 插值方法 克里金插值 高精度 降雨量 通行能力
在线阅读 下载PDF
基于HurberM极大似然估计求解的CKF滤波定位算法
10
作者 戴卿 《控制工程》 CSCD 北大核心 2018年第9期1760-1764,共5页
针对传统移动机器人定位算法精度欠佳的问题,设计了一种基于无线传感器网络滤波算法的移动机器人定位算法(HM_CKF)。该算法利用HurberM极大似然估计代价函数,将线性化后CKF观测矩阵求出,从而解决CKF滤波算法在未知非高斯白噪声干扰下... 针对传统移动机器人定位算法精度欠佳的问题,设计了一种基于无线传感器网络滤波算法的移动机器人定位算法(HM_CKF)。该算法利用HurberM极大似然估计代价函数,将线性化后CKF观测矩阵求出,从而解决CKF滤波算法在未知非高斯白噪声干扰下估计精度欠佳的问题。然后在体育馆中基于WSNs网络构建了移动机器人定位实验环境,融合移动机器人动力学模型,对比了HM_CKF、CKF算法的定位精度。实验结果表明,在不含噪声干扰和含未知噪声干扰的两种情况下,HM_CKF算法定位精度依次比CKF算法提高了7%和15%。 展开更多
关键词 无线传感器网络 HurberM代价函数 CKF滤波 移动机器人 定位
在线阅读 下载PDF
利用迁移学习的机载激光雷达点云分类 被引量:31
11
作者 赵传 张保明 +2 位作者 余东行 郭海涛 卢俊 《光学精密工程》 EI CAS CSCD 北大核心 2019年第7期1601-1612,共12页
为解决现有机载激光雷达点云分类方法存在难以在获得较高精度的点云分类结果的同时降低分类过程所需时间等问题,提出了一种利用迁移学习的机载激光雷达点云分类方法。首先,计算点云的归一化高程、强度和法向量三个特征,通过设置不同邻... 为解决现有机载激光雷达点云分类方法存在难以在获得较高精度的点云分类结果的同时降低分类过程所需时间等问题,提出了一种利用迁移学习的机载激光雷达点云分类方法。首先,计算点云的归一化高程、强度和法向量三个特征,通过设置不同邻域大小,利用所提出的点云特征图生成策略生成多尺度点云特征图;然后,利用预训练的深度残差网络从每个点的多尺度点云特征图提取其多尺度深度特征;最后,为了实现快速地训练,构建仅包含两层全连接神经网络模型,再利用训练好的模型对点云进行分类。两组ISPRS提供的标准点云数据集的试验结果表明:提出的方法所需训练时间少,分类结果的整体精度为89.6%,较ISPRS官网上所报道的最佳点云方法分类精度高4.4%。分类结果可为机载激光雷达点云的后续处理与应用提供可靠的信息。 展开更多
关键词 迁移学习 特征图 分类 机载激光雷达点云
在线阅读 下载PDF
基于BP神经网络的多特征融合变化检测方法 被引量:11
12
作者 郝睿 徐俊峰 +1 位作者 王庆宝 张保明 《海洋测绘》 CSCD 2016年第1期79-82,共4页
当前,随着遥感影像空间分辨率的提高,传统的变化检测方法已经无法满足实际应用的需要。针对这一问题,引入BP神经网络,提出了一种多特征融合的遥感影像变化检测方法。在提取影像光谱特征、纹理特征和空间特征的基础上,构建了用于进行二... 当前,随着遥感影像空间分辨率的提高,传统的变化检测方法已经无法满足实际应用的需要。针对这一问题,引入BP神经网络,提出了一种多特征融合的遥感影像变化检测方法。在提取影像光谱特征、纹理特征和空间特征的基础上,构建了用于进行二类分类的BP神经网络模型,并在一定数量的样本支持下,有效完成了变化信息的提取。实验表明,文中方法相比于传统的影像代数法和仅使用对象光谱特征的方法,在保持较高正确率的情况下,有效降低了漏检率和虚警率。 展开更多
关键词 变化检测 多特征融合 BP神经网络 高空间分辨率 多光谱影像
在线阅读 下载PDF
网络电子地图多尺度显示的分析与研究 被引量:8
13
作者 张昊 江南 曹亚妮 《测绘通报》 CSCD 北大核心 2014年第1期90-93,124,共5页
通过具体试验获取大量数据,并对比分析多个地图网站的多尺度模型,从不同网站的尺度级别设置、重点表达区域、载负量变化规律、同一区域在不同网站表达的详略程度等方面入手,探索网络电子地图载负量随尺度缩放的变化规律,并提出适宜人眼... 通过具体试验获取大量数据,并对比分析多个地图网站的多尺度模型,从不同网站的尺度级别设置、重点表达区域、载负量变化规律、同一区域在不同网站表达的详略程度等方面入手,探索网络电子地图载负量随尺度缩放的变化规律,并提出适宜人眼视觉感受的多尺度模型计算方法。 展开更多
关键词 电子地图 多尺度 载负量 关键比例尺
在线阅读 下载PDF
面向新闻推荐的用户兴趣模型构建与更新 被引量:8
14
作者 袁仁进 陈刚 李锋 《计算机应用研究》 CSCD 北大核心 2019年第12期3593-3596,共4页
针对新闻推荐系统中用户兴趣模型构建与用户兴趣漂移问题,提出了一种面向新闻推荐的用户兴趣模型构建与更新方法。首先采用向量空间模型与bisecting K-means聚类算法构建了原始用户兴趣模型;然后以艾宾浩斯遗忘曲线为基础构造了遗忘函数... 针对新闻推荐系统中用户兴趣模型构建与用户兴趣漂移问题,提出了一种面向新闻推荐的用户兴趣模型构建与更新方法。首先采用向量空间模型与bisecting K-means聚类算法构建了原始用户兴趣模型;然后以艾宾浩斯遗忘曲线为基础构造了遗忘函数,并以此对用户兴趣模型进行时间加权,从而达到对用户兴趣模型更新的目的。实验以基于用户的协同过滤推荐、基于物品的协同过滤推荐为baseline,实验结果表明所构建的原始用户兴趣模型推荐性能更优,在F值上提升了4%,更新后的模型与原始模型相比F值提高了1. 3%。 展开更多
关键词 个性化推荐 向量空间模型 用户兴趣模型 用户兴趣漂移 遗忘函数
在线阅读 下载PDF
卫星导航欺骗干扰装备效能评估方法 被引量:8
15
作者 王月 郝金明 刘伟平 《兵工学报》 EI CAS CSCD 北大核心 2020年第1期108-118,共11页
作为“导航战”重要工具之一的卫星导航欺骗干扰装备,其干扰效能直接关系到被保护目标的生存概率。为给出置信度较高的效能评估,有利于装备高效应用于实际,在合理的评估指标体系和评估试验平台基础上,通过建立指标的检测方法及计算模型... 作为“导航战”重要工具之一的卫星导航欺骗干扰装备,其干扰效能直接关系到被保护目标的生存概率。为给出置信度较高的效能评估,有利于装备高效应用于实际,在合理的评估指标体系和评估试验平台基础上,通过建立指标的检测方法及计算模型,得到理论评估方法;采用模糊综合评价思路、半定性和半定量的等级量化定权方法及加权积的算法,得到理论评估阈值;以全数字、半实物仿真和全实物外场试验,分析效能指标对装备干扰效果的具体影响并得到结论。研究结果表明,缩短信号捕获时间,适宜增加干扰信号功率等方式都可全面提高装备干扰效果。 展开更多
关键词 导航战 欺骗式干扰 指标体系 模糊综合评价 干扰效能测试
在线阅读 下载PDF
基于RegEM算法的GPS坐标时间序列插值应用分析 被引量:9
16
作者 王方超 吕志平 +1 位作者 吕浩 邝英才 《大地测量与地球动力学》 CSCD 北大核心 2020年第1期45-50,共6页
将一种基于数据驱动的RegEM算法引入GPS坐标时间序列插值中,分别采用模拟不同比例连续缺失数据与实测含缺失数据,比较RegEM与拉格朗日方法、三次样条方法、正交多项式方法的插值效果与性能。结果表明,对于模拟不同比例连续缺失数据,RegE... 将一种基于数据驱动的RegEM算法引入GPS坐标时间序列插值中,分别采用模拟不同比例连续缺失数据与实测含缺失数据,比较RegEM与拉格朗日方法、三次样条方法、正交多项式方法的插值效果与性能。结果表明,对于模拟不同比例连续缺失数据,RegEM算法插值效果均优于传统方法,且在大量数据连续缺失的情况下效果最优;对于实测含缺失数据,RegEM方法插值所得序列保留方差最大化效果最好,约为正交多项式方法的1.17倍、三次样条方法的1.38倍。 展开更多
关键词 GPS坐标时间序列 插值 数据驱动模型 RegEM
在线阅读 下载PDF
基于路况数据的城市道路交通事件点检测 被引量:5
17
作者 孙梦婷 魏海平 +1 位作者 李星滢 徐立 《地理与地理信息科学》 CSCD 北大核心 2019年第6期9-14,共6页
城市道路交通事件检测是道路交通智慧应急系统的关键技术之一。针对当前难以自动检测交通事件问题,提出了一种利用路况数据发现造成交通拥堵的交通事件检测方法。首先,分析空间维和时间维下的路况异常模式,并利用交通事件在时间和空间... 城市道路交通事件检测是道路交通智慧应急系统的关键技术之一。针对当前难以自动检测交通事件问题,提出了一种利用路况数据发现造成交通拥堵的交通事件检测方法。首先,分析空间维和时间维下的路况异常模式,并利用交通事件在时间和空间上的相关性提取交通事件点;进一步计算交通事件的持续时间,根据持续时间的长短对交通事件划分预警等级,为交通管理部门的应急部署提供决策支持。运用郑州市252个拥堵路段的路况数据进行实验分析,检测交通事件点的完整率为84.85%。经验证,在检测高峰期城市快速路和主干道的交通事件点时,该算法具有较高的检测率,对高峰期和平峰期城市快速路的交通事件误报率均较低。 展开更多
关键词 交通事件检测 交通拥堵 路况时空序列 路况异常模式 预警分级
在线阅读 下载PDF
基于Hough变换的DInSAR环状地表形变区自动提取 被引量:3
18
作者 徐娇 靳国旺 +1 位作者 熊新 张红敏 《海洋测绘》 CSCD 2019年第3期45-49,共5页
针对当前DInSAR地表形变区检测依赖于人工的现状,提出基于Hough变换的环状地表形变区自动提取方法。该方法在外部DEM的辅助下,由干涉数据获取去地形效应的差分干涉图,在对差分干涉图进行边缘提取的基础上,利用Hough变换实现环状地表形... 针对当前DInSAR地表形变区检测依赖于人工的现状,提出基于Hough变换的环状地表形变区自动提取方法。该方法在外部DEM的辅助下,由干涉数据获取去地形效应的差分干涉图,在对差分干涉图进行边缘提取的基础上,利用Hough变换实现环状地表形变区的自动提取。采用青海门源震区Sentinel-1A数据和山东济宁矿区ALOS PALSAR数据分别进行环状地表形变区自动提取实验,验证了本文方法的有效性。 展开更多
关键词 合成孔径雷达差分干涉测量 地表形变 边缘提取 HOUGH变换 形变区提取
在线阅读 下载PDF
无人机影像中直线类电力线检测算法研究 被引量:9
19
作者 李朋月 郝向阳 +2 位作者 李建胜 杨戬峰 王高杰 《电子测量技术》 2019年第19期148-153,共6页
随着无人机技术的高速发展,无人机巡航环境中的障碍物检测成为热点。针对目前复杂环境下直线类障碍物检测准确率低、精度差的问题,提出了一种基于像素梯度和距离编组的直线类电力线检测算法。首先结合图像灰度特征提取图像直线边缘点,... 随着无人机技术的高速发展,无人机巡航环境中的障碍物检测成为热点。针对目前复杂环境下直线类障碍物检测准确率低、精度差的问题,提出了一种基于像素梯度和距离编组的直线类电力线检测算法。首先结合图像灰度特征提取图像直线边缘点,根据图像电力线特征构建直线支撑区域并对其进行调整,为改善图像电力线检测断裂的问题,最后根据直线间的距离度量完成直线段的编组和拟合,实现复杂背景环境中直线类电力线的检测。通过实验证明,算法有效解决了无人机复杂背景环境电力线图像中直线的误检、断裂与漏检问题,降低了光照变化对电力线的检测影响,提高了电力线检测的鲁棒性。 展开更多
关键词 电力线检测 灰度特征 边缘检测 直线编组 无人机影像
在线阅读 下载PDF
BDS三频非差周跳探测与修复 被引量:3
20
作者 王赛 刘长建 +2 位作者 田翌君 杜莹 许岭峰 《大地测量与地球动力学》 CSCD 北大核心 2016年第6期539-543,共5页
通过比较伪距/载波组合与相位组合阶段取整两种方法在考虑电离层延迟下的综合噪声,选择噪声最小的组合(1,-3,2)、系数之和不为0时噪声最小的组合(4,-3,-2)以及组合(0,1,-1)对北斗三频周跳进行探测与修复。使用IGMAS提供的BDS三频数据,针... 通过比较伪距/载波组合与相位组合阶段取整两种方法在考虑电离层延迟下的综合噪声,选择噪声最小的组合(1,-3,2)、系数之和不为0时噪声最小的组合(4,-3,-2)以及组合(0,1,-1)对北斗三频周跳进行探测与修复。使用IGMAS提供的BDS三频数据,针对BDS 3种不同星座分别进行分析比较,验证了此方法可实时准确探测出非差的各类大小周跳。 展开更多
关键词 BDS 三频 相位组合 伪距/载波组合 超宽巷
在线阅读 下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部