Creating a polyline from a start to end point, following points of the same elevation

[SimbarasheKapfumo]

How are you  people . I have an interesting question which i hope  u might help me with . 

I have  a DEM  and  two point shapefiles. I want  to draw  a line from point 1 in shapefile 1 to any point  in shapefile 2 which has the  same  elevation with point 1. The  line  should follow points of the  same elevation in the DEM hence the  line will likely be an winding line rather than a straight one.

 Thank  you  for your  time.

[DanPatterson_Retired · Accepted solution]

Contour list .... or even contour, then convert to polylines.

Finally, intersect with your points using Split at point and keep the bits you want

[SimbarasheKapfumo]

Thank you  . Its a very long process though  because i had to  do some  clean up to remove  the  other  unwanted contours ..  

[DanPatterson_Retired]

Or if the number is smaller, select the ones you want and make a new featureclass from them

[DanPatterson_Retired]

Do you have an image of the input and final result?

[SimbarasheKapfumo]

still writing code to select the wanted lines. there a lot of conditions the lines have to meet and there are a lot of lines.

[SimbarasheKapfumo]

is there any way of writing the statement below in arcpy 

arcpy.gp.RasterCalculator_sa("(\"%FlowAcc%\" > 5000)", Rstream)
i want to select values greater than 5000 in Flow Accumulation raster (FlowAcc) and save them to Rstream .
when i put it in python it returned an error : 'int' object has no attribute 'save

[DanPatterson_Retired]

You have to provide the raster's name... right now the "\Flo..." stuff is an integer I bet and Rstream needs to be a filename

[SimbarasheKapfumo]

Rstream is a file. the code is so long i just shorted it. But whole code is almost the same as the 1 i used and it worked well 

arcpy.gp.RasterCalculator_sa("Con(\"%Fill (2)%\" == \"%ZonalStast_Max%\",\"%Fill (2)%\")", MaxRasterPoints)