Pairwise shortest distance between points in a road network

At least one row of your roads is causing this error when SpatialLinesNetwork tries to get its length:

> SpatialLinesLengths(l[65403,])
Error in FUN(X[[i]], ...) : non-finite line lengths

This row appears to have zero length, which is messing up the length calculation when it is using lat-long coordinates:

> l[65403,]
class       : SpatialLinesDataFrame 
features    : 1 
extent      : -41.98302, -41.98302, -15.25437, -15.25437  (xmin, xmax, ymin, ymax)

I suspect a divide-by-zero somewhere when working out the distance on a sphere by lat-long.

The adjacent segments are fine:

> SpatialLinesLengths(l[65402,])
[1] 1.813751
> SpatialLinesLengths(l[65404,])
[1] 1.802447

You can trick it into not using lat-long coordinates - there is an option but I think having a CRS on the lines overrides this anyway. So set the CRS to nothing:

> proj4string(l) = ""

And then I can build a spatial lines network without errors:

> ln = SpatialLinesNetwork(l)
> 

But that will give you a spatial network based on euclidean lat-long distance which might not be very valid.

Or you could project to an appropriate cartesian grid system for your data and use that.

Alternatively you could try an eliminate these rows from your data, but in some cases it seems to be caused by zero-length segments within road segments, which is hard to track down.