城市表层土壤重金属污染分析

城市层土壤重金属污染分析(建模大赛) 1.Excel转为数据库,展点,建TIN 转数据库:将A套题中Excel(将“高程”改为H,不要单位)另存为“采样点坐标.dbf” (为了能更顺利地转为泰森多边形,将采样点坐标数据另存为TXT文件。) 展点:在arc map中加载“采样点坐标.TXT”——右键:display x,y date——X field:x,Y field:y ——不要坐标——OK。将展点后的文件输出,展点后文件右键——data——export data——自己文件夹下“采样点” 建TIN:工具栏3D Analyst——create /modfiy TIN——create TIN from features ——采样点——Hight source:H——自己文件夹“output TIN”——OK 2.依采样点划分功能区 划分功能区:泰森多边形V oroni Voroni:工具栏spatial analyst——distance——straight line——distance to:采样点—output cell size:50 —create allocation自己文件夹:vo2(多边形) 将vo2转为矢量:spatial analyst——convent——Raster to Features——input raster: vo2—output geometry type:polygon——自己文件夹“采样区”Save 点给采样区赋值(Join):采样区右键join——join data from another layer based on spatial location以位置赋值——整个属性表——自己文件夹“功能区片 ” 功能区融合(Dissolve):工具箱Data Management Tools——Generalization——Dissolve——input:功能区片output:功能区,lables标注功能区 TIN转矢量,剪切功能区:3D analyst——convert——TIN to Features——input:TIN conversion:interpolation zone to polygon——“tinfeatures1” 用tinfeatures1裁剪功能区clip:工具箱clip input:功能区clip:tinfeatures1——OK lables加标签 3.各功能区各污染物含量分析 1)“重金属.dbf”与“采样点”join:污染物去单位另存为“重金属.dbf”,arc map 中加载“重金属”,“采样点”右键join and relates——join—join attributes from table —编号重金属——OK,“采样点”export data “污染含量” 2)栅格内差,对污染物分析:spatial analyst——Interpolate to Raster——Inverse distance weighted反距离加权内差——input:污染含量——Z value field:(污染物) AS power:加权2 (经优化1.5较好)number:20 cell:50 output:As (将各污染物进行分析,将功能区设置透明值,污染物与功能区对照),可将污染物分析制成三维3D analyst——convert——Raster to TIN——OK 4.将各污染物重分类,计算,生成综合污染物分析 1)各污染物无量纲化,重分类:spatial analyst——reclassify—input raster:(污染物)Zn——reclass field:value——classify:5级——自己文件夹“Zn_rec”(“污染物_rec”)——save 2)各污染物叠加:spatial analyst——Raster Calulator——calculation=As_rec+…+Zn_rec 5.将污染区统计裁剪成统计范围形状(赋值为1,相乘) 1)TIN的矢量图将ID改为1:TIN to Feature“tinfeatures1”打开属性表——start editing——ID由0改为1 2)“tinfeatures1”矢量转栅格:spatial analyst——convert——Features to Raster Input: tinfeatures1—field:ID—output cell size:50——自己文件夹“wuranquRas”——OK 3)栅格图与综合污染calculation相乘(ID为1的地方保留,其余消除):spatial analyst——Raster Calulator——wuran=wuranquRas* calculation 6.将功能区栅格化,与污染物进行统计分析 1)将功能区由矢量转为栅格:spatial analyst——convert——Features to Raster Input:功能区—field:功能区—output cell size:50——自己文件夹“gongnengqu_1”——OK 2)统计分析:spatial analyst——Zonal Statisic——zone dataset:gongnengqu_1—value raster:calculation—chart statistic:Mean——自己文件夹“zonghewuranfenqu.dbf” 将每一项污染都进行统计分析eg: zone dataset:gongnengqu_1—value raster:As——“As_fenqutongji” 7.污染物与地形坡度、坡向统计分析 1)坡度:依TIN坡度分析:spatial analyst——surface analysis/3D Analyst——surface analysis——slope—input:TIN—out size:50—“tinslope”——OK 重分类(降低复杂度,便于统计分析):spatial analyst——reclassify—Input:tinslope —reclass field:不选—classify:5级—“slope_rec” 统计分析:spatial analyst——Zonal Statisic——zone dataset:slope_rec—value raster:calculation—chart statistic:Mean——自己文件夹“slope_wurantongji”—Save 2)坡向:依TIN坡向分析:spatial analyst——surface analysis/3D Analyst——surface analysis——Aspect—input:TIN—out size:50—“tinaspect”——OK 重分类:reclassify—Input:tinaspect—reclass field:不选—classify:10级—“Aspect_rec” 统计分析:“Aspect_wurantongji” 统计分析结果路径:spatial analyst——options——General——working自己文件夹 8.Arc toolbox中生成列表与统计分析图 1)Tabulate Area:平坦区:Raster Calulator——pingtanqu=slop<=5(0为非平坦区,1为平坦区);Arc toolbox——spatial analyst tools——Zonal——Tabulate Area —Input:wuran—Input:pingtanqu——OK;source中表,打开属性表可看统计完成0,1的值 Slope_rec:5级坡度的污染物列表分析,Arc toolbox——spatial analyst tools——Zonal——Tabulate Area—Input:wuran—Input:slope_rec——OK;source中打开属性表可看5级坡度统计结果。 2)统计分析图:Arc toolbox——spatial analyst tools——Zonal——Statistics as table —input:slope_rec—input:wuran—OK;在source中找到统计分析表,打开属性表options—Create Graph—value field:Mean—X field:value—Title:污染值坡度统计图——OK;“污染值坡度统计图”右键Save自己文件夹 9.分析图输出 比例尺:Data view:layer右键属性General——Map:centimeters Display:meters;layout view:Insert——scale bar 横向图幅:File——page and print setup—landscape;图右键——distribut——fit height to margins 标题,指北针,图例:Insert——Title,Noth arrow,legend 坐标:制图区——右键——属性——Grids——New Grids 输出:File——export map