活塞式柱状淤泥采样器:用mapgis结合excel作资源量估算图

用mapgis结合excel作资源量估算图
摘要：传统的资源量估算图的作法效率低而且容易出错，用mapgis结合excel进行资源量估算及制图，不仅在精度上很好的满足了要求，而且给地质人员减轻了负担。具有很好的实用性。
关键词：mapgis；资源储量核实；资源储量估算图；应用；
传统的资源储量估算图的作法效率低而且容易出错，借助mapgis及excel来完成资源量估算图，不失为一种新的选择。下面结合自己的经验，论述一下用mapgis结合excel作资源量估算图的过程。
一、图框的生成
以前手工绘图的时候，做一幅图一般是先打好网格，再逐网格来绘制，这样能够减少误差。现在有了电脑，一般是先用电脑生成需要的图框，再把图配准到图框中，再来处理。生成图框很重要。Mapgis提供了两大类图框的生成，一是经纬网图框，二是公里网图框。一般在资源储量核实中用得最多的是公里网图框。公里网图框的生成很简单，打开mapgis投影变换模块，打开键盘生成矩形图框，根据需要填写相应的参数即可生成。需要注意的是，坐标系一般填写国家坐标系，起始带号填写相应的带号，标注一般为公里值，因为一般图都不规则，矩形分幅方法选择任意公里矩形分幅，比例尺分母在图上都能找到，网起始值指的是公里网从哪儿开始，比如左下角坐标为2735232.33，35437248.11，1:5000的图，网起始值x填写为437.5，网起始值y填写为2735.5。对1:5000的图网间隔一般为0.5，而且将左下角平移为原点不要打钩，这样才方便后续工作（见图1）。

二、拐点投影及坐标转换
在作资源量估算图的过程中，把已知的坐标点投影到图上，并进行坐标变换经常要用到，比如矿权范围坐标，作为资源量估算来讲，矿权范围很重要，把拐点坐标很精确的投影到图上，以便确定面积，并根据需要进行转换，以确定经纬度位置等等。而mapgis软件就提供了这样的方便，打开mapgis，打开投影变换模块，无论用用户文件投影或单点投影变换，可以把点投影到相应比例尺的图上或进行坐标变换。需要注意的是，用户投影参数及结果投影参数一定要搞清楚，否则会出现错误的结果，如果仅仅是为了坐标变换，比如要结果为3度带或6度带的高斯坐标，结果投影参数的比例尺分母一定要设置为1，单位为米。而要经纬度坐标，单位则要定为或度、或分、或秒，虽然用DDDMMSS.SS也可以，但有些误差。在要求不是很严格的时候，也可以用。
当然，如果点很少的情况下，也可以直接在输入编辑模块中用定位点来投影（见图2），实际上在资源储量核实中的坐标，一般都是大比例尺的高斯平面直角坐标，在mapgis中的坐标表示公式为：

图2  1:5000的图定位点
横坐标x＝高斯坐标y（不要带号）/（比例尺分母/1000）
纵坐标y＝高斯坐标x/(比例尺分母/1000)
比如：高斯坐标为x＝2764650，y＝35400500，放到比例尺为1:5000的图上，图上坐标一般为：横坐标x＝80100，y＝552930。
这样做的前提条件是用mapgis生成的标准图框没有平移到原点。
三、图像配准
有了图框，投了拐点，在拷贝相应的图签，写上相应的图名，一幅的图的框架基本构成了，在这样一个框架的基础上，把各煤层扫描的储量图配准的这样一个统一的图框中，根据需要进行切割，图像处理模块发挥了重要作用。
打开mapgis图像分析或镶嵌配准模块，把扫描的tif图像转换为msi图像，进行配准，在这个过程中需要注意的是，配准所选的点一定要均匀分布，而且点越多越好，校正的时候可以选择用影像精校正，这样出来的图较准确。最后可以用标准图框生成的区来进行切割。以方便打印光栅。
四、矢量化
矢量化和其它的矢量化没有什么区别，需要注意的是，造区的时候，如果为煤矿，且有底板等高线，块段的边界一般以底板等高线和勘探线为界。另外每一块段最好造一个区，并且在区属性结构中增加一个块段号字段，填写相应的块段号作为属性，方便区属性的输出（见图3、图4）。

图3 编辑区属性结构

图4 输入区属性
五、量算面积及资源量估算
最后一步就是计算资源量，以前的老办法是用求积仪量算面积，再用计算器算出资源量，很麻烦，而且容易出错，现在通过mapgis的区属性，很容易读出面积，再调入excel中进行计算，不容易出错。在mapgis中，打开属性管理模块，输出相应的区属性（见图5），区属性输出时，一般输出为dbf的格式，方便调入excel中进行资源量计算。

图 5 属性输出
然后打开excel，读入区属性，因为区很多，怎样识别你所要的区呢，块段号字段发挥了重要作用。需要注意的是，区属性中的面积单位一般为平方毫米，要根据比例尺转换为平方米才能满足需要，即乘于常数即可，常数为比例尺分母除以1000后的平方值，然后再填写相应的参数进行计算即可。在计算的过程中，为了跟计算器算的结果一致，需要用到sum及round函数（见图6），sum函数用于求和，round函数的功能是返回某个数字按指定位数取整后的数字，语法是：
ROUND(number,num_digits)
Number    需要进行四舍五入的数字。
Num_digits    指定的位数，按此位数进行四舍五入。
最后再把excel算的结果表“复制”，在cad中“选择性粘贴”，修改完善后再转入mapgis放到资源估算图上。
以上过程，笔者已经在很多项目中实践过，笔者的体会是越大、块段数越多的项目，越能简化工作，通过mapgis结合excel来作资源量估算图，简化了烦琐的投点、求积仪量算面积，计算器算资源量等，给地质人员带来了很多方便，随着软件的发展和软件使用者的熟练，必将给地质人员带来越来越多的方便。

图6 资源量估算表
参  考  文  献
[1] 武汉中地信息工程有限公司.mapgis地理信息系统实用教程[M].2000
[2] 云南省国土资源厅.云南省固体矿产资源储量核实报告编写暂行规定.2004