KR101160959B1 - Method and apparatus of interpolating delay in ultrasound beamformer - Google Patents

Method and apparatus of interpolating delay in ultrasound beamformer Download PDF

Info

Publication number
KR101160959B1
KR101160959B1 KR1020100016437A KR20100016437A KR101160959B1 KR 101160959 B1 KR101160959 B1 KR 101160959B1 KR 1020100016437 A KR1020100016437 A KR 1020100016437A KR 20100016437 A KR20100016437 A KR 20100016437A KR 101160959 B1 KR101160959 B1 KR 101160959B1
Authority
KR
South Korea
Prior art keywords
delay time
delay
increase
interpolation
weight
Prior art date
Application number
KR1020100016437A
Other languages
Korean (ko)
Other versions
KR20110096903A (en
Inventor
유양모
송태경
강지운
윤창한
김기덕
Original Assignee
서강대학교산학협력단
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 서강대학교산학협력단 filed Critical 서강대학교산학협력단
Priority to KR1020100016437A priority Critical patent/KR101160959B1/en
Publication of KR20110096903A publication Critical patent/KR20110096903A/en
Application granted granted Critical
Publication of KR101160959B1 publication Critical patent/KR101160959B1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00Ultrasound therapy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4427Device being portable or laptop-like
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5284Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving retrospective matching to a physiological signal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00Ultrasound therapy
    • A61N7/02Localised ultrasound hyperthermia
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00Ultrasound therapy
    • A61N2007/0056Beam shaping elements
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00Ultrasound therapy
    • A61N2007/0086Beam steering

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

본 발명은 초음파 빔포머에서의 지연시간 보간 방법에 관한 것으로서 3 개의 영상점 P1, P2, P3의 지연시간을 레지스터에 저장하고, P2가 P1과 P3 사이에 위치할 때, P1과 P2 사이의 지연시간 증가분과 P1과 P3 사이의 지연시간 증가분을 산출한 다음, P1과 P2 사이의 지연시간 증가분과 P1과 P3 사이의 지연시간 증가분을 이용하여 P1과 P2 사이에 존재하는 영상점들의 지연시간을 보간하는 것을 특징으로 하며, 모든 채널에 대해서 하나의 지연 시간 계산기를 가지고 지연 시간을 계산할 수 있고, 동일한 하드웨어 복잡도를 갖는 경우에 지연시간 보간의 정확도를 높일 수 있다.The present invention relates to a delay time interpolation method in an ultrasonic beamformer, and stores delay times of three image points P1, P2, and P3 in a register, and delays between P1 and P2 when P2 is located between P1 and P3. Calculate the time increment and the delay increase between P1 and P3, and then interpolate the delay times of the image points existing between P1 and P2 using the delay increase between P1 and P2 and the delay increase between P1 and P3. The delay time can be calculated using a delay calculator for all channels, and the accuracy of delay interpolation can be increased when the hardware complexity is the same.

Description

초음파 빔포머에서의 지연시간 보간 방법 및 장치{Method and apparatus of interpolating delay in ultrasound beamformer}Method and apparatus for delay time interpolation in ultrasonic beamformer

본 발명은 초음파 빔포머에서의 지연시간 보간 방법에 관한 것으로서, 더욱 상세하게는 모든 채널에 대해서 하나의 지연 시간 계산기를 가지고 지연 시간을 계산하되 그 정확도를 더 향상시킬 수 있는 초음파 빔포머에서의 지연시간 보간 방법에 관한 것이다.The present invention relates to a delay time interpolation method in an ultrasonic beamformer, and more particularly, a delay in an ultrasonic beamformer that can further improve the accuracy of the delay time with one delay time calculator for all channels. Relates to a time interpolation method.

기존의 초음파 의료장비 시스템의 동적 수신 빔 집속기에서는 각 채널별 지연 시간 계산을 통하여 점 광원 신호에서 온다고 가정되는 원형 방사 신호를 하나의 신호로 다시 구성한다.In the dynamic receiving beam focuser of the conventional ultrasound medical equipment system, the circular radiation signal assumed to come from the point light source signal is reconstructed into one signal by calculating the delay time for each channel.

도 1은 일반적인 수신 동적 빔집속을 나타낸 개념도이다.1 is a conceptual diagram illustrating a general reception dynamic beam focusing.

도 1을 참조하면, 주사선 상의 화소점으로부터 온다고 가정되는 원형 방사 신호가 배열 변환자에 도달하는 것이 나타나 있다. 배열 변환자에 도달한 수신 신호는 각각 다른 지연시간을 갖고 매질을 통해 배열 변환자에 도달한 신호이므로, 가변 시간지연 과정을 통해 위상차를 보상한 신호를 집속하여 하나의 신호로 구성하게 됨을 알 수 있다.Referring to FIG. 1, it is shown that a circular radiated signal assumed to come from a pixel point on a scan line reaches an array transducer. Since the received signals reaching the array transducers have different delay times and reached the array transducers through the medium, the signals obtained by compensating the phase difference through the variable time delay process are composed of one signal. have.

초소형 초음파 의료장비 시스템에 적용하기 위하여는 지연 시간 계산기를 각 채널별로 두는 것이 아니라, 하나의 지연 시간 계산기를 가지고 모든 채널의 지연시간을 계산하게 되며 (채널수-1) 만큼의 지연시간이 계산되지 않은 영상점에 대해서는 선형 보간을 통하여 지연 시간을 얻어내므로, 더 적은 하드웨어 자원을 가지고 빔 집속을 수행할 수 있다.In order to apply to the ultra-small ultrasound medical equipment system, instead of having a delay calculator for each channel, the delay time of all channels is calculated with one delay calculator. Since the delay time is obtained through linear interpolation for the missing image point, beam focusing can be performed with less hardware resources.

물론, 이 경우에는 계산된 지연 시간값 사이, 즉 (채널수-1) 만큼의 영상점에 대해서는 직접 계산한 지연 시간값에 비교하여 오차를 갖는다. 이러한 오차는 초소형 의료장비 시스템의 성능을 개선하는 데 불리한 조건이 될 수 있으므로 기존 방식과 동일하게 모든 채널에 대해서 하나의 지연 시간 계산기를 가지고 지연 시간을 계산하되 그 정확도를 더 향상시킬 수 있는 보간 방식이 필요하다.Of course, in this case, there is an error compared to the delay time value calculated directly between the calculated delay time values, that is, the number of image points corresponding to (number of channels-1). This error can be an adverse condition for improving the performance of the micro medical device system. Therefore, as in the conventional method, an interpolation method can be used to calculate the delay time with one delay calculator for all channels but improve the accuracy more. This is necessary.

따라서, 본 발명이 해결하고자 하는 첫 번째 과제는 모든 채널에 대해서 하나의 지연 시간 계산기를 가지고 지연 시간을 계산하되 그 정확도를 더 향상시킬 수 있는 지연시간 보간 방법을 제공하는 것이다.Accordingly, the first problem to be solved by the present invention is to provide a delay time interpolation method that can calculate the delay time with one delay calculator for all channels, and further improve the accuracy.

본 발명이 해결하고자 하는 두 번째 과제는 모든 채널에 대해서 하나의 지연 시간 계산기를 가지고 지연 시간을 계산하되 그 정확도를 더 향상시킬 수 있는 지연시간 보간 장치를 제공하는 것이다.A second problem to be solved by the present invention is to provide a delay time interpolation device that can calculate the delay time with one delay time calculator for all channels and further improve its accuracy.

또한, 상기된 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록 매체를 제공하는데 있다.Further, the present invention provides a computer-readable recording medium having recorded thereon a program for executing the above method on a computer.

본 발명은 상기 첫 번째 과제를 달성하기 위하여, 3 개의 영상점 P1, P2, P3의 지연시간을 레지스터에 저장하는 단계; 상기 P2가 상기 P1과 상기 P3 사이에 위치할 때, 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 산출하는 단계; 및 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 이용하여 상기 P1과 상기 P2 사이에 존재하는 영상점들의 지연시간을 보간하는 단계를 포함하는 지연시간 보간 방법을 제공한다.The present invention comprises the steps of storing the delay time of the three image points P1, P2, P3 in the register to achieve the first object; Calculating an increase in delay time between P1 and P2 and an increase in delay time between P1 and P3 when P2 is located between P1 and P3; And interpolating the delay times of the image points existing between the P1 and the P2 by using the delay time increment between the P1 and the P2 and the delay time increment between the P1 and the P3. To provide.

본 발명의 일 실시예에 의하면, 상기 P1과 상기 P2 사이에 존재하는 보간할 영상점의 위치에 따라 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분에 적용하는 가중치를 다르게 적용함으로써, 상기 P1과 상기 P2 사이에 존재하는 영상점들의 지연시간을 보간할 수 있다.According to an embodiment of the present invention, the delay time increases between the P1 and the P2 and the delay time increases between the P1 and the P3 depending on the position of the video point to be interpolated between the P1 and the P2. By applying the weights differently, the delay times of the image points existing between the P1 and the P2 may be interpolated.

또한, 상기 P1과 상기 P2 사이에 존재하는 보간할 영상점이 상기 P1에 가까울수록 상기 P1과 상기 P2 사이의 지연시간 증가분의 가중치를 증가시키고, 상기 P3에 가까울수록 P1과 상기 P3 사이의 지연시간 증가분의 가중치를 증가시킬 수 있다.In addition, as the image point to be interpolated between P1 and P2 is closer to P1, the weight of the delay time increase between P1 and P2 increases, and the closer to P3, the delay time increases between P1 and P3. It can increase the weight of.

본 발명의 다른 실시예에 의하면, 상기 P2와 상기 P3 사이에 존재하는 보간할 영상점의 위치에 따라 상기 P2와 상기 P3 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분에 적용하는 가중치를 다르게 적용함으로써, 상기 P2와 상기 P3 사이에 존재하는 영상점들의 지연시간을 보간할 수 있다.According to another embodiment of the present invention, the delay time increases between the P2 and the P3 and the delay time increases between the P1 and the P3 according to the position of the video point to be interpolated between the P2 and the P3 By applying the weights differently, the delay times of the image points existing between the P2 and the P3 can be interpolated.

또한, 상기 P2와 상기 P3 사이에 존재하는 보간할 영상점이 상기 P2에 가까울수록 상기 P1와 상기 P3 사이의 지연시간 증가분의 가중치를 증가시키고, 상기 P3에 가까울수록 P2와 상기 P3 사이의 지연시간 증가분의 가중치를 증가시킬 수 있다.In addition, as the image point to be interpolated between the P2 and the P3 is closer to the P2, the weight of the delay time increment between the P1 and the P3 increases, and the closer to the P3, the delay time increases between the P2 and the P3. It can increase the weight of.

또한, 상기 P1과 상기 P2 사이의 지연시간 증가분은 상기 P1의 지연시간과 상기 P2의 지연시간 간의 차이를 채널 수로 나눔으로써 산출되는 것이 바람직하다.In addition, the increase in delay time between P1 and P2 is preferably calculated by dividing the difference between the delay time of P1 and the delay time of P2 by the number of channels.

본 발명은 상기 두 번째 과제를 달성하기 위하여, 3 개의 영상점 P1, P2, P3의 지연시간을 저장하는 레지스터; 상기 P2가 상기 P1과 상기 P3 사이에 위치할 때, 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 산출하는 증가분 산출부; 및 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 이용하여 상기 P1과 상기 P2 사이에 존재하는 영상점들의 지연시간을 보간하는 지연시간 보간부를 포함하는 지연시간 보간 장치를 제공한다.The present invention is a register for storing the delay time of the three image points P1, P2, P3 to achieve the second object; An increment calculator configured to calculate an increase in delay time between P1 and P2 and an increase in delay time between P1 and P3 when P2 is located between P1 and P3; And a delay time interpolation unit for interpolating the delay time of the image points existing between the P1 and the P2 by using the delay time increase between the P1 and the P2 and the delay time increase between the P1 and the P3. Provide a time interpolation device.

상기 다른 기술적 과제를 해결하기 위하여, 본 발명은 상기된 지연시간 보간 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록 매체를 제공한다. In order to solve the above other technical problem, the present invention provides a computer-readable recording medium that records a program for executing the above-described delay time interpolation method on a computer.

본 발명에 따르면, 모든 채널에 대해서 하나의 지연 시간 계산기를 가지고 지연 시간을 계산할 수 있다. 또한, 본 발명에 따르면, 동일한 하드웨어 복잡도를 갖는 경우에 지연시간 보간의 정확도를 높일 수 있다.According to the present invention, the delay time can be calculated with one delay calculator for all channels. In addition, according to the present invention, it is possible to increase the accuracy of delay time interpolation in the case of having the same hardware complexity.

도 1은 일반적인 수신 동적 빔집속을 나타낸 개념도이다.
도 2는 본 발명의 바람직한 일 실시예에 따른 초음파 빔집속 장치의 블록도이다.
도 3은 본 발명의 일 실시예에 따른 초음파 빔집속 장치에서 지연시간을 보간하는 방법을 나타낸 것이다.
도 4는 지연시간이 레지스터에 저장된 영상점들 사이에 존재하는 영상점들의 지연시간을 보간하는 방법을 개략적으로 나타낸 것이다.
도 5는 본 발명의 일 실시예에 따른 초음파 빔집속 장치에서 지연시간을 보간하는 방법을 보다 상세하게 나타낸 것이다.
도 6은 본 발명의 바람직한 일 실시예에 따른 지연시간 보간 방법의 흐름도이다.
도 7은 본 발명의 바람직한 일 실시예에 따른 지연시간 보간 장치의 블록도이다.
도 8은 채널축에 대한 누적 에러량을 나타낸 그래프와 선형 보간기와 이중 선형 보간기의 누적 에러량에 대한 비율을 나타낸 그래프를 도시한 것이다.
도 9는 영상점축에 대한 누적 에러량을 나타낸 그래프와 선형 보간기와 이중 선형 보간기의 누적 에러량에 대한 비율을 나타낸 그래프를 도시한 것이다.
1 is a conceptual diagram illustrating a general reception dynamic beam focusing.
2 is a block diagram of an ultrasonic beam focusing apparatus according to a preferred embodiment of the present invention.
3 illustrates a method of interpolating a delay time in an ultrasonic beam focusing apparatus according to an embodiment of the present invention.
4 schematically illustrates a method of interpolating delay times of image points existing between image points stored in a register.
5 illustrates in more detail a method of interpolating a delay time in an ultrasonic beam focusing apparatus according to an embodiment of the present invention.
6 is a flowchart of a delay time interpolation method according to an exemplary embodiment of the present invention.
7 is a block diagram of a delay time interpolation apparatus according to an exemplary embodiment of the present invention.
8 is a graph showing the cumulative error amount with respect to the channel axis, and a graph showing the ratio of the cumulative error amount between the linear interpolator and the dual linear interpolator.
9 is a graph showing the cumulative error amount with respect to the image axis and a graph showing the ratio of the cumulative error amount of the linear interpolator and the dual linear interpolator.

본 발명에 관한 구체적인 내용의 설명에 앞서 이해의 편의를 위해 본 발명이 해결하고자 하는 과제의 해결 방안의 개요 혹은 기술적 사상의 핵심을 우선 제시한다.Prior to the description of the specific contents of the present invention, for the convenience of understanding, the outline of the solution of the problem to be solved by the present invention or the core of the technical idea will be presented first.

본 발명의 일 실시예에 따른 초음파 빔포머에서의 지연시간 보간 방법은 3 개의 영상점 P1, P2, P3의 지연시간을 레지스터에 저장하고, 상기 P2가 상기 P1과 상기 P3 사이에 위치할 때, 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 산출한 다음, 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 이용하여 상기 P1과 상기 P2 사이에 존재하는 영상점들의 지연시간을 보간하는 단계를 포함한다.Delay time interpolation method in the ultrasonic beamformer according to an embodiment of the present invention stores the delay time of the three image points P1, P2, P3 in the register, when the P2 is located between the P1 and the P3, After calculating the delay time increase between P1 and P2 and the delay time increase between P1 and P3, the delay time increase between P1 and P2 and the delay time increase between P1 and P3 are used. And interpolating a delay time of the image points existing between the P1 and the P2.

이하, 바람직한 실시예를 들어 본 발명을 더욱 상세하게 설명한다. 그러나 이들 실시예는 본 발명을 보다 구체적으로 설명하기 위한 것으로, 본 발명의 범위가 이에 의하여 제한되지 않는다는 것은 당업계의 통상의 지식을 가진 자에게 자명할 것이다. 본 발명이 해결하고자 하는 과제의 해결 방안을 명확하게 하기 위한 발명의 구성을 본 발명의 바람직한 실시예에 근거하여 첨부 도면을 참조하여 상세히 설명하되, 도면의 구성요소들에 참조번호를 부여함에 있어서 동일 구성요소에 대해서는 비록 다른 도면상에 있더라도 동일 참조번호를 부여하였으며 당해 도면에 대한 설명시 필요한 경우 다른 도면의 구성요소를 인용할 수 있음을 미리 밝혀둔다. 아울러 본 발명과 관련된 공지 기능 혹은 구성에 대한 구체적인 설명 그리고 그 이외의 제반 사항이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우, 그 상세한 설명을 생략한다.Hereinafter, the present invention will be described in more detail with reference to preferred examples. However, these examples are intended to illustrate the present invention in more detail, it will be apparent to those skilled in the art that the scope of the present invention is not limited thereby. The configuration of the invention for clarifying the solution to the problem to be solved by the present invention will be described in detail with reference to the accompanying drawings based on the preferred embodiment of the present invention, the same in the reference numerals to the components of the drawings The same reference numerals are given to the components even though they are on different drawings, and it is to be noted that in the description of the drawings, components of other drawings may be cited if necessary. In the following description, a detailed description of known functions and configurations incorporated herein will be omitted when it may make the subject matter of the present invention rather unclear.

도 2는 본 발명의 바람직한 일 실시예에 따른 초음파 빔집속 장치의 블록도이다.2 is a block diagram of an ultrasonic beam focusing apparatus according to a preferred embodiment of the present invention.

도 2를 참조하면, 본 실시예에 따른 초음파 빔집속 장치는 초음파 배열 변환자(200), 메모리부(210), 레지스터 블록부(220), 소수지연 선택 합성부(230), 소수지연 필터부(240), 합성부(250), 확장구경 누산부(260), 및 지연시간 계산부(270)로 구성된다.Referring to FIG. 2, the ultrasonic beam focusing apparatus according to the present exemplary embodiment includes an ultrasonic array transducer 200, a memory unit 210, a register block unit 220, a fractional delay select synthesis unit 230, and a fractional delay filter unit. 240, a synthesizer 250, an expansion diameter accumulator 260, and a delay time calculator 270.

초음파 배열 변환자(200)는 집속 송신한 신호가 매질을 지나 조직에서 반사되어 온 초음파 신호를 수신한다.The ultrasound array transducer 200 receives an ultrasound signal in which the focused signal is reflected from tissue through a medium.

메모리부(210)는 초음파 배열 변환자(200)가 수신한 초음파 신호를 저장한다. 메모리부(210)에 포함된 메모리는 채널 수만큼 각 채널에 대응하여 존재하는 것이 바람직하다. 이때 메모리부(210)는 초음파 신호를 4f0로 샘플링한 데이터를 저장하고, 4f0로 샘플링된 데이터는 16f0의 지연 해상도(delay resolution)를 갖도록 보간될 것이다.The memory unit 210 stores the ultrasonic signal received by the ultrasonic array transducer 200. The memory included in the memory unit 210 may exist corresponding to each channel by the number of channels. In this case, the memory unit 210 may store data sampled at 4f 0 of the ultrasonic signal, and the data sampled at 4f 0 may be interpolated to have a delay resolution of 16f 0 .

레지스터 블록부(220)는 지연시간 계산부(270)에서 계산된 코스 지연(Coarse delay)에 따라 메모리부(210)에 저장된 데이터를 레지스터에 저장한다.The register block 220 stores data stored in the memory 210 in a register according to a coarse delay calculated by the delay time calculator 270.

레지스터 블록부(220)에 포함된 레지스터는 각 채널별로 존재하고, 소수지연 필터부(240)의 탭수와 같은 갯수의 레지스터가 각 채널별로 존재하는 것이 바람직하다. 도 2를 참조하면, 각 채널별 레지스터의 갯수는 8개이므로, 소수지연 필터부(240)의 탭수는 8탭이다.The registers included in the register block unit 220 exist for each channel, and the number of registers equal to the number of taps of the fractional delay filter unit 240 preferably exists for each channel. Referring to FIG. 2, since the number of registers for each channel is eight, the number of taps of the fractional delay filter unit 240 is eight taps.

소수지연 선택 합성부(230)는 지연시간 계산부(270)에서 계산된 파인 지연(fine delay)에 따라 레지스터 블록부(220)에 저장된 데이터를 각 채널별로 선택하고, 합성한 다음 소수지연 필터부(240)로 출력한다.The fractional delay selector 230 selects and stores data stored in the register block 220 for each channel according to a fine delay calculated by the delay time calculator 270, and then synthesizes the data. Output to 240.

소수지연 필터부(240)는 지연시간별로 하나씩 존재하며, 지연시간 D는 0.00, 0.25, 0.50, 0.75로 표시된다.The fractional delay filter unit 240 exists for each delay time, and the delay time D is expressed as 0.00, 0.25, 0.50, and 0.75.

합성부(250)는 소수지연 필터부(240)로부터 수신한 데이터들을 합성한 후 확장 구경 누산부(260)로 출력한다.The synthesizer 250 synthesizes data received from the fractional delay filter 240 and outputs the expanded aperture accumulator 260.

확장구경 누산부(260)는 합성부(250)에서 합성된 데이터를 누산한다. 누산된 데이터는 후처리된 후 디스플레이된다.The expansion diameter accumulator 260 accumulates the data synthesized by the combiner 250. The accumulated data is displayed after being processed.

지연시간 계산부(270)는 코스 지연(Coarse delay)과 파인 지연(fine delay)을 각각 메모리부(210)와 소수지연 선택 합성부(230)로 출력한다. 지연시간 계산부(270)가 코스 지연(Coarse delay)과 파인 지연(fine delay)을 보간하는 방법을 도 3을 이용하여 상세하게 살펴보기로 한다.The delay calculator 270 outputs a coarse delay and a fine delay to the memory 210 and the fractional delay select synthesis unit 230, respectively. A method of interpolating the coarse delay and the fine delay by the delay time calculator 270 will be described in detail with reference to FIG. 3.

도 3은 본 발명의 일 실시예에 따른 초음파 빔집속 장치에서 지연시간을 보간하는 방법을 나타낸 것이다.3 illustrates a method of interpolating a delay time in an ultrasonic beam focusing apparatus according to an embodiment of the present invention.

일반적으로 수신 빔 집속을 위해서는 모든 채널에 대해 지연시간 계산부가 필요하다. 하지만 이는 하드웨어 복잡도가 채널수에 비례하게 증가하기 때문에 본 발명의 일 실시예에서는 하나의 지연시간 계산부를 가지고 모든 채널의 지연시간을 계산한 이후에 각 채널에서 (채널수-1)만큼 계산되지 못한 영상점에 대해서는 본 발명의 일 실시예에 따른 지연시간 보간방법을 사용하여 보간을 수행함으로써 보다 적은 하드웨어 자원을 가지고 빔집속을 수행할 수 있게 된다. In general, a delay calculation unit is required for all channels for the reception beam focusing. However, since the hardware complexity increases in proportion to the number of channels, in one embodiment of the present invention, after calculating the delay times of all channels with one delay calculation unit, the channel complexity cannot be calculated for each channel (number of channels-1). For image points, beam interpolation can be performed with less hardware resources by performing interpolation using a delay time interpolation method according to an embodiment of the present invention.

도 3을 참조하면, 채널에 따라 돌아가면서 한 영상점씩 지연시간이 저장되고 각 채널에서 저장되지 않은 영상점의 지연시간 부분에 대해서는 보간을 통해 지연시간을 구한다.Referring to FIG. 3, the delay time is stored by one image point as the channel is rotated, and the delay time is obtained through interpolation for the delay time portion of the image point which is not stored in each channel.

도 4는 지연시간이 레지스터에 저장된 영상점들 사이에 존재하는 영상점들의 지연시간을 보간하는 방법을 개략적으로 나타낸 것이다.4 schematically illustrates a method of interpolating delay times of image points existing between image points stored in a register.

l과 l+1 사이의 지연시간 보간은 (l)-(l+1) 지연시간 보간선과 (l)-(l+2) 지연시간 보간선을 이용한다. l에 가까운 영상점들의 지연시간은 (l)-(l+1) 지연시간 보간선의 가중치를 (l)-(l+2) 지연시간 보간선의 가중치보다 많이 부여하고, l+1에 가까운 영상점들의 지연시간은 (l)-(l+2) 지연시간 보간선의 가중치를 (l)-(l+1) 지연시간 보간선의 가중치보다 많이 부여한다. 이와 같이 가중치를 부여함으로써, 실제의 지연시간에 더 가깝게 보간될 수 있을 것이다.Delay interpolation between l and l + 1 uses (l)-(l + 1) delay interpolation lines and (l)-(l + 2) delay time interpolation lines. Delay times of image points close to l give the weight of the (l)-(l + 1) delay time interpolation line more than the weight of the (l)-(l + 2) delay time interpolation line and close to l + 1. The delay time of the image points gives more weight of the (l)-(l + 2) delay time interpolation line than the weight of the (l)-(l + 1) delay time interpolation line. By weighting in this way, it can be interpolated closer to the actual delay time.

또한, l+1과 l+2 사이의 지연시간 보간은 (l+1)-(l+2) 지연시간 보간선과 (l)-(l+2) 지연시간 보간선을 이용한다. l+1에 가까운 영상점들의 지연시간은 (l)-(l+2) 지연시간 보간선의 가중치를 (l+1)-(l+2) 지연시간 보간선의 가중치보다 많이 부여하고, l+2에 가까운 영상점들의 지연시간은 (l+1)-(l+2) 지연시간 보간선의 가중치를 (l)-(l+2) 지연시간 보간선의 가중치보다 많이 부여한다.In addition, the delay interpolation between l + 1 and l + 2 uses a (l + 1)-(l + 2) delay time interpolation line and a (l)-(l + 2) delay time interpolation line. The delay time of image points near l + 1 gives the weight of (l)-(l + 2) delay time interpolation line more than the weight of (l + 1)-(l + 2) delay time interpolation line, and l Delay times of image points near +2 give more weight of (l + 1)-(l + 2) delay time interpolation line than weights of (l)-(l + 2) delay time interpolation line.

도 5는 본 발명의 일 실시예에 따른 초음파 빔집속 장치에서 지연시간을 보간하는 방법을 보다 상세하게 나타낸 것이다.5 illustrates in more detail a method of interpolating a delay time in an ultrasonic beam focusing apparatus according to an embodiment of the present invention.

먼저 지연시간을 보간하고자 하는 영상점의 인덱스를 n(n은 자연수)이라고 하고, 채널 수를 N(N은 자연수)으로 놓으면 지연시간 계산부(260)는 각 채널의 지연 시간을 N마다 한번씩 돌아가면서 계산하게 되고 본 발명의 일 실시예에 따른 초음파 빔집속 장치에서는 N-1개의 계산되지 않은 지연시간을 보간한다. 지연시간을 보간하기 위하여 종래에는 2개의 영상점의 지연시간이 필요하였으나, 본 발명의 일 실시예에 따른 지연시간 보간방법에서는 3개의 영상점의 지연시간이 필요하다. First, if the index of an image point to which the delay time is to be interpolated is n (n is a natural number), and the number of channels is set to N (N is a natural number), the delay time calculation unit 260 returns the delay time of each channel once every N times. The ultrasound beam focusing apparatus according to an embodiment of the present invention interpolates N-1 uncalculated delay times. In order to interpolate a delay time, a delay time of two image points is conventionally required, but a delay time of three image points is required in a delay time interpolation method according to an embodiment of the present invention.

2개의 영상점의 지연시간을 이용하여 보간하는 방법을 선형 지연시간 보간방법이라고 하고, 3개의 영상점의 지연시간을 이용하여 보간하는 방법을 이중 선형 지연시간 보간방법이라고 하기로 한다. 선형 지연시간 보간방법는 2개의 영상점의 지연시간 차이를 채널 수로 나누어 단위 증가량을 구한 후 보간할 영상점마다 더해가는 보간 방식인 반면, 이중 선형 지연시간 보간방법의 경우에는 3개의 영상점의 지연시간을 가지고 보간을 수행하게 된다.The interpolation method using the delay time of two image points is called linear delay time interpolation method, and the interpolation method using the delay time of three image points is called dual linear delay time interpolation method. Linear delay interpolation is a method of dividing the difference between the delay times of two video points by the number of channels to obtain a unit increment, and then adds it to each video point to be interpolated.In the case of the dual linear delay interpolation method, the delay time of three video points is applied. Interpolation is performed with.

먼저 3개의 영상점의 지연시간을 각각 들어온 순서대로 P1, P2, P3라고 하면 (P2-P1)/N으로 표현되는 slopea와 (P3-P1)/2N으로 나타나는 slopec를 구한다. 이 두 단위 증가량을 가지고 P1에서부터는 slopea로 증가시켜나가 각 영상점마다 slopea에 의한 지연시간 Su n(n은 영상점의 인덱스)를 구성하고, P2에서부터 N*slopec를 뺀 값에다가 영상점이 진행됨에 따라 다시 slopec를 더해나가 slopec에 의한 지연시간 Sd n(n은 영상점의 인덱스)를 구성한다. 이 두 값의 특성은 구성하고자 하는 영상점이 P1에 가까울 때에는 구하고자 하는 이상적인 지연시간에 Su n가 더 가깝고, 영상점이 진행되어 계산된 값인 P2에 가까워지게 되면 Sd n에 가까워지게 된다. 따라서 지연시간을 보간하고자하는 영상점이 보간 구간에서 어디까지 왔는지의 비율을 이용하여 Su n와 Sd n를 비율을 다르게 하여 더함으로써, 선형 지연시간 보간기보다 더 이상적인 계산값에 가까운 결과를 갖는 이중 선형 지연시간 보간기의 출력 Tn를 구성할 수 있게 된다.First, if the delay times of the three image points are P1, P2, and P3 in the order of arrival, the slope a represented by (P2-P1) / N and the slope c represented by (P3-P1) / 2N are obtained. With these two unit increments, increase the slope a from P1 to the slope a , and configure the delay time S u n (n is the index of the image point) by slope a for each image point, and subtract N * slope c from P2. As the image point progresses, slope c is added again to form a delay time S d n (n is an index of the image point) by slope c . The characteristics of these two values are closer to S d n when the image point to be configured is closer to P1 and S u n is closer to the ideal delay time to be obtained, and closer to P2, the calculated value. Therefore, by adding different ratios of S u n and S d n by using the ratio of how far the image point to be interpolated from the interpolation interval has, the result is closer to the calculated value than the linear delay interpolator. The output T n of the dual linear delay interpolator can be configured.

도 6은 본 발명의 바람직한 일 실시예에 따른 지연시간 보간 방법의 흐름도이다.6 is a flowchart of a delay time interpolation method according to an exemplary embodiment of the present invention.

도 6을 참조하면, 본 실시예에 따른 지연시간 보간 방법은 도 2에 도시된 지연시간 계산부(270)에서 시계열적으로 처리되는 단계들로 구성된다. 따라서, 이하 생략된 내용이라 하더라도 도 2에 도시된 지연시간 계산부(270)에 관하여 이상에서 기술된 내용은 본 실시예에 따른 지연시간 보간 방법에도 적용된다. Referring to FIG. 6, the delay time interpolation method according to the present exemplary embodiment includes steps processed in time series by the delay time calculator 270 illustrated in FIG. 2. Therefore, even if omitted below, the above description of the delay time calculator 270 shown in FIG. 2 is also applied to the delay time interpolation method according to the present embodiment.

600 단계에서 지연시간 계산부는 지연시간 레지스터0, 지연시간 레지스터1, 지연시간 레지스터2에 영상점 P1, P2, P3의 지연시간을 계산하여 각각 저장한다.In step 600, the delay calculation unit calculates and stores the delay times of the image points P1, P2, and P3 in the delay time register 0, the delay time register 1, and the delay time register 2, respectively.

610 단계에서 지연시간 계산부는 P1과 P2 사이의 지연시간 증가분과 P1와 P3 사이의 지연시간 증가분을 산출한다. P1과 P2 사이의 지연시간 증가분 slopea와 P1와 P3 사이의 지연시간 증가분 slopec은 다음의 수학식 1과 같이 산출될 수 있다.In step 610, the delay calculation unit calculates an increase in delay time between P1 and P2 and an increase in delay time between P1 and P3. The delay increase slope a between P1 and P2 and the delay increase slope c between P1 and P3 may be calculated by Equation 1 below.

Figure 112010011916700-pat00001
Figure 112010011916700-pat00001

여기서, P1은 지연시간 레지스터0에 저장된 지연시간, P2는 지연시간 레지스터1에 저장된 지연시간, P3은 지연시간 레지스터2에 저장된 지연시간을 의미하며, intp_level은 보간 레벨로서, 채널 수 N과 같은 값을 갖는 것이 바람직하다.Here, P1 is a delay time stored in delay register 0, P2 is a delay time stored in delay register 1, P3 is a delay time stored in delay register 2, intp_level is an interpolation level, and the same value as the number of channels N. It is preferable to have.

620 단계에서 지연시간 계산부는 P1과 P2 사이의 지연시간 증가분을 이용한 지연시간 보간선과 P1와 P3 사이의 지연시간 증가분을 이용한 지연시간 보간선을 산출한다.In step 620, the delay calculation unit calculates a delay time interpolation line using an increase in delay time between P1 and P2 and a delay time interpolation line using an increase in delay time between P1 and P3.

P1과 P2 사이의 지연시간 증가분을 이용한 지연시간 보간선을 intp_addr_upper라고 하고, P1와 P3 사이의 지연시간 증가분을 이용한 지연시간 보간선을 intp_addr_lower라고 하면, 다음의 수학식 2와 같이 산출될 수 있다.If the delay time interpolation line using the delay time increase between P1 and P2 is called intp_addr_upper, and the delay time interpolation line using the delay time increase between P1 and P3 is called intp_addr_lower, the following equation (2) may be calculated.

Figure 112010011916700-pat00002
Figure 112010011916700-pat00002

여기서, index는 지연시간을 보간하고자 하는 영상점의 인덱스 n을 의미하고, intp_level은 보간 레벨로서, 채널 수 N과 같은 값을 갖는 것이 바람직하다.Here, index means index n of an image point to which the delay time is to be interpolated, and intp_level is an interpolation level and preferably has a value equal to the number of channels N.

630 단계에서 지연시간 계산부는 P1과 P2 사이의 지연시간 증가분을 이용하여 산출된 지연시간 보간선과 P1와 P3 사이의 지연시간 증가분을 이용하여 산출된 지연시간 보간선 각각에 가중치를 부여하고 합산하여 가중치합을 산출한다. 산출된 가중치합은 P1과 P2 사이의 영상점들의 보간된 지연시간을 나타낸다. 이때, P1과 P2 사이의 영상점들의 보간된 지연시간 intp_address는 다음의 수학식 3과 같이 나타낼 수 있다. In step 630, the delay calculation unit weights and sums each of the delay time interpolation lines calculated using the delay time increase between P1 and P2 and the delay time interpolation lines calculated using the delay time increase between P1 and P3. Calculate the sum. The calculated weight sum represents the interpolated delay time of the image points between P1 and P2. In this case, the interpolated delay time intp_address of the image points between P1 and P2 may be expressed by Equation 3 below.

Figure 112010011916700-pat00003
Figure 112010011916700-pat00003

여기서,

Figure 112010011916700-pat00004
이다.here,
Figure 112010011916700-pat00004
to be.

이상에서는 P1과 P2 사이에 존재하는 영상점의 지연시간을 보간하는 것을 살펴보았으며, P2와 P3 사이에 존재하는 영상점의 지연시간은 P1과 P3 사이의 지연시간 증가분과 P2와 P3 사이의 지연시간 증가분을 이용하여 산출할 수 있을 것이다.In the above, the interpolation of the delay time of the image point existing between P1 and P2 has been examined, and the delay time of the image point existing between P2 and P3 is the increase of the delay time between P1 and P3 and the delay between P2 and P3. It can be calculated using the time increment.

도 7은 본 발명의 바람직한 일 실시예에 따른 지연시간 보간 장치의 블록도이다.7 is a block diagram of a delay time interpolation apparatus according to an exemplary embodiment of the present invention.

도 7을 참조하면, 본 실시예에 따른 지연시간 보간 장치는 레지스터(710), 증가분 산출부(720), 및 지연시간 보간부(730)로 구성된다.Referring to FIG. 7, the delay interpolation apparatus according to the present exemplary embodiment includes a register 710, an increment calculator 720, and a delay time interpolator 730.

레지스터(710)는 지연시간 레지스터0, 지연시간 레지스터1, 지연시간 레지스터2에 영상점 P1, P2, P3의 지연시간을 각각 저장한다.The register 710 stores the delay times of the image points P1, P2, and P3 in the delay time register 0, the delay time register 1, and the delay time register 2, respectively.

증가분 산출부(720)는 P1과 P2 사이의 지연시간 증가분과 P1와 P3 사이의 지연시간 증가분을 산출한다. 산출된 지연시간 증가분은 P1과 P2 사이에 있는 영상점의 지연시간을 보간하는데 이용된다. 반면, P2와 P3 사이에 있는 영상점을 보간하기 위해서는 P1과 P3 사이의 지연시간 증가분과 P2와 P3 사이의 지연시간 증가분을 산출하여 이용한다.The increment calculator 720 calculates the delay time increase between P1 and P2 and the delay time increase between P1 and P3. The calculated delay time increment is used to interpolate the delay time of the image point between P1 and P2. On the other hand, in order to interpolate the image point between P2 and P3, the delay time increase between P1 and P3 and the delay time increase between P2 and P3 are calculated and used.

지연시간 보간부(730)는 P1과 P2 사이의 지연시간 증가분을 이용하여 산출된 지연시간 보간선과 P1와 P3 사이의 지연시간 증가분을 이용하여 산출된 지연시간 보간선 각각에 가중치를 부여하고 합산하여 가중치합을 산출한다. The delay time interpolation unit 730 weights and sums each of the delay time interpolation lines calculated using the delay time increase between P1 and P2 and the delay time interpolation lines calculated using the delay time increase between P1 and P3. Calculate the sum of weights.

본 발명의 일 실시예에 따른 지연시간 보간방법을 이용한 초음파 빔집속 장치를 VHDL을 이용하여 하드웨어로 구현을 할 때 소요되는 하드웨어 비용은 기존의 선형 지연시간 보간 방법에 비하여 2배가량 하드웨어 자원이 많이 사용된다.Hardware cost of implementing the ultrasonic beam focusing apparatus using the VHDL in hardware using the delay time interpolation method according to an embodiment of the present invention is twice as much hardware resources as the conventional linear delay time interpolation method. Used.

그러나 실제로 동적 수신 빔집속기 상에서 가장 하드웨어 자원이 많이 필요한 부분은 초음파 배열 변환자를 통하여 받아들인 집속되지 않은 데이터를 담아둘 메모리가 대부분을 차지하게 되며 실제 지연시간 보간기를 포함한 지연시간 계산부는 전체의 1% 가량을 사용한다. 따라서 전체 시스템의 입장에서 보았을 때 이중 지연시간 보간기를 사용하는 것은 2% 정도의 자원만을 차지하므로 전체에 미치는 하드웨어 비용은 미비하다. 그러나 이하에서 살펴볼 채널축 누적 에러량과 영상점축 누적 에러량에서 볼 수 있듯이 보다 정확하게 지연시간을 보간할 수 있다.In practice, however, the most hardware resource-intensive part of the dynamic receive beam focuser is the memory that will hold the unfocused data received through the ultrasonic array transducer, and the delay calculation part including the actual delay time interpolator is 1% of the total. Use about. Therefore, from the standpoint of the whole system, using the double latency interpolator takes only about 2% of the resources, so the hardware cost for the whole is insignificant. However, the delay time can be interpolated more accurately as can be seen in the channel axis cumulative error amount and the image point accumulation error amount, which will be described below.

도 8은 채널축에 대한 누적 에러량을 나타낸 그래프와 선형 보간기와 이중 선형 보간기의 누적 에러량에 대한 비율을 나타낸 그래프를 도시한 것이다.8 is a graph showing the cumulative error amount with respect to the channel axis, and a graph showing the ratio of the cumulative error amount between the linear interpolator and the dual linear interpolator.

붉은색 선은 이중 선형 보간기, 파란 선은 선형 보간기의 채널축에 대한 누적 에러량이다.The red line is the double linear interpolator, and the blue line is the cumulative error amount along the channel axis of the linear interpolator.

도 8을 참조하면, 채널축에 대한 누적 에러량에 대하여 이중 선형 지연시간 보간기가 기존의 선형 보간기보다 좋은 성능을 지닌 것을 알 수 있다. 좌측의 누적에러량 그림에서 더 적은 에러량과 더 빠른 변동량의 감소를 보이는 것을 볼 수 있다. 32채널의 경우 앞의 200 영상점부터는 큰 차이를 보이지 않고 두 방식이 모두 비슷한 누적 에러량을 갖게 되는 것을 볼 수 있으나 채널 수가 증가함에 따라 에러의 비율은 더 깊은 영상점에 대해서도 이중 선형 보간기의 성능이 더 우월한 것을 볼 수 있다. 최대 5배의 에러 감소를 보이는데 1024개의 영상점이 샘플링 주파수가 40MHz일 때의 영상에서는 2cm를 의미하기 때문에 32채널일 때와는 달리 64 혹은 128 채널의 시스템을 구성하였을 때에는 1cm 이상의 영상 영역에 잘못된 지연시간 계산값을 반영하게 된다.Referring to FIG. 8, it can be seen that the dual linear delay time interpolator has better performance than the conventional linear interpolator with respect to the cumulative error amount on the channel axis. On the left, the cumulative error amount plot shows less error amount and a faster decrease in the amount of change. In the case of 32 channels, the two previous methods do not show a big difference, and both methods have similar cumulative error amount. However, as the number of channels increases, the error rate increases even for the deeper image points. It can be seen that performance is superior. There is a maximum error reduction of 5 times, which means 2cm for 1024 image points when the sampling frequency is 40MHz. Therefore, when 32 or 64 channels system is configured, an error of 1 cm or more is delayed. This will reflect the time calculation.

도 9는 영상점축에 대한 누적 에러량을 나타낸 그래프와 선형 보간기와 이중 선형 보간기의 누적 에러량에 대한 비율을 나타낸 그래프를 도시한 것이다.9 is a graph showing the cumulative error amount with respect to the image axis and a graph showing the ratio of the cumulative error amount of the linear interpolator and the dual linear interpolator.

도 9를 참조하여, 영상점축 누적 에러량이 채널수가 증가함에 따라 바깥쪽 채널에서 에러의 양이 어떻게 변화되는지 살펴보기로 한다.Referring to FIG. 9, it will be described how the amount of error in the outer channel changes as the number of accumulated images error is increased.

32채널일 경우 비슷한 수준을 유지하던 누적 에러 비율이 채널수가 증가하게 되면 바깥쪽으로 갈수록 급격하게 감소함을 알 수 있다. 중앙 채널에서 0.8~0.9의 비율을 갖던 누적 에러 비율이 바깥쪽 채널로 갈수록 64채널일 경우 0.2~0.3, 128채널일 경우 0.1~0.2로 낮아지는 것을 볼 수 있는데, 결국 바깥쪽 채널에 대해 더 정확한 보간을 통해 에러를 줄일 수 있다.In the case of 32 channels, the cumulative error rate, which remained at the same level, decreases rapidly as the number of channels increases. The cumulative error rate from 0.8 to 0.9 in the center channel decreases to 0.2 to 0.3 for 64 channels and 0.1 to 0.2 for 128 channels, which is more accurate for the outer channel. Interpolation can reduce errors.

본 발명의 실시예들은 다양한 컴퓨터 수단을 통하여 수행될 수 있는 프로그램 명령 형태로 구현되어 컴퓨터 판독 가능 매체에 기록될 수 있다. 상기 컴퓨터 판독 가능 매체는 프로그램 명령, 데이터 파일, 데이터 구조 등을 단독으로 또는 조합하여 포함할 수 있다. 상기 매체에 기록되는 프로그램 명령은 본 발명을 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 컴퓨터 판독 가능 기록 매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media), 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. 프로그램 명령의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 코드뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 코드를 포함한다. 상기된 하드웨어 장치는 본 발명의 동작을 수행하기 위해 하나 이상의 소프트웨어 모듈로서 작동하도록 구성될 수 있으며, 그 역도 마찬가지이다.Embodiments of the present invention may be implemented in the form of program instructions that can be executed by various computer means and recorded in a computer readable medium. The computer readable medium may include program instructions, data files, data structures, etc. alone or in combination. The program instructions recorded on the medium may be those specially designed and constructed for the present invention or may be available to those skilled in the art of computer software. Examples of computer-readable recording media include magnetic media such as hard disks, floppy disks, and magnetic tape, optical media such as CD-ROMs, DVDs, and magnetic disks, such as floppy disks. Magneto-optical media, and hardware devices specifically configured to store and execute program instructions, such as ROM, RAM, flash memory, and the like. Examples of program instructions include not only machine code generated by a compiler, but also high-level language code that can be executed by a computer using an interpreter or the like. The hardware device described above may be configured to operate as one or more software modules to perform the operations of the present invention, and vice versa.

