檢測，得到多個相應尺度的位場邊 緣；
[0021] 采用形態(tài)學骨骼算法將計算得到的各尺度的位場邊緣分別細化為單像素寬度，得 到多個相應尺度的構造格架圖。
[0022] 優(yōu)選地，該方法進一步包括將計算得到的所述多個相應尺度的構造格架圖疊置生 成綜合構造格架圖。
[0023] 優(yōu)選地，該方法進一步包括將各尺度構造格架圖上的每一邊緣點處的梯度的模作 為該尺度構造格架圖中該邊緣點處的強度值，得到多個相應尺度的構造強度格架圖。
[0024] 優(yōu)選地，該方法進一步包括將所述多個相應尺度的構造強度格架圖疊置生成綜合 構造強度格架圖。
[0025] 優(yōu)選地，所述位場數據為重力位場數據或磁法位場數據，所述預處理進一步包括 對磁法數據進行化極計算得到化極磁異?；蜻M行偽重力計算得到偽重力異常；或對重力數 據進行預處理得到布格重力異常。
[0026] 優(yōu)選地，針對每一尺度的位場數據進行多方向邊緣檢測，包括以下步驟：
[0027] 設尺度參數s = z/Z(l，且z>Z(l，Ztl代表測量高度，z代表向上延拓后的高度，定義高 度為零的位置（X，y)處的重力異?；虼女惓閒〇(x，y)，
[0028] 尺度s的平滑函數定義為：
[0029]
【主權項】
1. 一種位場構造格架自動提取方法，包括以下步驟： 對來自待研究區(qū)域的位場數據進行預處理； 將經預處理的位場數據向上延拓多個預定高度得到多個相應尺度的位場數據； 分別針對每一尺度的位場數據進行多方向邊緣檢測，得到多個相應尺度的位場邊緣； 采用形態(tài)學骨骼算法將計算得到的各尺度的位場邊緣分別細化為單像素寬度，得到多 個相應尺度的構造格架圖。
2. 如權利要求1所述的位場構造格架自動提取方法，其特征在于，該方法進一步包括 將計算得到的所述多個相應尺度的構造格架圖疊置生成綜合構造格架圖。
3. 如權利要求1所述的位場構造格架自動提取方法，其特征在于，該方法進一步包括： 將各尺度構造格架圖上的每一邊緣點處的梯度的模作為該尺度構造格架圖中該邊緣 點處的強度值，得到多個相應尺度的構造強度格架圖。
4. 如權利要求3所述的位場構造格架自動提取方法，其特征在于，該方法進一步包括 將所述多個相應尺度的構造強度格架圖疊置生成綜合構造強度格架圖。
5. 如權利要求1所述的位場構造格架自動提取方法，其特征在于，所述位場數據為重 力位場數據或磁法位場數據，所述預處理進一步包括 對重力數據進行預處理得到布格重力異常；或 對磁法數據進行化極計算得到化極磁異?；蜻M行偽重力計算得到偽重力異常。
6. 如權利要求1所述的位場構造格架自動提取方法，其特征在于， 針對每一尺度的位場數據進行多方向邊緣檢測，包括以下步驟： 設尺度s = z/Z(l，且Z>Z(I，Ztl代表測量高度，z代表向上延拓后的高度，定義高度為零的 位置（X，y)處的重力異常或磁異常為fQ(x，y)， 尺度s的平滑函數定義為：
其中L
k(x, y, z)為格林函數， 在方向α的小波函數定義為：
其中，D表不一階導數； 在尺度s和位置（x，y)的情況下，重力異常或磁異常L(Xj)在方向α的小波變換定 義為：
其中，*表示卷積運算， 由位場向上延拓公式已知.V) = ./n fz(x，y)為f〇(x，y)從高度零上延到高度Z = 82(|處的重力異?；虼女惓?，通過將在測 量高度Ztl測量得到的重力異?；虼女惓Ｏ蛏涎油馗叨萙-Ztl獲得， 因此，[f。] (X，y, s) = sD°fz(x，y) =(z/z〇)Dafz(x, y)； 進一步，在尺度s和位置（x，y)，重力異?；虼女惓 (x，y)在方向《+f的小波變換 定義為：
貝ij，LUy)的二維方向小波變換用梯度表示為：
w[./〇](s,a) = {z! Zn)V/； (λ-, v,a) 其中V為二維梯度， 對于高度z，定義梯度的模為：
該梯度相應的沿水平方向的輻角為：
邊緣點為模M[fz] (X，y, a )沿福角方向Afz(x, y, a )有局部極大值的點， 針對每一方向a，梯度的模的局部極大值點沿梯度的垂直方向連接得到的曲線構成邊 緣， 針對同一高度，以多個不同的方向a計算邊緣，對計算得到的各邊緣求并集得到相應 尺度的位場邊緣。
7. 如權利要求6所述的位場構造格架自動提取方法，其特征在于，針對同一高度以 多個不同的方向a計算邊緣的步驟進一步包括，各方向a取值為其中k = 0, 1，2,…，（2η-1)，η為大于或等于2的整數，以完整覆蓋二維平面。
8. 如權利要求6所述的位場構造格架自動提取方法，其特征在于，位場向上延拓多個 預定高度后提取的邊緣對應于不同深度的構造，對得到的各尺度的構造格架圖進行疊加得 到反映不同深度信息的綜合構造格架圖。
9. 如權利要求2或8所述的位場構造格架自動提取方法，其特征在于，以漸變顏色表征 高度的高低，來形成所述綜合構造格架圖。
10. 如權利要求4所述的位場構造格架自動提取方法，其特征在于，以漸變顏色表征強 度值的大小，來形成所述綜合構造強度格架圖。
【專利摘要】本發(fā)明提供一種位場構造格架自動提取方法。該方法包括以下步驟：對來自待研究區(qū)域的重力位場數據或磁法位場數據進行預處理；對經預處理的重力位場數據或磁法位場數據進行多個尺度下多方向的邊緣檢測，分別得到各尺度的邊緣；采用形態(tài)學骨骼算法將計算得到的各尺度的邊緣細化為單像素寬度，并且每個點具有深度和強度屬性，得到重力或磁綜合構造格架圖和綜合構造強度格架圖。根據本發(fā)明的方法，可以實現(xiàn)對控制礦床形成的地質構造的識別和定性解釋，根據研究區(qū)先驗知識確定潛在的礦床類型和控制礦床形成的構造屬性，對不同類型構造格架進行篩選，從而實現(xiàn)金屬礦床的靶區(qū)定位。
【IPC分類】G06T17-00, G06T7-00
【公開號】CN104658037
【申請?zhí)枴緾N201510096981
【發(fā)明人】曹殿華
【申請人】中國地質科學院礦產資源研究所
【公開日】2015年5月27日
【申請日】2015年3月4日