3D立体とC/C++との接続  

LibraryLink は2Dおよび3D画像を直接外部ライブラリに送ったり，外部ライブラリから送られたりすることをサポートするよう拡張されている．

次は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関数の呼出しがcolor_negate CC++プロシージャにリンクされる．

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]=