[Insight-users] Automatic Junction Detection--runtime error

[Pascal]

Dear all

I've compiled and built the 2D test version of the code using ITK 3.8. I used the IJ_324_Neuroninput.mha input data.

http://www.insight-journal.org/browse/publication/324


Upon debugging, I obtain a runtime error :  map/set iterator not incrementable.

When I execute the **** lines below..... it fails on the next iteration.


 while( !distQueue.empty() )
          {
            IndexType oneIndex = distQueue.front();
            distQueue.pop();
            for(dsIt = distSet.begin(); dsIt!=distSet.end(); ++dsIt)
           
 {
              IndexType anotherIndex = *dsIt;
              OffsetType offset = anotherIndex - oneIndex;
              typename OffsetType::OffsetValueType offsetSumSquared = 0;
              for(int j=0; j<ImageDimension; j++)
                offsetSumSquared += offset[j]*offset[j];
              if( offsetSumSquared <= 3 )
              {
       ****        
 distQueue.push(anotherIndex);
       ****         distSet.erase(dsIt);
                
              }
            }
          }



Any clues on this would be appreciated.

regards

Pascal



      __________________________________________________________________
Connect with friends from any web browser - no download required. Try the new Yahoo! Canada Messenger for the Web BETA at http://ca.messenger.yahoo.com/webmessengerpromo.php
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20090722/2d48fbf8/attachment-0001.htm>