一种汽车运输过程粮食质量安全评估方法及系统与流程

本发明涉及粮食质量控制领域，更具体地，涉及一种汽车运输过程粮食质量安全评估方法及系统。

背景技术：

在散粮运输途中，可能经过汽车、火车、船舶多次转运，多次装卸，运输方式的变换、远距离外部环境的变化等都会对粮食质量产生影响，很难获得运输过程中的信息，难以评估粮食质量变化风险。而汽车运输是粮食运输最为常用的运输手段，而因为汽车运输的移动不确定性和运输特殊性，目前仍无法很好的汽车运输过程中的粮食质量安全进行评估。

技术实现要素：

为了解决背景技术存在的对汽车运输过程中的粮食质量安全进行评估存在困难的问题，本发明提供了一种汽车运输过程粮食质量安全评估方法及系统；所述方法及系统通过对运载时间、粮食损耗以及行车定位三个不同维度中汽运过程中可能发生粮食质量风险进行评估，并可根据实际环境的设置不同的风险权重比例，获得汽车运输过程中的粮食质量安全评估结果；所述一种汽车运输过程粮食质量安全评估方法，包括：

采集汽车运输过程中的预设种类的数据；所述数据包括运输该批粮食的每一辆车的发运时间和到货时间、每一辆车的装车质量和卸车质量、以及每一辆车的行驶轨迹数据；

根据所述发运时间和到货时间计算获得汽运时间风险系数；

根据所述装车质量与所述卸车质量计算获得损耗风险系数；

根据所述行驶轨迹数据计算获得位置风险系数；

根据所述时间风险系数、损耗风险系数、位置风险系数以及预设的权重计算获得汽运风险系数，根据所述汽运风险系数生成粮食质量安全评估结果。

进一步的，根据所述发运时间和到货时间计算获得汽运时间风险系数，包括：

根据所述每辆车的发运时间以及到货时间，计算获得每辆车的在途时间；

计算每辆车的在途时间所形成的在途时间队列的方差；

所述时间风险系数rlogtime的计算公式为：

其中，tmax为在途时间队列里的最大在途时间；为在途时间队列的方差。

进一步的，所述根据所述装车质量与所述卸车质量计算获得损耗风险系数，包括：

根据所述每辆车的装车质量与所述卸车质量，计算获得每辆车的损耗值；

计算所有车损耗值的平均值，并统计低于所述平均值的车辆数量；

所述损耗风险系数的计算公式为：

其中，rwaste为损耗风险系数。

进一步的，所述根据所述行驶轨迹数据计算获得位置风险系数，包括：

预先设置标准行驶路线，并根据预设时间间隔以及预设车速，在地图上提取标准横坐标数列以及标准纵坐标数列；

根据所述每辆车的行驶轨迹数据，生成每辆车对应的横坐标数列以及纵坐标数列；所述采集每辆车的行驶轨迹数据，包括在地图上每辆车在行使过程中，按预设时间间隔提取的横坐标以及纵坐标；

计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数；

根据同一辆车的横坐标相关性系数与纵坐标相关性系数通过预设规则计算获得该辆车对应的风险系数；

根据所有车辆的风险系数的平均值，作为位置风险系数。

进一步的，计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，以及计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数的公式分别为：

其中，ikhorizontal为第k辆车的横坐标相关性系数，ikvertical为第k辆车的纵坐标相关性系数；ak为第k辆车的横坐标数列，bk为第k辆车的纵坐标数列；x为标准横坐标数列，y为标准纵坐标数列。

进一步的，根据同一辆车的横坐标相关性系数与纵坐标相关性系数通过预设规则计算获得该辆车对应的风险系数，包括：

根据每一辆车的横坐标相关性系数的绝对值的倒数，获得该辆车的横坐标风险系数；

根据每一辆车的纵坐标相关性系数的绝对值的倒数，获得该辆车的纵坐标风险系数；

根据每辆车的横坐标风险系数与纵坐标风险系数的平均值，获得该辆车对应的风险系数。

所述一种汽车运输过程粮食质量安全评估系统包括：

数据采集单元，所述数据采集单元用于采集汽车运输过程中的预设种类的数据；所述数据包括运输该批粮食的每一辆车的发运时间和到货时间、每一辆车的装车质量和卸车质量、以及每一辆车的行驶轨迹数据；

