見積競争の公告 国立大学法人筑波大学において 次のとおり請負契約について見積競争を実施します 1. 見積競争に付する事項 (1) 件 名 BNCT 用治療計画システムへの患者位置追従機能追加 (2) 規格及び数量 仕様書のとおり (3) 納入期限 平成 31 年 3 月 15 日 (4) 納入場所 仕様書のとおり 2. 仕様書の交付場所等 (1) 場 所 茨城県つくば市天王台 1 丁目 1 番地の1 国立大学法人筑波大学財務部契約課 (2) 連 絡 先 担当折戸 電話番号 029-853-5099 3. 技術審査申請書等の提出先及び提出期限この見積競争に参加を希望する者は 見積書提出の注意事項 8 により技術審査に合格した者だけが見積競争に参加できるものとする (1) 提出先上記 2 の仕様書の交付場所と同じ (2) 提出期限平成 31 年 1 月 24 日 14 時 00 分 4. 見積書等提出場所及び提出期限 (1) 提出先上記 2 の仕様書の交付場所と同じ (2) 提出期限平成 31 年 1 月 31 日 14 時 00 分見積競争結果については 電話等により行う 5. 見積の方法 (1) 国立大学法人筑波大学契約事務取扱細則を熟知し 仕様書及び契約条項を承諾のうえ 見積るものとする (2) 契約決定に当たっては 見積金額に 8 パーセントに相当する額を加算した金額 ( 当該金額に 1 円未満の端数があるときは その端数金額を切り捨てるものとする ) をもって契約金額とするので 見積者は 消費税及び地方消費税に係る課税事業者であるか免税事業者であるかを問わず 契約金額の 108 分の 100 に相当する金額で見積るものとする 6. 見積競争に参加する者に必要な資格 (1) 国立大学法人筑波大学財務規則施行規程 ( 以下 規程 という ) 第 46 条の規定に該当しない者であること なお 未成年者 被保佐人又は被補助人であって 契約締結のために必要な同意を得ている者は 同条中 特別の理由がある場合に該当する (2) 規程第 47 条の規定に該当しない者であること (3) 平成 30 年度の全省庁統一資格又は本学の競争参加資格において 関東 甲信越地域の 役務の提供 で有効な資格を有している者であること 又は当該資格を有しない者であって 過去 1 年以内に本学との取引実績を有する者であること (4) 契約担当役から取引停止の措置を受けている期間中の者でないこと 7. 契約の方式 (1) 最低価格の見積書を提出した者及び次順位者を契約予定者として 価格交渉を行う (2) 契約予定者との価格交渉により 契約相手方及び契約金額を決定する
平成 31 年 1 月 17 日 以 上 国立大学法人筑波大学契約担当役財務担当副学長石野利和
見積書提出の注意事項 1 技術審査書類及び参考見積書提出日時期限 平成 31 年 1 月 24 日 14 時 00 分 場所 305-8577 茨城県つくば市天王台 1-1-1 国立大学法人筑波大学本部棟 3 階財務部契約課 担当折戸 TEL:029-853-5099 2 見積書提出日時期限平成 31 年 1 月 31 日 14 時 00 分場所上記 1 の参考見積書提出場所と同じ 3 見積書作成の注意 (1) 見積金額は算用数字を用いて明確に記入すること (2) 住所氏名を記入し押印すること (3) 日付を必ず記入すること 4 上記注意事項に適合しない見積書は無効とすることがある 5 契約決定に当たっては 見積書に記載された金額に当該金額の 8 パーセントに相当する額を加算した金額 ( 当該金額に 1 円未満の端数があるときは その端数金額を切り捨てるものとする ) をもって契約金額とするので 見積者は 消費税及び地方消費税に係る課税事業者であるか免税事業者であるかを問わず 見積もった契約金額の 108 分の 100 に相当する金額を見積書に記載すること 6 いったん提出された見積書は引換え 変更 取消しをすることができない 7 仕様書 その他の書類は見積書提出の当日返却すること 8 技術審査申請書類等提出書類この見積競争に参加を希望する者は 下記の書類を参加者負担において作成し期日までに提出すること 技術審査に合格した者だけが見積競争に参加できるものとする (1) 平成 30 年度の資格審査結果通知書 ( 全省庁統一資格又は国立大学法人筑波大学の競争参加資格 ) の写し (2) 上記の資格 8(2) を有しない場合は過去 1 年以内の本学との契約実績 (3) 参考見積書 ( 工数等を明記 ) (4) 作業単価表 (5) CT MRI 等医療画像 DICOM フォーマット及び DICOM-RT フォーマットのデータを扱った治療計画システム関連のソフトウェア開発の経験があること 特に計算コードにモンテカルロ法を扱うことから モンテカルロベースの医療用線量評価システムの開発経験を有すること (6) 汎用モンテカルロ計算コード PHITS を用いた線量計算の経験があり PHITS のインプットファイル及び出力結果のフォーマットを全て理解できること (7) 医療画像データを基に構築した人体の 3 次元モデルを 3 次元描画させることから このような画像処理技術の制作を実施できること (8) Java によるシステム開発を行うため システム開発担当者は Oracle 認定の Java プログラマーの資格を有していること (9) 個人情報である医療データを取り扱う技術の開発を行うことから情報セキュリティマネ
ジメントシステム (ISO27001/ISMS) 認証を取得していること また 情報セキュリティマニュアルを提出できること
仕様書 1. 件名 BNCT 用治療計画システムへの患者位置追従機能追加 2. 内容別紙のとおり 3. 期限平成 31 年 3 月 15 日 4. 場所筑波大学陽子線医学利用研究センター 3 階 306 号室 5. 実施作業に係る条件本作業は医療分野と工学分野の境界領域 ( 医学物理分野 ) の技術に関する製作作業であり また モンテカルロ法による線量計算も実施するシステムの製作であり原子力学 放射線工学の分野のノウハウも必要となる 従って本作業を納期までに確実に実施するため 本作業実施においては以下に示す条件を満たすこと (1)CT MRI 等医療画像 DICOM フォーマット及び DICOM-RT フォーマットのデータを扱った治療計画システム関連のソフトウェア開発の経験があること 特に計算コードにモンテカルロ法を扱うことから モンテカルロベースの医療用線量評価システムの開発経験を有すること (2) 汎用モンテカルロ計算コード PHITS を用いた線量計算の経験があり PHITS のインプットファイル及び出力結果のフォーマットを全て理解できること (3) 医療画像データを基に構築した人体の 3 次元モデルを 3 次元描画させることから このような画像処理技術の制作を実施できること (4)Java によるシステム開発を行うため システム開発担当者は Oracle 認定の Java プログラマーの資格を有していること (5) 個人情報である医療データを取り扱う技術の開発を行うことから情報セキュリティマネジメントシステム (ISO27001/ISMS) 認証を取得していること また 情報セキュリティマニュアルを提出できること 6. 成果物書類名 提出時期 提出部数 (1) 工程表 着手前 3 (2) 打ち合せ議事録 打ち合せ後すみやかに 3 (3) 製作報告書 すべての作業終了後速やかに 3 (4) 制作したプログラムのソース すべての作業終了後速やかに 1 (5) 操作マニュアル すべての作業終了後速やかに 1 CD-ROM 等による電子媒体での提出 7. 検収製作したソフトウェアが別紙業務内容 4. に示す動作確認試験に合格し 6. に定める全ての成果物とともに製作したソフトウェアが提出され 貸与したデータ等が返却されたことをもって検収とする
8. 支払支払は 検査終了後適法な請求書を受理した日から起算して 40 日以内に支払うものとする 9. その他 (1) 本業務に伴って貸与する医療画像データ DICOM-RT データ 第三者への提示 コピー等を行わないこと 作業終了後すみやかに返却すること (2) 本業務遂行のために必要な図書類は貸与する 貸与した図書類は作業終了時に返却すること (3) 本業務の内容は筑波大学の承認なく第三者に開示しないこと (4) 本契約に必要な細目は 国立大学法人筑波大学契約事務取扱細則によるものとする (5) 他詳細については 本学担当職員の指示によるものとする
別紙 1. 概要筑波大学陽子線医学利用研究センターでは 加速器中性子による中性子捕捉療法 (BNCT) の治療を実現するため 治療実施に不可欠な治療計画システム ( 以下 TPS) や患者位置合わせ装置などの周辺装置の開発研究を行っている さらに 科研費 基盤研究 (B) の研究テーマ BNCT 照射中の患者の動きに追従してリアルタイム線量評価を行うシステムの開発 として BNCT 実施中の患者の位置変動に追従して リアルタイムで線量評価を実施するシステムの基盤技術研究を実施している 本制作作業は この科研費 基盤研究 B の開発研究の一環として これまでに開発した BNCT 用治療計画システム ( 開発コード : ツクバプラン ) に対し 照射中の患者の位置情報を取り込み ツクバプラン上の患者モデルの位置を実患者位置に一致させて その条件下での線量評価を実施できるシステム基盤を制作するものである 2. 開発条件 動作環境 O.S. : Windows 8, 10 系 (x64 Edition) もしくは Linux ハードウェア : Windows もしくは Linux が動作する PC 開発言語 : JAVA C 言語 扱うシステム : 筑波大学でこれまで開発したモンテカルロ治療計画システム 扱う画像データ : DICOM フォーマット 扱う医療データ : DICOM-RT フォーマット 扱う計算コード : 汎用モンテカルロコード PHITS 扱う核データ : JENDL 4.0 グラフィックカードなどのハードウェア OS のマイナーバージョンアップに依存しないこと 取り扱うデータの大規模化を踏まえ 64bit システムとして開発を行うこと 作業を的確かつ効率的に進めるため 作業に必要な打ち合わせを適宜行うとともに 打ち合わせ後速やかに議事録を作成すること また 製作物の不具合への対応 修正作業を速やかに実施できること 著作権はすべて筑波大学に譲渡するものとし ソフトウェアに対して Copyright 2018 University of Tsukuba. All Right Reserved. 及びバージョン番号等の記述を入れること 3. 製作仕様システムの全体概要これまでに筑波大学で開発した汎用モンテカルロ治療計画システム ( 以下 ツクバプラン )( 1) に対し 別途開発している患者位置合わせ 位置計測装置 ( 以下 患者位置計測装置 )( 2) から出力された患者の実際の位置情報 ( 3) をツクバプランに取り込ませ その情報を基にツクバプラン上の患者モデルのビーム孔に対する位置条件を実際の位置に合うように移動させ 実際の位置での線量評価を行う機能を追加する また この機能追加に合わせて 以下に詳述する関連する機能を追加するとともに 既存の機能を高度化する 追加 高度化するそれぞれの機能の詳細仕様を以下に示す また 詳細な仕様 要求される機能については適宜の打ち合わせによって決定するものとする 1 ツクバプラン : 患者の CT MRI などの医療画像データを取り込んで患者の 3 次元モデル ( 患者 3D モデル ) を構築し 関心領域や病巣領域を設定して放射線外部照射シミュレーションを
実施して 各関心領域 病巣領域に付与される線量率を算出する 線量計算には モンテカルロ輸送計算コード :PHITS を用いて粒子輸送計算を実施して患者 3 次元モデル内の線量分布を算出する 2 患者位置計測装置 : 複数のデジタルカメラを用いて照射中の患者を撮影し モーションキャプチャー技術を利用して患者上の複数の特徴点 ( 鼻尖 外眼核 外耳道等 ) や患者に設置したマーカーを認識して これらの点の 3 次元座標 ( ビーム孔の中心を基準座標 :( 0,0,0) とする ) をリアルタイムで計測する装置する 図 1 に患者位置計測装置の概略を示す 3 患者の実際の位置情報 : 患者位置計測装置からは照射室内の患者上 及び 設備上の各点の 3 次元座標が DICOM フォーマットのデータで逐次出力されるものとする 患者位置計測用デジタルカメラ 図 1 患者位置計測装置の概略 患者位置計測装置 (1) ツクバプランへの患者の実際の位置情報の取り込み機能の追加現状のツクバプランに対して 患者の位置情報を取り込む機能を追加する 患者の位置情報データは 患者位置計測装置から出力されるデータであり 患者上の特徴点やマーカー 設備上の特徴点 マーカー点の 3 次元座標情報が記録されているものとする 患者上の特徴点としては目の外眼角 内眼角 鼻根部 鼻尖 外耳孔などである マーカーは 直径約 1 cm程度の球形状のもの ( 材質はプラスチックなど ) を複数設置し 色を赤 青 黄色などとして カメラで撮像した画像上で患者人体とは明らかに異なる部位として認識できる点とする これらの部位を設定できる個数は 位置承認精度を担保するため 指定する部位 マーカーを異なるアングルのカメラから最低 3 点以上認識させる必要があるため 患者モデル上では 6 点以上の点を設定できるようにしておくこと 3 次元座標情報は ビーム孔の表面中心の位置を基準点 :(0,0,0) として 中性子ビームの進行方向を Z 軸座標系 設備の水平軸を X 軸座標系 垂直方向を Y 軸座標系として記録されているものとする また これらの患者位置情報データは DICOM フォーマットで記述されているものとする (2) ツクバプラン上の患者の位置再現機能の追加 (1) で記載した患者の位置情報取り込み機能によってツクバプランに取り込んだ複数の各点の 3 次元座標を基に ツクバプラン上の仮想ビーム孔に対する患者 3D モデルの位置を再現する 患者 3D モデルの位置の再現方法は 患者 3D モデル上の複数点 ( 患者位置計測装置より出力される点 )
の 3 点以上の重心 G 0 を算出し さらにこの重心点 G 0 からある点 :P 00 への方向ベクトル P 0 を算出しておく 次に取り込んだ実際の患者の位置情報での重心点 G 1 及びある点 :P 10 への方向ベクトル :P 1 を算出し 元の患者モデルの重心 G 0 と方向ベクトル P 0 を G 1 と P 1 に移動させることで一致させるものとする 再現した位置での PHITS 計算用インプットを作成し その位置条件での線量分布計算を実行できるようにする 本機能の概略を図 2 に示す ここで 重心を求める点の数が多くなるほど重心位置の算出精度を向上できるアルゴリズムとすること 患者を位置を再現するタイミングは 将来的には適宜 ( 最短で 1 分毎 ) で行えることを考慮して機能制作を行う 本制作では ユーザーが任意のタイミングの位置情報を指定して その位置での再現を実行するものとし 再現処理は 1 分以内で完了してその位置条件での線量計算を速やかに開始できること 患者モデル側の重心の座標を患者位置計測装置で算出した実際の患者の位置の重心に移動 患者モデル側の重心 ツクバプラン ツクバプラン上の患者モデル 実際の患者側のベクトルに一致するように患者モデルのベクトルを移動 図 2 ツクバプラン上のモデルを患者位置計測装置で計測した実際の患者位置に移動させる方法 (3) PHITS Version 3.0.2 への対応線量計算を実施するツクバプラン用モンテカルロ計算エンジン :MC-Tsukuba は 汎用モンテカルロ計算コード :PHITS をベースに構築している 現在の PHITS の最新版は Version3.0.2 であり 種々の新機能追加とバグフィックスが行われている よって MC-Tsukuba を PHITS の Version 3.0.2 をベースに再構築し MC-Tsukuba のバージョンを 1.2 に変更する この MC-Tsukuba の変更に伴い ツクバプランが作成するインプットファイルも MC-Tsukuba Ver.1.2 に対応に関するインプットを出力できるように修正 改良する MC-Tsukuba Ver.1.2 に対応するための記述の変更点は筑波大学側から指示するものとする (4) PHITS 計算用ボクセルモデルの構築法の変更ツクバプランはモンテカルロコード :PHITS をベースにした専用モンテカルロコード :MC- Tsukuba による輸送計算によって線量計算を実行するため 患者の医療画像データを基に製作した 3 次元モデルからボクセルモデルを作成する このボクセルモデルは PHITS の Lattice 機能を用いて まず 2D スライスのボクセルを複数作成する 1 スライスは最大で CT の解像度と一致するため 512 512 個のボクセルで構成される 場合によっては計算速度を優先し 256 256 個のボクセル 128 128 個のボクセルのスライスに変更するものとする この 2D ボクセルスライスを順番に並べて 3D のボクセルモデルを作成している スライス数は 2 枚以上で 100 枚以上のスライスを読み込めることとする この方法に対して 最初の 2D ボクセルスライスは CT データの各スライスに一致した 2D ボクセルモデルを作成するように変更する ( 4) したがって原則として 2D ボクセルスライスは CT のスライス数分作成されることになる そして 3D ボクセ
ルモデルは この CT スライスに対応した 2D ボクセルスライスを並べて 3D モデルを構築することとする Tally 定義も原則としてこの 2D ボクセルスライスに一致するように定義する ( 4) したがって計算結果は 1 スライス当たり最大で 512 512=262,144 個の計算結果が出力されることになり 各 CT スライス毎の線量 2 次元分布を算出できることになる この処理をスライス毎に行うことで人体内の 3 次元の線量分布を算出することとなる 具体的な処理方法や ユーザーインターフェイスの詳細は筑波大学との打合せによって決定するものとする 4 現状のツクバプランでは 2 スライスもしくは 3 スライスを 1 つにまとめた複数スライスでの 2D ボクセルスライスを作成する機能がある 今回の仕様変更においても複数スライスをまとめた 2D ボクセルスライスを作成できるものとし この場合 元の CT スライスデータの複数で 1 つの 2D ボクセルスライスを作成するものとする Tally 定義も同様に 複数の CT スライスにまたがって Tally を定義できるものとし 計算された 2 次元線量分布は 内挿補間処理によって各ピクセル毎に適切な値が割り当てられるものとする (5) 包括的なバグ検索と改修作業潜在しているバグを洗い出し すべてのバグに対する改修作業を実施する 4. 動作検査 3. で示す各機能を追加したシステムを用いた動作試験を実施し 製作物が仕様どおりの機能を有することを確認する 本確認作業は筑波大学において実施する