ArcGIS中的空間插值和面積計(jì)算

1、 #說(shuō)明:本文闡述了空間插值和污染面積估算的方法，供群內(nèi)交流學(xué)習(xí)用，若要用于商業(yè)用途或轉(zhuǎn)載，請(qǐng)與原作者聯(lián)系。本文若有不正確之處，敬請(qǐng)指出！一、空間插值插值方法種類很多，每種插值方法里參數(shù)也很多，至于哪種最好，沒(méi)有定論，只能根據(jù)需求以及制圖的效果來(lái)選定。建議：插值效果圖與網(wǎng)格圖進(jìn)行對(duì)比，哪種效果最接近網(wǎng)格圖（能體現(xiàn)局部）而且又能反映整體趨勢(shì)就取哪種。 1.1、 3a口Editw日Validationphape日gfieldn日nUseMODATAvlue日日日0申|l：5(1,288口phapeOil,156B35D色PCYDGeostatisticalWizard:ChooseInputDat

2、aandlethodgridOSgridl6admin_cntyFrojett口-:eldmin_region_Proiect口一：=Lilninjirov_Projectri_EllE4liYl*tBialaiu-hnIpiars鼬iwiz】鼻TinflorHi3pjjduHWbgrr:4Q少Hbvirh：I八UK4(rwiMhW：jbowtTttghb：可以，但精度會(huì)受到影響，看平滑后的效果來(lái)決定是否進(jìn)行平滑處理。建議不做431K633MSLTDE.iZW.-K7SS4I瓏歆柩rffiSi110.稍稱CcoTti.tiflirTixsEd:IDTInrtcrpoli.tion:StcjiI

3、n-f2Srtraiuctrcsj1*?請(qǐng)問(wèn)：此處有可選smooth，可以做進(jìn)行平滑”處理嗎?9|ii-i-Zfc1HUll.fMia&lTe-71+&DbLk?Mririi-hniDL3.c?njarLvlIdsi3nivivilukdlF；arii申|MfekinktfaFaiivrClva-hQi険怦口_巧.b智til-rhFr*!聲Cfvi.LviN&Lkitii1CrivHGmLmFrArftpKuLtflJ!BqKaiklaaahihiji口i*i11:fTrMhvRj|Tip口QpYiEliMi&TijriHpBjf.*T-clx+：QtnluLhilLtKlAoadyiL1Q|

4、.iTrfi*r!tfBrasarinT-ioliQImKiiiTad.1iL!)m.tiribTPbJinalwTlull丼比WaBKrUr*BKT/開(kāi)始倉(cāng)FC-BD*r“id電ArcMa.InverseDistaneeWeightimFileEditViewInsertSele匚tionTooIejVirLdowHelpLayer1ine_regiorL_FrojectMethods:Blinejirov_Froject日l(shuí)int：jirovs_FrojectDatasetInputdata:Attribute:Yfield:日city_cnty_FrojectoBcity_region_

5、Froject日l(shuí)ine_cnt7_Froject叩njkxiiJTii丁?:童卓*口A險(xiǎn)fliflxtitiTtird細(xì)訂E”山佔(zhàn)GlobalPolynomialInterpolationLocalPolynomidlInterpolationRadialBasisFun匚tion生KrigingCokriqinqFC_YH_Freject日Layers日尊污染面積計(jì)算1-E0.面積計(jì)算- #3、擴(kuò)展研究區(qū)域 EilNilYf#tflalaiu-ko.TjiiarS鼬hdgtbnjorjilpCnajvritE.54wci?th帖些|1xIid!LInLbcdSsrinirjr|DidW41

6、llIlE-dewMMNA口回昨jcs2?HngvWma?iin曲r-112s；-323$STIH耐甜師IA：I3o旳yrufaiKrl-faC-hllE.IUidiEA.Edalin|idTa:I.XhkFI|BT*ll114JFfdTrada-ilian陽(yáng)亡訂*rcd.1li11*4XOLBSHII*Till創(chuàng).4HlhKI-THLTEulitrvpsyppEC4ftictatfi|*1|MrwiqF!*?-f #J?3J.2T4繪5、這是采用“反距離權(quán)重法，1次方”來(lái)插值的。=_云可選用“局部多項(xiàng)式”或“普通克里格插值”方法來(lái)試試，看哪種和網(wǎng)格分級(jí)圖更接近些。但無(wú)論哪種方法聚類誤差可能都

7、較大，一部分高值可能被掩蓋。二、下面轉(zhuǎn)成柵格圖層再進(jìn)行分層設(shè)色圖制作，這樣精度較高，且圖層可用來(lái)進(jìn)行面積估算2.1、導(dǎo)出成柵格圖層11ne_regiorL_rroject13Vgrti匚呂”pmHorizontal:_FrojectOutputraster:E:5OILbackupp_cd1ine_prov_FrojectCar（2.3和2.4均非必要步驟，只是為了另外的處理或制圖的美觀性。如果是為了制圖的美觀性有可能這兩個(gè)步驟會(huì)弄巧成拙，是否須要請(qǐng)根據(jù)具體需要和效果來(lái)定）2.3、并可對(duì)柵格圖層重分類，生成新的柵格圖層如（ah_cd）面積計(jì)算-Arclap-Arclnfo-ciy_Frojme

8、t-city_region_Froject-1ine_cnty_ProjectNeighborhijod占tsLtiEt1ce.0.686731-1.3941491.394149-3.5607913.560791-10.196667FC_YD_FrojectLITl凰BI一因因I園alBIEIBIEB一晶一Q:4西因k髓黑國(guó)o-比OWN伽1護(hù)伽亦I2口由輕1#|t；MZ俯二田曰匕：鳥(niǎo)tp峪攀淞鮒警：*_u占二曲鬲汕Jkx.drri5OjOJTMI世“PfifkiaTri-la&lagdi_h+jtLlrcm*mi0.1押比片LicmnurLvlTdjtjjif|vkgpikuidgpd.idi

9、血“41FiilwCla*krrCiJqTHM.bm.taihPr*jCtfiai-lv令11&iLtiitrZ*LriiUCmbLk-住申D*fiTiiiPrtKliidTKaifctiantiiiyClilVMiEacrlwi-fr請(qǐng)問(wèn)：此處重分類又該如何設(shè)置？有無(wú)相應(yīng)依據(jù)？同上也面積計(jì)算-Arclap-ArcinfoFileEditViewInsertSelectionToolsWimlowHelp0.455T560.6867311.3941490.247642-0.3230560.323056-0.34T6T9-0.355719-0.380342-0.455756-0.686731-1

10、.394149-3.560791-0.3476T90.3557190.3803423.560T9110.19666T 2.4做平滑處理（可選）Typenlh*-1電；：Ln：中IMmi0TrnMLi-TtKnc*Hi-ill3FearLacLnanIBprlJ*j*riiCMiM昇R時(shí)如上UtttT-.TM9tTl.MUiTi:3?Wf：.7KJ琢那D誡Cl曲勰幻寒：國(guó)田：:11rUEWtiLvRflSabn&j|pliii!|制竊|叩克iHLeiirMEiIHcairirn9|止廿jj冠心i.dmuldnkiiiEtAdkireMefiljMG#ir$IMMLf#.|9d|常lifEndii

11、Td5lH4nbaFeatM*|sa）號(hào)出hfMdM/n沖1舟口wtt*AlunnwySatelst（4MM）$ufWtLentfh|3jjt.CeCCMiinHiWMri11111BUT店J1dMeHwwmerA:T血e0dB則CrwJtipkl(cctwbT申T曲柘GTE*IndKT4Tng;兩Sfitraiwrf?W仇他二|皿#j尊口k?Edit吃TSpati：1AtlisIvet亍Layer：p工Geo.EtiLtiwticsLLAil：=l17et亍日LayersGRIDCODE口口口口pc口同口liefinitionQueryGerLer：alSurn-ceShow:Feature

12、sCategoriesQuantities-0p_cdV：=lLusHigh:10.Low:0.0：El0InverseDist：Fredictio:pc_yd_Pi-0.01666T0.2476420.3230560.3476790.3557190.3803420.455T560.686T311.3941493.560T91B0FC_YD_FrojecGraduatedcolorsDispl：yLabelsSeiectionClassificationNaturalBreaksUenkslDrawquantitiesusingcolortoshowvalues.GraduatEdsymbol

13、sProportionalsymbolsDotdensityChartsMultipleAttributes確定取消I弓IElci17_cnt7_Froject2.6再掩膜（或切割），最后成圖掩膜可以提高作圖效率 #:趣冊(cè)擴(kuò)入箱式妙xAd痂0）口0:-衛(wèi)鼻厶_：_黑匙耳,m*jj_b.曲畫(huà)童上軀乂血一山冊(cè)堞釦.陸.妻一二科”玉蘭二戶g*昌.A茅二畫(huà)H邑葺基*i三氓參M|_L3I一1I?1I2I3it4fia9!lD71iiDi耳1&!屬|(zhì)麗116|!宿|亦打|越|曲丨24侮|孫靜心估|卻対|穴曲|如第|3$山|如到占呦切輻丨2ArcL磁令樨嚇？1-ar*iHicrosoftLiindonsfP

14、?.5.3.2GOSJ(C)阪權(quán)所也19&5-2PQ1MicrosoftCorp.C:DociJHnt&andSettingBzMdarcF1niHii$F1119E1DE匚rFar4nIMrE3*n.rHTwin.T(JUN19B7).FARSection52227-19(JUN19B7).and/arFARSection12.211/12.212CoiviiarciBlTvctinicBlBiB/CcuHputoTSoftMaro1)andDFARSSection252.227-7915(NQU1995)TechnicBlDataand/orDFSRSSection227.7202Conp

15、uterSoftware,asapplicBble.ContrBGtor/MBnufBGturerisEnuiroiniaentBlS-ysteiiBResearchInstitute,Inc,3B0NewVorkStreet,Redlands,CA92373-0199,USA. # # # #再輸入ARC:&WO:SILBACKUP占劇叱丸13?H1H15/ir.甌球再輸入ARC:&runtjgdldwrmj.aml*U盤(pán)侮-.丈件圮】aiE現(xiàn)國(guó)iftflAa?林式小:賓m*ai簾助陽(yáng)、T命專1亂翱口巧盤(pán)如快4- # # # #巴J也-事&岑覽畫(huà)昱j堇這-）三三奉Liseofworkspac

16、esdtlocation：E：50lLBRCltUPAudilablworkspacesRrc：ue：5oilbackupEARNING：N$ulocationisnotaworkspace.flrc：1mThissoftwareisduplication,andrMtriclionsassetforthinFARSection52.227-14AlternateIII(g)(3)hhNONEhhRrc：赴runCjgdldwrmj.amiprovideduiehrestrictedandlimitedrights.u.disclourbytheU.S.Gouernmenl:aresubjec

17、tto兀命專捱喬的KTtIT16:05：PDJZ00&) # # # #UseOUTPUTtodctuictheoutputpiojccTionandEND墜削婦哎自螢匱血”、口o峑上丿二蟲(chóng)：勒=、吾圣、工、丄*再運(yùn)彳丁SzPollueteArea.exe（需要.NET2.0支扌寸，如無(wú)，請(qǐng)安裝）$ # #網(wǎng)格冊(cè)一.口Kii-rciSo.I帝專埠. # WCI*&皆棋v平乜梓疸匯屜最后得到污染面積匯總表HIS #117IS20按*udt-ft.-審jhwtH將用地類型數(shù)據(jù)先轉(zhuǎn)成Database格式，然后采用下面命令整合成六大類，并進(jìn)行投影轉(zhuǎn)換Updatehldisebtmpcl=left(LD

