一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法
【專利摘要】本發(fā)明屬于物化探數(shù)據(jù)處理領(lǐng)域,具體涉及一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法。目的是在利用物化探數(shù)據(jù)編制等值線圖時,消除由于數(shù)據(jù)的測量年份不同、測量方法不同或測量批次不同所導(dǎo)致的背景差異。該方法包括:將整個圖區(qū)內(nèi)不同年份或者不同方法測量的物化探數(shù)據(jù)分成若干個子圖區(qū),計算各子圖區(qū)平均值和標(biāo)準(zhǔn)方差;取其中面積最大圖幅子圖區(qū)的平均值作為基準(zhǔn)背景值,將其他子圖區(qū)數(shù)據(jù)根據(jù)所述基準(zhǔn)背景值進(jìn)行調(diào)整;計算調(diào)整后的各子區(qū)全部數(shù)據(jù)點(diǎn)作為整體求取平均值,重新編制平面等值圖。該方法可以將不同年份、不同儀器或不同批次所測量的數(shù)據(jù)調(diào)節(jié)到一個區(qū)域背景下編制等值圖,使全區(qū)具有統(tǒng)一的背景值,突出了低值片區(qū)的異常。
【專利說明】一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于物化探數(shù)據(jù)處理領(lǐng)域,具體涉及一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,適用于利用物化探數(shù)據(jù)在編制等值線圖時,由于數(shù)據(jù)的測量年份不同、測量方法不同或測量批次不同所導(dǎo)致的背景差異的消除。
【背景技術(shù)】
[0002]物化探數(shù)據(jù)取得后需要編制等值圖對數(shù)據(jù)進(jìn)行分析利用。但對于較大區(qū)域的數(shù)據(jù)取得,由于區(qū)域范圍大,相鄰區(qū)域數(shù)據(jù)往往會分階段或分年份取得,也可能由靈敏度不同的測量儀器或者不同的實(shí)驗(yàn)室來分析、測量取得數(shù)據(jù),那么如果不加調(diào)整直接用來編制等值圖,勢必會造成圖面由于背景值的差異所造成的“臺階”狀出現(xiàn),同時,也壓抑了某些平均值較低區(qū)域的異常,造成失真,因此,有必要對數(shù)據(jù)進(jìn)行調(diào)整。
【發(fā)明內(nèi)容】
[0003]本發(fā)明的目的是提供一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,在利用物化探數(shù)據(jù)編制等值線圖時,消除由于數(shù)據(jù)的測量年份不同、測量方法不同或測量批次不同所導(dǎo)致的背景差異。
[0004]本發(fā)明的技術(shù)方案如下所述:
[0005]一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,包括如下步驟:
[0006]步驟S1:將整個圖區(qū)內(nèi)不同年份或者不同方法測量的物化探數(shù)據(jù)分成若干個子圖區(qū),按子圖區(qū)分別計算子圖區(qū)平均值和標(biāo)準(zhǔn)方差;
[0007]步驟S2:取其中面積最大圖幅子圖區(qū)的平均值或者已存在的相鄰區(qū)區(qū)域背景作為基準(zhǔn)背景值Xb,將其他子圖區(qū)數(shù)據(jù)根據(jù)所述基準(zhǔn)背景值進(jìn)行調(diào)整;
[0008]調(diào)整的方法是:假設(shè)某一子圖區(qū)平均值是Xi,如果X1-Xb < 0,則該子圖區(qū)各點(diǎn)數(shù)值加上X1-Xb的絕對值;如果X1-Xb > 0,則該子圖區(qū)各點(diǎn)數(shù)值減去X1-Xb的絕對值;
[0009]步驟S3:計算調(diào)整后的各子區(qū)全部數(shù)據(jù)點(diǎn)作為整體求取平均值,重新編制平面等值圖。
[0010]如上所述的一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,其中:所述步驟I中,在計算各子圖區(qū)平均值和標(biāo)準(zhǔn)方差時,需要反復(fù)循環(huán)計算剔除大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn),直到無大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)存在為止。
[0011]如上所述的一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,其中:所述步驟S3中計算平均值時,反復(fù)循環(huán)計算剔除大于正n倍標(biāo)準(zhǔn)方差和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn),直到無大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)存在為止。
[0012]本發(fā)明的有益效果為:
[0013](I)本發(fā)明用于物化探數(shù)據(jù)等值圖制作時人為因素所導(dǎo)致的子區(qū)背景差異的調(diào)節(jié),可以將不同年份、不同儀器或不同批次所測量的數(shù)據(jù)調(diào)節(jié)到一個區(qū)域背景下編制等值圖,從而使全區(qū)具有統(tǒng)一的背景值,突出了低值片區(qū)的異常,也保證了整幅圖件的尺度統(tǒng)
O
[0014](2)方法操作簡便,可操作性強(qiáng)。
[0015](3)對于具有類似問題的物化探數(shù)據(jù)處理具有普遍性,適用性強(qiáng)。
【專利附圖】
【附圖說明】
[0016]圖1為本發(fā)明提供的一種用于等值圖制作時人為因素所導(dǎo)致的子區(qū)背景差異的調(diào)節(jié)方法的流程圖;
[0017]圖2為未經(jīng)調(diào)節(jié)的數(shù)據(jù)繪制的等值圖;
[0018]圖3為經(jīng)本發(fā)明提供的方法調(diào)節(jié)后的數(shù)據(jù)繪制的等值圖;
[0019]圖中①、②、③、④、⑤片是不同時間、不同測量儀器所測物探數(shù)據(jù),可以看出①、②之間,①、③之間和①、③、④之間有明顯的區(qū)域界線。從圖2可以看出,各片數(shù)據(jù)之間的系統(tǒng)差異已經(jīng)消除,整個圖面具有一致的數(shù)據(jù)尺度。
【具體實(shí)施方式】
[0020]下面結(jié)合附圖和實(shí)施例對本發(fā)明提供的一種用于等值圖制作時人為因素所導(dǎo)致的子區(qū)背景差異的調(diào)節(jié)方法進(jìn)行詳細(xì)說明。
[0021]如圖1所示,一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,包括如下步驟:
[0022]步驟S1:將整個圖區(qū)內(nèi)不同年份或者不同方法測量的物化探數(shù)據(jù)分成若干個子圖區(qū),按子圖區(qū)分別計算子圖區(qū)平均值和標(biāo)準(zhǔn)方差。為使計算準(zhǔn)確,便于后續(xù)使用,在計算各子圖區(qū)平均值和標(biāo)準(zhǔn)方差時,需要反復(fù)循環(huán)計算剔除大于正n(n的大小視測區(qū)范圍和比例尺而定)倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn),直到無大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)存在為止。
[0023]步驟S2:取其中面積最大圖幅子圖區(qū)的平均值或者已存在的相鄰區(qū)(指整個圖區(qū)的相鄰區(qū))區(qū)域背景作為基準(zhǔn)背景值Xb,將其他子圖區(qū)數(shù)據(jù)根據(jù)所述基準(zhǔn)背景值進(jìn)行調(diào)整;
[0024]調(diào)整的方法是:假設(shè)某一子圖區(qū)平均值是Xi,如果X1-Xb < 0,則該子圖區(qū)各點(diǎn)數(shù)值加上X1-Xb的絕對值;如果X1-Xb > 0,則該子圖區(qū)各點(diǎn)數(shù)值減去X1-Xb的絕對值;
[0025]步驟S3:計算調(diào)整后的各子區(qū)全部數(shù)據(jù)點(diǎn)作為整體求取平均值,重新編制平面等值圖。為優(yōu)化計算,計算平均值時,反復(fù)循環(huán)計算剔除大于正n倍標(biāo)準(zhǔn)方差和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn),直到無大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)存在為止。
[0026]試驗(yàn)時,以新疆伊犁盆地南緣氡測量數(shù)據(jù)處理為例,進(jìn)行如下操作:
[0027]步驟I)將整個圖區(qū)內(nèi)不同年份和不同方法測量的物化探數(shù)據(jù)分成五個子圖區(qū),按子圖區(qū)分別計算平均值和標(biāo)準(zhǔn)方差。做此計算時,需要反復(fù)循環(huán)計算剔除大于正3.5倍標(biāo)準(zhǔn)方差和小于負(fù)3.5倍標(biāo)準(zhǔn)方差的異值,直到無大于正3.5倍標(biāo)準(zhǔn)方差和小于負(fù)3.5倍標(biāo)準(zhǔn)方差的異值存在為止,此時各子圖區(qū)的平均值作為該子圖區(qū)的背景值;
[0028]步驟2)取其中面積最大子圖區(qū)平均值作為基準(zhǔn)背景值Xb,如圖2所示,本實(shí)施例中取第④號子圖反復(fù)剔除正負(fù)3.5倍標(biāo)準(zhǔn)方差后的平均值作為基準(zhǔn)背景值Xb,其他①、②、③、⑤子圖區(qū)數(shù)據(jù)據(jù)此進(jìn)行調(diào)整,調(diào)整的具體方法是:①號子圖區(qū)背景值是XI,因?yàn)閄l-Xb< O,因此把①號子圖區(qū)各點(diǎn)數(shù)值加上Xl-Xb的絕對值,(當(dāng)然如果Xl-Xb > O,則該子圖區(qū)各點(diǎn)數(shù)值減去X1-Xb的絕對值),然后按此方法對②、③、⑤號子圖分別做如上處理,從而將由于分析年份不同或分析批次不同所造成的背景值差異消除,將整個編圖區(qū)歸結(jié)到一個區(qū)域背景上;
[0029]步驟3)將背景值歸一后的①、②、③、④、⑤號子圖區(qū)數(shù)據(jù)作為一個整體數(shù)據(jù)區(qū)重新反復(fù)剔除正負(fù)3.5倍標(biāo)準(zhǔn)方差后,取該平均值作為整個圖區(qū)的區(qū)域背景值編制平面等值圖。
[0030]處理后的等值圖如圖3所示,從圖3可以看出,各片數(shù)據(jù)之間的系統(tǒng)差異已經(jīng)消除,整個圖面具有一致的數(shù)據(jù)尺度。
【權(quán)利要求】
1.一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,包括如下步驟: 步驟S1:將整個圖區(qū)內(nèi)不同年份或者不同方法測量的物化探數(shù)據(jù)分成若干個子圖區(qū),按子圖區(qū)分別計算子圖區(qū)平均值和標(biāo)準(zhǔn)方差; 步驟S2:取其中面積最大圖幅子圖區(qū)的平均值或者已存在的相鄰區(qū)區(qū)域背景作為基準(zhǔn)背景值Xb,將其他子圖區(qū)數(shù)據(jù)根據(jù)所述基準(zhǔn)背景值進(jìn)行調(diào)整; 調(diào)整的方法是:假設(shè)某一子圖區(qū)平均值是Xi,如果X1-Xb < O,則該子圖區(qū)各點(diǎn)數(shù)值加上X1-Xb的絕對值;如果X1-Xb > O,則該子圖區(qū)各點(diǎn)數(shù)值減去X1-Xb的絕對值; 步驟S3:計算調(diào)整后的各子區(qū)全部數(shù)據(jù)點(diǎn)作為整體求取平均值,重新編制平面等值圖。
2.根據(jù)權(quán)利要求1所述的一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,其特征在于:所述步驟I中,在計算各子圖區(qū)平均值和標(biāo)準(zhǔn)方差時,需要反復(fù)循環(huán)計算剔除大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn),直到無大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)存在為止。
3.根據(jù)權(quán)利要求2所述的一種物化探數(shù)據(jù)的子區(qū)背景差異的調(diào)節(jié)方法,其特征在于:所述步驟S3中計算平均值時,反復(fù)循環(huán)計算剔除大于正n倍標(biāo)準(zhǔn)方差和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn),直到無大于正n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)和小于負(fù)n倍標(biāo)準(zhǔn)方差的數(shù)據(jù)點(diǎn)存在為止。
【文檔編號】G06T7/00GK103778638SQ201410042869
【公開日】2014年5月7日 申請日期:2014年1月29日 優(yōu)先權(quán)日:2014年1月29日
【發(fā)明者】李林強(qiáng), 趙寧博, 寧斯同 申請人:核工業(yè)北京地質(zhì)研究院