이상과 같이 본 발명에서는 구체적인 구성 요소 등과 같은 특정 사항들과 한정된 실시예 및 도면에 의해 설명되었으나 이는 본 발명의 보다 전반적인 이해를 돕기 위해서 제공된 것일 뿐, 본 발명은 상기의 실시예에 한정되는 것은 아니며, 본 발명이 속하는 분야에서 통상적인 지식을 가진 자라면 이러한 기재로부터 다양한 수정 및 변형이 가능하다. 따라서, 본 발명의 사상은 설명된 실시예에 국한되어 정해져서는 아니되며, 후술하는 특허청구범위뿐 아니라 이 특허청구범위와 균등하거나 등가적 변형이 있는 모든 것들은 본 발명 사상의 범주에 속한다고 할 것이다.In the present invention as described above has been described by the specific embodiments, such as specific components and limited embodiments and drawings, but this is provided to help a more general understanding of the present invention, the present invention is not limited to the above embodiments. For those skilled in the art, various modifications and variations are possible from these descriptions. Therefore, the spirit of the present invention should not be limited to the described embodiments, and all of the equivalents or equivalents of the claims as well as the claims to be described later will belong to the scope of the present invention. .

Claims (13)

3 개의 영상점 P1, P2, P3의 지연시간을 레지스터에 저장하는 단계;
상기 P2가 상기 P1과 상기 P3 사이에 위치할 때, 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 산출하는 단계; 및
상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 이용하여 상기 P1과 상기 P2 사이에 존재하는 영상점들의 지연시간을 보간하는 단계를 포함하는 것을 특징으로 하는 지연시간 보간 방법.
Storing delay times of three image points P1, P2, and P3 in a register;
Calculating an increase in delay time between P1 and P2 and an increase in delay time between P1 and P3 when P2 is located between P1 and P3; And
Interpolating the delay time of the image points existing between the P1 and the P2 by using the delay time increase between the P1 and the P2 and the delay time increase between the P1 and the P3. Time interpolation method.
제 1 항에 있어서,
상기 P1과 상기 P2 사이에 존재하는 영상점들의 지연시간을 보간하는 단계는,
상기 P1과 상기 P2 사이에 존재하는 보간할 영상점의 위치에 따라 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분에 적용하는 가중치를 다르게 적용하는 것을 특징으로 하는 지연시간 보간 방법.
The method of claim 1,
Interpolating the delay time of the image points existing between the P1 and the P2,
According to the position of the video point to be interpolated between the P1 and the P2, the delay time increase between the P1 and the P2 and the weight applied to the delay time increase between the P1 and the P3 are applied differently. Latency interpolation method.
제 2 항에 있어서,
상기 P1과 상기 P2 사이에 존재하는 보간할 영상점이 상기 P1에 가까울수록 상기 P1과 상기 P2 사이의 지연시간 증가분의 가중치를 증가시키고, 상기 P3에 가까울수록 P1과 상기 P3 사이의 지연시간 증가분의 가중치를 증가시키는 것을 특징으로 하는 지연시간 보간 방법.
The method of claim 2,
As the image point to be interpolated between P1 and P2 is closer to P1, the weight of the delay increase between P1 and P2 increases, and the closer to P3, the weight of delay increase between P1 and P3 Delay interpolation method characterized in that to increase the.
제 1 항에 있어서,
상기 P2와 상기 P3 사이에 존재하는 영상점들의 지연시간을 보간하는 단계는,
상기 P2와 상기 P3 사이에 존재하는 보간할 영상점의 위치에 따라 상기 P2와 상기 P3 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분에 적용하는 가중치를 다르게 적용하는 것을 특징으로 하는 지연시간 보간 방법.
The method of claim 1,
Interpolating the delay time of the image points existing between the P2 and the P3,
According to the position of the video point to be interpolated between the P2 and the P3, the delay time increase between the P2 and the P3 and the weight applied to the delay time increase between the P1 and the P3 is applied differently. Latency interpolation method.
제 4 항에 있어서,
상기 P2와 상기 P3 사이에 존재하는 보간할 영상점이 상기 P2에 가까울수록 상기 P1와 상기 P3 사이의 지연시간 증가분의 가중치를 증가시키고, 상기 P3에 가까울수록 P2와 상기 P3 사이의 지연시간 증가분의 가중치를 증가시키는 것을 특징으로 하는 지연시간 보간 방법.
The method of claim 4, wherein
As the image point to be interpolated between P2 and P3 is closer to P2, the weight of the delay increase between P1 and P3 increases, and the closer to P3, the weight of delay increase between P2 and P3 Delay interpolation method characterized in that to increase the.
제 1 항에 있어서,
상기 P1과 상기 P2 사이의 지연시간 증가분은 상기 P1의 지연시간과 상기 P2의 지연시간 간의 차이를 채널 수로 나눔으로써 산출되는 것을 특징으로 하는 지연시간 보간 방법.
The method of claim 1,
The delay time interpolation method between the P1 and the P2 is calculated by dividing the difference between the delay time of the P1 and the delay time of the P2 by the number of channels.
3 개의 영상점 P1, P2, P3의 지연시간을 저장하는 레지스터;
상기 P2가 상기 P1과 상기 P3 사이에 위치할 때, 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 산출하는 증가분 산출부; 및
상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분을 이용하여 상기 P1과 상기 P2 사이에 존재하는 영상점들의 지연시간을 보간하는 지연시간 보간부를 포함하는 것을 특징으로 하는 지연시간 보간 장치.
A register for storing delay times of three image points P1, P2, and P3;
An increment calculator configured to calculate an increase in delay time between P1 and P2 and an increase in delay time between P1 and P3 when P2 is located between P1 and P3; And
And a delay time interpolation unit for interpolating the delay times of the image points existing between the P1 and the P2 by using the delay time increase between the P1 and the P2 and the delay time increase between the P1 and the P3. Delay time interpolation device.
제 7 항에 있어서,
상기 지연시간 보간부는
상기 P1과 상기 P2 사이에 존재하는 보간할 영상점의 위치에 따라 상기 P1과 상기 P2 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분에 적용하는 가중치를 다르게 적용하는 것을 특징으로 하는 지연시간 보간 장치.
The method of claim 7, wherein
The delay time interpolation unit
According to the position of the video point to be interpolated between the P1 and the P2, the delay time increase between the P1 and the P2 and the weight applied to the delay time increase between the P1 and the P3 are applied differently. Latency interpolation device.
제 8 항에 있어서,
상기 지연시간 보간부는
상기 P1과 상기 P2 사이에 존재하는 보간할 영상점이 상기 P1에 가까울수록 상기 P1과 상기 P2 사이의 지연시간 증가분의 가중치를 증가시키고, 상기 P3에 가까울수록 P1과 상기 P3 사이의 지연시간 증가분의 가중치를 증가시키는 것을 특징으로 하는 지연시간 보간 장치.
The method of claim 8,
The delay time interpolation unit
As the image point to be interpolated between P1 and P2 is closer to P1, the weight of the delay increase between P1 and P2 increases, and the closer to P3, the weight of delay increase between P1 and P3 Delay interpolation apparatus, characterized in that for increasing.
제 7 항에 있어서,
상기 지연시간 보간부는
상기 P2와 상기 P3 사이에 존재하는 보간할 영상점의 위치에 따라 상기 P2와 상기 P3 사이의 지연시간 증가분과 상기 P1과 상기 P3 사이의 지연시간 증가분에 적용하는 가중치를 다르게 적용하는 것을 특징으로 하는 지연시간 보간 장치.
The method of claim 7, wherein
The delay time interpolation unit
According to the position of the video point to be interpolated between the P2 and the P3, the delay time increase between the P2 and the P3 and the weight applied to the delay time increase between the P1 and the P3 is applied differently. Latency interpolation device.
제 10 항에 있어서,
상기 P2와 상기 P3 사이에 존재하는 보간할 영상점이 상기 P2에 가까울수록 상기 P1와 상기 P3 사이의 지연시간 증가분의 가중치를 증가시키고, 상기 P3에 가까울수록 P2와 상기 P3 사이의 지연시간 증가분의 가중치를 증가시키는 것을 특징으로 하는 지연시간 보간 장치.
11. The method of claim 10,
As the image point to be interpolated between P2 and P3 is closer to P2, the weight of the delay increase between P1 and P3 increases, and the closer to P3, the weight of delay increase between P2 and P3 Delay interpolation apparatus, characterized in that for increasing.
제 7 항에 있어서,
상기 P1과 상기 P2 사이의 지연시간 증가분은 상기 P1의 지연시간과 상기 P2의 지연시간 간의 차이를 채널 수로 나눔으로써 산출되는 것을 특징으로 하는 지연시간 보간 장치.
The method of claim 7, wherein
The delay time interpolation apparatus of claim 1, wherein the delay time increase between P1 and P2 is calculated by dividing the difference between the delay time of P1 and the delay time of P2 by the number of channels.
제 1 항 내지 제 6 항 중에 어느 한 항의 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체. A non-transitory computer-readable recording medium having recorded thereon a program for executing the method of claim 1.
KR1020100016437A 2010-02-23 2010-02-23 Method and apparatus of interpolating delay in ultrasound beamformer KR101160959B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020100016437A KR101160959B1 (en) 2010-02-23 2010-02-23 Method and apparatus of interpolating delay in ultrasound beamformer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100016437A KR101160959B1 (en) 2010-02-23 2010-02-23 Method and apparatus of interpolating delay in ultrasound beamformer

Publications (2)

Publication Number Publication Date
KR20110096903A KR20110096903A (en) 2011-08-31
KR101160959B1 true KR101160959B1 (en) 2012-06-29

Family

ID=44932000

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100016437A KR101160959B1 (en) 2010-02-23 2010-02-23 Method and apparatus of interpolating delay in ultrasound beamformer

Country Status (1)

Country Link
KR (1) KR101160959B1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101442221B1 (en) * 2012-12-28 2014-09-19 서강대학교산학협력단 Multi beam former apparatus and method using polyphase interpolation filter
KR20200100469A (en) * 2019-02-18 2020-08-26 삼성메디슨 주식회사 Analog beam former

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101590481B1 (en) * 2014-07-31 2016-02-01 서강대학교산학협력단 Delay generating method and apparatus based on LUT for dynamic receive beamforming
KR102087269B1 (en) * 2017-12-13 2020-03-10 서강대학교산학협력단 Delay apparatus and method for dynamic receive beamforming

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070232906A1 (en) 2006-03-30 2007-10-04 Aloka Co., Ltd. Delay controller for ultrasound receive beamformer
KR100901787B1 (en) 2006-12-15 2009-06-11 서강대학교기술지주 주식회사 Fractional delay filter-based beamformer apparatus using post filtering
KR100911879B1 (en) 2007-12-03 2009-08-11 서강대학교기술지주 주식회사 Ultrasound Synthetic Aperture Beamformer and Ultrasound image apparatus using the beamformer

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070232906A1 (en) 2006-03-30 2007-10-04 Aloka Co., Ltd. Delay controller for ultrasound receive beamformer
KR100901787B1 (en) 2006-12-15 2009-06-11 서강대학교기술지주 주식회사 Fractional delay filter-based beamformer apparatus using post filtering
KR100911879B1 (en) 2007-12-03 2009-08-11 서강대학교기술지주 주식회사 Ultrasound Synthetic Aperture Beamformer and Ultrasound image apparatus using the beamformer

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
J.Biomed.Eng Vol.30 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101442221B1 (en) * 2012-12-28 2014-09-19 서강대학교산학협력단 Multi beam former apparatus and method using polyphase interpolation filter
KR20200100469A (en) * 2019-02-18 2020-08-26 삼성메디슨 주식회사 Analog beam former
KR102708716B1 (en) 2019-02-18 2024-09-23 삼성메디슨 주식회사 Analog beam former

