The code would look like the code below. Notice the following changes:
- Line 6 holds the interval, I changed it from 1 to 3, so now it picks points every 3 meters on each contour line.
- Line 14 contains a where clause that filters the contour lines with Contour > 1000
- Line 21 applies the where clause
import arcpy, os
fc_in = r"C:\Forum\Pasture\gdb\Contours.gdb\contours2"
fc_out = r"C:\Forum\Pasture\gdb\Contours.gdb\contour_pnts03"
fld_contour = "Contour"
interval = 3 # pick points every x meter on every contour line
# create empty output featureclass
sr = arcpy.Describe(fc_in).spatialReference
out_ws, out_name = os.path.split(fc_out)
arcpy.CreateFeatureclass_management(out_ws, out_name, "POINT", fc_in, "DISABLED", "DISABLED", sr)
# create a where clause
where = "{0} > {1}".format(arcpy.AddFieldDelimiters(fc_in, fld_contour), 1000)
# add field with contour value to output fc
arcpy.AddField_management(fc_out, fld_contour, "DOUBLE")
cnt = 0
with arcpy.da.InsertCursor(fc_out, ("SHAPE@", fld_contour)) as curs_out:
with arcpy.da.SearchCursor(fc_in, ("SHAPE@", fld_contour), where_clause=where) as curs_in:
for row_in in curs_in:
polyline = row_in[0]
contour = row_in[1]
max_len = polyline.length
d = 0
while d < max_len:
pnt = polyline.positionAlongLine(d, False)
curs_out.insertRow((pnt.firstPoint, contour, ))
cnt += 1
if cnt % 100 == 0:
print "Processing point: {0}".format(cnt)
# pnt.firstPoint
d += interval
curs_out.insertRow((polyline.lastPoint, contour, ))
cnt += 1
