本申請(qǐng)涉及數(shù)據(jù)分析領(lǐng)域，具體涉及一種用于污染物溯源的設(shè)備和方法。

背景技術(shù)：

目前，大氣、水源等污染狀況日益嚴(yán)重。污染的形成受排放物、擴(kuò)散條件、地理環(huán)境等因素影響，其成因復(fù)雜多樣，這為準(zhǔn)確地分析某地區(qū)污染的分布和擴(kuò)散趨勢(shì)帶來(lái)了困難。

現(xiàn)有的方法主要基于高斯模型建立大氣污染擴(kuò)散模型。然而，該方法僅適用于對(duì)較大區(qū)域(市區(qū)與城鎮(zhèn))進(jìn)行較長(zhǎng)時(shí)間且尺度較大(以月份為單位)的污染擴(kuò)散分析，無(wú)法用于短時(shí)間、細(xì)時(shí)空粒度的污染物溯源。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明提出一種對(duì)不同監(jiān)測(cè)站點(diǎn)所監(jiān)測(cè)到的污染物進(jìn)行溯源的技術(shù)方案。主要思想是：通過(guò)迭代計(jì)算不同監(jiān)測(cè)點(diǎn)污染濃度在相鄰時(shí)間窗的相關(guān)性，找出污染源頭。另外，本發(fā)明的技術(shù)方案同樣也適用于針對(duì)與大氣污染具有相似時(shí)序特征數(shù)據(jù)(如水源污染、交通擁塞流等)的事故源定位追蹤。

根據(jù)本發(fā)明的一個(gè)方面，提供了一種用于污染物溯源的設(shè)備，包括：獲取單元，被配置為獲取監(jiān)測(cè)點(diǎn)的監(jiān)測(cè)數(shù)據(jù)；計(jì)算單元，被配置為從監(jiān)測(cè)數(shù)據(jù)中選取參數(shù)，并對(duì)每個(gè)參數(shù)均勻劃分級(jí)別，以及設(shè)定時(shí)間窗和時(shí)間間隔，計(jì)算監(jiān)測(cè)點(diǎn)之間的激活力矩陣；以及確定單元，被配置為根據(jù)激活力矩陣中的數(shù)值大小，在距離約束下依次找出前一時(shí)間窗的激活力最強(qiáng)的監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。

在一個(gè)實(shí)施例中，計(jì)算單元被配置為：對(duì)選取的參數(shù)進(jìn)行分組；以及通過(guò)對(duì)分組后的參數(shù)進(jìn)行直方圖均衡來(lái)劃分級(jí)別。

在一個(gè)實(shí)施例中，計(jì)算單元被配置為：計(jì)算監(jiān)測(cè)點(diǎn)i到達(dá)共現(xiàn)級(jí)別的次數(shù)fi；計(jì)算監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的次數(shù)fj；計(jì)算監(jiān)測(cè)點(diǎn)i在特定時(shí)間窗內(nèi)先于或同時(shí)和監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的次數(shù)fij；計(jì)算監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的距離dij；以及根據(jù)fi、fj、fij和dij來(lái)計(jì)算監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的激活力。

在一個(gè)實(shí)施例中，計(jì)算單元被配置為根據(jù)以下參數(shù)來(lái)計(jì)算監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的距離djj：監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的地理距離、實(shí)時(shí)風(fēng)力級(jí)別、風(fēng)向與監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j的連線之間的夾角、監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j的海拔高度、以及風(fēng)速監(jiān)測(cè)點(diǎn)的海拔高度。

在一個(gè)實(shí)施例中，確定單元被配置為：設(shè)定污染源數(shù)n以及距離約束，其中n為大于1的正整數(shù)；在目標(biāo)時(shí)間窗的激活力矩陣中找出滿足距離約束的值最大的前n對(duì)監(jiān)測(cè)點(diǎn)；以及以找到的前n對(duì)監(jiān)測(cè)點(diǎn)為目標(biāo)，依次在前一時(shí)間窗的激活力矩陣中找出滿足距離約束的值最大的n對(duì)監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。