时间风险系数计算单元，所述时间风险系数计算单元用于根据所述发运时间和到货时间计算获得汽运时间风险系数；

损耗风险系数计算单元，所述损耗风险系数计算单元用于根据所述装车质量与所述卸车质量计算获得损耗风险系数；

位置风险系数计算单元，所述位置风险系数计算单元用于根据所述行驶轨迹数据计算获得位置风险系数；

质量安全评估单元，所述质量安全评估单元用于根据所述时间风险系数、损耗风险系数、位置风险系数以及预设的权重计算获得汽运风险系数，根据所述汽运风险系数生成粮食质量安全评估结果。

进一步的，所述时间风险系数计算单元用于根据所述每辆车的发运时间以及到货时间，计算获得每辆车的在途时间；

所述时间风险系数计算单元用于计算每辆车的在途时间所形成的在途时间队列的方差；

所述时间风险系数rlogtime的计算公式为：

其中，tmax为在途时间队列里的最大在途时间；为在途时间队列的方差。

进一步的，所述损耗风险系数计算单元用于根据所述每辆车的装车质量与所述卸车质量，计算获得每辆车的损耗值；

所述损耗风险系数计算单元用于计算所有车损耗值的平均值，并统计低于所述平均值的车辆数量；

所述损耗风险系数的计算公式为：

其中，rwaste为损耗风险系数。

进一步的，所述位置风险系数计算单元用于预先设置标准行驶路线，并根据预设时间间隔以及预设车速，在地图上提取标准横坐标数列以及标准纵坐标数列；

所述位置风险系数计算单元用于根据所述每辆车的行驶轨迹数据，生成每辆车对应的横坐标数列以及纵坐标数列；所述采集每辆车的行驶轨迹数据，包括在地图上每辆车在行使过程中，按预设时间间隔提取的横坐标以及纵坐标；

所述位置风险系数计算单元用于计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数；

所述位置风险系数计算单元用于根据同一辆车的横坐标相关性系数与纵坐标相关性系数通过预设规则计算获得该辆车对应的风险系数；

所述位置风险系数计算单元用于根据所有车辆的风险系数的平均值，作为位置风险系数。

进一步的，所述位置风险系数计算单元用于计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，以及计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数的公式分别为：

其中，ikhorizontal为第k辆车的横坐标相关性系数，ikvertical为第k辆车的纵坐标相关性系数；ak为第k辆车的横坐标数列，bk为第k辆车的纵坐标数列；x为标准横坐标数列，y为标准纵坐标数列。

进一步的，所述位置风险系数计算单元用于根据每一辆车的横坐标相关性系数的绝对值的倒数，获得该辆车的横坐标风险系数；

所述位置风险系数计算单元用于根据每一辆车的纵坐标相关性系数的绝对值的倒数，获得该辆车的纵坐标风险系数；

所述位置风险系数计算单元用于根据每辆车的横坐标风险系数与纵坐标风险系数的平均值，获得该辆车对应的风险系数。

本发明的有益效果为：本发明的技术方案，给出了一种汽车运输过程粮食质量安全评估方法及系统；所述方法及系统通过对运载时间、粮食损耗以及行车定位三个不同维度中汽运过程中可能发生粮食质量风险进行评估，并可根据实际环境的设置不同的风险权重比例，获得汽车运输过程中的粮食质量安全评估结果；所述方法及系统针对粮食在汽车运输过程中可能存在的风险隐患进行量化的评估。管理者根据不同维度的风险指数情况，可以指导管理工作，或者重新审核/变更供应商，或者对粮食重新进行卫生检查，或者重新规划物流线路、流程。

附图说明

通过参考下面的附图，可以更为完整地理解本发明的示例性实施方式：

图1为本发明具体实施方式的一种汽车运输过程粮食质量安全评估方法的流程图；

图2为本发明具体实施方式的一种汽车运输过程粮食质量安全评估系统的结构图。

具体实施方式

现在参考附图介绍本发明的示例性实施方式，然而，本发明可以用许多不同的形式来实施，并且不局限于此处描述的实施例，提供这些实施例是为了详尽地且完全地公开本发明，并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中，相同的单元/元件使用相同的附图标记。

除非另有说明，此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外，可以理解的是，以通常使用的词典限定的术语，应当被理解为与其相关领域的语境具有一致的含义，而不应该被理解为理想化的或过于正式的意义。

