一種基于樣點(diǎn)時(shí)空代表性的npp監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)方法
【專利摘要】本發(fā)明涉及一種基于樣點(diǎn)時(shí)空代表性的NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)方法，該方法根據(jù)監(jiān)測(cè)區(qū)域的大小和監(jiān)測(cè)目標(biāo)，利用相應(yīng)的時(shí)間序列NPP遙感產(chǎn)品，在對(duì)其時(shí)空重建的基礎(chǔ)上，逐像元計(jì)算平均相對(duì)偏差及其標(biāo)準(zhǔn)差，分析各像元的時(shí)空代表性，根據(jù)可接受的平均相對(duì)偏差和可投入建立的監(jiān)測(cè)點(diǎn)規(guī)模最終確定監(jiān)測(cè)網(wǎng)的樣點(diǎn)分布。本發(fā)明簡(jiǎn)單實(shí)用、考慮空間關(guān)系、不需要任何假設(shè)且能充分利用現(xiàn)有多種遙感數(shù)據(jù)的NPP監(jiān)測(cè)網(wǎng)采樣，實(shí)現(xiàn)了充分利用了容易獲取的NPP遙感產(chǎn)品，可快速、低成本的設(shè)計(jì)NPP長(zhǎng)期監(jiān)測(cè)網(wǎng)，并可用于已有觀測(cè)網(wǎng)的冗余檢查與優(yōu)化。
【專利說明】—種基于樣點(diǎn)時(shí)空代表性的NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)方法
[0001]
【技術(shù)領(lǐng)域】
[0002]本發(fā)明涉及地理學(xué)與生態(tài)科學(xué)【技術(shù)領(lǐng)域】，具體說是一種通過分析樣點(diǎn)的時(shí)空代表性來確定NPP監(jiān)測(cè)點(diǎn)位置和數(shù)量的方法，該方法也可用于優(yōu)化已有監(jiān)測(cè)網(wǎng)，以減少觀測(cè)冗余。
【背景技術(shù)】
[0003]植被凈第一性生產(chǎn)力(NPP)是指綠色植物在單位面積、單位時(shí)間內(nèi)所累積的有機(jī)物數(shù)量，是由光合作用所產(chǎn)生的有機(jī)質(zhì)總量(GPP)中扣除自養(yǎng)呼吸后的剩余部分，它直接反映植物群落在自然環(huán)境條件下的生產(chǎn)能力。NPP在全球碳循環(huán)研究中扮演著重要的角色，是評(píng)價(jià)生態(tài)系統(tǒng)結(jié)構(gòu)與功能協(xié)調(diào)性、人類承載力的重要指標(biāo)，是進(jìn)行生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動(dòng)研究的重要基礎(chǔ)。
[0004]在全球變化和人類活動(dòng)加劇的背景下，認(rèn)識(shí)和了解某一區(qū)域NPP的實(shí)際變化，需要建立長(zhǎng)期NPP監(jiān)測(cè)網(wǎng)，如何科學(xué)合理的設(shè)計(jì)監(jiān)測(cè)網(wǎng)，使其觀測(cè)值能長(zhǎng)期的反映該區(qū)域的NPP變化，即確定監(jiān)測(cè)網(wǎng)中觀測(cè)點(diǎn)的位置和數(shù)量，是這一問題的關(guān)鍵。
[0005]目前，常用的采樣設(shè)計(jì)方法包括經(jīng)典采樣方法和基于空間統(tǒng)計(jì)的采樣方法兩大類。經(jīng)典采樣方法如隨機(jī)采樣、系統(tǒng)采樣和分層隨機(jī)采樣等，執(zhí)行簡(jiǎn)單、需要的數(shù)據(jù)少，但其不考慮采樣要素的空間關(guān)系，用其設(shè)計(jì)的觀測(cè)網(wǎng)一般需要較多的觀測(cè)點(diǎn)。基于空間統(tǒng)計(jì)是目前應(yīng)用和研究較多的一類方法，該方法考慮采樣要素的空間自相關(guān)性，設(shè)計(jì)的監(jiān)測(cè)網(wǎng)可得到最優(yōu)無偏的NPP估計(jì)結(jié)果，但該方法一般需要較多的先驗(yàn)數(shù)據(jù)，且需要滿足二階平穩(wěn)假設(shè)，而這個(gè)假設(shè)在實(shí)際中很難滿足。
[0006]隨著地球觀測(cè)系統(tǒng)的建立，NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)中可用的數(shù)據(jù)已經(jīng)越來越多，多種分辨率的時(shí)間序列NPP遙感產(chǎn)品可以容易低廉的獲取，雖然這些遙感數(shù)據(jù)的精度還有待進(jìn)一步驗(yàn)證，但其空間覆蓋能力是空前的，且能很好反映NPP的空間結(jié)構(gòu)，這些數(shù)據(jù)可為地面監(jiān)測(cè)網(wǎng)的設(shè)計(jì)提供豐富的先驗(yàn)知識(shí)。因此，設(shè)計(jì)一種簡(jiǎn)單實(shí)用、不需要任何假設(shè)、且能充分利用豐富的先驗(yàn)知識(shí)的NPP監(jiān)測(cè)網(wǎng)設(shè)計(jì)方法。(Vachaud G, Silans APD, Balabanis P，et al.Temporal stability of spatially measured soil water probability densityfunction.Soil Science Society of America, 1985,49: 822-28)等提出的空穩(wěn)定性的概念為發(fā)展這樣的方法提供了可能。
[0007]時(shí)空穩(wěn)定性的概念是指空間采樣點(diǎn)與環(huán)境屬性統(tǒng)計(jì)參數(shù)之間關(guān)系的時(shí)間不變性。Kachanoski (Kachanoski, R.G., & De Jong, E.Scale dependence and the temporalpersistence of spatial patterns of soil waters storage.Water ResourcesResearch, 1988，24，85 - 91)等又對(duì)這一概念進(jìn)行了發(fā)展，他們認(rèn)為時(shí)空穩(wěn)定性是與空間尺度有關(guān)的空間分布結(jié)構(gòu)或模式在時(shí)間上的持續(xù)性?？梢?，時(shí)空穩(wěn)定性描述的是環(huán)境屬性空間結(jié)構(gòu)在時(shí)間上的持續(xù)性，可用于表達(dá)某一采樣點(diǎn)相對(duì)于整個(gè)區(qū)域平均值的時(shí)間穩(wěn)定性，用于評(píng)價(jià)某一觀測(cè)點(diǎn)估計(jì)區(qū)域平均的能力。時(shí)空穩(wěn)定性概念在土壤水分觀測(cè)數(shù)據(jù)的分析中得到了較多的應(yīng)用，研究者認(rèn)為:在點(diǎn)觀測(cè)向區(qū)域平均的尺度轉(zhuǎn)換，時(shí)空穩(wěn)定性有利于減少尺度轉(zhuǎn)換的偏差。如果換一個(gè)應(yīng)用對(duì)象和角度，利用時(shí)空穩(wěn)定性分析方法分類先驗(yàn)數(shù)據(jù)的每一像元的時(shí)空代表性，則可將其變?yōu)橐环NNPP監(jiān)測(cè)網(wǎng)采樣方法。

