Loop through list with python and get geometry

The listBuilder does not perserver geometry, but the Aggregator does, and can build the corresponding list.

Hi,

Although I can see the aggregate in data inspector, I still have the same result in the python caller. The list is None with this line of code

my_list = feature.getAttribute("_list{}")

I think I miss something from your explanation

regards,

Sylvain

@david_r

@jdh

Hi, thank for your reply,

Although I can now see all point aggregated by partition in Inspector, I still have the same result with this line

my_list = feature.getAttribute("_list{}") --> NoneType

I can't retrieve the geometry

Maybe I missed something in your explanation

Regards,

Sylvain

You won't find the geometries inside the attribute list (it only contains attributes), you'll have to iterate over the aggregate parts and extract the geometry from each part, probably something like this:

from fmeobjects import *

if feature.getGeometryType() == FME_GEOM_AGGREGATE:
    aggregate_parts = feature.splitAggregate(True)
elif feature.getGeometryType() == FME_GEOM_DONUT:
    aggregate_parts = feature.getDonutParts()
else:
    aggregate_parts = [feature,]

for n, part in enumerate(aggregate_parts):
    part_geom = part.getGeometry()
    part_objectid = feature.getAttribute('_list{%s}.OBJECTID' % n)
    # Do something exciting with "part_geom" and "part_objectid" here...

You'll find more information about the FMEGeometry object in the API documentation:

def input(self, feature):

    # flush the list
    self.fList = feature.splitAggregate(False)

    self.logger.logMessageString("There is {} point : ".format(len(self.fList)), fmeobjects.FME_ERROR)

    self.logger.logMessageString("Starting DBScan: ", fmeobjects.FME_ERROR)
    self.DBScan(self.fList, self.eps, self.minPts)


def DBScan(self, D, eps, minPts):

    corePoints = [f for f in D if self.getNbOfNeighbours(D,f,eps) >= minPts ]

    self.logger.logMessageString("Outputing Core Point: ", fmeobjects.FME_ERROR)

    for f in corePoints:
        self.logger.logMessageString("Outputing point {}: ".format(str(f)), fmeobjects.FME_ERROR)
        self.pyoutput(f)

def distance3D(self, featureA, featureB):

    coorA = featureA.getAllCoordinates()
    coorB = featureB.getAllCoordinates()

    # [node number of current feature][coord component x, y, z]

    xa = coorA[0][0]
    ya = coorA[0][1]
    za = coorA[0][2]

    xb = coorB[0][0]
    yb = coorB[0][1]
    zb = coorB[0][2]

    d = math.sqrt((xa - xb) ** 2 + (ya - yb) ** 2 + (za - zb) ** 2)
    return d