基于ArcGIS_10.0的地表高程序貫高斯模擬和地質(zhì)統(tǒng)計建模1.數(shù)據(jù)需求：總共需要四個文件夾：即兩個輸入型文件夾和兩個輸出型文件夾，其中輸入型文件夾包括一個存放原始數(shù)據(jù)的文件夾和一個存放shp格式數(shù)據(jù)的文件夾；而輸出型文件夾包含一個存放簡單克里金估計結(jié)果（dbf格式）的文件夾和一個存放序貫高斯模擬結(jié)果（dbf格式）的文件夾。以及mjgdb數(shù)據(jù)庫：即一個存放簡單克里金輸出結(jié)果（Table格式）的數(shù)據(jù)庫，一個存放100次序貫高斯模擬的直接云圖結(jié)果的數(shù)據(jù)庫以及一個存放序貫高斯模擬輸出結(jié)果（Table格式）的數(shù)據(jù)庫（之所以需要給出簡單克里金插值結(jié)果，是因為序貫高斯模擬只能在簡單克里金插值結(jié)果的基礎(chǔ)上進行計算和分析）2.數(shù)據(jù)關(guān)系1）Import_Surface_Data&Background_Grid_xlsx&txt文件夾：①Import_Surface_Data.xlsx（原始數(shù)據(jù)，包含位置坐標(biāo)和屬性值，供建立shp用）②Import_Surface_Data.txt（原始數(shù)據(jù)，包含位置坐標(biāo)和屬性值，供Surfer插值生成背景網(wǎng)格用以及后期和ArcGIS插值結(jié)果作對比）③Background_Grid.xlsx（只包含待插值點的位置坐標(biāo)，一般來說插值高程的話為x，y；插值屬性的話可以為x，y，z；）Import_Surface_Data.grdOutput_grid_of_kriging_standard_deviation_from_surfer.grdOut_import_Surface_Data.datOut_Transform.datOut_Convert.dat）Import_Map_shp文件夾：Import_Surface_Data.dbfBackground_Grid.dbfKriging_Estimation_Result_Second.dbfKriging_Estimation_Result_Third.dbf3）Export_Kriging_Result_shp文件夾：Export_Output_Second.dbfExport_Output_Third.dbf）Extract_Simulation_Results_shp文件夾：?100次的模擬結(jié)果（dbf格式），在此不詳細列出Sequential_Gaussian_Simulation.gdb數(shù)據(jù)庫:1?100次的模擬結(jié)果，在此不詳細列出Export_Kriging_Result.gdb數(shù)據(jù)庫：

Export_Output_SecondExport_Output_ThirdExtract_Simulation_Results.gdb數(shù)據(jù)庫：1?100次的模擬結(jié)果（Table格式），在此不詳細列出.實施方案本文以廣州龍頭山表面高程數(shù)據(jù)為例，本次序貫高斯模擬總共含有3大步驟:A.首先需要建立基于龍頭山表面高程原始數(shù)據(jù)的shp文件和利用Surfer軟件生成的待插值背景網(wǎng)格節(jié)點的shp文件，其中需要通過Surfer軟件先生成10m*10m的插值結(jié)果grd文件，在把其轉(zhuǎn)換成dat格式文件，方便得到背景網(wǎng)格的shp文件；B.然后通過簡單克里金對原始數(shù)據(jù)進行插值，在得到克里金插值結(jié)果的同時結(jié)合背景網(wǎng)格文件生成一個基于驗證和預(yù)測的克里金插值結(jié)果shp文件，從中可以看到每個網(wǎng)格節(jié)點上的預(yù)測結(jié)果值與預(yù)測標(biāo)準差；C.最后利用上述克里金插值結(jié)果來獲得100次序貫高斯隨機模擬的結(jié)果，其中包含每個網(wǎng)格節(jié)點上的最大值、最小值、均值和標(biāo)準差等等，以及每一次模擬所生成的云圖，并通過提取功能得到模擬的輸出結(jié)果以便于后期外部處理和利用。.具體技術(shù)路線ReEdk?Ga EuisiseWMf島同瞄他ReEdk?Ga EuisiseWMf島同瞄他9IBK 軍口 .L4MCmr3 EwaL>Jna&fircl￡S^arxxLa i,m_31s*Ls'*i, .“apigs?IX)田Iripw二還電力口?二日25_^^>=*_口^必加西HiEtEFEwd_￡ind口!；*LtTD5he?2iTT|49lfipcn_5irTVCv_ELiLi|lrsport_Sir<wc?_Ori>5?禺品肅一匚3出啟動ArcMap導(dǎo)入之前生成的shp文件，如下圖所示:『■如BLsntwi,.二.,L.：drOdDlll-rumNlajE0fl￡1t?tAiLayimiDokpfeJ.La^$15rs:?15（業(yè)ldajn￡*d OpduK..UHgDomtnLalHlj匕船mu?an_CwrrrtFsrtwoECi'Bpl'kanh-UA?廿?rap?rivi..S3-1■■+:;ad-=[[■■由十c:[[esJImcK13二二一一iHHEiiaa；Hi一-U444U;￡■■*■“HTtf-t,-!1混(:￡q;i..f■■?.！.i.h.Bi^Hl儲拄已措卷三強乏■拄笠FmiT《覆"中騁事U4UC-三二5:一USI二3laTrJnrg5'5T=:一U4UL一ntrnnffH修改標(biāo)記，方便查看和論文成圖，如下圖所示:>.Bji￡rf'Lr:Hla,-+liHri1租用 4ttt*注曲：N升I4UI如果第一次安裝ArcGIS,地統(tǒng)計學(xué)工具條并未可用，則需要點擊菜單欄中的Customize下的Extensions按鈕，勾選其中的GeostatisticalAnalyst前置方框，然后Close，地統(tǒng)計學(xué)工具條就可以用了，如下圖所示：(C事?|匚 二二?隼加?一一一.RiEifeMwwUrr—ki口巴過色ft*丁小比舊― v>iOIwaH9(C事?|匚 二二?隼加?一一一.RiEifeMwwUrr—ki口巴過色ft*丁小比舊― v>iOIwaH9UnQ]Bad^iard.C-ndEapiBT#DriJiQSSB- -JUrkHaRiEifeMwwtsc4*r*kioe■昌g?穌NFiQ曾XT^it-accrtra*kdij廳QSSB- -JUrkHaRiEifeMwwtsc4*r*kioe■昌g?穌NFiQ曾XT^it-accrtra*kdij廳-Q]Hi^ginnd.Cnd|lfF?Mi/drl3e*jaxJ點擊已經(jīng)變綠的地質(zhì)統(tǒng)計學(xué)分析工具條，對Import金統(tǒng)計分析（這里盡量使用簡單克里金估計類型，下圖所示: 理i-miu必皿"10?！?toV>_Surface_Data進行簡單克里

為了方便后期模擬使用），如/■*1uvUFxUtiMihnsn￡*aaaiMMBlAi?Hudclrt??jK4ai?iEhHmminiHY?,iXMJ如修修5ft口bitSdeHctfr?*,w?herwlrwitc"vxmetoEridH口wrgtrtrEmti■mRnlc/Ei^nManHWrElEiiiBiHilEPu?>atotrtw?TKkt<EBlJkfcRTh SE^dO1^UHPiAKt^tTflT(^ i .「什HTICff'rrRBT 愴 卜H>3StffcXTFT^d'OO概IwiFRiSP"TIknJ^CTRFTlW,BriF hjtBdMX11注意選擇正態(tài)變換和去趨勢，因為一般數(shù)據(jù)通常不符合插值所需的標(biāo)準條件:

