Skip to main content
  • EN

I have the following issue:

A bridge and the road below, standard 2D lines, both have nodes with identical coordinates - i.e. they share a "node location". When I want to output this dataset to OSM XML, the writer creates a connection in this point, even though they are not connected in realit". I've tried to separate them by Z-coordinate, but that had no effect. Therefore, my guess is that it creates a "node cache" with the coordinates as key, which means they will share a node because they share a coordinate.

With this in mind, my solution would be to slightly adjust the position of these internal nodes that overlap. Probably with just a minor change(cm? mm?), the nodes will not be cached to the same location. However, I've got no idea how to A) access all non-start or non-end nodes of a line and subsequently 😎 slightly and randomly alter their position.

Hi @fhilding, For A see the CoordinateExtractor, mode Specify Coordinate, Coordinate Index =-1 and for B you can use the RandomNumberGenerator to generate a value between 0.001 and 0.010 (if your data is in meters) to add to the last coordinate via the Offsetter transformer.

Hope this helps.

An alternative approach would be to use the CoincidentPointOffsetter transformer from the FME Hub. Since it works on points only you'll have to deconstruct your lines and keep track of both a line id and a point id so you can later on reconstruct the lines using the PointConnector.

offset-points.fmw

Hope this helps.

I went the Python route on this, inspired by redgeographics. The following snippet did the trick:

(Although it corrects ALL internal nodes, the offset is so small that it's acceptable for the data provided.)

    line = feature.getGeometry()
    ptTupleList = line.getPoints()
    ptTupleListNoEnds = ptTupleLists1:len(ptTupleList)-1]
    randomDifference = 0
    while -0.01 < randomDifference < 0.01:
        randomDifference = random.uniform(-0.10,0.10)
    correctedPoints = r]
    for point in ptTupleListNoEnds:
        point.setXYZ(point.getXYZ()b0]+randomDifference, point.getXYZ()L1], point.getXYZ()Â2])
        correctedPoints.append(point)
    for i in range(len(correctedPoints)):
        line.setPointAt(i+1, correctedPointsni].getXYZ())
    feature.setGeometry(line)
    #feature.setAttribute("corr",randomDifference)
    return feature

Reply


Community Stats

30,126
Posts
114,368
Replies
37,764
Members
Terms & ConditionsCookie settings
A product of
Safe Software respectfully acknowledges that we live, learn and work on the traditional and unceded territories of the Kwantlen, Katzie, and Semiahmoo First Nations.