Mapgis结合excel做资源量估算图

　　关键词:资源量估算、投影变换、区属性、矢量化。

　　Abstract: The traditional resource estimation graph practice efficiency is low and prone to error, combine with the actual work discusses mapgis frame to generate the inflection point projection, area calculation and excel to calculate the amount of resources the process, to improve the accuracy and resource estimation the efficiency of a certain reference.

　　Key words: resource estimation; projection transformation; district property; vectorization

　　中图分类号：U41 文献标识码：A 文章编号：

　　传统的资源量估算图的做法效率低且容易出错，借助Mapgis及excel来完成资源量估算图，不失为一种新的选择，下面结合自己的经验，论述一下作图过程。

　　1、图框自动生成

　　以前手工绘图的时候，一般是先打好网格，再逐格网来绘制，这样能减少误差，现在有了电脑，一般是先用电脑生成需要的图框，再把图配准到图框中，Mapgis提供了两大类图框的生成，一是经纬网图框，二是公里网图框，在资源量估算中用的最多的是公里网图框，公里网图框生成很简单，利用Mapgis“投影变换”模块，键盘生成矩形图框，根据需要填写相应参数即可生成。需要注意的是，坐标系一般填写国家坐标系，起始带号填写相应的带号，标注为公里值，Mapgis坐标系与大地坐标系的x、y互换，因为一般图都不规则，分幅方法选择任意公里矩形分幅，网起始值是指公里网从哪里开始，比如左下角坐标为4545900,40551100，比例尺1:1000的图，网起始值x填写为551.1，y填写为4545.9，网间隔为0.1。而且“将左下角平移为原点”，“底边旋转为水平”不要打勾。以方便以后作图。

　　2、拐点投影

　　在做资源量估算图的过程中，把已知的坐标点投影到图上，并进行坐标变换经常要用到，比如矿权范围坐标，作为资源量估算图来讲，要把矿权范围很精确地投影到图上，以便确定面积，而Mapgis软件就提供了这样的方便，打开Mapgis实用服务—投影变换模块，利用用户文件投影转换，将矿权坐标和其他带有大地坐标的点位，投影到相应比例尺的图上，需要注意的是，用户投影参数和结果投影参数的设置要正确，否则，会出现错误。

　　用户投影参数：坐标系类型为大地坐标系，椭球参数为一般为“北京54/克拉索夫斯基【1940椭球】”或“西安80/1975年I.U.G.G推荐椭球”，投影类型为高斯——克吕格【横切椭圆柱等角】投影，比例尺分母设为1，坐标单位设为米，其他参数按实际填写。

　　结果投影参数：坐标系类型为投影平面直角，椭球参数为一般为“北京54/克拉索夫斯基【1940椭球】”或“西安80/1975年I.U.G.G推荐椭球”，投影类型为高斯——克吕格【横切椭圆柱等角】投影，比例尺分母设为所作图的比例尺分母，单位设为豪米。其他参数按实际填写。

　　用户投影参数和结果投影参数设置后，设置其他选项，即可将坐标点投影到图框内。

　　如果要投影的点较少，也可以直接在输入编辑中用定位点来投影。作图时，坐标一般都是高斯坐标，在Mapgis中坐标的表示公式为：

　　横坐标x=高斯坐标y(不要带号)/(比例尺分母/1000)

　　纵坐标y=高斯坐标x/(比例尺分母/1000)

　　比如：高斯坐标为X=4545900、Y=40551100的点放到比例尺为1:1000的图上，图上坐标为：横坐标X=551100、 纵坐标Y=4545900，

　　这样做的前提条件是用Mapgis生成的标准图框没有平移到原点。

　　3、图像配准

　　有了图框，投了拐点，再做一些修饰，一幅图的框架基本构成了，在此基础上把扫描的储量图配准到图框上，根据需要进行切割，图形分析模块发挥了重要作用。

　　打开Mapgis“图像分析”模块，把扫描的图像转换为msi图像，进行配准，在这个过程中，需要注意的是：配准所选的点一定要均匀分布，而且点越多越好，校正时可以选择影像精校正，这样出来的图较准确。

　　4、矢量化

　　矢量化要时注意，造区时如果为煤矿，且有底板等高线，块段边界一般以底板等高线和勘探线为界，每个块段造一个区，并且在区属性结构中增加一个块段号字段，填写相应的块段号作为属性值。方便属性输出。

　　如果等高线较多，矢量化量大，可以使用交互矢量化功能，可以提高矢量化速度。步骤如下：

　　4.1首先使用Photoshop软件将图像进行处理，

　　(1)用Photoshop打开扫描的图像， 转换为灰度图，

　　(2)用“图像—调整—色阶”将线条调正清晰，

　　(3)用“图像—模式—位图”将图像转换为位图，

　　(4)用“图像—调整—反相”将位图(二制图)反相，存为tiff格式，

　　4.2使用MAPGIS输入编辑模块，装入反相二制图，进行等高线矢量化(不要进行图像转换和配准)，矢量化参数设置为：

　　节点搜索半径=0.000000000001

　　Buffer分析半径=1

　　插密光滑半径=0.005

　　坐标点间最小距离=0.000000000001

　　抽稀因子=0.03，用“线编辑—光滑线”光滑等高线，光滑类型选择“三次Bizer插值”，插密距离选择“0.005，然后用“线编辑—抽稀线” 抽稀等高线，抽稀因子设置为0.03，接下来，再矢量化其他图上要素。

　　4.3最后使用MAPGIS误差校正模块，将矢量化的图校正到标准图框上，步骤如下：

　　(1)使用MAPGIS误差校正模块，将要校正的文件装入工作区，

　　(2)设置采集控制点参数，数据采集类型为实际值，采集搜索范围为5，采集选择项勾选搜索T型相交线、显示搜索信息、显示控制点号、设置询问操作。

　　(3)选择采集文件，也就是选择要采集实际控制点的文件，一般为图框线文件，至少采集4个控制点，

　　(4)采集理论控制点，数据采集类型为理论值，其他参数不变，选择采集理论控制点的文件，一般为图框线文件，采集的4个控制点与实际控制点要一一对应。

　　4.4文件校正，将要校正的文件全部选中，注意不要选择采集理论控制点的文件。至此，交互矢量化完成。

　　五、量算面积及资源量估算

　　通过投影变换模块，将属性输出为文本文件，然后打开excel读入区属性，因为区很多，如何识别所要的区，块段号字段发挥了作用，需要注意的是，区属性中的面积单位一般为平方毫米，要根据比例尺转换为平方米才能满足需要，即乘以一个常数即可，常数为图的比例尺分母除以1000后的平方值。然后在填写相应的参数进行计算。在计算过程中，为了跟计算器的结果一致，需要用到sun及round函数。sun函数用于求和，round函数的功能是返回某个数字按指定位数取整后的数字，

　　ROUND 函数可将某个数字四舍五入为指定的位数。例如，如果单元格 A1 含有 23.7825 并且希望将该数字四舍五入为小数点后两位，则可以使用以下公式：

　　=ROUND(A1, 2)

　　此函数的结果为 23.78。

　　语法是：

　　Round(number,num_digits) number需要进行四舍五入的数字，

　　num_digits指定的位数，按此位数对number指定的参数进行四舍五入。

　　最后，把excel计算结果转换为Mapgis表文件(*.wb)。将Mapgis表文件(*.wb)与资源量估算图的区属性挂接。这样可以方便地在Mapgis中查看每个块段的储量信息。

　　以上过程，块段数越多的项目，越能简化工作，通过Mapgis结合excel来做资源量估算图，给地质人员带来极大方便，随着软件的更新和使用者熟练程度的提高，必将带来更大的方便。