Commit 0e660263 authored by Michiel Cottaar's avatar Michiel Cottaar
Browse files

BUG: don't divide target density by vertex area twice

parent f92f6af2
Pipeline #4762 passed with stage
in 24 minutes and 10 seconds
......@@ -75,10 +75,7 @@ surf_mask = nib.load(args.surf_mask).darrays[0].data != 0
white_surf = CorticalMesh.read(args.white)[surf_mask]
if args.target_density is None:
target_density = (
Cortex([white_surf, CorticalMesh.read(args.pial)[surf_mask]]).wedge_volume() /
white_surf.size_faces()
)
target_density = Cortex([white_surf, CorticalMesh.read(args.pial)[surf_mask]]).wedge_volume()
else:
target_density_vertex = nib.load(args.target_density).darrays[0].data[surf_mask]
target_density = white_surf.graph_connection_point().T.dot(target_density_vertex) / 3.
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment