Volume 3D e ligação com C/C++ 

A LibraryLink é estendida para suportar imagens 2D e 3D sendo diretamente enviadas de/para bibliotecas externas.

Aqui temos um exemplo de inversão dos valores dos pixels de um objeto Image ou Image3D. A função color_negate funciona com imagens planares ou intercaladas com qualquer número de canais ou qualquer tipo de dados. Veja o exemplo completo aqui.

O pacote CCompilerDriver permite que você crie arquivos executáveis​​, bibliotecas e arquivos de objeto a partir de código fonte C/C++ utilizando compiladores instalados em seu computador.

In[1]:=

X Needs["CCompilerDriver`"];

mostre o input completo de Wolfram Languageoculte o input

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; ... } }";

Crie uma biblioteca de ligação dinâmica a partir do código acima.

In[3]:=

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

Aqui a biblioteca é carregada e a chamada da função colorNegateCpp é ligada à função C++ color_negate.

In[4]:=

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

Out[4]=

Importe um volume.

In[5]:=

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

Out[5]=

Segmente o fundo ruidoso usando crescimento de regiões a partir dos oito cantos.

In[6]:=

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

Out[6]=

Aplique colorNegateCpp da biblioteca C++ para inverter o fundo e criar uma máscara no primeiro plano.

In[7]:=

X mask = colorNegateCpp[background]

Out[7]=

Aplique uma máscara no volume para remover o fundo ruidoso.

In[8]:=

X clean = ImageMultiply[vol, mask]

Out[8]=