Lowwiitik"如心DOEtWias■熱如沿tev^wi軻3呵&Lowwiitik"如心DOEtWias■熱如沿tev^wi軻3呵&?kfirw如HdreJwiBuntaCw[ImWa|-<iu_]m：>]]rwi |-ciKriGiMiMrUrciJ-n-3drd七力炯cupJaf5L“1制TJXCliaSKfE^3HUflKiTJd(n^na.□gUfiTLHK^rnhsfitvM '-tarn￡<?■Oo^jite'brt?r舊曲??比?flbr》DB1厘IWjsflfiW WW2*301a)h3/7"*百—‘叁?mv注意漂移階數(shù)，一定要選擇合適的漂移階數(shù)使得去趨勢和轉(zhuǎn)換后的數(shù)據(jù)盡可能的呈現(xiàn)出正態(tài)分布的特征，如下圖所示（例如本次模擬選擇2階多項式漂移）：協(xié)方差函數(shù)和變異函數(shù)模型的擬合曲線如下圖所示:選擇合適的模型，并計算其基臺值，偏基臺值，塊金值以及變程等參數(shù):

GeneralExaminebivariatedistributionFalse

VariableSemivariogramOptimizemodelModelNugget

EnableTrue

CalculateNuggetTrue

Nugget0.0071928488230167

MeasurementError100Model#1

TypeGaussian

MajorRange255.65093131342041

AnisotropyFalse

CalculatePartialSillTrue

PartialSill1.2536437163695078Model#2

Type<none>

Model#3

Type<none>

Lag

LagSize29.37115822112131NumberofLags12得到最終結(jié)果，如下圖所示：預(yù)測值和測量值有很高的相關(guān)性，如下圖所示:預(yù)測值和測量值有很高的相關(guān)性，如下圖所示:Predicted-10一，0.5S50.4271.370.8990.T421BS515281Predicted-10一，0.5S50.4271.370.8990.T421BS515281213105B1S420.271130.459 O,B05 1.15 1,496 1.842Meazured-102as3S“wnrd工卻購?at?>Cro:d'/jl川m■h，htnin>o?das3S“wnrd工卻購?at?>Cro:d'/jl川m■h，htnin>o?dSivutfcicfrraiEXEn■F2172”，Ini區(qū)樹皿―1中聯(lián)J"Li.U73uii-.a...工35COL.，3黃759.(n31^3JLI1,.3...工FH]?.?■a.IWLZi??f..E3MEd”2MB故對，濟班由一看照5西宗“TM取胃A..上皿中團一*X：lLkKH._inme-XuEl..]...Z.3H1EM..E.EE11349n..EH.5HJldia...M5Al直畤“TU峰原?MCUi3Lli-.廂也..i"nne-包31-.卻出，一fnM■i..4.WI3i]5U..伯21LZJLI9..IA3WLraHL+..*5j9lH?-I,-fM八。3&JS..樣..ZlU2nn..2?i33_.4.L^X'MCT..F231ELFHLUWSIgr咐213Hll.i1初次四溺TM*ILif..LaL-nSOdE..在MERV的..(nM.>=3LIT-.211B3XVM-.ab?*?FHT._E3aB3L]]jiQ..Z3?L1*W^i.ssamswns四6HI5<01H蟒吟“N<iif.]L1E..0一.2.U3n?!3..i"n3113LJ3,.a...上3axM.，皿4M■?：>□.-Yn￡LEar.iT..■a..ZKca?-H..伯UE2JLI1-..LSDM^nb]I.XM73139m..*5MSHi2,b-1-啊jL命熟5麗?.?iSmMi3...上櫛占U7C.I.：L-fMu73￡1.-Tvine-35.M,.a...z-imKi..Mfr，,陞-i..3.WH1JSqWWMflM*5融工2lHm?.DTMH.I]LH..3...U3>￡?aK3-.j.MPKM..hl3121_..LZMJrH]..[4D3K^na.E?.e-B.M..Q..LPM9M1nesML9電?/網(wǎng)網(wǎng)■ft]<ZUEEW513M<a女組wa百分k■cum口皿FIXARTas?w算Ki海H湎也門HKmsgg的認3皿口事HHE0aWMlM9mTJfl在八國嬌SB陰甚0.7M?4C(Zl?￡i?a-mFjJMMWgJE!2E黃則京凌財弋+S0EEK13<i^6^5]U刑%4.MX4:U3^G?i-Nx￡>1142?【口1出加O-WJWL'JSllHJE!*蝴第4刀加泗認出旺必的港隹■D.39LIU33K7ZM94.IML!31<I4CC<T1Z3■ftqjBMiUSnKUH?5期號:3期21l泗swgmLR?5：HEMiSR以臥EgI文"呼也B弟聞向好的遜I」也*力的宣中ErQ303MJLMTWIW-T。序毋M國附MJiimTPirdWX^LE<ONkhffita|rnti LWKui堂聰君??P回考"EH干"E工1 h_y=M5?八根據(jù)本文數(shù)據(jù)特征，這里可以改變漂移次數(shù)，生成兩種漂移階數(shù)的結(jié)果：