根據(jù)本發(fā)明的另一個(gè)方面，提供了一種用于污染物溯源的方法，包括：獲取監(jiān)測(cè)點(diǎn)的監(jiān)測(cè)數(shù)據(jù)；從監(jiān)測(cè)數(shù)據(jù)中選取參數(shù)，并對(duì)每個(gè)參數(shù)均勻劃分級(jí)別，以及設(shè)定時(shí)間窗和時(shí)間間隔，計(jì)算監(jiān)測(cè)點(diǎn)之間的激活力矩陣；以及根據(jù)激活力矩陣中的數(shù)值大小，在距離約束下依次找出前一時(shí)間窗的激活力最強(qiáng)的監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。

在一個(gè)實(shí)施例中，對(duì)選取的參數(shù)進(jìn)行分組；以及通過(guò)對(duì)分組后的參數(shù)進(jìn)行直方圖均衡來(lái)劃分級(jí)別。

在一個(gè)實(shí)施例中，計(jì)算監(jiān)測(cè)點(diǎn)i到達(dá)共現(xiàn)級(jí)別的次數(shù)fi；計(jì)算監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的次數(shù)fj；計(jì)算監(jiān)測(cè)點(diǎn)i在特定時(shí)間窗內(nèi)先于或同時(shí)和監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的次數(shù)fij；計(jì)算監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的距離dij；以及根據(jù)fi、fj、fij和dij來(lái)計(jì)算監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的激活力。

在一個(gè)實(shí)施例中，根據(jù)以下參數(shù)來(lái)計(jì)算監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的距離djj：監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的地理距離、實(shí)時(shí)風(fēng)力級(jí)別、風(fēng)向與監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j的連線之間的夾角、監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j的海拔高度、以及風(fēng)速監(jiān)測(cè)點(diǎn)的海拔高度。

在一個(gè)實(shí)施例中，設(shè)定污染源數(shù)n以及距離約束，其中n為大于1的正整數(shù)；在目標(biāo)時(shí)間窗的激活力矩陣中找出滿足距離約束的值最大的前n對(duì)監(jiān)測(cè)點(diǎn)；以及以找到的前n對(duì)監(jiān)測(cè)點(diǎn)為目標(biāo)，依次在前一時(shí)間窗的激活力矩陣中找出滿足距離約束的值最大的n對(duì)監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。

采用本發(fā)明的技術(shù)方案，能夠?qū)崿F(xiàn)較短時(shí)間(例如數(shù)小時(shí))、細(xì)時(shí)空粒度(例如1km*1km，每15分鐘)上的污染物溯源。

附圖說(shuō)明

通過(guò)下文結(jié)合附圖的詳細(xì)描述，本發(fā)明的上述和其它特征將會(huì)變得更加明顯，其中：

圖1是示出了根據(jù)本發(fā)明的用于污染物溯源的設(shè)備的框圖。

圖2是示出了根據(jù)本發(fā)明的用于污染物溯源的方法的流程圖。

圖3-13是示出了根據(jù)本發(fā)明一個(gè)具體示例的數(shù)據(jù)計(jì)算結(jié)果的示意圖。

具體實(shí)施方式

下面，通過(guò)結(jié)合附圖對(duì)本發(fā)明的具體實(shí)施例的描述，本發(fā)明的原理和實(shí)現(xiàn)將會(huì)變得明顯。應(yīng)當(dāng)注意的是，本發(fā)明不應(yīng)局限于下文所述的具體實(shí)施例。另外，為了簡(jiǎn)便起見，省略了與本發(fā)明無(wú)關(guān)的公知技術(shù)的詳細(xì)描述。

圖1是示出了根據(jù)本發(fā)明一個(gè)實(shí)施例的用于污染物溯源的設(shè)備的框圖。如圖1所示，設(shè)備10包括獲取單元110、計(jì)算單元120和確定單元130。下面，詳細(xì)描述設(shè)備10中的各個(gè)單元的操作。

獲取單元110被配置為獲取監(jiān)測(cè)點(diǎn)的監(jiān)測(cè)數(shù)據(jù)，該監(jiān)測(cè)數(shù)據(jù)例如可以是時(shí)空數(shù)據(jù)。在本申請(qǐng)中，“時(shí)空數(shù)據(jù)”是指同時(shí)具有時(shí)間和空間維度的數(shù)據(jù)，例如大氣污染監(jiān)測(cè)數(shù)據(jù)、交通流量數(shù)據(jù)等等。

在一個(gè)示例中，對(duì)于大氣污染監(jiān)測(cè)站點(diǎn)來(lái)說(shuō)，監(jiān)測(cè)點(diǎn)的時(shí)空數(shù)據(jù)可以包括6種主要大氣污染物(pm2.5、pm10、so2、no2、co、o3)的濃度及其對(duì)應(yīng)空氣質(zhì)量指數(shù)(iaqi)值。

計(jì)算單元120被配置為從監(jiān)測(cè)數(shù)據(jù)中選取參數(shù)，并對(duì)每個(gè)參數(shù)均勻劃分級(jí)別，使得每個(gè)級(jí)別的監(jiān)測(cè)數(shù)據(jù)的數(shù)量大致相同。計(jì)算單元120還設(shè)定時(shí)間窗和時(shí)間間隔，并計(jì)算監(jiān)測(cè)點(diǎn)之間的激活力矩陣。

確定單元130被配置為根據(jù)激活力矩陣中的數(shù)值大小，在距離約束下依次找出前一時(shí)間窗的激活力最強(qiáng)的監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。從而，最終確定污染源頭。

下面，通過(guò)具體實(shí)施例來(lái)描述圖1所示的設(shè)備10的操作。

在本實(shí)施例中，用激活力矩陣來(lái)表示監(jiān)測(cè)點(diǎn)之間的相關(guān)性。在本申請(qǐng)中，激活力是一種可以表述網(wǎng)絡(luò)節(jié)點(diǎn)關(guān)系的統(tǒng)計(jì)量，能夠用來(lái)構(gòu)建復(fù)雜網(wǎng)絡(luò)。

獲取單元110獲取多個(gè)監(jiān)測(cè)點(diǎn)的監(jiān)測(cè)數(shù)據(jù)。就大氣污染監(jiān)測(cè)站點(diǎn)來(lái)說(shuō)，獲取單元110可以獲取6種主要大氣污染物(pm2.5，pm10，so2，no2，co，o3)的濃度及其對(duì)應(yīng)空氣質(zhì)量指數(shù)(iaqi)值。

計(jì)算單元120從監(jiān)測(cè)數(shù)據(jù)中選取參數(shù)，并對(duì)每個(gè)參數(shù)均勻劃分級(jí)別，使得每個(gè)級(jí)別的監(jiān)測(cè)數(shù)據(jù)的數(shù)量大致相同。優(yōu)選地，計(jì)算單元120對(duì)選取的參數(shù)進(jìn)行分組，以及通過(guò)對(duì)分組后的參數(shù)進(jìn)行直方圖均衡來(lái)劃分級(jí)別。例如，以分組后的大氣污染物(iaqi1，iaqi2，iaqi3)為例，對(duì)于不同季節(jié)/時(shí)間段，直方圖均衡后的所劃分的等級(jí)數(shù)目和等級(jí)之間的界限不盡相同，可以視實(shí)際需求而定。以2013.12-2014.2期間大氣污染監(jiān)測(cè)數(shù)據(jù)為例，可以將其分為三組，每組10個(gè)污染等級(jí)，如圖3所示。

