ARCGIS纵投影图的制作(以S2为例)(second)

ARCGIS垂直纵投影图的制作及储量估算流程

How to make vertical projection map and estimate reserve by ARCGIS

——以沙沟S2为例（Take S2 IN YING Project For Example） 1数据整理（data preparation）

1.1钻孔样品及样线数据（Drilling and Tunnel samples）

新建一数据库，沙沟矿区命名为 shagou.mdb，在该数据库中包含以下6个基础表格；

1.1.1钻孔表（collar）

1.1.2测量表（survey）

1.1.3分析表（assay）

1.1.4矿脉特征（vein）

1.1.5样品坐标（samples_coordinate）

由Surpac中导出的样品坐标

1.1.6各矿脉的投影参数（projection_parameter）

1.2数据检查（Data check ）

1.2.1在数据库中检查（Data check in the database）

通过一些列查询操作检查数据的完整性；例如Collar表中的钻孔在Assay、Survey表中有无对应的数据。

1.2.2在Surpac中显示数据（Display the drillhole in Surpac）

在Surpac中显示钻孔数据，可以检查数据的一致性，例如；采样位置超过了钻孔的最大深度或样品的起始位置互相重叠。

1.3投影变换及样品组合（projection transformation and sample combination ）

1.3.1新建一个名为Arcgis_shagou.mdb的数据库（Create a new name Arcgis_shagou. MDB database）把shagou.mdb中的6个表全部链接过来。

1.3.2样品组合（sample combination ）

创建一查询，命名为averagecompositesample，计算出组合样品的品位、长度及坐标。

设计视图（design view）

SQL视图（SQL view）

SELECT assay.hole_id, assay.vein, Avg(samples_coordinate.x0) AS x, Avg(samples_coordinate.y0) AS y, Avg(samples_coordinate.z0) AS z,

Sum([assay]![ag]*([assay]![depth_to]-[assay]![depth_from]))/Sum([assay]![depth_to]-[assay]![ depth_from]) AS ag1,

Sum([assay]![pb]*([assay]![depth_to]-[assay]![depth_from]))/Sum([assay]![depth_to]-[assay]![ depth_from]) AS pb1,

Sum([assay]![zn]*([assay]![depth_to]-[assay]![depth_from]))/Sum([assay]![depth_to]-[assay]![ depth_from]) AS zn1, Sum([assay]![depth_to]-[assay]![depth_from]) AS length,

Count(assay.samp_id) AS [count]

FROM assay INNER JOIN samples_coordinate ON assay.samp_id = samples_coordinate.sample_id WHERE (((assay.cal)="Y"))

GROUP BY assay.hole_id, assay.vein;

数据视图（Data view）

hole_id vein x y z ag1 pb1 zn1 ZKY6-7 S21 524,071.62 3,782,846.72 507.04 0.01 0.01 0.01 ZKY6-7 S7 524,069.05 3,782,809.94 566.03 35.39 3.16 5.22 ZKY6-5 S7 524,069.04 3,782,809.86 579.95 0.01 0.03 0.08 ZKY6-5 S21 524,071.99 3,782,851.99 551.46 1,741.09 10.01 24.66 ZKY6-3 S7 524,061.17 3,782,796.08 584.98 0.01 0.06 0.03 ZKY6-3 S21 524,019.90 3,782,802.62 554.62 211.90 6.15 7.50 ZKY6-1 S7 524,064.82 3,782,801.25 584.90 13.20 0.29 0.79 ZKY6-1 S21 524,049.28 3,782,831.74 560.04 0.01 0.01 0.06 ZKY5-3 S7_2 524,745.83 3,783,467.98 504.79 12.83 1.38 0.19 ZKY5-2 S8 524,722.37 3,783,436.63 399.95 0.01 0.20 0.12

1.3．2投影变换（projection transformation）

创建一查询，命名为coordinator_of_projection，依据矿脉的产状和钻孔的测斜数据计算出组合样品的投影坐标及水平厚度。

设计视图（design view）

SQL 视图（SQL view）

SELECT averagecompositesample.hole_id, averagecompositesample.vein,

IIf([projection_parameter]![group]="14",([averagecompositesample]![y]-[projection_parameter]![y1])/0.906307787+([averagecompositesample]![x]-[projection_parameter]![x1])*0.422618262-([averagecompositesample]![y]-[projection_parameter]![y1])/2.14450692*0.422618262+[projection_parameter]![base_line],[Y]) AS X2, [averagecompositesample]!z AS Y2, Round([ag1],2) AS ag, Round([pb1],2) AS pb, Round([zn1],2) AS zn,

IIf([collar]![hole_type]="DRILL",Abs([averagecompositesample]![length]*Cos(0-[survey_0]![dip ]/57.3)*Sin([survey_0]![azi]/57.3-[vein]![trend_radians])+[averagecompositesample]![length]*Sin(0-[survey_0]![dip]/57.3)/Tan([vein]![dip_radians])),[averagecompositesample]![length]) AS width, averagecompositesample.length, collar.hole_type,

[averagecompositesample]![ag1]*0.92+[averagecompositesample]![pb1]*36.69+[averagecompositesa mple]![zn1]*28.84 AS ageq, Round([ag1]*[width],2) AS agm, Round([pb1]*[width],2) AS pbm, Round([zn1]*[width],2) AS znm,

([averagecompositesample]![ag1]*0.92+[averagecompositesample]![pb1]*36.69+[averagecomposites ample]![zn1]*28.84)*[width] AS ageqm

FROM (((averagecompositesample LEFT JOIN projection_parameter ON averagecompositesample.vein = projection_parameter.vein) LEFT JOIN vein ON averagecompositesample.vein = vein.vein) INNER JOIN survey_0 ON averagecompositesample.hole_id = survey_0.hole_id) INNER JOIN collar ON averagecompositesample.hole_id = collar.hole_id; 数据视图（Data view）

1.3.3 各条矿脉的地质样品表（samples for every vein） 为每一条矿脉创建一个生成表查询，生成各个矿脉的地质组合样品表