Python FIMBULの例

コード例 #1

ファイル: distance_stack_fimbul.py プロジェクト: kmunve/processgpr

'''
Created on 16. nov. 2010

@author: Karsten Mueller
'''

import os

import matplotlib.pyplot as plt

from processgpr.ext.fimbul import FIMBUL


wdir = r'C:\hubradata\Fimbul091230'
os.chdir(wdir)
file = 'data1.h5'
pro = FIMBUL()
pro.readH5(file)
print "Total Distance: ", pro.coords.total_distance, pro.traces
pro.distanceStack(10.0)
print "Total Distance: ", pro.coords.total_distance, pro.traces
pro.saveH5(os.path.splitext(file)[0] + "_ds10")
plt.scatter(pro.coords.x, pro.coords.y, c=pro.coords.elev, edgecolors='none')
plt.show()

コード例 #2

ファイル: convert_fimbul.py プロジェクト: kmunve/processgpr

'''
Created on 20.10.2010

@author: Karsten Mueller
'''
import os

import matplotlib.pyplot as plt

from processgpr.ext.fimbul import FIMBUL


wdir = r'C:\hubradata\Fimbul091230'
file = 'data1'
pro = FIMBUL()

pro.convertHDF(os.path.join(wdir, file), wdir,
               6000, combine=True, stack=3, filerange="all-all",
               envelope=True, gpsfile=os.path.join(wdir, file + '.cor'),
               projection='PS71')

print "Total Distance: ", pro.coords.total_distance, pro.traces
plt.scatter(pro.coords.x, pro.coords.y, c=pro.coords.elev, edgecolors='none')
plt.show()