y23 = gaussian_filter1d(y22, sigma)
z23 = gaussian_filter1d(z22, sigma)
ax.plot(z13,y13,x13,color = 'b')
ax.plot(z23,y23,x23,color = 'b')
plt.show()
# This entire script is used mainly for debugging and visualizing the cross sections to make
# that everything is as it should be. Outputs multiple images as visual representation
if __name__ == "__main__":
parser = ap.ArgumentParser()
parser.add_argument('im')
parser.add_argument('crossSections')
parser.add_argument('centerPoints')
args = parser.parse_args()
centerPoints = EdgeFinder.getCenterPoints(args.centerPoints)
img = cv2.imread(args.im)
crossSections = []
with open(args.crossSections) as f:
for line in f:
crossSection = line.split()
xMid = float(crossSection[0])
yMid = float(crossSection[1])
dist = float(crossSection[2])
angle = float(crossSection[3])
height = float(crossSection[4])
crossSections.append((xMid,yMid,dist,angle,height))
displayCrossSections(crossSections,img,"interpolatedCrossSections1.txt")
plotTrack(crossSections,centerPoints)