【發(fā)明內(nèi)容】

[0008]基于對(duì)問題的上述思考，本發(fā)明的目的在于提供一種基于樣點(diǎn)時(shí)空代表性的NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)方法。該方法根據(jù)監(jiān)測(cè)區(qū)域的大小和監(jiān)測(cè)目標(biāo)，利用相應(yīng)的時(shí)間序列NPP遙感產(chǎn)品，在質(zhì)量控制的基礎(chǔ)上，逐像元分析其時(shí)空穩(wěn)定性，根據(jù)時(shí)空穩(wěn)定性分析結(jié)果和監(jiān)測(cè)網(wǎng)規(guī)模確定監(jiān)測(cè)網(wǎng)樣點(diǎn)的數(shù)量和分布。
[0009]本發(fā)明的目的是通過以下技術(shù)方案來實(shí)現(xiàn):
一種基于樣點(diǎn)時(shí)空代表性的NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)方法，其步驟是:
第一步:先驗(yàn)數(shù)據(jù)準(zhǔn)備
根據(jù)監(jiān)測(cè)區(qū)域的大小，綜合考慮NPP觀測(cè)儀器本身的足跡(Footprint)，選用不同空間分辨率的時(shí)間序列NPP遙感產(chǎn)品；
第二步:時(shí)空重建
對(duì)獲取的時(shí)間序列NPP遙感產(chǎn)品進(jìn)行時(shí)空重建，得到時(shí)空連續(xù)的NPP遙感產(chǎn)品； 第三步:時(shí)空穩(wěn)定性計(jì)算
計(jì)算時(shí)空穩(wěn)定性的兩個(gè)指標(biāo)平均相對(duì)偏差及其標(biāo)準(zhǔn)差；
第四步:時(shí)空穩(wěn)定性分析
根據(jù)平均相對(duì)偏差對(duì)所有像元進(jìn)行排序并標(biāo)出相應(yīng)的標(biāo)準(zhǔn)差，分析并初步選擇可作為采樣點(diǎn)的像元；
第五步:確定采樣點(diǎn)數(shù)量和位置
根據(jù)擬投入建立的觀測(cè)點(diǎn)數(shù)量和精度要求，從平均相對(duì)偏差排序中選擇時(shí)空代表性最好的像元作為采樣點(diǎn)。
[0010]本發(fā)明的優(yōu)點(diǎn):
1、本發(fā)明簡(jiǎn)單實(shí)用、可充分利用已有的先驗(yàn)知識(shí)、考慮了時(shí)間和空間關(guān)系，且不需要任何假設(shè)；
2、本發(fā)明可快速、低成本的設(shè)計(jì)NPP監(jiān)測(cè)網(wǎng)，確定了監(jiān)測(cè)點(diǎn)的數(shù)量和位置，并可用于已有觀測(cè)網(wǎng)的冗余檢查與優(yōu)化。
【專利附圖】

【附圖說明】
[0011]圖1基于樣點(diǎn)時(shí)空代表性的NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)流程圖。
[0012]圖2樣點(diǎn)時(shí)空代表性分析。
【具體實(shí)施方式】
[0013]綜合考慮監(jiān)測(cè)區(qū)域的大小、地面NPP觀測(cè)方式等信息，選用不同時(shí)空分辨率的NPP遙感產(chǎn)品，并利用時(shí)空濾波的方法對(duì)其進(jìn)行時(shí)間序列重建，利用重建后時(shí)空完整的NPP遙感產(chǎn)品，逐像元做時(shí)空穩(wěn)定性分析，即計(jì)算每個(gè)像元的平均相對(duì)偏差及其標(biāo)準(zhǔn)差，最后根據(jù)擬投入的觀測(cè)成本和精度要求，確定觀測(cè)點(diǎn)的數(shù)量和位置，即選擇平均相對(duì)偏差接近0且其標(biāo)準(zhǔn)差小的像元建立觀測(cè)點(diǎn)。
[0014]下面，就具體的方法表述如下:
一種基于樣點(diǎn)時(shí)空代表性的NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)方法，其步驟是:
第一步:獲取時(shí)間序列NPP遙感廣品； 根據(jù)監(jiān)測(cè)區(qū)域的大小，綜合考慮擬采用的NPP地面觀測(cè)儀器本身的空間代表性，即足跡(Footprint),選用不同空間分辨率的時(shí)間序列NPP遙感產(chǎn)品。如果監(jiān)測(cè)區(qū)域較大,且采用渦動(dòng)相關(guān)儀(EC)觀測(cè)，則利用MODIS的250米NPP產(chǎn)品即可。NPP遙感產(chǎn)品的時(shí)間分辨率可根據(jù)遙感產(chǎn)品的可獲取性和對(duì)時(shí)間代表性的要求，一般在1-20天為宜，最大不宜超過30天?？蛇x的時(shí)間序列NPP遙感產(chǎn)品可來源于以下幾個(gè)方面:
UMODIS遙感產(chǎn)品，空間分辨率從250米到I公里，時(shí)間分辨率逐日；
2、基于TM/ETM+估算NPP，空間分辨率為30米，時(shí)間分辨率為16天；
3、環(huán)境(HJ)衛(wèi)星估算NPP，空間分辨率為30米，時(shí)間分辨率可達(dá)到2天。
[0015]第二步:時(shí)間序列重建；
獲取所需要的時(shí)間序列NPP遙感產(chǎn)品，需要評(píng)估產(chǎn)品的時(shí)空完整性。對(duì)于時(shí)空不完整的產(chǎn)品需要做時(shí)間序列重建處理。本發(fā)明中，缺失像元的NPP值采用Garcia, (2010)發(fā)展的基于離散余弦變換(Discrete Cosine Transforms, DCT)的懲罰最小二乘回歸算法計(jì)算得到，該算法可明確的使用時(shí)間序列信息估計(jì)缺失值，該算法的matlab函數(shù)(smoothn)可在 http://www.biomecardi0.com/matlab/smoothn.html 下載，下面是該網(wǎng)站提供了一個(gè)計(jì)算示例:
X = Iinspace(0, 100, 2'8);
y = cos (x/10) + (x/50).'2 + randn(size(x))/5;
y([70 75 80]) = [5.5 5 6];
[z, s] = smoothn(y) ; %如果時(shí)間分辨率較粗,則建議使用普通平滑。
[0016]zr = smoothn (y, ’ robust’)； ％如果時(shí)間分辨率高,如逐日的時(shí)間分辨率,貝U建議使用強(qiáng)烈的平滑，即增加Robust參數(shù)。
[0017]subplot (121), plot (x, y, ' r.' , x, z)
title (’Regular smoothing’)
subplot (122), plot (x, y, ’ r.’，x, zr)
title (’Robust smoothing’)
第三步:時(shí)空穩(wěn)定性計(jì)算；
時(shí)空穩(wěn)定性用平均相對(duì)偏差及其標(biāo)準(zhǔn)差來評(píng)價(jià)，是Vachaud等(1985)提出的，在這里，我們用其評(píng)價(jià)每個(gè)遙感像元NPP的時(shí)空代表性。
[0018]平均相對(duì)偏差的定義如下:
si=~2^—f—'(3)
1J-1
(3)式中，j表示時(shí)間序列NPP產(chǎn)品的具體時(shí)間，i表示像元編號(hào)，則表示整個(gè)監(jiān)測(cè)區(qū)內(nèi)第n個(gè)像元j時(shí)相的NPP值。.?表示j時(shí)相監(jiān)測(cè)區(qū)內(nèi)所有像元的平均NPP值。平均相對(duì)偏差表示某一像元的NPP相對(duì)于整個(gè)監(jiān)測(cè)區(qū)域NPP平均值的狀態(tài)，正值表示該像元NPP相對(duì)較高，負(fù)值表示該像元NPP相對(duì)較低。
[0019]任意一監(jiān)測(cè)點(diǎn)平均相對(duì)偏差的標(biāo)準(zhǔn)差定義為:
【權(quán)利要求】
1.一種基于樣點(diǎn)時(shí)空代表性的NPP監(jiān)測(cè)網(wǎng)采樣設(shè)計(jì)方法，其步驟是: 第一步:先驗(yàn)數(shù)據(jù)準(zhǔn)備 根據(jù)監(jiān)測(cè)區(qū)域的大小，綜合考慮NPP觀測(cè)儀器本身的足跡，選用不同空間分辨率的時(shí)間序列NPP遙感廣品； 第二步:時(shí)空重建 對(duì)獲取的時(shí)間序列NPP遙感產(chǎn)品進(jìn)行時(shí)空重建，得到時(shí)空連續(xù)的NPP遙感產(chǎn)品； 第三步:時(shí)空穩(wěn)定性計(jì)算 計(jì)算時(shí)空穩(wěn)定性的兩個(gè)指標(biāo)平均相對(duì)偏差及其標(biāo)準(zhǔn)差； 第四步:時(shí)空穩(wěn)定性分析 根據(jù)平均相對(duì)偏差對(duì)所有像元進(jìn)行排序并標(biāo)出相應(yīng)的標(biāo)準(zhǔn)差，分析并初步選擇可作為采樣點(diǎn)的像元； 第五步:確定采樣點(diǎn)數(shù)量和位置 根據(jù)擬投入建立的觀測(cè)點(diǎn)數(shù)量和精度要求，從平均相對(duì)偏差排序中選擇時(shí)空代表性最好的像元作為采樣點(diǎn)。
【文檔編號(hào)】G06F17/50GK103631999SQ201310612246
【公開日】2014年3月12日 申請(qǐng)日期:2013年11月28日 優(yōu)先權(quán)日:2013年11月28日
【發(fā)明者】冉有華 申請(qǐng)人:中國(guó)科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所