Volúmenes en 3D y enlace a C/C++  

LibraryLink ha sido ampliada para dar soporte a imágenes en 2D y 3D que son directamente enviadas a/de bibliotecas externas.

Aquí hay un ejemplo para negar los valores de píxel de un objeto Image o Image3D. La función color_negate trabaja con imágenes intercaladas o planas con cualquier número de canales y cualquier tipo de datos. Vea el ejemplo completo aquí.

El paquete CCompilerDriver le permite construir archivos de ejecutables, bibliotecas y objetos desde un código fuente C/C++ utilizando compiladores instalados en su equipo.

In[1]:=

X Needs["CCompilerDriver`"];

muestre la entrada completa de Wolfram Languageoculte la entrada

In[2]:=

X code = " #include \"WolframLibrary.h\" #include \"WolframImageLibrary.h\" template <> char maxValue<char>() { return 1; } template <> raw_t_ubit8 maxValue<raw_t_ubit8>() { return 255; } ... template <typename T> static void icolor_negate(void *out0, const void *in0, mint length) { mint ii; T *out = reinterpret_cast<T *>(out0); const T *in = reinterpret_cast<const T *>(in0); for (ii = 0; ii < length; ii++) { out[ii] = maxValue<T>() - in[ii]; } } /* Negate image colors */ EXTERN_C DLLEXPORT int color_negate(WolframLibraryData libData, mint Argc, MArgument *Args, MArgument res) { ... WolframImageLibrary_Functions imgFuns = \ libData->imageLibraryFunctions; ... type = imgFuns->MImage_getDataType(image_in); length = imgFuns->MImage_getFlattenedLength(image_in); data_in = imgFuns->MImage_getRawData(image_in); data_out = imgFuns->MImage_getRawData(image_out); if (data_in == NULL || data_out == NULL) goto cleanup; switch (type) { case MImage_Type_Bit: icolor_negate<char>(data_out, data_in, length); break; case MImage_Type_Bit8: icolor_negate<raw_t_ubit8>(data_out, data_in, length); break; ... } }";

Cree un enlace dinámico de biblioteca a partir del código anterior.

In[3]:=

X colornegateLibrary = CreateLibrary[code, "libcolornegate", "Language" -> "C++"];

Aquí, la biblioteca es descargada y la función llamada colorNegateCpp está enlazada al proceso de color_negate C++.

In[4]:=

X colorNegateCpp = LibraryFunctionLoad[colornegateLibrary, "color_negate", {LibraryDataType[Image | Image3D]}, LibraryDataType[Image | Image3D]]

Out[4]=

Importe un volumen.

In[5]:=

X vol = ImageRotate[ ExampleData[{"TestImage3D", "MRbrain"}], {\[Pi]/2, {1, 0, 0}}]

Out[5]=

Segmente el fondo ruidoso usando la región de crecimiento, a partir de las ocho esquinas.

In[6]:=

X corners = Tuples@Thread[{1, ImageDimensions[vol]}]; background = RegionBinarize[vol, corners, 0.1]

Out[6]=

Aplique colorNegateCpp desde una biblioteca C++ para invertir el fondo y crear una máscara de primer plano.

In[7]:=

X mask = colorNegateCpp[background]

Out[7]=

Enmascare el volumen para remover el fondo ruidoso.

In[8]:=

X clean = ImageMultiply[vol, mask]

Out[8]=