PolygonRoi trimFloatPolygon(PolygonRoi roi, double length) {
FloatPolygon poly = roi.getFloatPolygon();
float[] x = poly.xpoints;
float[] y = poly.ypoints;
int n = poly.npoints;
x = smooth(x, n);
y = smooth(y, n);
float[] curvature = getCurvature(x, y, n);
double threshold = rodbard(length);
// IJ.log("trim: "+length+" "+threshold);
double distance = Math.sqrt((x[1] - x[0]) * (x[1] - x[0]) + (y[1] - y[0]) * (y[1] - y[0]));
int i2 = 1;
double x1, y1, x2 = 0, y2 = 0;
for (int i = 1; i < n - 1; i++) {
x1 = x[i];
y1 = y[i];
x2 = x[i + 1];
y2 = y[i + 1];
distance += Math.sqrt((x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1)) + 1;
distance += curvature[i] * 2;
if (distance >= threshold) {
x[i2] = (float) x2;
y[i2] = (float) y2;
i2++;
distance = 0.0;
}
}
int type = roi.getType() == Roi.FREELINE ? Roi.POLYLINE : Roi.POLYGON;
if (type == Roi.POLYLINE && distance > 0.0) {
x[i2] = (float) x2;
y[i2] = (float) y2;
i2++;
}
PolygonRoi p = new PolygonRoi(x, y, i2, type);
if (roi.getStroke() != null) p.setStrokeWidth(roi.getStrokeWidth());
p.setStrokeColor(roi.getStrokeColor());
p.setDrawOffset(roi.getDrawOffset());
p.setName(roi.getName());
imp.setRoi(p);
return p;
}
