Alpha shape

alpha.py

def get_alpha_shape(infile,radius):
    command = "%s -A -aa %s -r -m10000000 -oN -oFpoints < %s" % (hull_path, str(radius), infile)
    print >> sys.stderr, "Running command: %s" % command
    retcode = subprocess.call(command, shell=True)
    results_file = open("points-alf")
    results_file.next() 
    results_indices = [[int(i) for i in line.rstrip().split()] for line in results_file]
    results_file.close()
    return results_indices