Added more plots to psf-4-modified.
authormkass <andy.kass@geo.au.dk>
Sun, 3 Sep 2017 07:27:19 +0000 (09:27 +0200)
committermkass <andy.kass@geo.au.dk>
Sun, 3 Sep 2017 07:27:19 +0000 (09:27 +0200)
EAGE_Abstract/Lemma/psf/psf-4-modified.py

index 8b4c5d5..a20bf3c 100644 (file)
@@ -71,12 +71,37 @@ for a in ax[1:]:
     
 ax[0].set_ylabel('Depth [m]')
 cb = plt.colorbar(im[0], axcb )
-cb.set_label(r'$\left| \mathcal{V}_0 \right|$~(nV)')
+#cb.set_label(r'$\left| \mathcal{V}_0 \right|$~(nV)')
 
-#plt.savefig('RES_SVD.pdf',bbox_inches='tight')
-#plt.show()
+plt.savefig('RES_SVD.pdf',bbox_inches='tight')
+plt.show()
 plt.close()
 
+#Plot combined SVD
+fig = plt.figure( figsize=(pc2in(80),pc2in(40)))
+ax1 = fig.add_axes( [0.091, .182, .3182, .6364] )
+ax2 = fig.add_axes( [0.455, .182, .3182, .6364], sharey=ax1,
+    sharex=ax1)
+ax2cb = fig.add_axes( [0.818, .182, .0455, .6364] )
+im = ax1.pcolormesh(X,Y,PS[4], cmap='viridis')#, vmin=0, vmax=3100)
+im2 = ax2.pcolormesh(X,Y,PS[0],cmap='viridis')#, vmin=0,vmax=3100)
+ax1.set_ylim( ifaces[-1],ifaces[0] )
+ax1.set_xlim( ifaces[0], ifaces[-1] )
+ax2.set_ylim( ifaces[-1],ifaces[0] )
+ax2.set_xlim( ifaces[0], ifaces[-1] )
+
+im.set_edgecolor('face')
+im2.set_edgecolor('face')
+ax1.set_title('Combined')
+ax2.set_title('Coincident')
+plt.setp(ax2.get_yticklabels(), visible=False)
+cb = plt.colorbar(im,ax2cb)
+plt.savefig('RES_COMB_SVD',bbox_inches='tight')
+ax1.set_xlabel('Depth [m]')
+ax1.set_ylabel('Depth [m]')
+ax2.set_xlabel('Depth [m]')
+plt.show()
+plt.close()
 
 #Plot explicit resolution 
 print(R[0].shape)