計(jì)算單元120設(shè)定時(shí)間窗和時(shí)間間隔，計(jì)算監(jiān)測(cè)點(diǎn)之間的激活力矩陣。例如，計(jì)算單元120設(shè)定時(shí)間窗itv，其由以t時(shí)刻為終點(diǎn)時(shí)刻的m個(gè)時(shí)間間隔(間隔時(shí)長(zhǎng)為l2)構(gòu)成，某個(gè)指定時(shí)間段(時(shí)長(zhǎng)為l)的時(shí)間窗(時(shí)長(zhǎng)為l1)的個(gè)數(shù)為：

n＝(l-l1)/l2+1

相應(yīng)地，第i和j監(jiān)測(cè)點(diǎn)在時(shí)間窗itv里的激活力可表示為：

激活力afij描述污染監(jiān)測(cè)點(diǎn)間的激活關(guān)聯(lián)關(guān)系，即其各自變化的關(guān)聯(lián) 關(guān)系。它是一個(gè)有向的度量，即afij不等于afji。fij為監(jiān)測(cè)點(diǎn)i在時(shí)間窗itv內(nèi)先于或同時(shí)和監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的的次數(shù)，fi為監(jiān)測(cè)點(diǎn)i到達(dá)共現(xiàn)級(jí)別的次數(shù)，fj為監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的次數(shù)，距離dij定義為受地理距離、風(fēng)向和地勢(shì)影響的距離：

其中，d0為監(jiān)測(cè)點(diǎn)i與j之間的地理距離，s為實(shí)時(shí)風(fēng)力級(jí)別，θ為風(fēng)向與i與j之間連線的夾角(其取值在[0，π])，h為監(jiān)測(cè)點(diǎn)i和j所處的平均海拔高度，h0為風(fēng)速監(jiān)測(cè)點(diǎn)的海拔高度。由此可見，若i與j之間為順風(fēng)，則dij相對(duì)減小。若風(fēng)向與i與j之間連線方向重合，則風(fēng)速影響大。此外，根據(jù)大氣運(yùn)動(dòng)理論，高度越大，風(fēng)速越大，則dij也相對(duì)減小。

例如，計(jì)算單元120計(jì)算監(jiān)測(cè)點(diǎn)i在時(shí)間窗itv(m＝5，共6個(gè)時(shí)刻：t-5，t-4，t-3，t-2，t-1，t)內(nèi)的級(jí)別依次為{5，1，2，2，3，4}，而監(jiān)測(cè)點(diǎn)j在時(shí)間窗itv內(nèi)的級(jí)別依次為{3，2，2，2，2，4}，則兩監(jiān)測(cè)點(diǎn)同時(shí)出現(xiàn)在同一級(jí)別的時(shí)刻為t，t-2，和t-3，其共現(xiàn)級(jí)別為{2，4}，則fij＝1，fji＝2，fi＝3，fj＝4。

確定單元130根據(jù)激活力矩陣中的數(shù)值大小，在距離約束下依次找出前一時(shí)間窗的激活力最強(qiáng)的監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。

例如，確定單元130設(shè)定污染源數(shù)n以及距離約束d，其中n為大于1的正整數(shù)，距離約束d是指監(jiān)測(cè)點(diǎn)i和監(jiān)測(cè)點(diǎn)j之間的距離滿足dij＜d。

然后，確定單元130在目標(biāo)時(shí)間窗的激活力矩陣中找出滿足距離約束的值最大的前n對(duì)監(jiān)測(cè)點(diǎn)。例如，設(shè)m個(gè)監(jiān)測(cè)點(diǎn)在目標(biāo)時(shí)間窗(即以目標(biāo)時(shí)刻t為終點(diǎn)時(shí)刻的時(shí)間窗)的激活力矩陣為：

