醫用圖像診斷輔助方法、裝置和圖像處理程序的製作方法

2023-06-02 07:40:41

專利名稱：醫用圖像診斷輔助方法、裝置和圖像處理程序的製作方法
技術領域：
本發明涉及高精度顯示肺氣腫的進展程度的醫用圖像診斷輔助方法和裝置。
本申請以基於日本專利法的專利申請特願第2004-222711號為基礎要求巴黎優先權，由於享受特願第2004-222711號的利益，參照引用其內容。
背景技術：
肺氣腫是一旦發病就無法停止其進展的疾病。為此，向醫生提供可以實現肺氣腫的早期發現、進展程度的鑑別的圖像診斷信息，成為重要的技術。
目前，就用於對肺氣腫的進展程度進行圖像診斷的技術而言，根據閾值從包括整個肺野區域的斷層圖像提取肺野區域和肺氣腫區域，計算肺氣腫區域相對於整個肺野區域的比例，根據該比例的大小改變肺氣腫區域的輪廓的顏色而進行顯示，以容易直觀了解肺氣腫的進展程度的方式進行顯示(例如，專利文件1)。
專利文件1特開2003-10171號公報但是，關於專利文件1中公開的肺氣腫區域的提取方法，提取的區域比原本的肺氣腫區域大數個像素，擔心過度提取。即，在專利文件1中並沒有提高肺氣腫區域的提取精度的考慮。

發明內容
因此，本發明的目的在於，提供可以提高病變候選區域的提取精度的醫用圖像診斷輔助方法、裝置和圖像處理程序。
本發明的醫用圖像診斷輔助方法，包括部位區域提取工序，得到由醫用圖像診斷裝置攝像的斷層圖像，並從該得到的斷層圖像提取規定部位區域；第一區域提取工序，根據由該規定部位區域提取工序提取的部位區域的像素值，從該部位區域提取第一病變候選區域；第二區域提取工序，根據由上述規定部位區域提取工序提取的部位區域的像素值的分布，從該部位區域提取第二病變候選區域；和區域校正工序，使用由該第二區域提取工序提取的第二病變候選區域，對由上述第一區域提取工序提取的第一病變候選區域進行校正。
在這裡，規定部位是指肺，病變候選適合為肺氣腫；部位為腸管等管腔臟器，病變候選可以用於檢測癌細胞，所以分別稱之為「規定部位」、「病變候選」。
本發明的醫用圖像診斷輔助裝置，包括部位區域提取機構，得到由醫用圖像診斷裝置攝像的斷層圖像，並從該得到的斷層圖像提取規定部位區域；第一區域提取機構，根據由該規定部位區域提取機構提取的部位區域的像素值，從該部位區域提取第一病變候選區域；第二區域提取機構，根據由上述規定部位區域提取機構提取的部位區域的像素值的分布，從該部位區域提取第二病變候選區域；和區域校正機構，使用由該第二區域提取機構提取的第二病變候選區域，對由上述第一區域提取機構提取的第一病變候選區域進行校正。
本發明的圖像處理程序，使計算機執行如下所述的工序，即讀入由醫用圖像診斷裝置攝像的斷層圖像的工序；從該得到的斷層圖像提取規定部位區域的部位區域提取工序；根據由該規定部位區域提取工序提取的部位區域的像素值，從該部位區域提取第一病變候選區域的第一區域提取工序；根據由上述規定部位區域提取工序提取的部位區域的像素值的分布，從該部位區域提取第二病變候選區域的第二區域提取工序；使用由該第二區域提取工序提取的第二病變候選區域，對由上述第一區域提取工序提取的第一病變候選區域進行校正的區域校正工序；和顯示上述第一病變候選區域的顯示工序。
通過本發明，可以提高病變候選區域的抽取精度。


圖1是表示本發明的醫用圖像診斷輔助裝置的構成的框圖。
圖2是表示圖1的醫用圖像診斷輔助裝置的動作順序的一例的流程圖。
圖3是表示圖2的步驟21的肺野區域提取處理的流程圖。
圖4是表示圖2的步驟22的氣管/支氣管區域刪除處理的詳細情況的圖。
圖5是表示圖2的步驟23的腸管區域刪除處理的詳細情況的圖。
圖6是表示圖2的步驟24的肺野區域分割處理的詳細情況的圖。
圖7是表示圖2的步驟25的肺氣腫區域分割處理的詳細情況的圖。
圖8是表示圖2的步驟26的分析處理的一部分的局部CT值偏差計算處理的說明圖。
圖9是測量結果顯示畫面的例子。
圖10是與圖9不同的測量結果和分析結果顯示畫面的例子。
圖11是表示基於支氣管的肺野區域的分割例子的模式圖，(a)是表示外側面的圖，(b)是表示內側面的圖，(c)是表示膈面的圖。
圖12是測量結果顯示畫面的例子。
圖13是表示圖象處理程序的構成的框圖。
圖中中央處理裝置(CPU)，11-醫用斷層圖像攝像裝置，12-區域網(LAN)，13-磁碟，14-主存儲器，15-控制器，16-滑鼠，17-鍵盤，18-顯示存儲器，19-顯示器，20-印表機。
具體實施例方式
醫用圖像診斷輔助裝置是由中央處理裝置(CPU)10、醫用斷層圖像攝像裝置11、磁碟13、主存儲器14、滑鼠15和滑鼠控制器16、鍵盤17、顯示存儲器18、CRT顯示器19、電連接印表機20的各構成要素的區域網(LAN)12以及公共總線12a構成。
CPU10對各構成要素的動作進行控制。醫用斷層圖像攝像裝置11是X線CT裝置、核磁共振成像裝置、超聲波診斷裝置、核醫學裝置等可以測量被檢查者的斷層圖像的醫用圖像診斷裝置。LAN12不單是指狹義的網絡，可以與網際網路、電話線路等各種通信網絡連接，與其他計算機或資料庫之間進行圖像數據的交換。磁碟13中儲存各患者的多個斷層圖像數據和動作程序等。主存儲器14中儲存醫用圖像診斷輔助裝置的控制程序。滑鼠15和滑鼠控制器16操作畫面上的軟開關(soft switch)。鍵盤17設定各種參數。顯示存儲器18臨時存儲用於顯示的圖像數據。CRT顯示器19根據來自顯示存儲器18的圖像數據顯示圖像。印表機20例如輸出圖像診斷的結果的報告。
在該實施方式中，作為主存儲器14以外的存儲裝置，只連接有磁碟13，但除此之外，還可以連接FDD、硬碟驅動、CD-ROM驅動、磁光碟(MO)驅動、ZIP驅動、PD驅動、DVD驅動等。
使用

圖1的圖像診斷輔助裝置的動作例子。圖2是表示圖1的醫用圖像診斷輔助裝置所執行的主流程的圖。圖1的CPU10按照該主流程進行動作。下面依照步驟順序詳細說明該主流程。
首先，在醫用圖像診斷輔助裝置的顯示器19上顯示被檢查者的ID輸入畫面，操作者輸入患者的ID編號。於是，從磁碟13，讀出由醫用斷層圖像攝像裝置11預先攝像的斷層圖像之中，與成為診斷對象的患者的ID編號相對應的斷層圖像。
CPU10對上述讀出的斷層圖像子程序調用肺野區域提取處理而執行。圖3是表示步驟21的肺野區域提取處理的子程序的流程圖。下面依照步驟順序詳細說明肺野區域提取處理。
對讀出的斷層圖像進行閾值處理，作成二進位化圖像。在這裡的閾值預先設定計算對象區域明確分離那樣的數值。另外，可以以顯示上述斷層圖像的CT值範圍內最多頻度的CT值為中央值，運算其CT值範圍內的下限值/上限值，將該運算的閾值作為CT值範圍內的下限值/上限值。
關於在步驟S30中作成的二進位化圖像，進行個別識別的標記處理，作成二維標記圖像。
用由步驟S32二維作成的標記圖像，提取最大面積的標記區域作為計算對象區域。
對於在步驟S33中提取的計算對象區域，對連接區域中缺損的部分進行填補處理，作成體區域。
通過在步驟S33中提取的區域和在步驟S34中作成的體區域的判斷辨別處理提取肺野區域。因此，移行至作為主程序的圖2的步驟22。
接著，子程序調用氣管/支氣管區域刪除處理並執行。圖4是表示步驟22的氣管/支氣管區域刪除處理的子程序的流程圖。下面依照步驟順序詳細說明該氣管/支氣管區域刪除處理。
對於在步驟S21中提取的肺野區域內的氣管/支氣管區域，利用其特徵性形狀刪除上述氣管/支氣管區域。作為該處理，在每一張斷層圖像的肺野區域，可以刪除其重心附近孤立的小區域。
刪除在步驟S21中提取的肺野區域和在步驟S41中提取的氣管/支氣管區域中重疊的區域。由此，移行至作為主程序的圖2的步驟23。
接著，子程序調用腸管刪除處理並執行。就該處理而言，在肺野和膈接觸的部分存在腸的影響，所以該處理刪除該影響。作為該處理，以在步驟S21和步驟S22中提取的肺野區域的重心為開始點，採用圖5的處理。圖5是表示步驟23的腸管刪除處理的子程序的流程圖。下面依照步驟順序詳細說明該腸管區域刪除處理。
在由步驟S21和S22提取的肺野區域中，判斷在上下片層(slice)相同位置處有無肺野區域。如果有肺野區域，進行至下一個片層。若無，進行至步驟S52。
刪除上下沒有重疊的區域。由此，移行至作為主程序的圖2的步驟24。
另外，代替圖5，可以採用在基於日本專利法的專利第2845995號公報中公開的區域提取方法。
接著，子程序調用肺野區域分割處理並執行。在本處理中，將肺野分割成右上/中/下肺層面野、右上/下肺野，但在各斷層圖像中，可以分割成肺野內部和外部。圖6是表示步驟24的肺野區域分割處理的子程序的流程圖。下面依照步驟順序詳細說明該肺野區域分割處理。
在由步驟S21～S23提取的肺野區域中，識別各斷層圖像中提取的肺野區域的邊緣部分。
設定以由步驟S61識別的邊緣的像素為重心的圓形區域作為屏蔽(mask)區域，區域半徑的初始值設為2cm，但操作者可以任意設定該值。
進行由步驟S21～S23提取的肺野區域和由步驟S62設定的屏蔽區域的判斷辨別處理。如果肺野區域和屏蔽區域重疊，肺野外部區域，如果不重疊，設定成肺野內部區域。由此，移行至作為主程序的圖2的步驟25。
接著，子程序調用肺氣腫區域檢測處理並執行。圖7是表示步驟25的肺氣腫區域檢測處理的子程序的流程圖。下面依照步驟順序詳細說明該肺氣腫區域檢測處理。
在由步驟S21～S23提取的肺野區域中，進行閾值處理，作成二進位化圖像。這裡的閾值設定成比健康者的肺野內的平均CT值低的-910。另外，可以使用判別分析法，根據其分析結果設定閾值。
在由步驟S71提取的肺氣腫區域中，刪除被看作肺泡的直徑不到3mm的區域。

在由步驟S21～S23提取的肺野區域中，計算出圖8的關心區域81內的CT值分布，根據其分布信息進行閾值處理。
上述CT值的分布可以通過上述規定部位區域和上述病變候選區域的面積比、均值比、標準差比的至少一個計算出來。
就面積比而言，只需計數各區域的像素數就可以得到，所以是最容易求得的值，均值比和標準差比使用區域的像素值，因此與面積比相比可以提高病變候選區域檢測的精度。如果組合面積比、均值比和標準差比，可以進一步提高病變候選區域檢測的精度。
另外，就上述CT值的分布而言，將肺野區域至少2分割，根據該至少2分割的規定部位區域和上述肺氣腫區域的面積比、均值比、標準差比的至少一個，可以計算出來。由此，根據分割區域的局部像素值的分布，求得面積比、均值比和標準差比，所以與上述部位區域整體相比可以實現極為細小的病變候選區域檢測的精度的調整。
僅提取由步驟S71/S72提取的區域和由步驟S73提取的區域相重疊的區域。由此，移行至作為主程序的圖2的步驟26。
接著，子程序調用分析處理並執行。根據由步驟S21～S23提取的肺野區域、由步驟S24分割的肺野區域和由步驟S25提取的肺氣腫區域的信息，以各斷層圖像單位，以各對應值成對的形式計算出肺野區域的面積、均值、標準差、均值的95％可靠區間的上限值、下限值，和肺氣腫區域的面積、比例、均值、標準差、均值的95％可靠區間的上限值、下限值。另外，如此計算出各全肺野/各肺野單位的分析結果。在這裡，將均值的可靠區域初始設定為95％，但操作者可以任意設定。
接著，執行測量結果顯示處理。在步驟S26中，以將全肺野作為對象的分析結果為基礎，切換顯示肺氣腫區域的顯示顏色。圖9表示此時的畫面顯示例。畫面91顯示讀入的斷層圖像。畫面92表示測量結果的顯示圖像。另外，在圖9中，顯示畫面93上顯示拍攝上述斷層圖像時的攝像條件，但顯示畫面93不是必須的，可以僅顯示畫面91和92。在這裡，提示了用2個畫面顯示畫面91和92的例子，但可以用1個畫面顯示斷層圖像和肺氣腫區域重疊的圖像。
另外，在圖10中顯示表示分析結果的畫面顯示例子。畫面101表示讀入的斷層圖像，畫面102表示測量結果的顯示圖像，以及畫面103表示分析結果。
圖10的畫面103顯示每個肺氣腫區域的實測面積值「30000cm2、1000cm2」、佔全肺野區域的比例「3.33％」，作為分析結果。由此，可以明示精度好的病變候選區域，並提供給醫生。
另外，上述顯示可以重疊顯示上述斷層圖像和上述病變候選區域。
由此，可以直觀診斷斷層圖像和病變候選區域的位置關係。
另外，上述顯示可以是將上述病變候選區域與上述斷層圖像鄰接而顯示。
由此，可以分析通過上述重疊顯示而隱藏起來的斷層圖像的像素值。
另外，上述顯示可以將表示上述病變候選區域的顯示/不顯示的軟開關顯示於顯示畫面，能夠通過指示設備進行選擇。由此，選擇是僅顯示斷層圖像來進行圖像診斷，還是顯示斷層圖像和病變候選區域兩者。
另外，可以列表顯示如圖10的103所示的各斷層圖像的分析結果。此外，向印表機20輸出該列表顯示的結果，並提供給患者，可以有助於戒菸指導等。另外，當用滑鼠16指定病變候選區域時，可以顯示彈出窗口，在該彈出窗口中顯示病變候選區域的特徵量和位置信息。根據圖11和12說明該顯示例子。
圖11顯示根據支氣管分割肺野區域的例子。左下斜線表示上葉區域，橫線表示中葉區域，右下斜線表示下葉區域。另外，「S1」、「S2」…是識別主要支氣管的識別編號。「S1a」、「S2b」…是表示以主要支氣管為中心的區域。圖12是顯示畫面例子，畫面121表示斷層圖像，畫面122表示測量結果的顯示圖像，畫面123表示彈出窗口。
醫用圖像診斷輔助裝置預先在磁碟23中儲存以圖11所示的支氣管為中心的區域信息。此外，根據由S21提取的肺區域和以圖11所示的支氣管為中心的區域信息，判斷肺氣腫區域位於以支氣管為中心的區域「S1a」、「S2b」…的哪一部位，記錄該結果作為肺氣腫區域的位置信息。其中，肺氣腫區域的位置信息不限於以支氣管為中心的區域，還可以是顯示屬於上葉、中葉、下葉區域的哪個區域的信息。
操作者用滑鼠16點擊顯示在顯示器19上的肺氣腫區域，此時，顯示圖12的彈出窗口123。在彈出窗口123顯示用滑鼠16點擊的肺氣腫區域的特徵量「30000cm2、1000cm2、3.33％」、和肺氣腫的位置「S9b區域」。在彈出窗口123顯示的信息可以是肺氣腫區域的特徵量、位置信息的任一方，還可以顯示其他信息。
由此，選擇是僅顯示斷層圖像進行圖像診斷，還是顯示斷層圖像和病變候選區域兩者。
在本實施方式中，包括得到由醫用斷層圖像攝像裝置11攝像的斷層圖像並從該得到的斷層圖像提取規定部位區域的工序21、根據由該規定部位區域提取工序提取的部位區域的像素值從該部位區域提取第一病變候選區域的工序22～24、根據上述所提取的部位區域的像素值的分布從該部位區域提取第二病變候選區域的工序25、使用上述所提取的第二病變候選區域並對由上述第一區域提取工序所提取的第一病變候選區域進行校正的工序26。
由此，可以校正由被檢查者的體型等個體差異等主要原因產生的病變候選區域的提取誤差。
圖13是表示圖像處理程序的框圖。
圖像處理程序包括從醫用斷層圖像攝像裝置11或磁碟、或藉助LAN12連接的圖像資料庫，讀入斷層圖像的讀入部10a；從所讀入的斷層圖像提取規定部位(例如肺區域)的部位區域提取部10b；根據該規定部位區域的像素值從該規定部位區域提取第一病變候選區域的第一區域提取部10c；根據該規定部位區域的像素值從該規定部位區域提取第二病變候選區域的第二區域提取部10d；使用第二病變候選區域校正第一病變候選區域的區域校正部10e；以及顯示上述第一病變候選區域的顯示控制部10f。顯示控制部10f按照來自滑鼠15或鍵盤17的輸入指示，進行顯示圖10的分析結果的畫面103或圖12的彈出窗口123的顯示控制。醫用圖像診斷輔助裝置在磁碟23中儲存圖13所示的圖像處理程序，CPU10適當加載在主存儲器14上來執行。由此。執行上述實施方式的動作處理。
在上述實施方式中，讀入斷層圖像進行各種處理，但不限於斷層圖像，還可以根據由X線診斷裝置或DR裝置拍攝的投影圖像進行同樣的處理。
工業上的可利用性本發明並不限於在基於醫用圖像的圖像處理中使用，還可以用於從讀入的圖像提取第一像素值區域和第二像素值區域並根據第二像素值區域校正第一像素值區域的圖像處理常用技術。
權利要求
1.一種醫用圖像診斷輔助方法，其特徵在於，包括部位區域提取工序，得到由醫用圖像診斷裝置攝像的斷層圖像，並從該得到的斷層圖像提取規定部位區域；第一區域提取工序，根據由該規定部位區域提取工序提取的部位區域的像素值，從該部位區域提取第一病變候選區域；第二區域提取工序，根據由所述規定部位區域提取工序提取的部位區域的像素值的分布，從該部位區域提取第二病變候選區域；和區域校正工序，使用由該第二區域提取工序提取的第二病變候選區域，對由所述第一區域提取工序提取的第一病變候選區域進行校正。
2.如權利要求1所述的醫用圖像診斷輔助方法，其特徵在於，所述規定部位區域的像素值的分布包括根據所述規定部位區域和所述病變候補區域的面積比、均值比、標準差比的至少一個計算出來。
3.如權利要求1所述的醫用圖像診斷輔助方法，其特徵在於，所述規定部位區域的像素值的分布包括將所述提取的規定部位區域至少2分割，根據該至少2分割的規定部位區域和所述病變候補區域的面積比、均值比、標準差比的至少一個計算出來。
4.如權利要求3所述的醫用圖像診斷輔助方法，其特徵在於，所述規定部位區域的分割包括分割成與所述規定部位區域相對應的分割數量。
5.如權利要求1～4中任意一項所述的醫用圖像診斷輔助方法，其特徵在於，還包括進展程度計算工序，計算出通過所述區域校正工序校正的病變候選區域的進展程度；和顯示工序，根據由該進展程度計算出的工序計算出來的病變候選區域的進展程度，改變所述校正的病變候選區域的顏色，使該顏色變化的病變候選區域與所述斷層圖像對應而顯示。
6.如權利要求5所述的醫用圖像診斷輔助方法，其特徵在於，所述顯示工序重疊所述斷層圖像和所述病變候選區域來顯示。
7.如權利要求5所述的醫用圖像診斷輔助方法，其特徵在於，所述顯示工序將所述病變候選區域鄰接於所述斷層圖像來顯示。
8.如權利要求5所述的醫用圖像診斷輔助方法，其特徵在於，所述顯示工序包括選擇所述病變候選區域的顯示/不顯示的選擇工序。
9.一種醫用圖像診斷輔助裝置，其特徵在於包括部位區域提取機構，得到由醫用圖像診斷裝置攝像的斷層圖像，並從該得到的斷層圖像提取規定部位區域；第一區域提取機構，根據由該規定部位區域提取機構提取的部位區域的像素值，從該部位區域提取第一病變候選區域；第二區域提取機構，根據由所述規定部位區域提取機構提取的部位區域的像素值的分布，從該部位區域提取第二病變候選區域；和區域校正機構，使用由該第二區域提取機構提取的第二病變候選區域，對由所述第一區域提取機構提取的第一病變候選區域進行校正。
10.如權利要求9所述的醫用圖像診斷輔助裝置，其特徵在於，所述規定部位區域的像素值的分布包括根據所述規定部位區域和所述病變候補區域的面積比、均值比、標準差比的至少一個計算出來。
11.如權利要求9所述的醫用圖像診斷輔助裝置，其特徵在於，所述規定部位區域的像素值的分布包括將所述提取的規定部位區域至少2分割，根據該至少2分割的規定部位區域和所述病變候補區域的面積比、均值比、標準差比的至少一個計算出來。
12.如權利要求11所述的醫用圖像診斷輔助裝置，其特徵在於，所述規定部位區域的分割包括分割成與所述規定部位區域相對應的分割數量。
13.如權利要求9～12中任意一項所述的醫用圖像診斷輔助裝置，其特徵在於，還包括進展程度計算機構，計算出通過所述區域校正機構校正的病變候選區域的進展程度；和顯示機構，根據由該進展程度計算機構計算出來的病變候選區域的進展程度，改變所述校正的病變候選區域的顏色，使該顏色變化的病變候選區域與所述斷層圖像對應而顯示。
14.如權利要求13所述的醫用圖像診斷輔助裝置，其特徵在於，所述顯示機構重疊所述斷層圖像和所述病變候選區域來顯示。
15.如權利要求13所述的醫用圖像診斷輔助裝置，其特徵在於，所述顯示機構將所述病變候選區域鄰接於所述斷層圖像來顯示。
16.如權利要求13所述的醫用圖像診斷輔助裝置，其特徵在於，所述顯示機構包括選擇所述病變候選區域的顯示/不顯示的選擇機構。
17.如權利要求13所述的醫用圖像診斷輔助裝置，其特徵在於，還包括用於指定所述顯示機構顯示的畫面上的位置坐標的位置坐標輸入機構，當由所述位置坐標輸入機構指定所述病變候選區域時，所述顯示機構顯示與已經顯示的畫面不同的其他畫面，在該其他畫面中顯示所述病變候選區域的特徵量、以及所述病變候選區域的位置信息的至少一個。
18.如權利要求9～17中任意一項所述的醫用圖像診斷輔助裝置，其特徵在於，所述斷層圖像是對被檢查者的肺野區域進行攝像而得到的圖像，所述病變候選區域是所述肺野區域中的肺氣腫區域。
19.如權利要求18所述的醫用圖像診斷輔助裝置，其特徵在於，所述第一區域提取機構提取肺氣腫區域，所述第二提取機構提取肺泡區域，所述區域校正機構進行從所述肺氣腫區域除去所述肺泡區域的校正。
20.一種圖像處理程序，其特徵在於，使計算機執行下述的工序讀入由醫用圖像診斷裝置攝像的斷層圖像的工序；從該得到的斷層圖像提取規定部位區域的部位區域提取工序；根據由該規定部位區域提取工序提取的部位區域的像素值，從該部位區域提取第一病變候選區域的第一區域提取工序；根據由所述規定部位區域提取工序提取的部位區域的像素值的分布，從該部位區域提取第二病變候選區域的第二區域提取工序；使用由該第二區域提取工序提取的第二病變候選區域，對由所述第一區域提取工序提取的第一病變候選區域進行校正的區域校正工序；和顯示所述第一病變候選區域的顯示工序。
全文摘要
本發明提供一種醫用圖像診斷輔助方法，包括部位區域提取工序，得到由醫用圖像診斷裝置攝像的斷層圖像，並從該得到的斷層圖像提取規定部位區域；第一區域提取工序，根據由該規定部位區域提取工序提取的部位區域的像素值，從該部位區域提取第一病變候選區域；第二區域提取工序，根據由所述規定部位區域提取工序提取的部位區域的像素值的分布，從該部位區域提取第二病變候選區域；和區域校正工序，使用由該第二區域提取工序提取的第二病變候選區域，對由所述第一區域提取工序提取的第一病變候選區域進行校正。
文檔編號G06T1/00GK1988850SQ20058002475
公開日2007年6月27日 申請日期2005年7月28日 優先權日2004年7月30日
發明者中島邦佳, 後藤良洋 申請人:株式會社日立醫藥