Also Published As

Publication number Publication date
KR20110096903A (en) 2011-08-31

Similar Documents

Publication Publication Date Title
KR101118515B1 (en) The apparatus of beamforming the ultrasound signal and the method using it
KR101160959B1 (en) Method and apparatus of interpolating delay in ultrasound beamformer
US9538987B2 (en) System and method for ultrasound imaging
JP4466658B2 (en) Signal processing apparatus, signal processing method, and program
KR101312309B1 (en) Apparatus and method of forming beams adaptively in ultrasound imaging
JP2017526950A (en) Audio signal resampling by interpolation for low-delay encoding / decoding
KR20160145177A (en) Ultrasonic imaging compression methods and apparatus
JPWO2009028366A1 (en) Ultrasonic imaging device
JP4982471B2 (en) Signal processing method, apparatus and program
KR20120095730A (en) Method of ultrasonic beamforming and apparatus therefor
US20040234165A1 (en) Image interpolation apparatus and method
CN102755174B (en) Beam-former, medical image system and the method for showing diagnostic image
US20130109969A1 (en) Sampling method, apparatus, probe, reception beamforming apparatus, and medical imaging system performing the sampling method
JP6261839B1 (en) Synthetic aperture radar signal processor
KR100890376B1 (en) Ultrasound diagnostic system and method for forming multiple receiving scan lines
US10743842B2 (en) Ultrasound signal processor, ultrasound signal processing method, and ultrasound diagnostic device
US20080129557A1 (en) Sampling Rate Converter With No Timing Drift and Bounded Amplitude Error
JP5450060B2 (en) Integer display of relative timing between desired output sample and corresponding input sample
JP5766896B1 (en) Peak frequency detection apparatus, method and program
KR101442221B1 (en) Multi beam former apparatus and method using polyphase interpolation filter
KR20110096443A (en) Method and apparatus of estimating ultrasound speed
KR100886932B1 (en) Multi Channel Beamformer using Single Interpolation
KR101135065B1 (en) Method and apparatus of forming multiple scanning lines
KR101334381B1 (en) A synthetic aperture beamforming system of using phase rotation beamformers
KR20160070604A (en) Phase rotation beamformer and beamforming method for envelope error compensating

Legal Events

Date Code Title Description
A201 Request for examination
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20150623

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20160602

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20170608

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20180521

Year of fee payment: 7

FPAY Annual fee payment

Payment date: 20190325

Year of fee payment: 8