图1为本发明具体实施方式的一种汽车运输过程粮食质量安全评估方法的流程图；如图1所示，所述方法包括：

步骤110，采集汽车运输过程中的预设种类的数据；所述数据包括运输该批粮食的每一辆车的发运时间和到货时间、每一辆车的装车质量和卸车质量、以及每一辆车的行驶轨迹数据；

在运输过程中遗撒损耗属于正常情况，因此有个别公司或司机个人会在运输过程中私扣少量货物牟利，同时也无法判断是否给粮食造成污染，影响到粮食质量安全。直接体现在某一辆车的发运时间和到货时间，以及装车质量和卸车质量，以及运输过程中的车辆地理位置信息。

通过对运载时间、粮食损耗以及行车定位三个不同维度中汽运过程中可能发生粮食质量风险进行评估，并根据需要评估的内容采集汽车运输过程中的如上数据。

步骤120，根据所述发运时间和到货时间计算获得汽运时间风险系数；

根据所述发运时间和到货时间计算获得汽运时间风险系数，包括：

根据所述每辆车的发运时间以及到货时间，计算获得每辆车的在途时间；

计算每辆车的在途时间所形成的在途时间队列的方差；

所述时间风险系数rlogtime的计算公式为：

其中，tmax为在途时间队列里的最大在途时间；为在途时间队列的方差。

步骤130，根据所述装车质量与所述卸车质量计算获得损耗风险系数；

所述根据所述装车质量与所述卸车质量计算获得损耗风险系数，包括：

根据所述每辆车的装车质量与所述卸车质量，计算获得每辆车的损耗值；

计算所有车损耗值的平均值，并统计低于所述平均值的车辆数量；

所述损耗风险系数的计算公式为：

其中，rwaste为损耗风险系数。

步骤140，根据所述行驶轨迹数据计算获得位置风险系数；

所述根据所述行驶轨迹数据计算获得位置风险系数，包括：

预先设置标准行驶路线，并根据预设时间间隔以及预设车速，在地图上提取标准横坐标数列以及标准纵坐标数列；

根据所述每辆车的行驶轨迹数据，生成每辆车对应的横坐标数列以及纵坐标数列；所述采集每辆车的行驶轨迹数据，包括在地图上每辆车在行使过程中，按预设时间间隔提取的横坐标以及纵坐标；

计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数；

进一步的，计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，以及计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数的公式分别为：

其中，ikhorizontal为第k辆车的横坐标相关性系数，ikvertical为第k辆车的纵坐标相关性系数；ak为第k辆车的横坐标数列，bk为第k辆车的纵坐标数列；x为标准横坐标数列，y为标准纵坐标数列。

根据同一辆车的横坐标相关性系数与纵坐标相关性系数通过预设规则计算获得该辆车对应的风险系数；

具体的，所述预设规则为：根据每一辆车的横坐标相关性系数的绝对值的倒数，获得该辆车的横坐标风险系数；根据每一辆车的纵坐标相关性系数的绝对值的倒数，获得该辆车的纵坐标风险系数；根据每辆车的横坐标风险系数与纵坐标风险系数的平均值，获得该辆车对应的风险系数。

根据所有车辆的风险系数的平均值，作为位置风险系数。

步骤150，根据所述时间风险系数、损耗风险系数、位置风险系数以及预设的权重计算获得汽运风险系数，根据所述汽运风险系数生成粮食质量安全评估结果。

本实施例所述的各项风险系数，均为0到1区间内的值，当所述风险系数越接近1时，则风险越小；反之则风险越大。

若所述汽运风险系数低于某一预设阈值时，说明此次汽车运输过程中粮食质量安全存储风险较大；横向对比其他批次汽车运输粮食中时间风险系数、损耗风险系数以及位置风险系数的波动，对于本次异常运输明显低于平均的风险系数，可能是存在风险的部分，根据此生成安全风险追溯结果。

图2为本发明具体实施方式的一种汽车运输过程粮食质量安全评估系统的结构图。如图2所示，所述系统包括：

数据采集单元210，所述数据采集单元210用于采集汽车运输过程中的预设种类的数据；所述数据包括运输该批粮食的每一辆车的发运时间和到货时间、每一辆车的装车质量和卸车质量、以及每一辆车的行驶轨迹数据；

