コード例 #1
0
ファイル: inverabel.py プロジェクト: huletlab/apparatus3-py
def inverAbel(path,show=1):
	try:
		data = np.loadtxt(path)
	except:
		print "File does not exist"
		return
		
	x =  data[:,0]
	y =  data[:,1]
	invx,invy = abel.inverAbel(x,y)
	x2 = np.append(x[:0:-1],x)
	y2 = np.append(y[:0:-1],y)
	invx2,invy2 = abel_scikit.inverAbel(x2,y2)
	
	fig = plt.figure()
	ax = fig.add_subplot(111)
	ax2 = ax.twinx()
	lns1 = ax.plot(x,y,label="Original Azimuthal Data")
	ax.set_xlabel("Pixels")
	ax.set_ylabel("Azimuthal Ave")
	lns2 = ax2.plot(invx,invy,"r--",label="Inv Abel Trans by direct Integral")
	lns3 = ax2.plot(invx2[len(invx2)/2+1:],invy2[len(invx2)/2+1:],"g--",label="Inv Abel Trans by Radon Trans")
	ax2.set_ylabel("Azimuthal Ave Inverse Abel Transfrom")
	lns = lns1+lns2+lns3
	labs = [l.get_label() for l in lns]
	ax.legend(lns, labs, loc=1)
	plt.savefig(path.replace(".AZASCII","_invAbel.png"))
	if (show):
		plt.show()
コード例 #2
0
ファイル: inverabel.py プロジェクト: huletlab/apparatus3-py
def inverAbelDat(x,y):
        # First one uses the direct integral
	invx,invy = abel.inverAbel(x,y)
        invx = np.array(invx)
        invy = np.array(invy)

        # Second one uses the radon transform
	x2 = np.append(x[:0:-1],x)
	y2 = np.append(y[:0:-1],y)
	invx2,invy2 = abel_scikit.inverAbel(x2,y2)
        skipn = len(invx2)/2+1
        invx2 = np.array( invx2[skipn:] )
        invy2 = np.array( invy2[skipn:] )

        return (invx, invy),(invx2,invy2)