Kriging_Estimation_Result_Second.dbf

Kriging_Estimation_Result_Third.dbf通過對比平均誤差，均方根誤差以及平均標(biāo)準誤差，發(fā)現(xiàn)本文數(shù)據(jù)利用三次漂移

多項式所產(chǎn)生的結(jié)果更精確。然后通過對簡單克里金估計結(jié)果進行驗證和預(yù)測，導(dǎo)出克里金估計結(jié)果shp文件，其中需要用到背景網(wǎng)格文件，如下圖所示：ChJftQi!'QflJMiC035同皿,ChJftQi!'QflJMiC035同皿,5iaM^jEf*O*注意下面這個對話框中的第二項Inputpointobservationlocations應(yīng)該選擇背景網(wǎng)格（背景網(wǎng)格是之前建立用來展示克里金估計結(jié)果）

nuiuu'.Uee.It—cj汨■Iknuiuu'.Uee.It—cj汨■Ik■2L94Mti39?黑KEE/』HE』"?98.55itau:cruBrK|ciriH>xhX?十匚匕匚于:3<ll-?■J3-M33LH■^>2244191-穌■LMLUd■H用“￡TT2.55293I7-73.IE5E3K?丁幽HE2?斡U2U5■電UL”d5-LLP/ini打■LLT.L!?in-UI7.CM3E&■L?7jCK.IlI5-出；最后開始準備執(zhí)行序貫高斯模擬，如下圖所示:在ArcCatalog里面新建立一個Sequential_Gaussian_Simulation.gdb數(shù)據(jù)庫文件，如下圖所示：

