Count polylines crossing each raster cell

Question

I'm looking for the best way to count polylines crossing each raster cell and assign that count value to the raster cell using arcpy.

On input I have a shapefile with hundred thousands of polylines and a raster with all cells 0 values.

The result should be a raster which cells value should be a number of crossing polylines from the shapefile.

Which ArcPy tools should I use to perform this analysis?

Hornbydd · Accepted Answer · 2023-06-16 13:54:03Z

2

This could easily be built in model builder. I would:

Convert your raster to points
Convert your points back into a raster but using FID as the value column, this will create unique raster cells
convert this new raster to a polygon dataset
Run the intersect tool, this will give your count of lines per individual pixel.
Join count to polygon and rasterize if required.

answered Jun 16, 2023 at 13:54

Hornbydd

44.2k5 gold badges41 silver badges82 bronze badges

Great. Regarding point 4: Intersect tool I run like this: arcpy.analysis.Intersect([lines_file, out_grid], out_intersect, "", "", "LINE") It gives me single lines cut to pixel borders. I still need to count lines in the pixels before performing join. Of course now it is easier to count, because I can count using pixels Id. Or is there another way to get count of the lines per pixel?
– Andžej Miloš
Commented Jun 21, 2023 at 7:02
1

You are correct, you use the ID and run that through the summary statistics tool, that should be step "4.5". :)
– Hornbydd
Commented Jun 21, 2023 at 9:24

Add a comment |

1 Answer 1