確定單元130找出該激活力矩陣中滿足距離約束d且值最大的前n對(duì)監(jiān)測(cè)點(diǎn)，每對(duì)監(jiān)測(cè)點(diǎn)中，對(duì)應(yīng)較早時(shí)刻的監(jiān)測(cè)點(diǎn)為另一監(jiān)測(cè)點(diǎn)污染物的源 頭。例如，監(jiān)測(cè)點(diǎn)a和b監(jiān)測(cè)點(diǎn)間的相關(guān)值較大，在該對(duì)監(jiān)測(cè)點(diǎn)中，監(jiān)測(cè)點(diǎn)b對(duì)應(yīng)時(shí)刻t-1早于監(jiān)測(cè)點(diǎn)a的時(shí)刻t，則監(jiān)測(cè)點(diǎn)b為監(jiān)測(cè)點(diǎn)a處污染的源頭。

接下來(lái)，確定單元130以找到的前n對(duì)監(jiān)測(cè)點(diǎn)為目標(biāo)，依次在前一時(shí)間窗的激活力矩陣中找出滿足距離約束的值最大的n對(duì)監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。

下面，以一個(gè)具體示例來(lái)詳細(xì)說(shuō)明本實(shí)施例中的設(shè)備10的操作。

假設(shè)獲取單元110獲取1號(hào)至6號(hào)監(jiān)測(cè)站點(diǎn)2014年7月1日9：00-7月2日11：00(共26小時(shí))的pm2.5濃度數(shù)據(jù)。

計(jì)算單元120設(shè)定時(shí)間窗長(zhǎng)度為25小時(shí)，每時(shí)間窗包含99個(gè)時(shí)間間隔，每個(gè)間隔長(zhǎng)度為15分鐘，共包含100個(gè)時(shí)刻。時(shí)間窗的個(gè)數(shù)n計(jì)算過(guò)程如下：

n＝(l-l1)/l2+1

＝(26*60-100*15)/15+1

＝5

即，2014年7月1日9：00-7月2日11：00(共26小時(shí))共有5個(gè)時(shí)間窗，如圖4所示。

計(jì)算單元120計(jì)算監(jiān)測(cè)點(diǎn)間的距離dij，如圖5所示。

