Changing coordinate of point shapefile in ArcGIS using ArcPy? [closed]

Question

I'm working on coding a function in ArcGIS trying to measure the number of pixels between a point (feature class) and its nearest building horizontally within a UDEM raster, then return the number.

Here is my idea:

1. Execute Zonal Statistics as Table and create a cursor to see whether the value in output table equals 0(the elevation of ground in my UDEM are set to 0 while elevation of building remain the real value)
2. If so(value=0), make x coordinate plus one, which means move the point 1 meter to the right(the resolution is 1×1) and refresh the coordinate of the point.
3. Iterate over step2 until the value doesn't equal 0,thus showing the location of point has moved to the "boundary" and I'll know the distance.

My problem now is having no idea how to perform Step2.

I tried some code to change its location but doesn't work.

Answer 1

To move a point in a feature class you can use the da.UpdateCursor and SHAPE@XY token. Select the point/points you want to move and then:

import arcpy
in_features = 'Pointlayer' #Change to match your layer name
x_shift = 1
y_shift = 1
with arcpy.da.UpdateCursor(in_features, ['SHAPE@XY']) as cursor:
    for row in cursor:
        cursor.updateRow([[row[0][0] + (x_shift),
                           row[0][1] + (y_shift)]])

With no selection all points will be moved. Backup your data before you try the code.

Answer 2

This is the case when raster solution is way faster and much more efficient. To implement it, convert your buildings to raster using any field, and:

arcpy.gp.RasterCalculator_sa('Con(IsNull("BUILDING_RASTER"),1,0)', "../Flow_Dir")
arcpy.gp.FlowLength_sa("Flow_Dir", "../FlowLength", "DOWNSTREAM", "")
arcpy.gp.EucDistance_sa("EASTERN BORDER", "../distance")

Sample both rasters to points table and compute difference between values in Distance and FlowLength fields. Points with 0 difference have no buildings to the East of them:

To compute distance West from points, change parameter in flow direction to UPSTREAM and replace DISTANCE by max(Distance)-Distance. Credits go to whuber, most certainly this was in one of his solutions, that I struggle to find.

Answer 3

To modify coordinates of a point, you just need to set its X or Y properties:

pointGeometry.X += 1

or, create a new point object:

newPoint = arcpy.Point(pointGeometry.X + 1, pointGeometry.Y)

See the relevant documentation here: https://pro.arcgis.com/en/pro-app/arcpy/classes/point.htm