ArcGIS中的空间插值和面积计算

1、 #说明:本文阐述了空间插值和污染面积估算的方法，供群内交流学习用，若要用于商业用途或转载，请与原作者联系。本文若有不正确之处，敬请指出！一、空间插值插值方法种类很多，每种插值方法里参数也很多，至于哪种最好，没有定论，只能根据需求以及制图的效果来选定。建议：插值效果图与网格图进行对比，哪种效果最接近网格图（能体现局部）而且又能反映整体趋势就取哪种。 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：可以，但精度会受到影响，看平滑后的效果来决定是否进行平滑处理。建议不做431K633MSLTDE.iZW.-K7SS4I瓏歆柩rffiSi110.稍稱CcoTti.tiflirTixsEd:IDTInrtcrpoli.tion:StcjiI

3、n-f2Srtraiuctrcsj1*?请问：此处有可选smooth，可以做进行平滑”处理吗?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/开始倉FC-BD*r“id电ArcMa.InverseDistaneeWeightimFileEditViewInsertSele匚tionTooIejVirLdowHelpLayer1ine_regiorL_FrojectMethods:Blinejirov_Froject日lint：jirovs_FrojectDatasetInputdata:Attribute:Yfield:日city_cnty_FrojectoBcity_region_

5、Froject日line_cnt7_Froject叩njkxiiJTii丁?:童卓*口A险fliflxtitiTtird細订E”山佔GlobalPolynomialInterpolationLocalPolynomidlInterpolationRadialBasisFun匚tion生KrigingCokriqinqFC_YH_Freject日Layers日尊污染面积计算1-E0.面积计算- #3、扩展研究区域 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阳亡订*rcd.1li11*4XOLBSHII*Till创.4HlhKI-THLTEulitrvpsyppEC4ftictatfi|*1|MrwiqF!*?-f #J?3J.2T4绘5、这是采用“反距离权重法，1次方”来插值的。=_云可选用“局部多项式”或“普通克里格插值”方法来试试，看哪种和网格分级图更接近些。但无论哪种方法聚类误差可能都

7、较大，一部分高值可能被掩盖。二、下面转成栅格图层再进行分层设色图制作，这样精度较高，且图层可用来进行面积估算2.1、导出成栅格图层11ne_regiorL_rroject13Vgrti匚吕”pmHorizontal:_FrojectOutputraster:E:5OILbackupp_cd1ine_prov_FrojectCar（2.3和2.4均非必要步骤，只是为了另外的处理或制图的美观性。如果是为了制图的美观性有可能这两个步骤会弄巧成拙，是否须要请根据具体需要和效果来定）2.3、并可对栅格图层重分类，生成新的栅格图层如（ah_cd）面积计算-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髓黑国o-比OWN伽1护伽亦I2口由輕1#|t；MZ俯二田曰匕：鳥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请问：此处重分类又该如何设置？有无相应依据？同上也面积计算-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时如上UtttT-.TM9tTl.MUiTi:3?Wf：.7KJ琢那D诫Cl曲勰幻寒：国田：:11rUEWtiLvRflSabn&j|pliii!|制窃|叩克iHLeiirMEiIHcairirn9|止廿jj冠心i.dmuldnkiiiEtAdkireMefiljMG#ir$IMMLf#.|9d|常lifEndii

11、Td5lH4nbaFeatM*|sa）号出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再掩膜（或切割），最后成图掩膜可以提高作图效率 #:趣册擴入箱式妙xAd痂0）口0:-卫鼻厶_：_黑匙耳,m*jj_b.曲画童上躯乂血一山冊堞釦.陆.妻一二科”玉兰二戶g*昌.A茅二画H邑葺基*i三氓參M|_L3I一1I?1I2I3it4fia9!lD71iiDi耳1&!屬|丽116|!宿|亦打|越|曲丨24侮|孙静心估|却対|穴曲|如第|3$山|如到占呦切輻丨2ArcL磁令樨吓？1-ar*iHicrosoftLiindonsfP

14、?.5.3.2GOSJ(C)阪权所也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盘侮-.丈件圮】aiE現国iftflAa?林式小:宾m*ai帘助阳、T命专1亂翱口巧盘如快4- # # # #巴J也-事&岑覽画昱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峑上丿二虫：勒=、吾圣、工、丄*再运彳丁SzPollueteArea.exe（需要.NET2.0支扌寸，如无，请安装）$ # #网格册一.口Kii-rciSo.I帝专埠. # WCI*&皆棋v平乜梓疸汇屜最后得到污染面积汇总表HIS #117IS20按*udt-ft.-审jhwtH将用地类型数据先转成Database格式，然后采用下面命令整合成六大类，并进行投影转换Updatehldisebtmpcl=left(LD

18、20Q6_lD1),曷（由于开始是直接在DOS下运行，不能拷贝，我自己根据做的在你的文档上修改的，可能大小写等有点不一样，只是表达了意思，这些也不需给作图人员，请你看看我的做法可对）3.1Cover坐标转换（如果转换不了，先不做跳过这个步骤，导出成栅格图层后在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 # #+- # # 下进行转换)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下增加一个字段tmpcl，用来存储六大类信息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