+28Hi,
I was just presented with a scenario where there are bunch of polygons, and some points, and requested to remove any polygon vertex that coincides with a point.
I initially thought about using the VertexCreator's Replace point at Index, but it turns out that the index can't be set to an attribute.
Is there a better alternative than chopping everything into points, filtering out the bad ones and then reconstructing the perimeter and rebuilding the area?
+17Hi @jdh, if you eventually need to remove the vertex, how will you go about it without breaking up the polygon?
+28Hi @jdh, if you eventually need to remove the vertex, how will you go about it without breaking up the polygon?
What I had hoped to do, was identify the index of the vertex to remove, use the VertexCreator Replace Point at Index to set it to NaN and then the GeometryValidator to remove the NaN.
+17What I had hoped to do, was identify the index of the vertex to remove, use the VertexCreator Replace Point at Index to set it to NaN and then the GeometryValidator to remove the NaN.
+28
Hi @jdh, I have no ideas except PythonCaller scripting.
import fme
import fmeobjects
class FeatureProcessor(object):
def __init__(self):
self.points = []
self.polygons = []
def input(self, feature):
geomType = feature.getGeometryType()
if geomType == fmeobjects.FME_GEOM_POINT:
self.points.append(feature)
elif geomType == fmeobjects.FME_GEOM_POLYGON:
self.polygons.append(feature)
def close(self):
points = set([p.getCoordinate(0) for p in self.points])
for polygon in self.polygons:
vertices = polygon.getAllCoordinates()
vertices = [p for p in vertices if p not in points]
if vertices[0] != vertices[-1]:
vertices.append(vertices[0])
boundary = fmeobjects.FMELine(vertices)
polygon.setGeometry(fmeobjects.FMEPolygon(boundary))
self.pyoutput(polygon)
+28Hi @jdh, I have no ideas except PythonCaller scripting.
import fme
import fmeobjects
class FeatureProcessor(object):
def __init__(self):
self.points = []
self.polygons = []
def input(self, feature):
geomType = feature.getGeometryType()
if geomType == fmeobjects.FME_GEOM_POINT:
self.points.append(feature)
elif geomType == fmeobjects.FME_GEOM_POLYGON:
self.polygons.append(feature)
def close(self):
points = set([p.getCoordinate(0) for p in self.points])
for polygon in self.polygons:
vertices = polygon.getAllCoordinates()
vertices = [p for p in vertices if p not in points]
if vertices[0] != vertices[-1]:
vertices.append(vertices[0])
boundary = fmeobjects.FMELine(vertices)
polygon.setGeometry(fmeobjects.FMEPolygon(boundary))
self.pyoutput(polygon)
You ca also use python factory in a python caller and be able to use the VertexCreatorFactory.
Please refer to the answer to this question to see how to do it.
Your factory string will look like this:
FACTORY_DEF * VertexCreatorFactory \
FACTORY_NAME VertexCreator \
INPUT FEATURE_TYPE Input_Port \
MODE REPLACE_AT_INDEX \
INDEX <index> \
CONTINUE_ON_ERROR YES \
XVAL "<x>" \
YVAL "<y>" \
ZVAL "<z>" \
OUTPUT OUTPUT FEATURE_TYPE VertexCreator_OUTPUT
You'll have to dynamically replace the <index> by the value of the attribute value holding your point index to replace and <x>, <y> and <z> by the coordinates values.
Regards,
Larry
+28You ca also use python factory in a python caller and be able to use the VertexCreatorFactory.
Please refer to the answer to this question to see how to do it.
Your factory string will look like this:
FACTORY_DEF * VertexCreatorFactory \
FACTORY_NAME VertexCreator \
INPUT FEATURE_TYPE Input_Port \
MODE REPLACE_AT_INDEX \
INDEX <index> \
CONTINUE_ON_ERROR YES \
XVAL "<x>" \
YVAL "<y>" \
ZVAL "<z>" \
OUTPUT OUTPUT FEATURE_TYPE VertexCreator_OUTPUT
You'll have to dynamically replace the <index> by the value of the attribute value holding your point index to replace and <x>, <y> and <z> by the coordinates values.
Regards,
Larry
That's a good workaround to the initial problem of the VertexCreator not accepting attributes. I was able to get it functioning with a little help from @DaveAtSafe. We'll have to run some tests to see if this solution is more or less efficient than @takashi's.
Enter your E-mail address. We'll send you an e-mail with instructions to reset your password.