时间风险系数计算单元220，所述时间风险系数计算单元220用于根据所述发运时间和到货时间计算获得汽运时间风险系数；

进一步的，所述时间风险系数计算单元220用于根据所述每辆车的发运时间以及到货时间，计算获得每辆车的在途时间；

所述时间风险系数计算单元220用于计算每辆车的在途时间所形成的在途时间队列的方差；

所述时间风险系数tlogtime的计算公式为：

其中，tmax为在途时间队列里的最大在途时间；为在途时间队列的方差。

损耗风险系数计算单元230，所述损耗风险系数计算单元230用于根据所述装车质量与所述卸车质量计算获得损耗风险系数；

进一步的，所述损耗风险系数计算单元230用于根据所述每辆车的装车质量与所述卸车质量，计算获得每辆车的损耗值；

所述损耗风险系数计算单元230用于计算所有车损耗值的平均值，并统计低于所述平均值的车辆数量；

所述损耗风险系数的计算公式为：

其中，rwaste为损耗风险系数。

位置风险系数计算单元240，所述位置风险系数计算单元240用于根据所述行驶轨迹数据计算获得位置风险系数；

进一步的，所述位置风险系数计算单元240用于预先设置标准行驶路线，并根据预设时间间隔以及预设车速，在地图上提取标准横坐标数列以及标准纵坐标数列；

所述位置风险系数计算单元240用于根据所述每辆车的行驶轨迹数据，生成每辆车对应的横坐标数列以及纵坐标数列；所述采集每辆车的行驶轨迹数据，包括在地图上每辆车在行使过程中，按预设时间间隔提取的横坐标以及纵坐标；

所述位置风险系数计算单元240用于计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数；

所述位置风险系数计算单元240用于根据同一辆车的横坐标相关性系数与纵坐标相关性系数通过预设规则计算获得该辆车对应的风险系数；

所述位置风险系数计算单元240用于根据所有车辆的风险系数的平均值，作为位置风险系数。

进一步的，所述位置风险系数计算单元240用于计算每一辆车横坐标数列与标准横坐标数列的横坐标相关性系数，以及计算每一辆车纵坐标数列与标准纵坐标数列的纵坐标相关性系数的公式分别为：

其中，ikhorizontal为第k辆车的横坐标相关性系数，ikvertical为第k辆车的纵坐标相关性系数；ak为第k辆车的横坐标数列，bk为第k辆车的纵坐标数列；x为标准横坐标数列，y为标准纵坐标数列。

进一步的，所述位置风险系数计算单元240用于根据每一辆车的横坐标相关性系数的绝对值的倒数，获得该辆车的横坐标风险系数；

所述位置风险系数计算单元240用于根据每一辆车的纵坐标相关性系数的绝对值的倒数，获得该辆车的纵坐标风险系数；

所述位置风险系数计算单元240用于根据每辆车的横坐标风险系数与纵坐标风险系数的平均值，获得该辆车对应的风险系数。

质量安全评估单元250，所述质量安全评估单元250用于根据所述时间风险系数、损耗风险系数、位置风险系数以及预设的权重计算获得汽运风险系数，根据所述汽运风险系数生成粮食质量安全评估结果。

在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本公开的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。

本领域那些技术人员可以理解，可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件，以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外，可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述，本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。本说明书中涉及到的步骤编号仅用于区别各步骤，而并不用于限制各步骤之间的时间或逻辑的关系，除非文中有明确的限定，否则各个步骤之间的关系包括各种可能的情况。

此外，本领域的技术人员能够理解，尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本公开的范围之内并且形成不同的实施例。例如，在权利要求书中所要求保护的实施例的任意之一都可以以任意的组合方式来使用。

本公开的各个部件实施例可以以硬件实现，或者以在一个或者多个处理器上运行的软件模块实现，或者以它们的组合实现。本公开还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如，计算机程序和计算机程序产品)。这样的实现本公开的程序可以存储在计算机可读介质上，或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到，或者在载体信号上提供，或者以任何其他形式提供。

应该注意的是上述实施例对本公开进行说明而不是对本公开进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本公开可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。

以上所述仅是本公开的具体实施方式，应当指出的是，对于本领域的普通技术人员来说，在不脱离本公开精神的前提下，可以作出若干改进、修改、和变形，这些改进、修改、和变形都应视为落在本申请的保护范围内。