def reportClustersAsTable(clusters,
allPoints,
XColumn='X',
YColumn='Y',
ZColumn='Z',
NRColumn='NR'):
'''
Report the clustered and unclustered points in the tables 'clusters' and 'unclustered'.
'''
rt = ResultsTable()
counter = 1
clusterCounter = 1
clusteredPoints = []
for c in clusters:
for dp in c.getPoints():
rt.incrementCounter()
p = dp.getPoint()
rt.addValue(NRColumn, counter)
rt.addValue(XColumn, p[0])
rt.addValue(YColumn, p[1])
rt.addValue(ZColumn, p[2])
rt.addValue("C", clusterCounter)
counter = counter + 1
clusteredPoints.append([p[0], p[1], p[2]])
clusterCounter = clusterCounter + 1
rt.show("clusters")
win = WindowManager.getWindow("Results")
rt = win.getResultsTable()
X, Y, Z = getColumns(XColumn, YColumn, ZColumn)
if not rt.columnExists(NRColumn):
for i in range(0, len(X)):
rt.setValue(NRColumn, i, i + 1)
rt.updateResults()
NR = getColumn(NRColumn)
unclusteredPoints = [
[point.getPoint()[0],
point.getPoint()[1],
point.getPoint()[2]] for point in allPoints
if [point.getPoint()[0],
point.getPoint()[1],
point.getPoint()[2]] not in clusteredPoints
]
counter = 1
rt = ResultsTable()
for p in unclusteredPoints:
rt.incrementCounter()
rt.addValue(NRColumn, counter)
rt.addValue(XColumn, p[0])
rt.addValue(YColumn, p[1])
rt.addValue(ZColumn, p[2])
counter = counter + 1
rt.show("unclustered")
WindowManager.setWindow(win)
