scikits-image/scikits-image — Issue #72: MCP_geometric.find_costs fails on large graphs
Author: ozak
Created: 2011-10-26 18:52:15+00:00
State: open

Hi,

I have been using graph.MCP_Geometric from the latest version of scikits.image 0.3 in order to find minimum travel costs with the find_costs function. On most cost surfaces it works without problems, but when I use large cost matrices it doesn't even try to find those costs. It looks as if the algorithm stopped before the second iteration, since all the starts have a cumulative cost of 0, while the rest of the matrix has inf. Currently I am trying to use it on a cost matrix of size 12837 x 43345. If I split it by half it works, but given what I want to do now, I cannot split the cost matrix and work on subsets anymore.

Is there a size limit? Any ideas as to what might be causing this?

I am working on a MacPro with 64GB memory, 2 CPU's, OS X Lion, Python 2.72, Numpy 1.6.1.

I appreciate any help for solving this.