网站建设要钱么,泰安房产网58同城,wordpress 导航栏图标,wordpress 附件清理DolphinDB 集高性能时序数据库与全面的分析功能为一体#xff0c;可用于海量结构化数据的存储、查询、分析、实时计算等#xff0c;在工业物联网场景中应用广泛。本文以纽约出租车行程时间预测为例#xff0c;介绍如何使用 DolphinDB 训练机器学习模型#xff0c;并进行实时…DolphinDB 集高性能时序数据库与全面的分析功能为一体可用于海量结构化数据的存储、查询、分析、实时计算等在工业物联网场景中应用广泛。本文以纽约出租车行程时间预测为例介绍如何使用 DolphinDB 训练机器学习模型并进行实时数据的预测为基于智能网联汽车的车联网企业提供基于机器学习方法的即时预测方案。 1. 概要
随着手机移动应用与网约车平台的迅速发展网约车出行逐渐成为城市生活中一种重要的出行方式。相较其他出行方式而言选择网约车的乘客对出行时效性有更高的要求本文将基于乘客上车时间及上下车地点等静态信息使用 DolphinDB 机器学习方法训练模型预测网约车行程时间。
在此基础上本文将介绍如何使用 DolphinDB 流数据处理系统对业务系统产生的持续增长的网约车订单动态数据进行实时的收集、清洗、统计、入库并实时展示行程时间预测结果。 行程时间实时预测流程
2. 数据介绍
2.1 数据来源及训练方法
本文训练和预测采用 Kaggle 提供的来自纽约出租车委员会的数据集训练方法参考了获奖者 beluga 的模型使用 DolphinDB 对原始数据进行数据预处理完成位置信息主成分分析PCA, Principal Component Analysis、位置信息聚类KMeans、新特征构建等工作并使用 DolphinDB XGBoost 插件完成模型训练及行程时间预测。
为对比 DolphinDB 在机器学习上的性能本文使用 Python Scikit-Learn 库及 XGBoost 在同一环境下进行了模型训练和预测DolphinDB 在训练耗时、模型精度等方面均有良好表现。
2.2 数据特征
该数据集预先分为训练数据集及测试数据集训练数据集共包含1458644条数据测试数据集共包括625134条数据训练数据集共包含以下11列信息。
列名列类型说明实例idSYMBOL行程的唯一标识id2875421vendor_idINT行程记录提供商代码2pickup_datetimeDATETIME出租车计价器开启时间2016/3/14 17:24:55dropoff_datetimeDATETIME出租车计价器关闭时间2016/3/14 17:32:30passenger_countINT乘客数量1pickup_longitudeDOUBLE出租车计价器开启位置经度-73.98215484619139pickup_latitudeDOUBLE出租车计价器开启位置纬度40.76793670654297dropoff_longitudeDOUBLE出租车计价器关闭位置经度-73.96463012695312dropoff_latitudeDOUBLE出租车计价器关闭位置纬度40.765602111816406store_and_fwd_flagCHAR标识来源是否为存储的历史数据Ntrip_durationINT行程时间按秒计455
行程时间预测的目标列为上表中 trip_duration 列即 dropoff_datetime 与 pickup_datetime 之差。测试数据集用于预测故其列信息不包括 dropoff_datetime 及 trip_duration 列测试数据集中行程标识、位置等列属性同上表。
上表的数据类型中SYMBOL 类型是 DolphinDB 中一种特殊的字符串类型在系统内部的存储结构为一个编码字典DATETIME 类型为包含了日期和时刻的时间类型。
DolphinDB 支持 loadText 方法读取 csv 等数据存储文件到内存表用户可以 schema 函数获取表的特征信息。DolphinDB 也支持使用 SQL 语句完成数据的查询。
train loadText(./taxidata/train.csv)
train.schema().colDefs
select count(*) from train
select top 5 * from train
2.3 数据存储
将数据加载到内存表后可以将训练数据与测试数据导入 DolphinDB 数据库中便于后续数据的读取与模型的训练数据导入分布式数据库的操作详见 database.md · dolphindb/Tutorials_CN - Gitee。
3. 模型构建
本节介绍行程时间预测模型的构建方法。
行程时间预测模型的构建分多个过程一是预处理原始数据对可能存在的空值进行转换并将字符等非数值型数据转换为可用于模型训练的数值型数据二是优化位置信息原始数据中的纬度经度信息集中在40.70 °N 至40.80 °N 及73.94 °W 至74.02 °W之间数据间位置特征差异不够显著使用主成分分析、聚类方法处理可以提取到特征更明显的信息三是新特征的构建位置信息和时间信息是订单数据的两个关键维度通过计算可以在位置信息基础上得到方位、距离信息提取更多空间特征而组合不同类别的位置信息和时间信息也可以得到更复杂的特征有利于模型学习深层次的时空规律。
3.1 数据预处理
在模型训练过程中首先需要检查数据集是否包含空值本训练数据集与测试数据集均不包含空值若存在缺失值还需要删除、插补等操作解决缺失数据问题。
其次需要检查数据集数据类型原始数据往往包含文本/字符数据由1.3节表可知本数据集中 store_and_fwd_flag 列为字符型数据pickup_datetime 及 dropoff_datetime 列为日期时间类型数据为充分利用这些信息训练模型需要将其转化为数值型数据。
此外考虑到该数据测试集评价指标为均方根对数误差Root Mean Squared Logarithmic Error, RMSLE同时最大行程时间接近1000小时离群值会影响模型训练效果对行程时间取对数作为预测值在评价时见3.6节可以直接使用均方根误差Root Mean Squared Error, RMSE指标。 RMSE
DolphinDB 提供多种计算函数可以帮助用户快速实现数据处理。DolphinDB 提供 isNull() 方法用于判断空值配合 sum() 等聚合函数使用可以快速完成整表数据的查询提供类似于条件运算符的 iif() 方法简化 if-else 语句date()、weekday()、hour() 等方法可以提取时间、日期数据的不同特征简洁高效类似于 Python 等编程语言DolphinDB 支持方括号([])索引简化了表的查找、更新和插入。
sum(isNull(train)) // 0不含空值
trainData[store_and_fwd_flag_int] iif(trainData[store_and_fwd_flag] N, int(0), int(1)) // 将字符N/Y转化为0/1值
trainData[pickup_date] date(trainData[pickup_datetime]) // 日期
trainData[pickup_weekday] weekday(trainData[pickup_datetime]) // 星期*
trainData[pickup_hour] hour(trainData[pickup_datetime]) // 小时
trainData[log_trip_duration] log(double(trainData[trip_duration]) 1)// 对行程时间取对数log(trip_duration1)
select max(trip_duration / 3600) from trainData // 训练集上最大行程时间为979h
3.2 位置信息主成分分析PCA
原始数据中的纬度经度信息集中在40.70 °N 至40.80 °N 及73.94 °W 至74.02 °W之间数据间位置特征差异不够显著使用 PCA 来转换经度和纬度坐标有助于 XGBoost 决策树的拆分DolphinDB PCA 函数使用详见 pca — DolphinDB 2.0 documentation。
DolphinDB PCA 返回的结果是一个字典包含 components、explainedVarianceRatio、singularValues 三个键分别代表对应大小为 size(colNames)*k 的主成分分析矩阵、前k个主成分每个特征的方差贡献率、主成分方差协方差矩阵特征值。可通过主成分分析矩阵转换待处理数据详见 Scikit-Learn PCA.transform()。
可从中取若干数据绘制经度-纬度散点图观察 PCA 结果。
经处理位置坐标分散在原点附近。 PCA前上客位置信息 PCA 后上客位置信息
pca() 接收一个或多个数据源为参数对指定列中的数据进行主成分分析用户可通过 table() 方法创建内存表用于 PCADolphinDB 也提供了 dot()、repmat() 等矩阵乘法、矩阵堆叠方法用户可使用内置函数快速完成矩阵运算处理位置信息。
PCApara table(trainData[pickup_latitude] as latitude, trainData[pickup_longitude] as longitude)
pca_model pca(sqlDS(select * from PCApara)) // 使用PCA计算数据集属性
pca_trainpick dot((matrix(trainPickPara) - repmat(matrix(avg(trainPickPara)), train_n, 1)), pca_model.components) // transform
trainData[pca_trainpick_0] flatten(pca_trainpick[:, 0])
DolphinDB 提供了plot函数供数据可视化。用户可通过chartType指定图表类型详见plot — DolphinDB 2.0 documentation。
x select top 1000 pca_trainpick_1 from trainData
y flatten(matrix(select top 1000 pca_trainpick_0 from trainData))
plot(x, y, chartTypeSCATTER)
3.3 位置信息聚类KMeans
原始数据位置数据规模庞大很难挖掘多条数据间的共同特征。KMeans 可以将经纬度相近的数据点归为同一个簇有助于更好地归纳组内数据特征。本模型指定要生成的聚类数为100质心最大迭代次数为100选择 KMeans 算法生成模型DolphinDB kmeans 可选参数及含义详见 kmeans — DolphinDB 2.0 documentation。
可以使用条形图观察聚类后的数据分布。 KMeans 聚类结果
kmeans() 接收一个表作为训练集。对于机器学习函数生成的模型DolphinDB 提供 saveModel 方法将模型保存到本地文件中用于后续预测用户可指定服务器端输出文件的绝对路径或相对路径DolphinDB 也提供了 predict 方法调用已训练的特定模型对相同表结构的测试集数据进行预测。
kmeans_set PCApara[rand(size(PCApara)-1, 500000)] // 随机选取500000数据用于聚类
kmeans_model kmeans(kmeans_set, 100, maxIter100, initk-means) // KMeans
saveModel(kmeans_model, ./taxidata/KMeans.model) // 保存模型训练结果
trainData[pickup_cluster] kmeans_model.predict(select pickup_latitude, pickup_longitude from trainData)
saveModel 和 predict 函数的使用方法可以参考
saveModel — DolphinDB 2.0 documentationpredict — DolphinDB 2.0 documentation
3.4 新特征构建
原始数据仅提供了经度及纬度位置信息可在此基础上增加位置特征如地球表面两经纬度点之间距离、两经纬度点之间的 Manhattan 距离、两个经纬度之间的方位信息等等。
地球表面两点间的距离可使用 haversine 公式精确得到而在本数据集中网约车实际驶过的往往是水平或竖直的街道所组成的路径Manhattan 距离也称城市街区距离标明两个点在标准坐标系上的绝对轴距总和可能能够更精确地反映实际的行驶距离。
在此基础上考虑到训练集包含了完整的时间信息还可以在训练集上添加速度特征。训练集上的速度特征无法直接用于测试集但在相同的位置聚类属性或相同的时间日期特征下行程所需时间及平均速度可能存在某些共性如郊区或凌晨时分网约车车速偏大而城区及早晚高峰时段网约车车速偏小可以将这种训练集上发现的经验应用到测试集上。可将数据按聚类属性或时间特征分组统计组内数据平均速度等特征组合生成新特征合并到测试集的相应分组中。
距离、方位的计算参数不同而方法相同DolphinDB 支持用户自定义函数通过独立的代码模块完成特定的计算任务。而对于不同类别聚类、时间内的特征可使用 groupby 方法在每个分组中计算需要的特征如平均值。groupby 接收三个参数将根据第三个参数指定的列进行分组取第一个参数为计算函数计算第二个参数对应列的特征并返回行数与分组数相等的表。用户可通过表连接操作将该组合特征合并入特征数据本文使用 fj(full join) 将特征表与 groupby 表合并fj() 指定第三个参数为连接列将前两个参数所传入的表合并。
// 两经纬度点距离、两个经纬度之间的 Manhattan 距离、两个经纬度之间的方位信息
trainData[distance_haversine] haversine_array(trainData[pickup_latitude], trainData[pickup_longitude], trainData[dropoff_latitude], trainData[dropoff_longitude])
trainData[distance_dummy_manhattan] dummy_manhattan_distance(trainData[pickup_latitude], trainData[pickup_longitude], trainData[dropoff_latitude], trainData[dropoff_longitude])
trainData[direction] bearing_array(trainData[pickup_latitude], trainData[pickup_longitude], trainData[dropoff_latitude], trainData[dropoff_longitude])
// 按时间、聚类等信息处理速度、行程时间产生新特征
for(gby_col in [pickup_hour, pickup_date, pickup_week_hour, pickup_cluster, dropoff_cluster]) {for(gby_para in [avg_speed_h, avg_speed_m, log_trip_duration]) {gby groupby(avg, trainData[gby_para], trainData[gby_col])gby.rename!(avg_ gby_para, gby_para _gby_ gby_col)trainData fj(trainData, gby, gby_col)testData fj(testData, gby, gby_col)}trainData.dropColumns!(gby gby_col)
}
3.5 模型训练XGBoost
在进行训练之前需要再一次检查训练集和测试集的数据需要剔除 ID、日期、字符等非数值型数据以及平均速度、行驶时间等仅在训练集上存在的数据保证训练数据与预测数据结构一致。
完成数据处理及特征构建后可以使用 XGBoost 等机器学习方法训练模型。为评价模型训练效果将训练数据集划分为训练集和验证集随机选取80%的数据作为训练集训练模型使用20%的数据作为验证集输出预测结果使用均方根误差指标计算验证集的预测值与真实值的偏差。最终可以在测试集上输出行程时间的预测结果。
本模型在验证集上的均方根误差为0.390可以绘制预测值-真值散点图定性分析模型预测效果。 验证集上的预测值与真实值
DolphinDB 提供了 XGBoost 插件实现模型训练及预测使用前需要下载插件到指定路径并加载 XGBoost 插件。DolphinDB XGBoost 插件使用详见 xgboost/README_CN.md · dolphindb/DolphinDBPlugin - Gitee。
xgb_pars {min_child_weight: 50, eta: 0.3, colsample_bytree: 0.3, max_depth: 10,subsample: 0.8, lambda: 1., nthread: 4, booster : gbtree, silent: 1,eval_metric: rmse, objective: reg:linear, nthread: 48} // xgb 参数设置
xgbModel xgboost::train(ytrain, train, xgb_pars, 60) // 训练模型
yvalid_ xgboost::predict(xgbModel, valid) // 使用模型进行预测
3.6 模型评价
为实现出租车行程时间的预测本文使用了三种机器学习方法。首先使用 PCA 对位置信息进行处理转换数据的经纬度特征使用 KMeans 对出租车上下客位置进行聚类将纽约市区划分为100个区域进行分析最后使用 XGBoost 对数据集特征进行训练。在验证集上模型均方根误差为0.390效果较好。
Python Scikit-Learn 也是主流的机器学习库之一本文在相同环境下使用 Python 对同一数据集进行训练PCA、KMeans、XGBoost 训练耗时如下表所示
模型DolphinDBPythonPCA0.325s0.396sKMeans45.711s104.568sXGBoost57.269s74.289s
DolphinDB 与 Python 训练模型在验证集上误差如下表所示
DolphinDBPythonRMSE0.3900.394
在本行程时间预测任务中在准确率上DolphinDB 与 Python 相近而在性能上DolphinDB 在 PCA、KMeans、XGBoost 上的性能均优于 Python。
4. 行程时间实时预测
本节结合现实场景介绍如何使用 DolphinDB 处理实时的订单流数据基于预测模型实时估计行程时间。
现实场景中网约车乘客对时效性要求高需要平台提供准确的行程时间估计而服务商也需要监控出行平台分析出行需求并完成资源调度仅仅使用预测模型无法高效处理实时数据难以完成即时预测任务无法满足乘客和服务商的实时需求。DolphinDB 流数据模块可以解决生产环境下实时数据的快速分析计算问题对服务商发送的即时数据DolphinDB 流数据引擎可高效完成数据预处理、信息提取、特征构建等工作使用预先训练的模型完成实时订单行程时间的快速准确预测为用户提供从模型训练、流数据注入到实时预测及在线监控的一站式解决方案。 DolphinDB 流数据处理框架
4.1 场景描述
DolphinDB 流数据模块采用发布-订阅-消费的模式流数据首先注入流数据表中通过流表来发布数据第三方应用可以通过 DolphinDB 脚本或 API 订阅及消费流数据。
为实现出租车行程时间的实时预测服务商可以创建 DolphinDB 流数据表订阅服务端消息获取乘客创建的行程信息使用离线训练完成的模型对行程时间进行实时预测最后可通过应用程序订阅预测数据并提供给乘客。
4.2 实时数据模拟及预测
为获取行程数据并使用机器学习模型预测行程时间用户需要创建三个流表实现实时预测一是创建订单信息表订阅乘客行程信息二是创建特征表完成对订单信息的特征提取三是创建预测表预测特征流数据表发送的行程特征信息输出预测结果。 用户可以使用 subscribeTable 完成流数据的订阅并通过 handler 指定处理订阅数据的方法详见subscribeTable — DolphinDB 2.0 documentation。在本例中特征表需订阅订单表完成原始信息的特征提取本模型定义 process 函数实现预测表需订阅特征表使用特征信息完成行程时间预测本模型定义 predictDuration 函数实现。函数实现详见6.2节所附代码。
为模拟实时数据使用 replay 函数回放历史数据。
// 订阅订单信息表数据从订单表流向特征表
subscribeTable(tableNameorderTable, actionNameorderProcess, offset0, handlerprocess{traitTable, hisData}, msgAsTabletrue, batchSize1, throttle1, hash0, reconnecttrue)
// 订阅特征表数据从特征表流向预测表
subscribeTable(tableNametraitTable, actionNamepredict, offset0, handlerpredictDuration{predictTable}, msgAsTabletrue, hash1, reconnecttrue)
// 回放历史数据模拟实时产生的生产数据
submitJob(replay, trade, replay{inputTablesdata, outputTablesorderTable, dateColumnpickup_datetime, timeColumnpickup_datetime, replayRate25, absoluteRatetrue, parallelLevel1})
4.3 Grafana 实时监控
服务商可通过第三方 API 连接 DolphinDB 数据库监控行程时间预测服务本文以 Grafana 为例简要介绍如何使用第三方应用程序动态展示实时数据。
Grafana 是一个用于时序数据动态可视化的数据展示工具DolphinDB 提供了 Grafana 的数据接口用户可在 Grafana 面板上编写查询脚本与 DolphinDB 进行交互实现 DolphinDB 时序数据的可视化并 Web 端进行实时数据分析详见 README.zh.md · dolphindb/grafana-datasource - Gitee。
添加 datasource 并新建 dashboard 后在 Query 中填写以下 DolphinDB 语句进行实时数据可视化
Query 1展示当日出行订单预估到达时间及预估行程时间
select id as ID, pickup_datetime as pickup_time, (pickup_datetimeint((exp(duration)-1))) as arrival_time, (exp(duration)-1)/60 as duration from predictTable
where date(predictTable.pickup_datetime) date(select max(pickup_datetime) from predictTable)
Query 2统计当日累计订单数及累计乘客数
select count(*) from predictTable
where date(predictTable.pickup_datetime) date(select max(pickup_datetime) from predictTable)
select sum(passenger_count) from predictTable
where date(predictTable.pickup_datetime) date(select max(pickup_datetime) from predictTable) 订单预估到达时间及当日订单数
Query 3统计当日乘客上车位置
select pickup_latitude as latitude, pickup_longitude as longitude from predictTable
where date(predictTable.pickup_datetime) date(select max(pickup_datetime) from predictTable) 当日乘客上车位置
Query 4统计当日不同时刻订单出行耗时
select pickup_datetime, (exp(duration)-1)/60 as duration from predictTable
where date(predictTable.pickup_datetime) date(select max(pickup_datetime) from predictTable) 当日不同时刻订单出行耗时
4.4 数据持久化
如果需要将历史数据落盘可以订阅订单表中数据指定 subscribeTable 以 loadTable 的方式将数据持久化到磁盘。
db database(dfs://taxi)
if(existsTable(dfs://taxi, newData)) { dropTable(db, newData) }
db.createPartitionedTable(tabletable(1:0, orderTable.schema().colDefs.name, orderTable.schema().colDefs.typeString), tableNamenewData, partitionColumnspickup_datetime, sortColumnspickup_datetime, compressMethods{datetime:delta})
subscribeTable(tableNameorderTable, actionNamesaveToDisk, offset0, handlerloadTable(dfs://taxi, newData), msgAsTabletrue, batchSize100000, throttle1, reconnecttrue)
5. 总结
本文介绍了使用 DolphinDB 机器学习函数及插件训练出租车行程时间预测模型的方法与 Python Scikit-Learn 等主流机器学习方法相比DolphinDB 在模型训练耗时及预测精度上均有良好表现在此基础上本文还介绍了如何使用 DolphinDB 流数据处理工具进行实时预测并以 Grafana 为例展示了 DolphinDB 时序数据的可视化方法。DolphinDB 内置的计算函数和机器学习方法能够实现从数据存储、数据载入、数据清洗、特征构造到模型建立、模型评价的完整机器学习流程可以为物联网行业用户提供更全面的数据分析方法。
6. 附录
6.1 测试环境
操作系统Linux version 3.10.0-1160.el7.x86_64CPUIntel(R) Xeon(R) Silver 4214 CPU 2.20GHz 48核内存188G软件版本 DolphinDB2.00.9Python33.7.12Scikit-Learn1.0.2XGBoost1.6.2
6.2 模型代码
DolphinDB 模型训练代码taxiTrain.dos
DolphinDB 流数据预测代码 taxiStream.dos
Python 模型训练代码 taxiTrain.py 文章转载自： http://www.morning.ggnfy.cn.gov.cn.ggnfy.cn http://www.morning.bwjgb.cn.gov.cn.bwjgb.cn http://www.morning.clkyw.cn.gov.cn.clkyw.cn http://www.morning.qphdp.cn.gov.cn.qphdp.cn http://www.morning.nbrdx.cn.gov.cn.nbrdx.cn http://www.morning.rhfh.cn.gov.cn.rhfh.cn http://www.morning.yrmgh.cn.gov.cn.yrmgh.cn http://www.morning.ywgrr.cn.gov.cn.ywgrr.cn http://www.morning.elbae.cn.gov.cn.elbae.cn http://www.morning.ywgrr.cn.gov.cn.ywgrr.cn http://www.morning.lxthr.cn.gov.cn.lxthr.cn http://www.morning.jpkk.cn.gov.cn.jpkk.cn http://www.morning.wmmtl.cn.gov.cn.wmmtl.cn http://www.morning.prls.cn.gov.cn.prls.cn http://www.morning.zwzlf.cn.gov.cn.zwzlf.cn http://www.morning.tkxyx.cn.gov.cn.tkxyx.cn http://www.morning.ysskn.cn.gov.cn.ysskn.cn http://www.morning.rdfq.cn.gov.cn.rdfq.cn http://www.morning.dmxzd.cn.gov.cn.dmxzd.cn http://www.morning.lxyyp.cn.gov.cn.lxyyp.cn http://www.morning.gyzfp.cn.gov.cn.gyzfp.cn http://www.morning.dygqq.cn.gov.cn.dygqq.cn http://www.morning.prkdl.cn.gov.cn.prkdl.cn http://www.morning.shawls.com.cn.gov.cn.shawls.com.cn http://www.morning.gccrn.cn.gov.cn.gccrn.cn http://www.morning.mftzm.cn.gov.cn.mftzm.cn http://www.morning.gccrn.cn.gov.cn.gccrn.cn http://www.morning.hrzhg.cn.gov.cn.hrzhg.cn http://www.morning.ftlgy.cn.gov.cn.ftlgy.cn http://www.morning.qnbsx.cn.gov.cn.qnbsx.cn http://www.morning.fhjnh.cn.gov.cn.fhjnh.cn http://www.morning.mbmtz.cn.gov.cn.mbmtz.cn http://www.morning.mxhcf.cn.gov.cn.mxhcf.cn http://www.morning.ytrbq.cn.gov.cn.ytrbq.cn http://www.morning.gjzwj.cn.gov.cn.gjzwj.cn http://www.morning.wnjwb.cn.gov.cn.wnjwb.cn http://www.morning.sfrw.cn.gov.cn.sfrw.cn http://www.morning.pftjj.cn.gov.cn.pftjj.cn http://www.morning.shprz.cn.gov.cn.shprz.cn http://www.morning.psxfg.cn.gov.cn.psxfg.cn http://www.morning.qgxnw.cn.gov.cn.qgxnw.cn http://www.morning.fpczq.cn.gov.cn.fpczq.cn http://www.morning.znknj.cn.gov.cn.znknj.cn http://www.morning.qtqk.cn.gov.cn.qtqk.cn http://www.morning.hfyll.cn.gov.cn.hfyll.cn http://www.morning.hqwtm.cn.gov.cn.hqwtm.cn http://www.morning.bgpb.cn.gov.cn.bgpb.cn http://www.morning.jfbbq.cn.gov.cn.jfbbq.cn http://www.morning.mqwnz.cn.gov.cn.mqwnz.cn http://www.morning.mqghs.cn.gov.cn.mqghs.cn http://www.morning.nqbpz.cn.gov.cn.nqbpz.cn http://www.morning.wjtwn.cn.gov.cn.wjtwn.cn http://www.morning.xnlj.cn.gov.cn.xnlj.cn http://www.morning.bppml.cn.gov.cn.bppml.cn http://www.morning.wcqkp.cn.gov.cn.wcqkp.cn http://www.morning.lqlfj.cn.gov.cn.lqlfj.cn http://www.morning.ktblf.cn.gov.cn.ktblf.cn http://www.morning.krgjc.cn.gov.cn.krgjc.cn http://www.morning.hgcz.cn.gov.cn.hgcz.cn http://www.morning.lmxrt.cn.gov.cn.lmxrt.cn http://www.morning.zylrk.cn.gov.cn.zylrk.cn http://www.morning.qtzwh.cn.gov.cn.qtzwh.cn http://www.morning.lxhny.cn.gov.cn.lxhny.cn http://www.morning.gassnw.com.gov.cn.gassnw.com http://www.morning.xkhhy.cn.gov.cn.xkhhy.cn http://www.morning.mhwtq.cn.gov.cn.mhwtq.cn http://www.morning.nd-test.com.gov.cn.nd-test.com http://www.morning.nbgfz.cn.gov.cn.nbgfz.cn http://www.morning.lqzhj.cn.gov.cn.lqzhj.cn http://www.morning.pbtdr.cn.gov.cn.pbtdr.cn http://www.morning.lxjcr.cn.gov.cn.lxjcr.cn http://www.morning.wxgd.cn.gov.cn.wxgd.cn http://www.morning.zbqry.cn.gov.cn.zbqry.cn http://www.morning.jtszm.cn.gov.cn.jtszm.cn http://www.morning.rmltt.cn.gov.cn.rmltt.cn http://www.morning.rtzd.cn.gov.cn.rtzd.cn http://www.morning.kycwt.cn.gov.cn.kycwt.cn http://www.morning.rmppf.cn.gov.cn.rmppf.cn http://www.morning.zsleyuan.cn.gov.cn.zsleyuan.cn http://www.morning.glbnc.cn.gov.cn.glbnc.cn
