ベクタをポリゴンで分割する SML
(Datum ワークステーション)

Ver.2.2 (UTF-8)

サンプルデータは、TNTmips/liteサンプルデータ(cb_data)中の、"土壌図"と"公的土地境界" を使いました。
前景の黄色いラインデータが "土地境界"(分割する側)で、後ろのハッチのデータが "土壌図"(分割される側) です。"土地境界"ポリゴンは、周囲の図郭線が欠けていますが、閉じたポリゴンは全部で36個あります。

「解析処理」>「ＳＭＬ」>「スクリプトの編集...」を選択します。
ダウンロードしたＳＭＬを開きます。
開いたら、ウィンドウ右下の[実行...]ボタンを押します。
「"Voperator"用のベクタオブジェクトを選択してください」と出ますので、分割する側(オペレータ)のベクタ(ポリゴン)を選択します。
次に「"Vsource"用のベクタオブジェクトを選択してください」と出ますので、分割される側(ソース)のベクタ(ポリゴン)を選択します。
「Please select Output RVC File !」と出ますので、分割後のオブジェクトを出力するＲＶＣファイル名を指定します。既定値は SPLIT.RVC です。

これまで分割後のベクタの名前は、"ソースのオブジェクト名"_ポリゴンの内部番号(Internal.ElemNum)でしたが、Internal.ElemNum ではなく、ユーザが作成したフィールドも選べるようにしました。
処理が開始します。１個１個の切り出しには、それなりの時間がかかります。今回のliteデータでも全部終わるまで約３分(Pentium4 2.4GHz, メモリ 504MB, WindowsXP Home)かかりました。コンソールウィンドウを見ると、何個目の切り出しをやっているか確認できます。
処理が終わると、「Finish !!」ウィンドウが出ます。
結果はこんな感じです。
分割したオブジェクトの名前は、自動的に"ソースのオブジェクト名"_"オペレータのポリゴン番号"と付けられます。今回のデータでは、CBSOILS_1から_36まであります。
次に、「20万分の１数値地質図」(産業技術総合研究所地質調査総合センター発行)を２次メッシュで分割してみました。
使用したデータは「東京」(1次メッシュ 5339)の図幅で、これを別途用意した２次メッシュで 8 x 8 = 64個に分割しました。
処理時間は同じパソコンで約５分でした。