O problema é o corte usando a vtkImpliticFunction em questão. Não consigo fazer funcionar direito. Hora de voltar pra prancheta.
Consegui fazer o corte corretamente seguindo o algoritmo de raytrace. (http://en.wikipedia.org/wiki/Point_in_polygon), mas está lerdo: no momento versão 1 do codigo, levo 0.0172 ms por voxel e 105828 ms no total. Esse tempo é inaceitável.
Eu realmente tenho que avaliar todos os pontos ou posso me concentrar em uma região?
Eu realmente tenho que percorrer toda a lista de pontos do contorno?
Eu estou alocando memória de maneira inteligente?
Eu ainda estou no single-thread pra facilitar a depuração. Num calculo simplório o tempo cairia pra 13000 se posto em 8 threads (provelmente cairia pra 20000 devido aos rituais dos threads e semáforos).