三维体积和 C/C++ 链接 

LibraryLink 扩展支持直接发送或来自外部程序库的二维和三维图像.

以下是对一个 Image 或 Image3D 对象的像素值取负的例子. color_negate 函数用于任何通道数或数据类型的交插或平面图像. 由此参见完整范例

CCompilerDriver 程序包可以创建可执行文件、程序库和利用编译器安装到计算机的 C/C++ 源代码对象文件.

In[1]:=

X Needs["CCompilerDriver`"];

显示完整的 Wolfram 语言输入隐藏输入

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

从以上编码创建动态链接程序库.

In[3]:=

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

在此，程序库已加载并且 colorNegateCpp 函数调用被链接到与 C++ 过程.

In[4]:=

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

Out[4]=

导入体积.

In[5]:=

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

Out[5]=

从8个角点开始，用区域增长分割噪声背景.

In[6]:=

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

Out[6]=

从 C++ 程序库运用 colorNegateCpp 来转换背景并创建前景掩模.

In[7]:=

X mask = colorNegateCpp[background]

Out[7]=

用掩模体积去除噪声背景.

In[8]:=

X clean = ImageMultiply[vol, mask]

Out[8]=