3D-Volumen und Anbindung an C/C++  

LibraryLink unterstützt nun die direkte Übertragung von 2D- und 3D-Bildern in und aus externen Bibliotheken.

Im folgenden Beispiel werden Pixelwerte eines Image- oder Image3D-Objekts negiert. Die color_negate-Funktion kann sowohl auf planare Bilder als auch auf Bilder, die über mehrere Ebenen verzahnt aufbereitet werden, angewendet werden, unabhängig von der Anzahl der Kanäle und des Datentyps. Das komplette Beispiel finden Sie hier.

Das CCompilerDriver-Package erlaubt es Ihnen, ausführbare Dateien, Bibliotheken und Objektdateien ausgehend von C/C++ Quelltext mit jenen Compilers zu erstellen, die auf Ihrem Rechner installiert sind.

In[1]:=

X Needs["CCompilerDriver`"];

Die gesamte Wolfram-Language Eingabe zeigenEingabe verbergen

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

Erstellen Sie aus dem obigen Code eine Bibliothek aus dynamischen Links.

In[3]:=

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

Die Bibliothek ist nun geladen und der Abruf der colorNegateCpp-Funktion ist mit dem color_negate C++ Vorgang verbunden.

In[4]:=

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

Out[4]=

Importieren Sie ein Volumen.

In[5]:=

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

Out[5]=

Segmentieren Sie den verrauschten Hintergrund mithilfe des Region Growing-Verfahrens ausgehend von den acht Ecken.

In[6]:=

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

Out[6]=

Wenden Sie colorNegateCpp aus der C++ Bibliothek an, um den Hintergrund zu invertieren und eine Maske des Vordergrunds zu erzeugen.

In[7]:=

X mask = colorNegateCpp[background]

Out[7]=

Maskieren Sie das Volumen, um den verrauschten Hintergrund zu entfernen.

In[8]:=

X clean = ImageMultiply[vol, mask]

Out[8]=