計(jì)算單元120在時(shí)間窗itv1(即以目標(biāo)t時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中，統(tǒng)計(jì)6個(gè)監(jiān)測(cè)點(diǎn)間的某監(jiān)測(cè)點(diǎn)i在該時(shí)間窗內(nèi)先于或同時(shí)和另一監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的次數(shù)fij，結(jié)果如圖6所示。

計(jì)算單元120在時(shí)間窗itv1(即以目標(biāo)t時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中，統(tǒng)計(jì)6個(gè)監(jiān)測(cè)點(diǎn)間的某監(jiān)測(cè)點(diǎn)i到達(dá)共現(xiàn)級(jí)別的次數(shù)fi，結(jié)果如圖7所示。然后，計(jì)算單元120在時(shí)間窗itv1(即以目標(biāo)t時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中，統(tǒng)計(jì)6個(gè)監(jiān)測(cè)點(diǎn)間的某監(jiān)測(cè)點(diǎn)j到達(dá)共現(xiàn)級(jí)別的次數(shù)fj，結(jié)果如圖8所示。

基于圖11和圖12的計(jì)算結(jié)果，計(jì)算單元120計(jì)算在時(shí)間窗itv1(即以目標(biāo)t時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中6個(gè)監(jiān)測(cè)點(diǎn)之間的激活力矩陣af^t(×10^-4)，結(jié)果如圖9所示。

以類似的方式，計(jì)算單元120計(jì)算在時(shí)間窗itv2(即以t-1時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中6個(gè)監(jiān)測(cè)點(diǎn)間的激活力矩陣結(jié)果如圖10所示。

以類似的方式，計(jì)算單元120計(jì)算在時(shí)間窗itv3(即以t-2時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中6個(gè)監(jiān)測(cè)點(diǎn)間的激活力矩陣結(jié)果如圖11所示。

以類似的方式，計(jì)算單元120計(jì)算在時(shí)間窗itv4(即以t-3時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中6個(gè)監(jiān)測(cè)點(diǎn)間的激活力矩陣結(jié)果如圖12所示。

以類似的方式，計(jì)算單元120計(jì)算在時(shí)間窗itv5(即以t-4時(shí)刻為終點(diǎn)時(shí)刻的時(shí)間窗)中6個(gè)監(jiān)測(cè)點(diǎn)間的激活力矩陣結(jié)果如圖13所示。

確定單元130設(shè)置n＝3，距離約束d＝2。

以監(jiān)測(cè)點(diǎn)6為例，確定單元130追溯其污染源。在目標(biāo)時(shí)間窗(即以目標(biāo)時(shí)刻t為終點(diǎn)時(shí)刻的時(shí)間窗)，激活力矩陣af^t的第6列為：(75.78854408，136.1259795，543.2346494，35.55249721，105.0967919，0)

確定單元130由此可以確定，與監(jiān)測(cè)點(diǎn)6相關(guān)性最強(qiáng)的3個(gè)點(diǎn)為監(jiān)測(cè)點(diǎn)3，2，5。同時(shí)，考慮距離約束d6，3，d6，2，d6，5＜d，所以監(jiān)測(cè)點(diǎn)6的最可能源頭為監(jiān)測(cè)點(diǎn)3。

接著，確定單元130以監(jiān)測(cè)點(diǎn)3為目標(biāo)，找到af^t-1的第3列：(259.0593409，546.6317198，0，37.39217989，14.82641246，4.256889997)

參考距離約束d后，確定單元130可以確定觀測(cè)點(diǎn)3的最可能源頭為觀測(cè)點(diǎn)2。

接著，確定單元130以監(jiān)測(cè)點(diǎn)2為目標(biāo)，找到af^t-2的第2列：(484.4417732，0，73.06454847，25.77477702，11.33722744，1.732755692)

參考距離約束d后，確定單元130可以確定觀測(cè)點(diǎn)2的最可能源頭為觀測(cè)點(diǎn)1。

接著，確定單元130以監(jiān)測(cè)點(diǎn)1為目標(biāo)，找到af^t-3的第1列：(0，72.16548944，24.59670026，13.84624713，1.360286145，4.393017483)

參考距離約束d后，確定單元130可以確定觀測(cè)點(diǎn)1的最可能源頭為觀測(cè)點(diǎn)2。

接著，確定單元130以監(jiān)測(cè)點(diǎn)2為目標(biāo)，找到af^t-4的第2列：(712.1877616，0，86.41876318，19.74760869，14.32023423，1.402560188)

參考距離約束d后，確定單元130可以確定觀測(cè)點(diǎn)2的最可能源頭為觀測(cè)點(diǎn)1。

至此，所有時(shí)間窗遍歷結(jié)束。確定單元130可以確定監(jiān)測(cè)點(diǎn)6的最可能的污染源頭的位置位于監(jiān)測(cè)點(diǎn)1。

采用本發(fā)明的技術(shù)方案，能夠?qū)崿F(xiàn)較短時(shí)間(例如數(shù)小時(shí))、細(xì)時(shí)空粒度(例如1km*1km，每15分鐘)上的污染物溯源。

圖2是示出了根據(jù)本發(fā)明一個(gè)實(shí)施例的用于污染物溯源的方法的流程圖。如圖2所示，方法20在步驟s210處開始。

在步驟s220，獲取監(jiān)測(cè)點(diǎn)的監(jiān)測(cè)數(shù)據(jù)。例如，該監(jiān)測(cè)數(shù)據(jù)可以包括大氣污染監(jiān)測(cè)數(shù)據(jù)或交通數(shù)據(jù)。

在步驟s230，從監(jiān)測(cè)數(shù)據(jù)中選取參數(shù)，并對(duì)每個(gè)參數(shù)均勻劃分級(jí)別，以及設(shè)定時(shí)間窗和時(shí)間間隔，計(jì)算監(jiān)測(cè)點(diǎn)之間的激活力矩陣。優(yōu)選地，對(duì)選取的參數(shù)進(jìn)行分組，以及通過(guò)對(duì)分組后的參數(shù)進(jìn)行直方圖均衡來(lái)劃分級(jí)別。

在步驟s240，根據(jù)激活力矩陣中的數(shù)值大小，在距離約束下依次找出前一時(shí)間窗的激活力最強(qiáng)的監(jiān)測(cè)點(diǎn)，直到所有時(shí)間窗被追溯。

最后，方法20在步驟s250處結(jié)束。

應(yīng)該理解，本發(fā)明的上述實(shí)施例可以通過(guò)軟件、硬件或者軟件和硬件兩者的結(jié)合來(lái)實(shí)現(xiàn)。例如，上述實(shí)施例中的系統(tǒng)內(nèi)的各種組件可以通過(guò)多種器件來(lái)實(shí)現(xiàn)，這些器件包括但不限于：模擬電路、數(shù)字電路、通用處理器、數(shù)字信號(hào)處理(dsp)電路、可編程處理器、專用集成電路(asic)、現(xiàn)場(chǎng)可編程門陣列(fpga)、可編程邏輯器件(cpld)，等等。

另外，本領(lǐng)域的技術(shù)人員可以理解，本發(fā)明實(shí)施例中描述的初始參數(shù)可以存儲(chǔ)在本地?cái)?shù)據(jù)庫(kù)中，也可以存儲(chǔ)在分布式數(shù)據(jù)庫(kù)中或者可以存儲(chǔ)在遠(yuǎn)程數(shù)據(jù)庫(kù)中。

此外，這里所公開的本發(fā)明的實(shí)施例可以在計(jì)算機(jī)程序產(chǎn)品上實(shí)現(xiàn)。更具體地，該計(jì)算機(jī)程序產(chǎn)品是如下的一種產(chǎn)品：具有計(jì)算機(jī)可讀介質(zhì)，計(jì)算機(jī)可讀介質(zhì)上編碼有計(jì)算機(jī)程序邏輯，當(dāng)在計(jì)算設(shè)備上執(zhí)行時(shí)，該計(jì)算機(jī)程序邏輯提供相關(guān)的操作以實(shí)現(xiàn)本發(fā)明的上述技術(shù)方案。當(dāng)在計(jì)算系統(tǒng)的至少一個(gè)處理器上執(zhí)行時(shí)，計(jì)算機(jī)程序邏輯使得處理器執(zhí)行本發(fā)明實(shí)施例所述的操作(方法)。本發(fā)明的這種設(shè)置典型地提供為設(shè)置或編碼在例如光介質(zhì)(例如cd-rom)、軟盤或硬盤等的計(jì)算機(jī)可讀介質(zhì)上的軟件、代碼和/或其他數(shù)據(jù)結(jié)構(gòu)、或者諸如一個(gè)或多個(gè)rom或ram或prom芯片上的固件或微代碼的其他介質(zhì)、或一個(gè)或多個(gè)模塊中的可下載的軟件圖像、共享數(shù)據(jù)庫(kù)等。軟件或固件或這種配置可安裝在計(jì)算設(shè)備上，以使得計(jì)算設(shè)備中的一個(gè)或多個(gè)處理器執(zhí)行本發(fā)明實(shí)施例所描述的技術(shù)方案。

盡管以上已經(jīng)結(jié)合本發(fā)明的優(yōu)選實(shí)施例示出了本發(fā)明，但是本領(lǐng)域的技術(shù)人員將會(huì)理解，在不脫離本發(fā)明的精神和范圍的情況下，可以對(duì)本發(fā)明進(jìn)行各種修改、替換和改變。因此，本發(fā)明不應(yīng)由上述實(shí)施例來(lái)限定，而應(yīng)由所附權(quán)利要求及其等價(jià)物來(lái)限定。