JP5643363B2 - Medical image display device and medical image display program - Google Patents
Medical image display device and medical image display program Download PDFInfo
- Publication number
- JP5643363B2 JP5643363B2 JP2013054649A JP2013054649A JP5643363B2 JP 5643363 B2 JP5643363 B2 JP 5643363B2 JP 2013054649 A JP2013054649 A JP 2013054649A JP 2013054649 A JP2013054649 A JP 2013054649A JP 5643363 B2 JP5643363 B2 JP 5643363B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- pixels
- display
- rows
- column direction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Liquid Crystal Display Device Control (AREA)
- Transforming Electric Information Into Light Information (AREA)
- Control Of Indicators Other Than Cathode Ray Tubes (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Description
本発明は、医療画像表示装置および医療画像表示プログラムに関し、更に詳しくは、医師がX線撮影画像,CT画像,MR画像などの医療画像を読影診断する際、その診断の効率化を図ることが出来る医療画像表示装置および医療画像表示プログラムに関する。 The present invention relates to a medical image display device and a medical image display program, and more particularly, to improve the efficiency of diagnosis when a doctor interprets a medical image such as an X-ray image, CT image, or MR image. The present invention relates to a medical image display device and a medical image display program.
従来、画素数の多い画像を、単純サンプリングまたは演算サンプリングにより縮退させて、画素数の少ない画像を生成する画像の縮退方法が知られている(特許文献1参照)。
他方、LCDの1つのピクセルを構成する3個のサブピクセルを独立に駆動するディスプレイ装置が知られている(特許文献2参照)。
On the other hand, a display device that independently drives three sub-pixels constituting one pixel of an LCD is known (see Patent Document 2).
例えば、画素数4096×3328のFPD(Flat Panel Detctor)を用いて撮像した医療画像を画素数1600×1200のLCDで表示したい場合がある。
この場合、特許文献1の画像の縮退方法を用いて、FPDで得た画素数4096×3328の医療画像を画素数1477×1200の画像に縮退させれば、画素数1600×1200のLCDで表示させることが可能である。なお、縮退した画像の短軸画素数が1200であるのはLCDの短軸画素数に合わせたためであり、縮退した画像の長軸画素数が1477であるのは原画像のアスペクト比を保存したためである。
しかし、原画像の短軸方向および長軸方向についてそれぞれ2.774(=4096/1477=3328/1200)個の画素がLCDの1個の画素に縮退させられているため、元の解像度が大きく損なわれる問題点がある。
具体例で説明すると、乳ガン診断においては、100μm以下の微小石灰化のような微小な病変を判別することが重要であり、FPDの画素ピッチは70μmであり、必要な解像度を有している。しかし、LCDでは、FPDの画素ピッチに換算した画素ピッチが194(=70*2.774)μmになるため、必要な解像度が得られなくなってしまう問題点がある。
For example, there is a case where a medical image captured using an FPD (Flat Panel Detctor) having 4096 × 3328 pixels is desired to be displayed on an LCD having 1600 × 1200 pixels.
In this case, if the medical image of 4096 × 3328 obtained by FPD is reduced to an image of 1477 × 1200 using the image reduction method of
However, since 2.774 (= 4096/1477 = 3328/1200) pixels are reduced to one pixel of the LCD in each of the short axis direction and the long axis direction of the original image, the original resolution is large. There is a problem that is damaged.
Specifically, in breast cancer diagnosis, it is important to discriminate minute lesions such as microcalcifications of 100 μm or less, and the pixel pitch of the FPD is 70 μm, which has a necessary resolution. However, in the LCD, since the pixel pitch converted to the pixel pitch of the FPD is 194 (= 70 * 2.774) μm, there is a problem that a necessary resolution cannot be obtained.
ここで、モノクロ液晶表示装置のLCDの1つのピクセルは、LCDの長軸方向に並ぶ3個のサブピクセルから構成されている。そこで、特許文献2のディスプレイ装置のようにLCDのサブピクセルを独立に駆動すれば、LCDの長軸方向の解像度を向上することが出来る。
しかし、サブピクセルの形状は長方形であり、正方形または円形形状のピクセルを前提としている特許文献1の画像の縮退方法を適用して表示すると、原画像のアスペクト比を保存できなくなる問題点がある。
Here, one pixel of the LCD of the monochrome liquid crystal display device is composed of three sub-pixels arranged in the major axis direction of the LCD. Therefore, if the LCD subpixels are driven independently as in the display device of
However, the shape of the subpixel is a rectangle, and there is a problem that the aspect ratio of the original image cannot be preserved when the image reduction method disclosed in
そこで、本発明の目的は、医師がX線撮影画像,CT画像,MR画像などの医療画像を読影診断する際、その診断の効率化を図ることが出来る医療画像表示装置および医療画像表示プログラムを提供することにある。 Accordingly, an object of the present invention is to provide a medical image display device and a medical image display program capable of improving the efficiency of diagnosis when a doctor interprets a medical image such as an X-ray image, CT image, or MR image. It is to provide.
第1の観点では、本発明は、m、n、qを自然数とし、m>qとするとき、モノクロ原医療画像におけるm行n列の原画像をq行(3*n*q/m)列の画素からなる表示画像に変換する画像変換手段と、モノクロLCDの行方向に3個並んで正方形の各画素を構成する長方形のサブピクセルを独立に駆動して前記表示画像を表示するモノクロ液晶表示装置とを具備した医療画像表示装置であって、前記画像変換手段が、前記原画像の全ての画素が表示画像に反映されるように前記原画像を列方向に加重平均演算してm行からq行に減少すると共に前記原画像のアスペクト比が維持されるように行方向に補間演算してn列から(3*n*q/m)列にすることを特徴とする医療画像表示装置を提供する。
上記第1の観点による医療画像表示装置では、m行n列の原画像を縮退してq行(3*n*q/m)列の画素からなる表示画像に変換する。原画像のアスペクト比(n/m)に応じた行方向の画素数は(n*q/m)であるから、表示画像の行方向の画素数は原画像のアスペクト比(n/m)に応じた画素数の3倍になっている。ところが、モノクロLCDの各画素では3個のサブピクセルが行方向に並んでおり、それらサブピクセルを独立に駆動して表示画像を表示するため、表示画像を表示するLCDの行方向の画素数は、(3*n*q/m)/3=(n*q/m)となる。よって、表示画像のアスペクト比は(n*q/m)/q=(n/m)となる。つまり、原画像のアスペクト比が保存されている。
そして、サブピクセルを独立に駆動しない場合に比べて行方向の解像度を3倍にすることが出来る。FPDの画素ピッチに換算したLCDのサブピクセル・ピッチは、FPDの画素ピッチのn/(3*n*q/m)=m/(3*q)倍になるので、例えばFPDの画素ピッチが70μm、m=3328、q=1200なら、FPDの画素ピッチに換算したLCDのサブピクセル・ピッチは65(=70*3328/(3*1200))μmになる。これにより、乳ガン診断において必要な解像度100μmに十分対応できることとなる。すなわち、医師がX線撮影画像,CT画像,MR画像などの医療画像を読影診断する際、その診断の効率化を図ることが出来る。
In the first aspect, the present invention relates to an original image of m rows and n columns in a monochrome original medical image as q rows (3 * n * q / m) where m, n, and q are natural numbers and m> q. Image conversion means for converting into a display image composed of pixels in a column, and monochrome liquid crystal for displaying the display image by independently driving three rectangular sub-pixels constituting each square pixel arranged in the row direction of the monochrome LCD A medical image display device comprising a display device, wherein the image conversion means performs a weighted average operation on the original image in a column direction so that all pixels of the original image are reflected in the display image, and m rows Medical image display apparatus characterized in that the number of columns is changed from n columns to (3 * n * q / m) by performing interpolation calculation in the row direction so that the aspect ratio of the original image is maintained while decreasing from n to q rows I will provide a.
In the medical image display apparatus according to the first aspect, the original image of m rows and n columns is degenerated and converted into a display image composed of pixels of q rows (3 * n * q / m) columns. Since the number of pixels in the row direction corresponding to the aspect ratio (n / m) of the original image is (n * q / m), the number of pixels in the row direction of the display image is equal to the aspect ratio (n / m) of the original image. The number of corresponding pixels is three times. However, in each pixel of the monochrome LCD, three subpixels are arranged in the row direction, and the subpixels are driven independently to display a display image. Therefore, the number of pixels in the row direction of the LCD that displays the display image is (3 * n * q / m) / 3 = (n * q / m). Therefore, the aspect ratio of the display image is (n * q / m) / q = (n / m). That is, the aspect ratio of the original image is stored.
Then, the resolution in the row direction can be tripled compared to the case where the subpixels are not driven independently. Since the LCD sub-pixel pitch converted to the FPD pixel pitch is n / (3 * n * q / m) = m / (3 * q) times the FPD pixel pitch, for example, the FPD pixel pitch is If 70 μm, m = 3328, and q = 1200, the subpixel pitch of the LCD converted to the pixel pitch of FPD is 65 (= 70 * 3328 / (3 * 1200)) μm. As a result, the resolution of 100 μm necessary for breast cancer diagnosis can be sufficiently handled. That is, when a doctor interprets a medical image such as an X-ray image, a CT image, or an MR image, the diagnosis can be made more efficient.
第2の観点では、本発明は、m、n、qを自然数とし、m>qとするとき、モノクロ原医療画像におけるm行n列の原画像をモノクロLCDの行方向に3個並んで正方形の各画素を構成する長方形のサブピクセルを独立に駆動してモノクロ液晶表示装置に表示するために、コンピュータを、前記原画像をq行(3*n*q/m)列の画素からなる表示画像に変換する画像変換手段、及び前記モノクロLCDのサブピクセルを独立に駆動して前記表示画像を表示させるべく前記表示画像を前記モノクロ液晶表示装置に出力する表示画像出力手段、として機能させるための医療画像表示プログラムであって、前記画像変換手段が、前記原画像の全ての画素が表示画像に反映されるように前記原画像を列方向に加重平均演算してm行からq行に減少すると共に前記原画像のアスペクト比が維持されるように行方向に補間演算してn列から(3*n*q/m)列にするステップを有することを特徴とする医療画像表示プログラムを提供する。
上記第2の観点による医療画像表示プログラムでは、m行n列の原画像を縮退してq行(3*n*q/m)列の画素からなる表示画像に変換する。原画像のアスペクト比(n/m)に応じた画素数は(n*q/m)であるから、表示画像の行方向の画素数は原画像のアスペクト比(n/m)に応じた画素数の3倍になっている。ところが、モノクロLCDの各画素では3個のサブピクセルが行方向に並んでおり、それらサブピクセルを独立に駆動して表示画像を表示するため、表示画像を表示する行方向の画素数は、(3*n*q/m)/3=(n*q/m)となる。よって、表示画像のアスペクト比は(n*q/m)/q=(n/m)となる。つまり、原画像のアスペクト比が保存されている。
そして、サブピクセルを独立に駆動しない場合に比べて行方向の解像度を3倍にすることが出来る。FPDの画素ピッチに換算したLCDのサブピクセル・ピッチは、FPDの画素ピッチのn/(3*n*q/m)=m/(3*q)倍になるので、例えばFPDの画素ピッチが70μm、m=3328、q=1200なら、FPDの画素ピッチに換算したLCDのサブピクセル・ピッチは65(=70*3328/(3*1200))μmになる。これにより、乳ガン診断において必要な解像度100μmに十分対応できることとなる。すなわち、医師がX線撮影画像,CT画像,MR画像などの医療画像を読影診断する際、その診断の効率化を図ることが出来る。
In a second aspect, the present invention provides a square in which three m rows and n columns of original images in a monochrome original medical image are arranged in the row direction of the monochrome LCD when m, n, and q are natural numbers and m> q. In order to independently drive the rectangular sub-pixels constituting each of the pixels and display them on the monochrome liquid crystal display device, the computer displays the original image consisting of pixels of q rows (3 * n * q / m) columns. Image conversion means for converting into an image, and display image output means for outputting the display image to the monochrome liquid crystal display device so as to display the display image by independently driving subpixels of the monochrome LCD A medical image display program, wherein the image conversion means reduces the original image from m rows to q rows by performing a weighted average operation in the column direction so that all pixels of the original image are reflected in the display image. And a medical image display program comprising a step of interpolating in the row direction so as to maintain the aspect ratio of the original image and changing from n columns to (3 * n * q / m) columns To do.
In the medical image display program according to the second aspect, the original image of m rows and n columns is degenerated and converted into a display image composed of pixels of q rows (3 * n * q / m) columns. Since the number of pixels corresponding to the aspect ratio (n / m) of the original image is (n * q / m), the number of pixels in the row direction of the display image is the pixel corresponding to the aspect ratio (n / m) of the original image. Three times the number. However, in each pixel of a monochrome LCD, three subpixels are arranged in the row direction, and these subpixels are driven independently to display a display image. Therefore, the number of pixels in the row direction for displaying the display image is ( 3 * n * q / m) / 3 = (n * q / m). Therefore, the aspect ratio of the display image is (n * q / m) / q = (n / m). That is, the aspect ratio of the original image is stored.
Then, the resolution in the row direction can be tripled compared to the case where the subpixels are not driven independently. Since the LCD sub-pixel pitch converted to the FPD pixel pitch is n / (3 * n * q / m) = m / (3 * q) times the FPD pixel pitch, for example, the FPD pixel pitch is If 70 μm, m = 3328, and q = 1200, the subpixel pitch of the LCD converted to the pixel pitch of FPD is 65 (= 70 * 3328 / (3 * 1200)) μm. As a result, the resolution of 100 μm necessary for breast cancer diagnosis can be sufficiently handled. That is, when a doctor interprets a medical image such as an X-ray image, a CT image, or an MR image, the diagnosis can be made more efficient.
本発明の画像表示装置および画像表示プログラムによれば、LCDの長軸方向の解像度を向上させることが出来ると共に原画像のアスペクト比を保存して縮退した画像を表示することが出来るので、医師がX線撮影画像,CT画像,MR画像などの医療画像を読影診断する際、その診断の効率化を図ることが出来る。 According to the image display device and the image display program of the present invention, the resolution in the long axis direction of the LCD can be improved, and the reduced image can be displayed by storing the aspect ratio of the original image. When a medical image such as an X-ray image, a CT image, or an MR image is subjected to interpretation diagnosis, the efficiency of the diagnosis can be improved.
以下、図に示す実施の形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
−実施例1−
図1は、実施例1に係る医療画像表示装置1を示す構成図である。
この医療画像表示装置1は、m、n、q、Pを自然数とし、m≦n、q<P、n/m<P/q、m>qとするとき、m行n列の画素からなるモノクロ医療画像データA(i,j)を入力し縮退してq行(3*n*q/m)列の画素からなる表示画像データM(k,L)に変換し出力する医療画像ビューワープログラム31を実行するパソコン3と、q行P列の画素からなるモノクロLCDを有すると共にモノクロLCDの行方向に並ぶサブピクセルを独立に駆動して表示画像データM(k,L)を表示するモノクロ液晶表示装置4とを具備してなる。
モノクロ医療画像データA(i,j)は、医療画像撮影装置Fのm行n列の画素からなるFPDで撮影され、医療画像データサーバDに記憶されたものである。
Hereinafter, the present invention will be described in more detail with reference to embodiments shown in the drawings. Note that the present invention is not limited thereby.
Example 1
FIG. 1 is a configuration diagram illustrating a medical
The medical
The monochrome medical image data A (i, j) is captured by an FPD composed of pixels of m rows and n columns of the medical image capturing device F and stored in the medical image data server D.
なお、m>nの場合、医療画像データA(i,j)の長軸方向が被検体の上下方向になることが多いため、LCDの長軸方向を上下方向にすることが多い。しかし、説明の都合上、医療画像データA(i,j)もLCDも横長の状態を想定する。 When m> n, the major axis direction of the medical image data A (i, j) is often the vertical direction of the subject, so the major axis direction of the LCD is often the vertical direction. However, for convenience of explanation, it is assumed that the medical image data A (i, j) and the LCD are horizontally long.
図2の(a)は、モノクロ医療画像データA(i,j)の数値例であり、m=3328、n=4096である。
図2の(b)は、LCDの数値例であり、q=1200、P=1600である。
FIG. 2A is a numerical example of monochrome medical image data A (i, j), where m = 3328 and n = 4096.
(B) of FIG. 2 is a numerical example of LCD, and is q = 1200 and P = 1600.
図3は、LCDと表示画像データM(k,L)の関係を示す説明図である。
n/m<P/qの関係があるため、LCDの短軸方向に表示画像データM(k,L)の短軸方向を一致させるようにモノクロ医療画像データA(i,j)を縮退させる。モノクロ医療画像データA(i,j)のアスペクト比n/mを保持すると、LCDの長軸方向の画素数Pより表示画像データM(k,L)の長軸方向の画素数pは短くなる。
例えば、m=3328、n=4096、q=1200、P=1600なら、p=4096*1200/3328=1477になる。
FIG. 3 is an explanatory diagram showing the relationship between the LCD and the display image data M (k, L).
Since there is a relationship of n / m <P / q, the monochrome medical image data A (i, j) is degenerated so that the minor axis direction of the display image data M (k, L) matches the minor axis direction of the LCD. . When the aspect ratio n / m of the monochrome medical image data A (i, j) is maintained, the number of pixels p in the major axis direction of the display image data M (k, L) is shorter than the number of pixels P in the major axis direction of the LCD. .
For example, if m = 3328, n = 4096, q = 1200, and P = 1600, then p = 4096 * 1200/3328 = 1477.
図4に示すように、LCDの1つの画素Uは、LCDの長軸方向に並ぶ3個のサブピクセルu1,u2,u3から構成されている。画素Uは正方形なので、サブピクセルu1,u2,u3の形状は長方形になる。 As shown in FIG. 4, one pixel U of the LCD is composed of three subpixels u1, u2, u3 arranged in the major axis direction of the LCD. Since the pixel U is a square, the shape of the subpixels u1, u2, and u3 is a rectangle.
図5,図6は、パソコン3の医療画像ビューワープログラム31による医療画像ビューワー処理の手順を示すフロー図である。
ステップS0では、m行n列の画素からなるモノクロ医療画像データA(i,j)を医療画像データサーバ2から読み込む。
ステップS1では、モノクロ医療画像データA(i,j)の短軸方向の画素数mとLCDの短軸方向の画素数qから短軸縮退率ry=m/qを求める。例えばm=3328、q=1200ならry=2.774になる。
FIG. 5 and FIG. 6 are flowcharts showing the procedure of the medical image viewer process by the medical
In
In step S1, the short axis degeneration rate ry = m / q is obtained from the number m of pixels in the short axis direction of the monochrome medical image data A (i, j) and the number of pixels q in the short axis direction of the LCD. For example, if m = 3328 and q = 1200, then ry = 2.774.
ステップS2では、短軸縮退率ryと調整値Bから短軸加重平均重みb=(ry−1)*Bを求める。例えばry=2.774、B=0.12ならb=0.2129になる。操作者が調整値Bを調整することにより、表示画像の見え方(画像ノイズなど)を調整することが出来る。 In step S2, a short axis weighted average weight b = (ry−1) * B is obtained from the short axis degeneration rate ry and the adjustment value B. For example, if ry = 2.774 and B = 0.12, then b = 0.129. By adjusting the adjustment value B, the operator can adjust the appearance of the display image (image noise or the like).
ステップS3では、医療画像データA(i,j)の全ての画素が表示画像に反映されるように、医療画像データA(i,j)と短軸加重平均重みbを用いて医療画像データA(i,j)を短軸方向に加重平均した短軸加重平均画像データA’(i,j)を求める。
図7に、医療画像データA(i,j)と短軸加重平均画像データA’(i,j)の概念を示す。
図8及び次式に示すように、隣接する3個の医療画像データA(i,j-1),A(i,j),A(i,j+1)を加重平均して短軸加重平均画像データA’(i,j)を求めているが、これは想定している短軸縮退率ry=2.774より大きい最小の整数が「3」だからである。
A’(i,j)=A(i,j-1)*b+A(i,j)*(1-2*b)+A(i,j+1)*b
但し、A’(i,0)=A(i,0)、A’(i,m-1)=A(i,m-1)とする。
i=0,1,…,n-1、j=0,1,…,m-1
一般的には、短軸縮退率ryより大きい整数を「2*g+1」とするとき、隣接する医療画像データA(i,j-g),…,A(i,j+g)を加重平均して短軸加重平均画像データA’(i,j)を求めればよい。短軸加重平均重みも「g」に応じて適宜調整すればよい。
In step S3, the medical image data A (i, j) and the short axis weighted average weight b are used so that all the pixels of the medical image data A (i, j) are reflected in the display image. Short-axis weighted average image data A ′ (i, j) is obtained by weighted average of (i, j) in the short-axis direction.
FIG. 7 shows the concept of the medical image data A (i, j) and the short axis weighted average image data A ′ (i, j).
As shown in FIG. 8 and the following equation, the three adjacent medical image data A (i, j-1), A (i, j), and A (i, j + 1) are weighted and short axis weighted. Average image data A ′ (i, j) is obtained because the smallest integer larger than the assumed short axis degeneration rate ry = 2.774 is “3”.
A '(i, j) = A (i, j-1) * b + A (i, j) * (1-2 * b) + A (i, j + 1) * b
However, it is assumed that A ′ (i, 0) = A (i, 0) and A ′ (i, m−1) = A (i, m−1).
i = 0,1, ..., n-1, j = 0,1, ..., m-1
In general, when an integer larger than the minor axis degeneration rate ry is “2 * g + 1”, the adjacent medical image data A (i, jg),..., A (i, j + g) are weighted averaged. The short axis weighted average image data A ′ (i, j) may be obtained. The short axis weighted average weight may be appropriately adjusted according to “g”.
ステップS4では、画像列座標i=0に初期化する。 In step S4, the image sequence coordinate i is initialized to 0.
ステップS5では、LCD短軸座標L=0に初期化する。 In step S5, the LCD minor axis coordinate L = 0 is initialized.
ステップS6では、LCD短軸座標Lに対応する画像基準短軸座標jq=int(L*ry)を求める。
ステップS7では、短軸補間重みWy(L)=(L*ry)−jqを求める。
In step S6, an image reference short axis coordinate jq = int (L * ry) corresponding to the LCD short axis coordinate L is obtained.
In step S7, the short axis interpolation weight Wy (L) = (L * ry) −jq is obtained.
ステップS8では、図9および次式に示すように、短軸加重平均画像データA’(i,j)と短軸補間重みWy(L)を用いて線形補間により短軸補間中間画像データC(i,L)を求める。
C(i,L)=A’(i,jq)*(1−Wy(L))+A’(i,jq+1)*Wy(L)
In step S8, as shown in FIG. 9 and the following equation, the short-axis interpolated intermediate image data C () is obtained by linear interpolation using the short-axis weighted average image data A ′ (i, j) and the short-axis interpolation weight Wy (L). i, L).
C (i, L) = A '(i, jq) * (1-Wy (L)) + A' (i, jq + 1) * Wy (L)
ステップS9及びS10では、L=q−1までステップS6〜S8を繰り返す。 In steps S9 and S10, steps S6 to S8 are repeated until L = q-1.
ステップS11及びS12では、i=n−1までステップS5〜S10を繰り返す。
以上により、図10に示すように、短軸加重平均画像データA’(i,j)が短軸補間中間画像データC(i,L)に変換される。
In steps S11 and S12, steps S5 to S10 are repeated until i = n-1.
As described above, as shown in FIG. 10, the short-axis weighted average image data A ′ (i, j) is converted into the short-axis interpolated intermediate image data C (i, L).
図6へ進み、ステップS13では、LCDの長軸方向の実際に表示に使用する画素数であるLCD長軸画素数p=n/ryを求める。例えばn=4096、ry=2.774なら、p=1477になる。
ステップS14では、LCDの長軸方向に画素1個当たり3個が並ぶサブピクセルを独立に駆動するため画素数pを3倍したサブピクセル数(p*3)を用いて長軸縮退率rx=n/(p*3)を求める。例えばn=4096、p=1477なら、rx=0.9244になる。
Proceeding to FIG. 6, in step S13, the LCD major axis pixel number p = n / ry, which is the number of pixels actually used for display in the major axis direction of the LCD, is obtained. For example, if n = 4096 and ry = 2.774, then p = 1477.
In step S14, in order to independently drive the subpixels in which three pixels are arranged in the major axis direction of the LCD, the number of subpixels (p * 3) obtained by multiplying the number of pixels p by three (p * 3) is used. Find n / (p * 3). For example, if n = 4096 and p = 1477, rx = 0.9244.
ステップS15では、LCD長軸座標k=0に初期化する。 In step S15, the LCD major axis coordinate k = 0 is initialized.
ステップS16では、LCD長軸座標kに対応する画像基準長軸座標ip=int(k*rx)を求める。
ステップS17では、長軸補間重みWx(k)=(k*rx)−ipを求める。
ステップS18では、LCD短軸座標L=0に初期化する。
In step S16, an image reference long axis coordinate ip = int (k * rx) corresponding to the LCD long axis coordinate k is obtained.
In step S17, the long axis interpolation weight Wx (k) = (k * rx) −ip is obtained.
In step S18, the LCD minor axis coordinate L = 0 is initialized.
ステップS19では、図11および次式に示すように、短軸補間中間画像データC(i,L)と長軸補間重みWx(k)を用いて線形補間により表示画像データM’(k,L)を求める。
M’(k,L)=C(ip,L)*(1−Wx(k))+C(ip+1,L)*Wx(k)
In step S19, as shown in FIG. 11 and the following equation, the display image data M ′ (k, L) is obtained by linear interpolation using the short-axis interpolation intermediate image data C (i, L) and the long-axis interpolation weight Wx (k). )
M ′ (k, L) = C (ip, L) * (1−Wx (k)) + C (ip + 1, L) * Wx (k)
ステップS20及びS21では、L=q−1までステップS19を繰り返し、ステップS22へ進む。 In steps S20 and S21, step S19 is repeated until L = q−1, and the process proceeds to step S22.
ステップS22及びS23では、k=p*3−1までステップS16〜S21を繰り返し、ステップS24へ進む。これにより、図12に示すように、表示画像データM’(k,L)が得られる。 In steps S22 and S23, steps S16 to S21 are repeated until k = p * 3-1, and the process proceeds to step S24. As a result, as shown in FIG. 12, display image data M ′ (k, L) is obtained.
ステップS24では、短軸方向の解像度向上のために、3点の重み付け計算による周波数強調処理を表示画像データM’(k,L)に施し、最終的な表示画像データM(k,L)を得る。これにより、図13に示すように、最終的な表示画像データM(k,L)が得られる。
M(k,L)=M’(k,L-1)*e+M’(k,L)*(1-2*e)+M’(k,L+1)*e
但し、M(k,0)=M’(k,0)、M(k,q-1)=M’(k,q-1)とする。
k=0,1,…,p*3-1、L=0,1,…,q-1
ここで、eは、周波数強調処理における重み付け係数であり、長軸方向とのバランスや表示された画像の印象を考慮して例えば−0.1〜−0.2の間で調整する。
In step S24, in order to improve the resolution in the minor axis direction, frequency enhancement processing by weighting calculation of three points is performed on the display image data M ′ (k, L), and the final display image data M (k, L) is obtained. obtain. As a result, final display image data M (k, L) is obtained as shown in FIG.
M (k, L) = M ′ (k, L−1) * e + M ′ (k, L) * (1-2 * e) + M ′ (k, L + 1) * e
However, it is assumed that M (k, 0) = M ′ (k, 0) and M (k, q−1) = M ′ (k, q−1).
k = 0,1, ..., p * 3-1, L = 0,1, ..., q-1
Here, e is a weighting coefficient in the frequency enhancement process, and is adjusted between, for example, −0.1 to −0.2 in consideration of the balance with the major axis direction and the impression of the displayed image.
ステップS25では、表示画像データM(k,L)をパソコン3からモノクロ液晶表示装置4へ出力する。
In step S25, the display image data M (k, L) is output from the
モノクロ液晶表示装置4では、LCDのサブピクセルを独立に駆動して、図14に示すように、表示画像データM(k,L)を表示する。
In the monochrome liquid
実施例1の医療画像表示装置1および医療画像ビューワープログラム31によれば、次の効果が得られる。
(1)LCDで実際に表示に使用する長軸方向のサブピクセル数は(p*3)個であるが、これを画素数に換算すると(n*q/m)個となる。よって、表示画像のアスペクト比は(n*q/m)/q=(n/m)となる。つまり、元の医療画像データA(i,j)のアスペクト比が保存されている。
(2)サブピクセルを独立に駆動しない場合に比べて長軸方向の解像度を3倍にすることが出来る。FPDの画素ピッチに換算したLCDのサブピクセル・ピッチは、FPDの画素ピッチのn/(3*n*q/m)=m/(3*q)倍になるので、例えばFPDの画素ピッチが70μm、m=3328、q=1200なら、FPDの画素ピッチに換算したLCDのサブピクセル・ピッチは65(=70*3328/(3*1200))μmになる。これにより、乳ガン診断において必要な解像度100μmに十分対応できることとなる。
(3)以上より、医師がX線撮影画像,CT画像,MR画像などの医療画像を読影診断する際、その診断の効率化を図ることが出来る。
According to the medical
(1) The number of sub-pixels in the major axis direction actually used for display on the LCD is (p * 3), but when converted into the number of pixels, it is (n * q / m). Therefore, the aspect ratio of the display image is (n * q / m) / q = (n / m). That is, the aspect ratio of the original medical image data A (i, j) is stored.
(2) The resolution in the major axis direction can be tripled compared to the case where the subpixels are not driven independently. Since the LCD sub-pixel pitch converted to the FPD pixel pitch is n / (3 * n * q / m) = m / (3 * q) times the FPD pixel pitch, for example, the FPD pixel pitch is If 70 μm, m = 3328, and q = 1200, the subpixel pitch of the LCD converted to the pixel pitch of FPD is 65 (= 70 * 3328 / (3 * 1200)) μm. As a result, the resolution of 100 μm necessary for breast cancer diagnosis can be sufficiently handled.
(3) From the above, when a doctor makes an interpretation diagnosis of a medical image such as an X-ray image, CT image, or MR image, the diagnosis can be made more efficient.
本発明の医療画像表示装置および医療画像表示プログラムは、医師がX線撮影画像,CT画像,MR画像などの医療画像を読影診断する際、その診断の効率化を図るのに利用できる。 The medical image display apparatus and medical image display program of the present invention can be used to improve the efficiency of diagnosis when a doctor interprets a medical image such as an X-ray image, CT image, or MR image.
1 医療画像表示装置
3 パソコン(パーソナルコンピュータ)
4 モノクロ液晶表示装置
31 医療画像ビューワープログラム
1 medical
4
Claims (6)
前記画像変換手段が、前記被変換画像の全ての画素が前記表示画像に反映されるように前記被変換画像を列方向にm行からq行に変換すると共に前記被変換画像のアスペクト比が維持されるように行方向にn列から(3*n*q/m)列に変換するに際し、列方向の変換は、列方向に隣接する複数の画素を参照画素として被変換画像における列方向の複数の画素から略同一数の画素を生成する加重平均演算処理と、列方向縮退補間演算にて列方向の画素数をq画素に削減し、行方向の変換は、列方向の演算にて得た行方向に隣接する複数の画素を参照画素とした行方向補間演算にて表示画像を生成する
ことを特徴とする医療画像表示装置。 When m, n, and q are natural numbers and m> q, the m-by-n converted image in the monochrome original medical image is converted into a display image composed of q rows (3 * n * q / m) columns of pixels. Medical image display comprising: an image converting means for displaying and a monochrome liquid crystal display device for independently driving the rectangular sub-pixels constituting each of the square pixels arranged in the row direction of the monochrome LCD to display the display image A device,
The image converting means, the aspect ratio of the object to be converted image and converts the q-line from m rows said to be converted image in the column direction so that all the pixels is reflected on the display image of the converted image is maintained As described above, when converting from n columns to (3 * n * q / m) columns in the row direction, the column direction conversion is performed by using a plurality of pixels adjacent in the column direction as reference pixels. The number of pixels in the column direction is reduced to q pixels by weighted average calculation processing for generating substantially the same number of pixels from a plurality of pixels and column direction degenerate interpolation calculation, and conversion in the row direction is obtained by calculation in the column direction. A medical image display device that generates a display image by a row direction interpolation calculation using a plurality of pixels adjacent in the row direction as reference pixels .
コンピュータを、
前記被変換画像をq行(3*n*q/m)列の画素からなる前記表示画像に変換する画像変換手段、及び
前記モノクロLCDのサブピクセルを独立に駆動して前記表示画像を表示させるべく前記表示画像を前記モノクロ液晶表示装置に出力する表示画像出力手段、
として機能させるための医療画像表示プログラムであって、
前記画像変換手段が、
前記被変換画像の全ての画素が前記表示画像に反映されるように前記被変換画像を列方向にm行からq行に変換すると共に前記被変換画像のアスペクト比が維持されるように行方向にn列から(3*n*q/m)列に変換するステップを有し、
列方向の変換は、列方向に隣接する複数の画素を参照画素として被変換画像における列方向の複数の画素から略同一数の画素を生成する加重平均演算処理と、列方向縮退補間演算にて列方向の画素数をq画素に削減し、行方向の変換は、列方向の演算にて得た行方向に隣接する複数の画素を参照画素とした行方向補間演算にて表示画像を生成する
ことを特徴とする医療画像表示プログラム。 When m, n, and q are natural numbers and m> q, three rectangular subpixels are formed by arranging three m rows and n columns of converted images in the monochrome original medical image in the row direction of the monochrome LCD. In order to drive the pixels independently and display them on a monochrome liquid crystal display device,
Computer
The displays the display image by driving the image conversion means for converting the target converted image on the display image composed of pixels q rows (3 * n * q / m ) column, and the monochrome LCD subpixels independently Display image output means for outputting the display image to the monochrome liquid crystal display device,
A medical image display program for functioning as
The image converting means
Row direction so that all the pixels aspect ratio of the object to be converted image and converts the q rows from m rows said to be converted image in the column direction as reflected in the display image is maintained in the object to be converted image comprising the step of converting the n column (3 * n * q / m ) column,
The conversion in the column direction is performed by weighted average calculation processing for generating substantially the same number of pixels from a plurality of pixels in the column direction in the converted image using a plurality of pixels adjacent in the column direction as reference pixels, and column direction degenerate interpolation calculation. The number of pixels in the column direction is reduced to q pixels, and the conversion in the row direction generates a display image by a row direction interpolation calculation using a plurality of pixels adjacent in the row direction obtained by the column direction calculation as reference pixels. <br/> A medical image display program characterized by the above.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013054649A JP5643363B2 (en) | 2013-03-18 | 2013-03-18 | Medical image display device and medical image display program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013054649A JP5643363B2 (en) | 2013-03-18 | 2013-03-18 | Medical image display device and medical image display program |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2011118553A Division JP2011235107A (en) | 2011-05-27 | 2011-05-27 | Medical image display device and medical image display program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2013179599A JP2013179599A (en) | 2013-09-09 |
| JP5643363B2 true JP5643363B2 (en) | 2014-12-17 |
Family
ID=49270819
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013054649A Active JP5643363B2 (en) | 2013-03-18 | 2013-03-18 | Medical image display device and medical image display program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5643363B2 (en) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03228488A (en) * | 1990-02-01 | 1991-10-09 | Nippon Television Network Corp | Television system |
| US6278434B1 (en) * | 1998-10-07 | 2001-08-21 | Microsoft Corporation | Non-square scaling of image data to be mapped to pixel sub-components |
| JP2004198792A (en) * | 2002-12-19 | 2004-07-15 | Konica Minolta Holdings Inc | Display device and display method |
| JP2005017559A (en) * | 2003-06-25 | 2005-01-20 | Sony Corp | Image display device |
-
2013
- 2013-03-18 JP JP2013054649A patent/JP5643363B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2013179599A (en) | 2013-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103034973B (en) | Based on the adaptive image scaling method of bicubic interpolation | |
| JP6523638B2 (en) | Display panel driver, display device, image processing apparatus and image processing method | |
| EP3522104A1 (en) | Real-time virtual reality acceleration method and device | |
| CN1871634A (en) | Systems and methods for image reconstruction and subpixel rendering to enable scaling for multimode displays | |
| CN109388448A (en) | Image display method, display system and computer readable storage medium | |
| JP2004085608A5 (en) | ||
| WO2015100987A1 (en) | Pixel array, driving method for pixel array, display panel and display device | |
| JP2019120749A (en) | Display controller, image projection system, control method and program | |
| JP2009053455A (en) | Image display device, image processing circuit and method thereof | |
| JP2017015751A5 (en) | Display panel driver and display panel driving method | |
| JP4773522B2 (en) | Image display device and image display program | |
| CN102067204A (en) | Image processing device, display device, image processing method, program, and recording medium | |
| TW201629928A (en) | Pixel driving method and associated display device | |
| Barshan et al. | 35.3: Resolution enhancement based on shifted superposition | |
| JP5643363B2 (en) | Medical image display device and medical image display program | |
| JP2011235107A (en) | Medical image display device and medical image display program | |
| JP2011040004A (en) | Image processing apparatus and image processing method | |
| JP2014126774A (en) | Image processor, image display device and image processing method | |
| CN105872306B (en) | Apparatus and method for image scaling | |
| Eo et al. | P‐192L: Late‐News Poster: A Novel Architecture and Algorithm for Real‐time Correction of Pincushion Distortion in Microdisplay Systems | |
| CN106652969B (en) | Driving method, driving circuit and the display device of pixel array | |
| EP2500899B1 (en) | Image displaying apparatus and image display program | |
| JP6326914B2 (en) | Interpolation apparatus and interpolation method | |
| TW201234856A (en) | Display device and display method | |
| EP2544145B1 (en) | Method, arrangement, computer programm and computer-readable storage medium for scaling two-dimensional structures |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20130718 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140128 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140225 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140428 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20141028 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141030 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5643363 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |