Hess = np.zeros([nParam * nLen, nParam * nLen, nKern], 'f')
print 'Starting to read data.'
# Read in data
for k in range(0, nKern - 2):
for n in range(0, nPos):
for m in range(0, mPos):
for i in range(0, nParam):
for j in range(0, nPerts):
# Read data
# PARAM AND PERT SWITCHED!!
toRead = '{0}-{1}-d{2}-{3:03d}.bin'.format(
kernel[k], param[j], perts[i], yPos[n])
temp = abp.read_data_2d(folder + toRead, nx, ny)
if nPos > mPos:
Hess[j * nLen:(j + 1) * nLen, i * nLen + n,
k] = temp[yPos, xPos]
else:
Hess[j * nLen:(j + 1) * nLen, i * nLen + m,
k] = temp[yPos, xPos]
## Remove diagonals
# if k != 3 and k != 4:
# for i in range(0,nParam*nLen):
# Hess[i,i,k] = 0.0
# if (i+nLen < nParam*nLen):
# Hess[i+nLen,i,k] = 0.0
# Hess[i,i+nLen,k] = 0.0
else:
return (x - x.min())
prefix = 'Gradient '
model = 'Grad5km/'
project = '/home/vegahag/msim-ext/SEG2018/'+model
save = 'SEG2018/'+model
folder = project
modelFolder = project
if not os.path.exists(save):
os.makedirs(save)
rh = abp.read_data_2d(modelFolder+"rh-true.bin",nx,ny);
vp = abp.read_data_2d(modelFolder+"vp-true.bin",nx,ny);
vs = abp.read_data_2d(modelFolder+"vs-true.bin",nx,ny);
model = np.zeros([ny,nx,3],'f')
model[:,:,0] = rh
model[:,:,1] = vp
model[:,:,2] = vs
param = [
#'rho-bg','vp-bg','vs-bg',
'rh-true','vp-true','vs-true',
#'rho-diff','vp-diff'
]
title = [
dt = 0.001
pclip = .1
tpow = 1.8
model = 'Gull0km/'
project = '/home/vegahag/msim-ext/SEG2018/' + model
save = 'SEG2018/' + model
folder = project
modelFolder = project
if not os.path.exists(save):
os.makedirs(save)
tru = abp.read_data_2d(modelFolder + "data-true.bin", nt, nx)
mod = abp.read_data_2d(modelFolder + "data-mod.bin", nt, nx)
res = abp.read_data_2d(modelFolder + "data-res.bin", nt, nx)
rec = np.zeros([nt, nx, 3], 'f')
rec[:, :, 0] = tru.transpose()
rec[:, :, 1] = mod.transpose()
rec[:, :, 2] = res.transpose()
param = [
'rec-tru',
'rec-mod',
'rec-res',
]
title = [
project = '/home/vegahag/msim-ext/SEG2018/'+model
save = 'SEG2018/'+model
xStart = nx/2-1;
xStop = xStart+1;
yStart = #ny/2-step-1;
yStop = #ny/2+step;
folder = project
modelFolder = project
if not os.path.exists(save):
os.makedirs(save)
rh = abp.read_data_2d(modelFolder+"model-rh-true.bin",nx,ny);
vp = abp.read_data_2d(modelFolder+"model-vp-true.bin",nx,ny);
vs = abp.read_data_2d(modelFolder+"model-vs-true.bin",nx,ny);
rhNorm = normalize(rh)
vpNorm = normalize(vp)
vsNorm = normalize(vs)
model = np.zeros([nx,ny,3],'f')
model[:,:,0] = rh
model[:,:,1] = vp
model[:,:,2] = vs
kernel = ['H1','H2','H3','H12','H123']
param = ['rh','vp','vs']
perts = ['rh','vp','vs']
prefix = 'Hessian '
model = 'Gull5km/'
project = '/home/vegahag/msim-ext/SEG2018/maur/' + model
save = 'SEG2018/' + model
folder = project
modelFolder = project
if not os.path.exists(save):
os.makedirs(save)
H1 = abp.read_data_2d(
modelFolder + "H1-" + par + "-d" + pert + "-" + '{0:03.0f}'.format(pos) +
".bin", nx, ny)
H2 = abp.read_data_2d(
modelFolder + "H2-" + par + "-d" + pert + "-" + '{0:03.0f}'.format(pos) +
".bin", nx, ny)
H3 = abp.read_data_2d(
modelFolder + "H3-" + par + "-d" + pert + "-" + '{0:03.0f}'.format(pos) +
".bin", nx, ny)
vp = abp.read_data_2d(modelFolder + "vp-true.bin", nx, ny)
Hess = np.zeros([ny, nx, 4], 'f')
Hess[:, :, 0] = H1
Hess[:, :, 1] = H2
Hess[:, :, 2] = H3
Hess[:, :, 3] = H1 + H2 + H3