摘要：傳統的資源量估算圖的做法效率低且容易出錯，結合實際工作，論述了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來做資源量估算圖，給地質人員帶來極大方便，隨著軟件的更新和使用者熟練程度的提高，必將帶來更大的方便。
