基于測點實測高程修正重力中區(qū)地改誤差評價方法的探討

耿濤,杜輝,馮治漢

(1.中國地質(zhì)調(diào)查局 西安地質(zhì)調(diào)查中心，陜西 西安 710054；2.西北地質(zhì)科技創(chuàng)新中心，陜西 西安 710054)

0 引言

現(xiàn)階段的重力調(diào)查工作中，隨著重力儀器精度的不斷提高和高精度差分GPS測地技術(shù)的廣泛使用，重力值觀測精度、布格改正精度、正常場改正精度均有大幅度提高，重力地形改正的精度對布格重力異?？偩鹊挠绊懺絹碓酱?，這一點大家已有共識。而在重力地形改正工作中，近區(qū)地形改正值一般通過野外實測獲得，在不斷引入新技術(shù)后，近區(qū)地形改正的效率和精度都有了極大的提高[1-3]；而遠區(qū)由于距離較遠，其地形影響對于一定范圍內(nèi)的重力異常已經(jīng)可視為常量[4]，對布格重力異常形態(tài)的影響較小，在研究局部異常特征時影響不大。因此，重力中區(qū)地形改正(下文簡稱中區(qū)地改)的精度是影響重力異常精度的關(guān)鍵。

幾十年來，眾多專家學(xué)者和專業(yè)技術(shù)工作者為提高地形改正精度持續(xù)不斷地做了大量工作。20世紀(jì)80年代以前，主要是進行改正方法方面的研究[5-8]，20世紀(jì)90年代后期，隨著計算機技術(shù)的發(fā)展，則在計算機算法研究[9-13]、不同算法的對比研究[14]、不同源數(shù)字高程模型(digital elevation model，DEM)適用性試驗研究[15-17]等方面取得了眾多成果，這些成果的應(yīng)用大大提高了中區(qū)地改的精度。然而，中區(qū)地改的誤差是客觀存在的，通過努力可以減小它，但不能消滅它。因此，在重力調(diào)查工作中，掌握中區(qū)地改誤差的實際大小，就顯得同樣重要，這有助于我們掌握重力異常的實際精度，對后續(xù)研究工作中的異常判別、認(rèn)識、反演、推斷、解釋都有幫助。誤差統(tǒng)計只是手段，真實掌握重力異常的精度才是目的，一個漂亮但不真實的誤差統(tǒng)計數(shù)據(jù)是沒有意義的。

具體到一個工作區(qū)，中區(qū)地改工作完成后，其精度到底能達到多少？這就涉及到中區(qū)地改的精度評價問題。

重力地形改正的目的就是要消除由于地形的起伏不平對測點重力觀測值的影響。因此，在不考慮地層密度變化的前提下，無論地形改正采用的地形資料(包括地形圖和DEM數(shù)據(jù))如何獲得，越逼近實際地形，地形改正的精度就越高。在實際工作中，所謂對實際地形的逼近包含兩個層面：首先是地形資料對實際地形的擬合程度(即地形資料的精度，如圖1a所示)，其次是在地形改正值計算過程中所采用的算法對地形資料的擬合程度(即計算方法的精度，如圖1b所示)。

圖1 重力地形改正中存在的兩種誤差示意

不難理解，任何算法的精度討論都是建立在所采用地形資料為真值的基礎(chǔ)之上的，如果所用地形資料本身對實際地形的描述差異就很大(即地形資料精度很低)，那么采用什么樣的算法地形改正精度也不可能高。因此，札喜旺登在1983年就指出，在影響重力地形改正精度的7個方面中，對地形改正精度影響最大的是地形圖的精度與地形模型單元大小的劃分，其他因素都是次要的[18]。

可見，在實際工作中，評價一個工作區(qū)中區(qū)地改精度時，不能回避所用地形資料精度對重力地形改正精度的影響。前人在方法研究過程中提出的一些精度評價方法，基本上是基于同一種地形資料對精度的評價方法進行研究[19-20]，或者是討論不同算法間的精度差異，或者是對比不同源DEM數(shù)據(jù)間的精度差異，而沒有考慮地形資料本身的精度問題，這在方法研究中無可非議，但把這些精度評價方法直接應(yīng)用在實際工作中是有問題的，因為這些精度評價方法沒有考慮地形資料精度對重力地形改正精度的影響。本文的目的就是探討通過具有可操作性的方法評價地形資料精度對重力地形改正精度的影響。

1 實際工作中中區(qū)地改精度評價常見問題

實際工作中，中區(qū)地改的精度評價都是依據(jù)各種重力規(guī)范中的要求進行的，首先了解一下規(guī)范中的精度評價要求。

1.1 規(guī)范對中區(qū)地改精度評價的要求

表1是目前地礦行業(yè)主要執(zhí)行的3種重力調(diào)查技術(shù)規(guī)范對中區(qū)地改精度評價的要求。為了能體現(xiàn)出技術(shù)要求的延革，表中列出了規(guī)范修訂前后的對照。

從表1可以看出，不同比例尺重力規(guī)范對中區(qū)地改的質(zhì)量檢查要求中，都是考慮到了中區(qū)地改所用地形圖或DEM的精度問題的，因此，在給出同精度質(zhì)量檢查方法的基礎(chǔ)上，建議有條件的地區(qū)采用更大比例尺的地形圖或DEM重新計算地形改正值，作為高精度檢查來評價中區(qū)地改精度。因為，一般來說，在同一地區(qū)，更大比例尺地形圖或DEM測繪精度更高，數(shù)據(jù)網(wǎng)度更密，描述的地形更逼近實際地形。

但是，規(guī)范給出的同精度檢查方式，由于采用的地形圖或DEM數(shù)據(jù)沒有變，實際上是忽略了地形圖或DEM本身的精度問題。

1.2 實際工作中常見問題

從表1可以看出，常用的重力規(guī)范中雖然都有建議采用更高精度的地形圖或DEM來檢查中區(qū)地改的精度，但因為這個建議不是強制性的，因此，在一般生產(chǎn)單位的實際工作中，鮮有按此建議進行中區(qū)地改精度檢查的。

表1 規(guī)范對中區(qū)地改精度評價的要求

由于沒有統(tǒng)一的要求，目前在實際工作中，中區(qū)地改的方法多種多樣。由于采用計算機計算的中區(qū)地改的方法已非常成熟，人工圓域量板讀圖改正法已被淘汰。在中區(qū)地改精度評價時，采用人工圓域量板讀圖法來檢查計算機計算中區(qū)地改結(jié)果，本身是一種低精度改正結(jié)果檢查高精度改正結(jié)果的方法[18]，因此也不再被采用，而主要采用變換高程數(shù)據(jù)節(jié)點位置法或采用不同算法重新計算。

但是，正如前文所述，這些方法都是基于同一種DEM數(shù)據(jù)進行的，回避了DEM本身的精度問題。對于一些地形相對平緩的工作區(qū)或采用的DEM數(shù)據(jù)節(jié)點本身較密的工作區(qū)，變換高程數(shù)據(jù)節(jié)點位置法或采用不同算法重新計算，其結(jié)果差異是非常小的，這也是經(jīng)常見到一些項目中區(qū)地改質(zhì)量檢查精度畸高的原因。

造成這一現(xiàn)象的原因，筆者認(rèn)為，并非是大家沒有意識到DEM精度對中區(qū)地改精度有影響，而是缺少對DEM精度進行評價的手段。規(guī)范中給出的建議出發(fā)點是好的，但在實際工作中操作性不強。原因如下：

首先，如果按規(guī)范的建議進行中區(qū)地改精度檢查，由于質(zhì)量檢查要求均勻分布的原則，則要求工作單位需準(zhǔn)備兩套DEM數(shù)據(jù)，一套符合規(guī)范要求的DEM數(shù)據(jù)進行中區(qū)地改工作，另一套精度更高的DEM數(shù)據(jù)專門用來進行質(zhì)量檢查。由于在同一地區(qū)，獲取更高精度的DEM數(shù)據(jù)往往需要花費更高的代價，甚至需要專門測繪、制作。出于項目經(jīng)費和工作成本考慮，除非是有特別的研究需要，沒有哪個單位愿意去準(zhǔn)備兩套數(shù)據(jù)。

其次，如規(guī)范中建議所言，有條件的地區(qū)(可理解為項目承擔(dān)單位本身就有工作區(qū)不同比例尺的兩套或以上DEM數(shù)據(jù))，那么，也不會采用規(guī)范中建議的方法進行中區(qū)地改精度檢查。因為，規(guī)范中規(guī)定的中區(qū)地改所應(yīng)采用的地形資料比例尺是精度下限，為了提高中區(qū)地改精度，在實際工作中一定是采用最高精度的DEM來完成中區(qū)地改工作，而不是僅僅用來做精度檢查。

因此，如何評價地形DEM的精度，是更真實地評價中區(qū)地改精度的重要內(nèi)容。

2 采用重力測點實測高程值評價DEM精度

2.1 方法的可行性

在當(dāng)前的重力調(diào)查工作中，隨著高精度差分GPS測地技術(shù)(包括GPS快速靜態(tài)、GPS-RTK、CORS系統(tǒng)等)的廣泛使用，中、小比例尺重力調(diào)查工作中，重力測點的三維坐標(biāo)測量精度已經(jīng)可以達到分米級，大比例尺重力勘探中重力測點的三維坐標(biāo)精度甚至可達到厘米級，并且，重力測點是均勻地分布于整個工作區(qū)的。

而用于進行中區(qū)地改的DEM數(shù)據(jù)主要是通過航空攝影測量法、航天遙感測量法、地形圖掃描矢量化法及數(shù)字線劃圖縮編法等方法生成數(shù)字線劃圖(DLG)[27]，在此基礎(chǔ)上再生成DEM[28]，其精度指標(biāo)分別按測繪行業(yè)標(biāo)準(zhǔn)CH/T9001.1-2013[29]和CH/T9001.2-2010[30]執(zhí)行。其中DEM的精度指標(biāo)見表2。

表2 數(shù)字高程模型(DEM)精度指標(biāo)

從表2可見，在實際工作中重力測點三維坐標(biāo)測量精度是優(yōu)于DEM精度的，利用重力測點實測高程對中區(qū)地改用DEM數(shù)據(jù)進行檢查，相當(dāng)于對其進行高精度檢查，是可行的。

2.2 檢查方法

在利用重力測點實測高程對中區(qū)地改用DEM數(shù)據(jù)進行檢查前，首先要保證兩種數(shù)據(jù)采用的是相同的坐標(biāo)系統(tǒng)，如果不是，則應(yīng)先進行坐標(biāo)轉(zhuǎn)換?，F(xiàn)在按規(guī)定均應(yīng)采用CGCS 2000國家坐標(biāo)系，1985國家高程基準(zhǔn)。

以實測重力測點的平面坐標(biāo)為準(zhǔn)，在DEM上讀取相應(yīng)點的高程；注意高程值的內(nèi)插方法應(yīng)與進行中區(qū)地改時選用的方法相同。

計算各個重力測點實測高程與DEM高程的差值，采用如下高精度質(zhì)量檢查均方誤差統(tǒng)計公式計算DEM的均方誤差：

(1)

式中：Δh為重力測點DEM讀取高程與實測高程之差；n為重力測點數(shù)。

目前在中、小比例尺重力調(diào)查工作中，中區(qū)地改采用的DEM比例尺一般是1∶10000和1∶50000，因此，計算DEM的均方誤差宜用高精度質(zhì)量檢查均方誤差統(tǒng)計公式，而對于大比例尺重力勘探工作，如果地形條件較好，DEM的精度較高，應(yīng)采用如下同精度質(zhì)量檢查均方誤差統(tǒng)計公式計算DEM的均方誤差：

(2)

式中各項含義同公式(1)。

如果實際工作中中區(qū)地改使用了其他來源的高程模型數(shù)據(jù)，檢查時采用何種均方誤差統(tǒng)計公式應(yīng)視所用數(shù)據(jù)具體情況確定。

3 DEM誤差引起的中區(qū)地改誤差的評價方法探討

用上述方法評價DEM的精度是符合抽樣檢查要求的，其結(jié)果可視為DEM的均方誤差。但是，我們的最終目的是評價中區(qū)地改的精度，因此，需要計算DEM的誤差引起的中區(qū)地改誤差。

然而，這是一個困難的問題。因為理論上講，只有用更高精度的DEM數(shù)據(jù)才能在更高精度標(biāo)準(zhǔn)上評價中區(qū)地改誤差，或是用同等精度但不同數(shù)據(jù)源的DEM數(shù)據(jù)做同精度的中區(qū)地改誤差評價。重力點實測高程雖然可以作為抽樣點檢查DEM高程精度，但要用來檢查中區(qū)地改的精度，實測高程點的密度(即對實際地形描述的精確度)是遠遠不夠的，因此只能通過別的途徑估算DEM的誤差造成的中區(qū)地改誤差。

如圖1a所示，DEM的誤差包括地形細節(jié)描述不準(zhǔn)確和地形高程不準(zhǔn)確兩方面，由地形改正的原理可知，其中地形高程不準(zhǔn)確造成的誤差是主要的，而地形細節(jié)描述不準(zhǔn)確造成的誤差相對要小很多。因此，可以考慮采用DEM讀取高程計算的中區(qū)地改值和實測高程計算的中區(qū)地改值之差，來近似替代DEM地形與實際地形不一致引起的中區(qū)地改誤差。

具體方法為，抽取一定比例的重力測點，采用DEM讀取高程和重力測點實測高程分別計算中區(qū)地改值，并計算其差值，對所有檢查點的差值做均方誤差統(tǒng)計。由于DEM與實際地形的誤差是非線性的，我們可以大致判斷實測重力點外圍的DEM網(wǎng)格點數(shù)據(jù)在近處與實際地形擬合情況，但無法判定其在遠處與實際地形的擬合程度，因此應(yīng)采用式(3)的同精度檢查公式進行誤差統(tǒng)計。

(3)

式中：δi為DEM高程和實測高程分別計算的中區(qū)地改值之差；n為檢查點數(shù)。

關(guān)于評價DEM誤差引起的中區(qū)地改誤差的方法，我們在討論時還提出了兩種方案：

一是計算在中區(qū)地改范圍內(nèi)、厚度為實測高程與DEM高程差值(hS-hD)的板狀模型引起的重力值，減去近區(qū)地改值后，近似視為該重力測點由DEM誤差引起的中區(qū)地改誤差，在此基礎(chǔ)上統(tǒng)計均方誤差。該方案還可進一步減化為利用全區(qū)實測高程與DEM高程差值的平均值或是均方誤差值為厚度計算板狀模型引起的重力值，將計算結(jié)果直接作為全區(qū)DEM誤差引起的中區(qū)地改均方誤差的近似值。

二是采用在DEM數(shù)據(jù)上增加噪聲的方法模擬計算由DEM誤差引起的中區(qū)地改誤差。方法是生成一組和DEM網(wǎng)格數(shù)據(jù)一致且滿足正態(tài)分布的隨機數(shù)，使其平均數(shù)μs和標(biāo)準(zhǔn)離差σs等于實測高程和DEM高程的誤差平均數(shù)μ和標(biāo)準(zhǔn)離差σ，將該隨機數(shù)與DEM數(shù)據(jù)相加(相當(dāng)于在DEM數(shù)據(jù)上增加了噪聲)，得到一個新的DEM數(shù)據(jù)；用新的DEM數(shù)據(jù)和原DEM數(shù)據(jù)分別計算中區(qū)地改值，用兩次計算結(jié)果差值的均方誤差模擬全區(qū)DEM誤差引起的中區(qū)地改均方誤差。

上述方法哪一種計算結(jié)果更接近實際誤差，需進行大量的對比實驗才能確定，本文先以第一種方法進行探討。

4 中區(qū)地改總的精度

上述方法只是評價了DEM的精度，并沒有改變DEM的數(shù)據(jù)網(wǎng)度，因此，如圖1所示的中區(qū)地改誤差中所包含的兩方面誤差——DEM地形與實際地形不一致引起的改正誤差和計算時選取的質(zhì)量模型擬合DEM不準(zhǔn)確引起的改正誤差——在中區(qū)地改精度評價時均應(yīng)予以統(tǒng)計。至于中區(qū)地改采用何種算法效果更好，前人已有很多研究[9-14]，本文不作探討。

4.1 質(zhì)量模型擬合DEM不準(zhǔn)確引起的改正誤差計算

質(zhì)量模型擬合DEM不準(zhǔn)確引起的改正誤差大小和DEM的數(shù)據(jù)網(wǎng)度有很大關(guān)系，網(wǎng)度越密，這部分誤差越小。在網(wǎng)度一定的情況下，可采用變換高程數(shù)據(jù)節(jié)點位置(包括平移、旋轉(zhuǎn))后重新計算中區(qū)地改值，利用其與原始計算值的差值進行均方誤差統(tǒng)計。計算公式應(yīng)采用式(3)的同精度檢查公式。該誤差也就是現(xiàn)階段重力工作中普遍作為中區(qū)地改總誤差來統(tǒng)計的那部分誤差，其計算統(tǒng)計方法已非常成熟，在此對其不多贅述。應(yīng)注意的是，這一過程中的前后兩次中區(qū)地改值計算，均應(yīng)采用DEM高程，而不應(yīng)采用實測高程。

4.2 中區(qū)地改總精度統(tǒng)計

DEM地形與實際地形不一致引起的改正誤差和計算時選取的質(zhì)量模型擬合DEM不準(zhǔn)確引起的改正誤差是同時且獨立存在的，因此，中區(qū)地改總的均方誤差應(yīng)采用式(4)進行計算：

(4)

式中：εgT為中區(qū)地改總的均方誤差；εgTD為DEM地形與實際地形不一致引起的改正誤差的均方差；εgTM為質(zhì)量模型擬合DEM不準(zhǔn)確引起的改正誤差的均方差。

以上兩部分誤差在計算時，如果所有測點都參與了，那么也可采用式(4)先求出每一個測點的中區(qū)地改誤差，然后用同精度檢查公式進行中區(qū)地改總的均方誤差統(tǒng)計。

5 應(yīng)用實例

以青海門源回族自治縣幅、山丹縣幅1∶25萬區(qū)域重力調(diào)查項目為例，對上述中區(qū)地改精度評價方法進行進一步說明。

5.1 工作區(qū)地形特征

工作區(qū)大部位于祁連山中部，東北部分為河西走廊和阿拉善地區(qū)，中間以龍首山相隔，地跨甘肅、青海、內(nèi)蒙古3省區(qū)。

工作區(qū)中部祁連山地勢較高峻，海拔4 000 m以上的山峰比比皆是；北段地勢不高，在龍首山脈以北地勢較低，海拔在1 500 m左右，龍首山脈為北端最高地段，海拔在3 000 m以上，最高點海拔3 616 m；在祁連山脈與龍首山之間，海拔在1 000～2 500 m之間，為著名的河西走廊。南段地勢較高，由一系列山脈及山間盆地構(gòu)成，其中與青海湖相接處海拔較低，大約3 200 m，最高峰在祁連山脈，海拔5 200 m以上，山間谷底海拔也在3 500 m左右。

工作區(qū)位于青藏高原邊部，地形整體高差大，山區(qū)地形陡峻，切割劇烈，山間盆地內(nèi)地形相對平坦。

5.2 工作區(qū)重力測點實測三維坐標(biāo)精度

重力測點三維坐標(biāo)采用GPS載波相位靜態(tài)差分測量方法進行測量。重力測點三維坐標(biāo)最終精度(包含控制網(wǎng)誤差、觀測誤差、坐標(biāo)轉(zhuǎn)換誤差)：平面均方誤差為±1.177 m，最大值為6.402 m；高程均方誤差為±0.332 m，最大值為0.927 m。可見，采用重力測點實測高程評價DEM高程精度屬高精度檢查。

5.3 工作區(qū)DEM精度評價

在開始評價前，首先將重力測點實測三維坐標(biāo)數(shù)據(jù)和DEM數(shù)據(jù)轉(zhuǎn)換為同一坐標(biāo)系統(tǒng)。在DEM中通過內(nèi)差讀取全區(qū)5 271個重力測點的DEM高程，采用式(1)計算得到DEM高程的均方誤差為7.49 m。

實測高程與DEM高程誤差分布直方圖見圖2。統(tǒng)計分析得出：實測高程與DEM高程差值(hS-hD)最大值64.846 m，最小值-59.009 m；差值絕對值在0～5 m間的測點數(shù)為3 281個，占比為 62.25%(其中差值絕對值在0～2 m間的測點數(shù)為1 675個，占比為31.78%)；絕對值在5～10 m間的測點數(shù)為1 203個，占比為22.82%；絕對值在10～15 m間的測點數(shù)為498個，占比為9.45%；絕對值在15～20 m間的測點數(shù)為171個，占比為3.24%；有118個測點的差值絕對值>20 m，占比為2.24%。

圖2 實測高程與DEM高程誤差分布直方圖

兩種高程誤差的分布與地形具有明顯的對應(yīng)關(guān)系，誤差較大的測點基本分布在山區(qū)地形起伏變化較大、切割劇烈的區(qū)域，而地形相對平坦的區(qū)域誤差均較小，其特征與表2所表述的一致。由超大誤差點(誤差值絕對值>20 m)分布與地形的關(guān)系圖(圖3)可以看出，在地形起伏、轉(zhuǎn)折變化劇烈的部位，DEM對地形的描述誤差顯然是比較大的。

圖3 超大誤差點(白色圓圈)分布與地形的關(guān)系

5.4 中區(qū)地改精度

首先，采用抽樣法評價質(zhì)量模型擬合DEM不準(zhǔn)確引起的改正誤差。全區(qū)共均勻分布抽取217個重力測點(占比4.12%)，將DEM旋轉(zhuǎn)45°后重新計算中區(qū)地改值，用式(3)統(tǒng)計其均方誤差，結(jié)果為：εgTM=±0.029×10-5m/s2。

其次，計算DEM地形與實際地形不一致引起的改正誤差。全區(qū)5 271個重力測點全部參與了計算，DEM高程與實測高程計算中區(qū)地改誤差分布直方圖見圖4。由于高程超大誤差的測點具有明顯的粗差，其計算的中區(qū)地改誤差對最終統(tǒng)計結(jié)果影響甚大，但這些點的誤差不能代表全區(qū)的真實誤差水平，因此首先予以剔除。本次統(tǒng)計剔除了高程誤差絕對值>20 m的118個測點，剩余5 153個測點用式(1)統(tǒng)計其均方誤差，最終統(tǒng)計結(jié)果為：εgTD=±0.125×10-5m/s2。

圖4 DEM高程與實測高程計算中區(qū)地改誤差分布直方圖

