Python: Yüzey alanını hesaplamak için bir voksel düzeneğini örter

Question

3D olarak bir taşı temsil eden bir voksel düzeneğim var. Taşımayı ve 0'ı başka yerde oluşturan voksellere 1'in atandığı ikili bir numpy 3D dizisidir.

1. üzerinde yüzey alanını hesaplamak ağ gibi saran yüzeyini
2. oluşturun: Ben istiyorum.

Ama nasıl?

Answer 1

o yararlı olabilir eğer. Yürüyen Küpler algoritması iyidir, ancak ikili dizilerde iyi çalışmaz. Böylece: 1) Gauss Filtre üçgen alanları toplamı) o (scikit-görüntü aracı) 3 mesh (scipy.filters) 3D dizisine 2) Bando küpleri algoritması uygulanır (scikit-görüntü aracı)

Answer 2

Uygulamanıza bağlı olarak birkaç seçenek görebiliyorum: 3D dışbükey gövde veya yürüyen küpler. senin dizi sen scipy.spatial.ConvexHull kullanabilirsiniz Daha sonra 1

eşit olduğu

3D dışbükey

Öncelikle, tüm voksel 3D bir dizi elde etmek gerekecektir koordine eder. Scipy yöntemleri, temel olarak qhull kitaplığı için bir sarmalayıcıdır. qhull kullanırken, çıkış seçeneklerinden (FA) biri, toplam alanı ve birimi hesaplamaktır.

Burada biraz daha bilgi bulabilirsiniz: Bando Küpler

Eğer içbükeyliliğe tutmak gerekiyorsa http://scipy.github.io/devdocs/generated/scipy.spatial.ConvexHull.html http://www.qhull.org/html/qconvex.htm

ardından önceki seçenek yeterli olmayacaktır. Yürüyen küpler iyi bir alternatif olacak. Scikit-image'de zaten bir python uygulaması var.

http://scikit-image.org/docs/dev/auto_examples/plot_marching_cubes.html

Yine önceki durumda olduğu gibi tüm noktaları elde etmek gerekecektir. Sen fiiller hesaplamak/ile karşı karşıya:

http://scikit-image.org/docs/dev/api/skimage.measure.html#marching-cubes

Ve alanı hesaplayabilirsiniz: Ben sorunu çözüldü

http://scikit-image.org/docs/dev/api/skimage.measure.html#mesh-surface-area