18、20Q6_lD1),曷（由于開(kāi)始是直接在DOS下運(yùn)行，不能拷貝，我自己根據(jù)做的在你的文檔上修改的，可能大小寫(xiě)等有點(diǎn)不一樣，只是表達(dá)了意思，這些也不需給作圖人員，請(qǐng)你看看我的做法可對(duì)）3.1Cover坐標(biāo)轉(zhuǎn)換（如果轉(zhuǎn)換不了，先不做跳過(guò)這個(gè)步驟，導(dǎo)出成柵格圖層后在toolbox-Il注M：=LrL：=LgementTools+Ilnta.Comp：ii_ison+Hatabase+毎iDiEcoruLectmd.EilitiTLg+磅jiDistr讓utmdGeoda+Dumains+iiFeaturmClass+彰iFeaturEE+每iFields+iiFileGeoda+Gener：elL+

19、Gener:=lLization+彰iIndexes+蒔iJoins # #+- # # 下進(jìn)行轉(zhuǎn)換)Arc:projectcoverld2006ld2006_new*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*vt*、t*、t*、t*、t*、t*、t*、t*、t*、t*、t*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*rT*rT*rT*rT*rT*rTxTheINPUTprojectionhasbeendefined.*、J、J、J、J、!、!、!k1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*

20、1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*1*UseOUTPUTtodefinetheoutputprojectionandENDtofinish.Project:outputProject:PROJECTIONALBERSProject:SpheroidKrasovsky000.000000.000:30.500000.000:33.500000.000:11700:0000.00000:5000000.00000:0Project:UNITSMETERSProject:PARAMETERS1ststandardparallel2ndstan

21、dardparallelcentralmeridianlatitudeofprojectionsoriginfalseeasting(meters)falsenorthing(meters)Project:endArc:buildld2006_newpoly3.1在workstation下增加一個(gè)字段tmpcl，用來(lái)存儲(chǔ)六大類信息Arc:wdsoilArc:ae/Arcedit:&woe:ahArcedit:editld2006_newTheeditcoverageisnowE:AHHLD2006_newArcedit:efpoly75061element(s)foreditfeaturePO

22、LYArcedit:additemtmpcl88iAddingTMPCLtoPAT.Arcedit:itemsCOLUMNITEMNAMEWIDTHOUTPUTTYPEN.DECALTERNATENAMEINDEXED?1AREA818F59PERIMETER818F517LD2006_new#21LD2006_NEW-ID25TMPCLArcedit:saveSavingchangesforE:AHLD2006_newSavingpolygons.ForPALofE:AHLD2006_new:75061polygonrecords(s)written.fromtheoriginal75061

23、,75061addedand75061deleted.Savingarcs.141234arc(s)writtentoE:AHLD2006_newfromtheoriginal141234,0addedand0deletedReopeningarcs.Savinglabels.75061polygonattributerecord(s)writtentoE:AHLD2006_new75060label(s)writtentoE:AHLD2006_newfromtheoriginal75060,0addedand0deletedReopeninglabels.Reopeningpolygons.

24、BNDreplacedintoE:AHLD2006_newSavingsettolerancestoTOLfile.Re-establishingeditfeaturePOLYGON3.2整合成六大類Arc:aeArcedit:editld2006_newTheeditcoverageisnowE:AHHLD2006Arcedit:efpoly75061element(s)foreditfeaturePOLYArcedit:selld2006_new-id1007870element(s)nowselectedArcedit:caltmpcl=1Arcedit:selld2006_new-id20andld2006_new-id304965element(s)nowselectedArcedit:caltmpcl=2Arcedit:saveSavingchangesforE:AHLD2006_newSavingpolygons.ForPALofE:AHLD2006_new:75061polygonrecords(s)written.fromtheoriginal75061,4965addedand4965deleted.Savingarcs.141234arc(s)writtentoE:AH