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)