brmrtg

doc: http://jenda.hrach.eu/brm/ct.pdf

parallel_beam.py patch:

30a31
> import matplotlib.pyplot as pl
383c384
< 	if (f_type.lower () == "none"):
---
> 	"""if (f_type.lower () == "none"):
387c388
< 			return filtered_proj
---
> 			return filtered_proj"""
415c416,417
< 		raise ValueError ("rho_filter: Unknown window type: %s" % f_type)
---
> 		#raise ValueError ("rho_filter: Unknown window type: %s" % f_type)
> 		window = 1
419c421
< 	
---
> 	print(filt)
428,429c430,432
< 	for i in range (proj_fft.shape [1]):
< 		proj_fft [:, i] *= filt
---
> 	if (f_type.lower () != "none"):
> 		for i in range (proj_fft.shape [1]):
> 			proj_fft [:, i] *= filt
435,436c438,442
< 	filtered_proj.resize ((proj.shape [0], filtered_proj.shape [1]))
< 	
---
> 	print(proj.shape [0], filtered_proj.shape [1])
> 	#filtered_proj.resize ((proj.shape [0], filtered_proj.shape [1]))
> 	arr=filtered_proj[:proj.shape [0]][:proj.shape [0]]
> 	print(arr.shape[0], arr.shape[1])
> 	filtered_proj=arr