ReEd±u—Go CufUMv^K Help也口型.?洸游出加◎因手：?曰扣，？1Laulim：C■9Hii：3iMK!gBsaL0-T￡miT[*J“jnL3RlfnmLinjShE*LElLR3KMLDocurriinl:ReEd±u—Go CufUMv^K Help也口型.?洸游出加◎因手：?曰扣，？1Laulim：C■9Hii：3iMK!gBsaL0-T￡miT[*J“jnL3RlfnmLinjShE*LElLR3KMLDocurriinl:廿^ropwriiai.jjHic-CwdrlalMM?IBrncrdiSiodal^aiirdoBnedDUElabk■XtdrHSbaUiKV,.Ecnpc-ifafAiftdriihhL&uui.謹慎進行序貫高斯模擬參數(shù)的選擇，序貫高斯模擬需要之前做的簡單克里金估計結(jié)果作為模板來進行條件隨機模擬（注意這里需要選擇條件數(shù)據(jù)，就是原始測量數(shù)據(jù)），如下圖所示：GMilrgaviE*h?whpErfdi■wr”改畫iirq3?山七1'GMilrgaviE*h?whpErfdi■wr”改畫iirq3?山七1'人*n.d.-a"tiQiJL-n庭Irnpari.iMHicl.QAiI喉wpe喇峰Toon口近既cTads|f^Tadsgl|rH￡I1LEiSB■AlMiuCjukigi -CA.UmYp1rlD■:丑 JSakJrHaMfHju?m?l上>uU-n」ii*i」?：qri^& s.“.itReEd*xgGdi CufbO^ise省描MH/i心口流工部排圾u0由國事口：>,?.”=LmHe:『i.lwi3!dpi!4M!：?sHF,T-ir-riC：''iAr<#D：-_JsfJ!ii*i<L_(Rmi _k_^+j121al.rm ?xmtwwmFj-tM-ifJLuuipLI4#R由RJgrc?ertw口與C^JJtfrnM tKuntri eO35_Le白UrfHrl.Maa.ihp屈白tr^sE一么QQillWlJ^hflL田陽lI□uipari_5inuIm￡o^J^Muhc[±ijt*Mriid.GfiBniTMn_^iiTULriKn?Td/Lk？e,由gWfcWUStiW?i%DuuJwdw出sc3毋CAi**wwT^pe三：flhiG?cd■山iiFIc■Gc^c*DbraeRliGabdmtaMriFk6wd?ifai??EleGZHAb出ehl?MiUi'Diul-Esbdiki3i?txwQsuc..OvtiL.Faswasias..MuirDiul-Ra岳'■8田,匚：1RliGabdmtaMriFk6wd?ifai??EleGZHAb出ehl?FIc■Gc^c*DbraeRliGabdmtaMriFk6wd?ifai??Fiek4C-jV?bH4txwQsuc..OvtiL.Faswasias..MuirDiul-Esbdiki3i?txwQsuc..OvtiL.Fasw&5l3Enhl?■^■PdK^MiFIc■Gc^c*DbraelihGiodiiutaiM*MuirDiul-Esbdiki3i?Z<wQjiiJiL.Fk6wd?ifai??EWmEh由用eOvtiL.Fasw0-5Ml?hl?■^■PdK^Mi*6ccd〃山mtMuirDiul-Fx山e口■口巴>RilGcd*山*￡*總IEW口11111[1_RaiibrOvtflL.FieGewiysbiMeFaiw6方！l匚R? ?dqdv>z-W丁整園千iE|"]㈣Q 口皿上=—屈”威嬴?叫|FkN5urft!!._D*?WltEpe?叫|FkN5urft!!._D*?WltEpe的 TMK叱e??mWMm他Ejima丫&-ToTaHePwttf^McToakitFmot^nTad*S??VfrTpghi■,V^piiai■■nd溶lorii

事Spiii^'SznnciiT?di■Ilr?ddrgArulaaTade再在ArcCatalog里面重新建立兩個gdb數(shù)據(jù)庫來存放克里金估計和序貫高斯模擬所導(dǎo)出的Table格式文件，即：Export_Kriging_Result.gdb；Extract_Simulation_Results.gdb再在ArcCatalog里面重新建立兩個文件夾來存放克里金估計和序貫高斯模擬所導(dǎo)出的shp格式的dbf文件，方便后期的excel處理，即：Export_Kriging_Result_shpExtract_Simulation_Results_shp導(dǎo)出克里金估計結(jié)果到Export_Kriging_Result.gdb數(shù)據(jù)庫，并得到Table格式數(shù)據(jù)，如下圖所示:取口>3.刈肥 AralyTi-1t■口>0i?5*ndflHaie-s□1dD^Cmbhih.，:二kuWWeSM即而中U<*卬他jqjwU.TghiJsi^LxLem:0-JW55772■1L7157IE?-1J?J■i靖m電恭.i$Ji.取口>3.刈肥 AralyTi-1t■口>0i?5*ndflHaie-s□1dD^Cmbhih.，:二kuWWeSM即而中U<*卬他jqjwU.TghiJsi^LxLem:0-JW55772■1L7157IE?-1J?J■i靖m電恭.i$Ji.1n ari=￡_口nigwwi?IS-ELMO?招l(wèi)■士1口必，豺?黑2?斑郎佻曲?總之ZCW?「。際*Z?nFcMihfTAIUJM0，r-J￡一IPFt上i呈Ikwin_$urlace_DM4S!不月IMgiOJM55772國用加一Viht也■班？?扎川國理?上l皿印*-瓦*?:ifl-lMLIM-融雪imm-7國%如父■wseitk立晶也乖■Si112345-117.167337■AL7ii7l5T-147JDUILK■L-<7DMLE5--LM-Z再從之前建立的Table中導(dǎo)出克里金估計結(jié)果到Export_Kriging_Result_shp文件夾，并得到dbf格式數(shù)據(jù)，方便excel直接打開處理，如下圖所示：ak呵咚）rjLM5：L53-WUUM3■如L2加ak呵咚）rjLM5：L53-WUUM3■如L2加■LL7.1SUJ7BlLLTWU3-UMU1I■1bHr^LI：■血3□Sep.AduOG>vkc_口中呷IYDu一5口1dL553iaiJlureurd注意這里文件格式應(yīng)選為dBASETable：irjl-H3jL53-U.33314：-?方溫L234E-LL7d注意這里文件格式應(yīng)選為dBASETable：irjl-H3jL53-U.33314：-?方溫L234E-LL7d￡lU7■LLT.Li？ll?■團口”二第■KT3W-EMJ_□it'nFg-53.552^67■U.553i997-從序貫高斯模擬結(jié)果中導(dǎo)出云圖到Extract_Simulation_Results.gdb數(shù)據(jù)庫并提取值到Table中方便進行編輯（這里需要從ArcCatalog中把所有1?100次模擬結(jié)果全部拖入ArcMap左邊框中，并在ArcMap中一個一個導(dǎo)出成相應(yīng)的數(shù)據(jù)格式，有點麻煩，但沒有辦法），如下圖所示：RiErifeVww&g^r<?kihwrt￡?Udw"De■含ft1葉由，耳二小守曾X事?*9gifOfgEWS flMr七回奇目idLjm.LL5LL2E 1axSQ『4 同1環(huán)f|1岷$衿3；-RIITc.S即.Hi6iimuLm:-UKL34E=El的i值視3:-NIL口x(小11BM07Lw:-LUL2EFiBSBWI IWriferwiiM>Jp-叱4*：'B1fl.：，垣1t飆psWiiliMiJ嗎4*卡比出iQ^C^hp-Afarr^-G^hc?eris'iArtGIS*AxGK__Cb'Srr^J!r,?NL&w>?n_?"■h^i■"ndbNiE3:嗔*向gGdyraoH匚耳CmCixtq+Wd—H*aaiiBtMiilfEgggB用季?口AW?L:<&bc^-'^i"1■liLrW'!k?iiVAElikiiiriCrkJr!?G35jKH?：■9]jGnruiuJinlBiit.,JEjlrirt flsTJh-f^lr?ZkfWiriJTJj^rrlrdCrfA+fci4rrmq|fajiaraini?nhUrw-■QUUAhl,WiW>師uMhlS,HI5*L4,UExk?MM,MltallTyp?FJh*Z*=dauhflL+FLiubiDiui—H?■□wduuiiJiw-RiudrDbim_口.GaMdluKl—kKlrirDlLM—a?SjffliJKiiMM-Rafler口?吃?口十一fcWtrUi-AcrjCQEkAcJMX-Dikrh：Ai 口 RflHirDvlo.Ar！amMMmmJ."Drtan-.■^i'ZwodntME*rAmrDnn_FM?5ceU0U3■口iwDaui—FJh*^KdautaflM-FLlUilDiUi—MGwdlCiJUL*RjUdf口BUM—口.GaMdlUiU—kKlrirDlLM—FijftiiJKiiMse陽h口nw-口總?曲陽MiM!"皿0Di—Ae<Knckta4M>rft?lE!Dikrs=Ai RdiHirDvlo.He