最后，采用式(4)計算出中區(qū)地改總的均方誤差為：εgT=±0.128×10-5m/s2。

從以上統(tǒng)計結(jié)果可見，考慮了由于DEM地形與實際地形不一致引起的中區(qū)地改誤差以后，得到的中區(qū)地改總均方誤差明顯大于只考慮質(zhì)量模型擬合DEM不準(zhǔn)確一項因素得出的中區(qū)地改總均方誤差，但這個結(jié)果更接近實際誤差，同時也說明，DEM地形與實際地形不一致引起的中區(qū)地改誤差是重力中區(qū)地形改正的主要誤差來源。

上述方法還有可改進的地方。比如，現(xiàn)在計算DEM地形與實際地形不一致引起的中區(qū)地改誤差時，采用的還都是相同的DEM數(shù)據(jù)，只是測點高程采用了不同的值，如果能采用某種算法對測點附近的DEM數(shù)據(jù)節(jié)點進行某種程度的修正，也許計算的中區(qū)地改誤差值更接近實際誤差值，這還需要進一步深入研究。

6 檢查結(jié)果的拓展應(yīng)用

中區(qū)地改中兩部分誤差的統(tǒng)計均可按規(guī)范要求進行5%～10%的抽樣檢查。但筆者建議全工作區(qū)所有測點均參與檢查，這可以為我們提供許多信息。例如，在區(qū)域重力調(diào)查工作中，中區(qū)地改所用DEM數(shù)據(jù)可能具有不同的來源，精度也不一致；可以通過繪制高程誤差分布圖、中區(qū)地改誤差分布圖及超差測點點位圖等方法，來分析工作區(qū)內(nèi)不同區(qū)域DEM的精度情況、中區(qū)地改誤差的分布情況及超差點的分布情況。如果出現(xiàn)區(qū)域性連片的超差情況，則應(yīng)考慮獲取其他精度更高的DEM[15-16]來重新完成這一區(qū)域的中區(qū)地改。即便無法解決，也可以做到心中有數(shù)，在后期的數(shù)據(jù)處理及解釋推斷過程中，對這一區(qū)域的重力異常有正確的認(rèn)識。

另外，現(xiàn)在進行中小比例尺重力調(diào)查時，一般遠區(qū)的重力地形改正也都和中區(qū)地改一樣采用DEM數(shù)據(jù)進行計算機改正，因此，該精度評價方法可以擴展至遠區(qū)。

7 結(jié)論

重力中區(qū)地改所用DEM數(shù)據(jù)是前人測量、加工形成的數(shù)據(jù)，因此，其在測量、加工過程中已經(jīng)包含了誤差。通過對DEM數(shù)據(jù)誤差的來源及大小分析，并采用實際資料對比可見，在山區(qū)等地形起伏變化較大的地方，這個誤差還是很大的，對重力中區(qū)地改精度的影響是不能忽視的。但是，如何評價這個影響，也是比較困難的。本文提出這種利用測點實測高程進行重力中區(qū)地形改正精度評價的方法，目的就是想通過具有可操作性的方法，來獲取更接近中區(qū)地改誤差真實水平的統(tǒng)計結(jié)果，為后續(xù)的研究工作打好基礎(chǔ)。但是，該方法畢竟是針對不具備DEM精度檢查條件時對中區(qū)地改精度的一種替代評價方法，有一定的局限性。這些方法都采用了等量代換的思路，哪一種方法計算的結(jié)果更接近實際誤差值，還需要在不同地形條件下進行大量的對比試驗，這是下一步研究的方向。

本文的目的主要是希望引起大家對現(xiàn)階段重力工作中中區(qū)地改精度評價不全面問題的重視，同時希望能拋磚引玉，大家共同探討，集思廣益，提出更好的解決方案，使重力中區(qū)地改的精度統(tǒng)計結(jié)果更符合實際。