I have a 3D (x,y,z) single point feature class which points are equally spaced, like this:
As the distance between each points are more than one meter, I was wondering if there's a possibility to automatically identify data by adding new points like this (maybe also in arcpy):
You can use Spatial Analyst extension Interpolate Tools (IDW) to create a Raster from the input point features. Then convert back to points if that is your needed format. Make sure to specify environment setting or specify cellSize to desired resolution (ie 1 = 1m) if you are in projected coordinate system
import arcpy
from arcpy.sa import *
arcpy.CheckOutExtension("Spatial")
inFeats = "C:\mydata.gdb\3dpoints"
cellsize= 1
zField = "zfld"
ras = Idw(inFeats, zField, cellsize)
outfeats = "C:\mydata.gdb\3dpoints_1m"
arcpy.RasterToPoint_conversion(ras, Outfeats,"Value")
If you have Geostatistical Analyst you can use the Inverse Distance Weighting tool to fill in your gaps.
I think IDW is the most appropriate for your use-case of the various interpolation tools available, though they do require Spatial Analyst license.
It will output a raster. Then, you can convert back to points at your desired resolution, then use Extract Values to Points to join back your z values. See: Assigning Raster Value to Point Data
