<div dir="ltr">Sorry, I forgot the attached files. </div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/5/21 Mathieu Rajerison <span dir="ltr"><<a href="mailto:mathieu.rajerison@gmail.com" target="_blank">mathieu.rajerison@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello,<div><br><div><br></div><div>I have a lines shapefile that I converted to a graph object, following the method in <a href="http://www.mail-archive.com/r-sig-geo@r-project.org/msg05836.html" target="_blank">http://www.mail-archive.com/r-sig-geo@r-project.org/msg05836.html</a></div>

<div><br></div><div>I have determined the farthest nodes of the graph.</div><div><br></div><div>Then, I tried to determine the diameter's path  but the result is odd as the attached file shows.</div>
<div><br></div><div>The problem may come from the SpatialLines to graph conversion..But I can't see where the problem is...</div><div><br></div><div>Any help would be greatly appreciated...</div>
<div><br></div><div><div>I've put the example data I used if you want to give it a try</div></div><div><br></div><div>Here is the code:</div><div><br></div><div># READ</div><div>
l <- readOGR("IN/testline.shp", "testline")<br></div><div>l.break <- gIntersection(l, l)<br></div><div>l.split <- SpatialLines(lapply(1:length(l.break@lines[[1]]@Lines), function(x) Lines(list(l.break@lines[[1]]@Lines[[x]]), ID=as.character(x))))<br>

</div><div><br></div><div>#GRAPH</div><div>edges = do.call(rbind,lapply(1:length(l.split),function(x){as.vector(t(l.split@lines[[x]]@Lines[[1]]@coords))}))<br></div><div><div>lengths = sqrt((edges[,1]-edges[,3])^2+(edges[,2]-edges[,4])^2)</div>

<div>  </div><div>froms = paste(edges[,1],edges[,2])</div><div>tos = paste(edges[,3],edges[,4])</div><div>g = graph.edgelist(cbind(froms,tos),directed=FALSE)</div><div>E(g)$weight=lengths</div><div><br></div><div># DIAMETER</div>

<div><div>sel <- get.diameter(g)</div><div>diam <- l.split[sel-1]</div></div><div>plot(diam)</div><span class="HOEnZb"><font color="#888888"><div><br></div><div><br></div></font></span></div><span class="HOEnZb"><font color="#888888"><div>
Mathieu</div></font></span></div></div>
</blockquote></div><br></div>