Skip to content

Instantly share code, notes, and snippets.

@ctrueden
Last active June 12, 2023 22:52
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save ctrueden/62f758807d2f0d6d6d1dafc3f001226d to your computer and use it in GitHub Desktop.
Save ctrueden/62f758807d2f0d6d6d1dafc3f001226d to your computer and use it in GitHub Desktop.
Convert a mesh to a binary mask using ImageJ Ops
Airplane dataset from:
https://people.math.sc.edu/Burkardt/data/ply/ply.html
import java.io.File;
import java.io.IOException;
import net.imagej.Dataset;
import net.imagej.ImageJ;
import net.imagej.mesh.Mesh;
import net.imagej.mesh.io.ply.PLYMeshIO;
import net.imagej.ops.OpService;
import net.imglib2.RandomAccessibleInterval;
import net.imglib2.converter.Converters;
import net.imglib2.type.logic.BitType;
import net.imglib2.type.numeric.integer.UnsignedByteType;
public class Voxelize {
public static void main(final String... args) throws IOException {
// Replace these paths with your desired ones.
String inPath = "/home/curtis/data/3d/burkardt/ply/airplane.ply";
String outPath = "/home/curtis/Desktop/airplane.tif";
System.out.println("Initializing ImageJ2...");
final ImageJ ij = new ImageJ();
final OpService ops = ij.op();
System.out.println("Opening mesh...");
//Object mesh = ij.io().open(path);
Mesh mesh = new PLYMeshIO().open(inPath);
int w = 128, h = 128, d = 128;
System.out.println("Voxelizing mesh at " + w + " x " + h + " x " + d + "...");
RandomAccessibleInterval<BitType> mask = ops.geom().voxelization(mesh, w, h, d);
// Convert image: [0, 1] -> [0, 255]
// Because ImageJ2/SCIFIO cannot save BitType images to TIFF. :-(
System.out.println("Wrapping voxelized mesh to 8-bit...");
RandomAccessibleInterval<UnsignedByteType> mask8bit = Converters.convert(mask, (in, out) -> out.setReal(in.get() ? 255 : 0), new UnsignedByteType());
System.out.println("Wrapping to ImageJ2 Dataset...");
// Convert RAI -> Dataset. Because ImageJ2/SCIFIO needs Dataset, not any old RAI. :-(
Dataset dataset = ij.dataset().create(mask8bit);
System.out.println("Saving result to TIFF file: " + outPath);
File outFile = new File(outPath);
if (outFile.exists()) outFile.delete();
ij.scifio().datasetIO().save(dataset, outPath);
System.out.println("All done! " + String.valueOf(Character.toChars(0x1f36d)));
}
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment