xy(z) table to point feature, varying coordinate-systems

Maps-Berlin

Hello everybody,

I was wondering if it is possible to create a point feature from an xy table which contains coordinates in different coordinate systems if the table has a column containing the epsg-code for each row. Same goes for z-value and a column containing information about the vertical reference. 

I know i can split the table by coordinate system, create point features for each and then merge but I was wondering if there is a possibility to to this in a single operation, basically choose the coordinate system by attribute.

Thanks is advance

 

JohannesLindner

Not out of the box, but it's possible with a little Python.

This is only for x and y, because I have no clue about vertical reference.

 

def xy_to_point_enhanced(in_table, x_field, y_field, epsg_field, out_fc, out_epsg):
    """Converts tabular coordinates into a point feature class, analogous to XY Table To Point. coordinates can be from different coordinate systems.

    in_table: str, path to the input table or table name in map
    x_field, y_field, epsg_field: str, names of the coordinate fields and the field that contains the epsg code
    out_fc: str, path of the output point feature class
    out_epsg: epsg code of the output feature class
    
    """
    from pathlib import Path
    import arcpy
    out_ws = str(Path(out_fc).parent)
    out_name = str(Path(out_fc).name)
    out_sr = arcpy.SpatialReference(out_epsg)
    in_fields = arcpy.ListFields(in_table)
    
    # create output feature class
    arcpy.management.CreateFeatureclass(out_ws, out_name, "POINT", spatial_reference=out_sr)
    for f in in_fields:
        if f.type not in ["OID", "GlobalID"]:
            arcpy.management.AddField(out_fc, f.name, f.type)
    # fill output feature class
    field_names = [f.name for f in in_fields]
    with arcpy.da.InsertCursor(out_fc, field_names + ["SHAPE@"]) as icur:
        with arcpy.da.SearchCursor(in_table, field_names) as scur:
            for row in scur:
                row_dict = dict(zip(field_names, row))
                # create a point geometry from the input coordinates
                in_sr = arcpy.SpatialReference(row_dict[epsg_field])
                in_point = arcpy.PointGeometry(arcpy.Point(row_dict[x_field], row_dict[y_field]), in_sr)
                # project the point geometry to the output coordinate system
                out_point = in_point.projectAs(out_sr)
                # insert the row values and the new geometry into the feature class
                icur.insertRow(list(row) + [out_point])


# and then you can call it like this
xy_to_point_enhanced("TestTable", "X", "Y", "EPSG", "memory/TestPoints", 25832)