def uti_plot2d1d(self, data, x_range, y_range, xc, yc, labels, _graphs_joined=True):
x0 = x_range[0]; x1 = x_range[1]; nx = x_range[2]
#y0 = x_range[0]; y1 = y_range[1]; ny = y_range[2]
y0 = y_range[0]; y1 = y_range[1]; ny = y_range[2] #OC090714
label2D = labels[0]; label1H = labels[1]; label1V = labels[2]
fig = None
if _graphs_joined:
#fig = _pl.figure(figsize=(12,5))
fig = self._pl.figure(figsize=(15,5.3))
self._plot_2D(data, x_range, y_range, label2D, fig, 131) #showing graphs in one panel
else: self.uti_plot2d(data, x_range, y_range, label2D)
xStep = 0
if nx > 1: xStep = (x1 - x0)/(nx - 1)
yStep = 0
if ny > 1: yStep = (y1 - y0)/(ny - 1)
inperpOrd = 1 #interpolation order to use (1 to 3)
#_plot_1D(data[iy,:],range_x,label1H,fig,132)
arCutX = array('d', [0]*nx)
xx = x0
for ix in range(nx):
arCutX[ix] = uti_math.interp_2d(xx, yc, x0, xStep, nx, y0, yStep, ny, data, inperpOrd, 1, 0)
xx += xStep
if _graphs_joined: self._plot_1D(arCutX, x_range, label1H, fig, 132) #OC150814
else: self.uti_plot1d(arCutX, x_range, label1H)
#_plot_1D(data[:,ix],range_y,label1V,fig,133)
arCutY = array('d', [0]*ny)
yy = y0
for iy in range(ny):
#arCutY[iy] = _interp_2d(xc, yy, x0, xStep, nx, y0, yStep, ny, data, inperpOrd, 1, 0)
arCutY[iy] = uti_math.interp_2d(xc, yy, x0, xStep, nx, y0, yStep, ny, data, inperpOrd, 1, 0)
yy += yStep
if _graphs_joined: self._plot_1D(arCutY, y_range, label1V, fig, 133)
else: self.uti_plot1d(arCutY, y_range, label1V)
if _graphs_joined: self._pl.tight_layout() #OC081115
return self._maybe_savefig(fig)
def uti_plot2d1d(self, data, x_range, y_range, xc, yc, labels, _graphs_joined=True):
x0 = x_range[0]; x1 = x_range[1]; nx = x_range[2]
#y0 = x_range[0]; y1 = y_range[1]; ny = y_range[2]
y0 = y_range[0]; y1 = y_range[1]; ny = y_range[2] #OC090714
label2D = labels[0]; label1H = labels[1]; label1V = labels[2]
fig = None
if _graphs_joined:
#fig = _pl.figure(figsize=(12,5))
fig = self._pl.figure(figsize=(15,5.3))
self._plot_2D(data, x_range, y_range, label2D, fig, 131) #showing graphs in one panel
else: self.uti_plot2d(data, x_range, y_range, label2D)
xStep = 0
if nx > 1: xStep = (x1 - x0)/(nx - 1)
yStep = 0
if ny > 1: yStep = (y1 - y0)/(ny - 1)
inperpOrd = 1 #interpolation order to use (1 to 3)
#_plot_1D(data[iy,:],range_x,label1H,fig,132)
arCutX = array('d', [0]*nx)
xx = x0
for ix in range(nx):
arCutX[ix] = uti_math.interp_2d(xx, yc, x0, xStep, nx, y0, yStep, ny, data, inperpOrd, 1, 0)
xx += xStep
if _graphs_joined: self._plot_1D(arCutX, x_range, label1H, fig, 132) #OC150814
else: self.uti_plot1d(arCutX, x_range, label1H)
#_plot_1D(data[:,ix],range_y,label1V,fig,133)
arCutY = array('d', [0]*ny)
yy = y0
for iy in range(ny):
#arCutY[iy] = _interp_2d(xc, yy, x0, xStep, nx, y0, yStep, ny, data, inperpOrd, 1, 0)
arCutY[iy] = uti_math.interp_2d(xc, yy, x0, xStep, nx, y0, yStep, ny, data, inperpOrd, 1, 0)
yy += yStep
if _graphs_joined: self._plot_1D(arCutY, y_range, label1V, fig, 133)
else: self.uti_plot1d(arCutY, y_range, label1V)
if _graphs_joined: self._pl.tight_layout() #OC081115
return self._maybe_savefig(fig)