CN101884544A - 图像获取方法及装置和x光拍片机 - Google Patents

图像获取方法及装置和x光拍片机 Download PDF

Info

Publication number
CN101884544A
CN101884544A CN2009101412101A CN200910141210A CN101884544A CN 101884544 A CN101884544 A CN 101884544A CN 2009101412101 A CN2009101412101 A CN 2009101412101A CN 200910141210 A CN200910141210 A CN 200910141210A CN 101884544 A CN101884544 A CN 101884544A
Authority
CN
China
Prior art keywords
interest
area
bulb
detector
image
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.)
Granted
Application number
CN2009101412101A
Other languages
English (en)
Other versions
CN101884544B (zh
Inventor
王德军
李寰中
张田田
倪先锋
K·贾布里
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to CN200910141210.1A priority Critical patent/CN101884544B/zh
Priority to EP10161544.1A priority patent/EP2250965B1/en
Priority to US12/777,414 priority patent/US20100290707A1/en
Publication of CN101884544A publication Critical patent/CN101884544A/zh
Application granted granted Critical
Publication of CN101884544B publication Critical patent/CN101884544B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
    • A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • A61B6/5241—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Human Computer Interaction (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种图像获取方法及装置和X光拍片机。所述图像获取方法包括，确定步骤：确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；计算步骤：根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；拍摄步骤：将所述球管和探测器移动到各个所述位置并在所述位置对所述感兴趣区域进行拍摄以得到若干图片；拼接步骤：将所述若干图片拼接成所述感兴趣区域的图像。所述图像获取装置包括：确定单元；计算单元；拍摄单元；拼接单元。所述X光拍片机包括所述图像获取装置。采用本发明的技术方案获得的图像更精确。

Description

图像获取方法及装置和X光拍片机
技术领域
本发明总体上涉及医用数字X光拍片机领域，尤其涉及一种图像获取方法及装置和X光的拍片机。
背景技术
在医疗领域，通常使用X光拍片机，例如X光机来对病人的感兴趣区域进行成像，医生再根据所得到的图像对病人进行诊断和医治。如图1所示，图示了一种X光机，其主要部件包括X光球管1，X光限束器2，患者固定器3，探测器4，其中X光球管1的主要功能是发射X光，X光限束器2的主要功能是限制由X光球管1所发射的X光的光野的照射范围，探测器4的功能在于接受X射线并且成像然后传到工作站进行进一步处理，患者固定器3的作用在于两点，第一点是将患者与探测器4隔离开来，从安全考虑，第二点是将患者固定从而使其在整个拍摄过程中尽量减少移动。
由该X光机可知，球管1发出X射线穿过感兴趣区域，然后到达X射线探测器4以获取感兴趣区域的图像。对于所得到的图像的尺寸而言，一般等于X射线探测器的尺寸。如果感兴趣区域的视域在X射线探测器的尺寸之内，那么在一幅图像中能够完全呈现出整个感兴趣区域。例如心脏、肺等感兴趣区域，其视域在X射线探测器的尺寸之内，因此，诸如心脏、肺之类的感兴趣区域就能够被完全呈现在一幅图片中。然后，有些感兴趣区域，其视域超出了X射线探测器的尺寸，诸如脊柱、大腿等，在一幅图片中不能将整个感兴趣区域都呈现出来，这就需要拍摄若干幅图片，然后再将这若干幅图片拼接起来以形成能够呈现出整个感兴趣区域的完整图像。
目前，对视域超出了X射线探测器尺寸的感兴趣区域而言，有几种方法来对其进行成像。其中，主要有两大类，一类是转角获取方法，其是通过转动球管，即改变球管的角度，来拍摄感兴趣区域的多幅图片，换句话说，就是当球管处于一定角度时拍摄一幅有关感兴趣区域的图片，然后再将球管转动到另一个角度再拍摄一幅有关感兴趣区域的图片，等等，直到感兴趣区域完全被包括在所有图片中为止，最后再将所有拍摄的图片拼接在一起形成感兴趣区域的图像。诸如本公司的美国专利US7177455，其采用了转动球管的角度的方式来获取感兴趣区域的图像。
该US7177455主要有下列缺点：
首先，由于需要对球管进行转角，因此其配备有球管转角定位器，该球管转角定位器的价格非常昂贵，使得采用该方法的拍片机成本高。
其次，该美国US7177455公开了所拍摄的第一图片和第二图片有重叠。如图2所示的情况，球管在球管平行移动的平面14移动并发出X射线照射病人，探测器处于探测器入射面12上，虽然第一图片和第二图片有重叠，但是在感兴趣区域平面上的感兴趣区域10并未有重叠，并且感兴趣区域平面10上的部分10并未被包括在任何图片中，这样最终得到的感兴趣区域的图像不精确。
另一类是采用球管与X光探测器并行移动的方法，即当球管与X光探测器处于第一位置时拍摄一幅图片，然后同时并行移动球管和X光探测器到第二位置，再拍摄一幅图片，等等依次并行移动球管和X光探测器直到感兴趣区域的结尾，最后再将所拍摄的图片拼接在一起形成该感兴趣区域的图像。这种图像拼接方法，是采用手动地并行移动球管和X光探测器的位置，即每拍摄完一幅图片后，由操作员手动地根据其经验将球管和X光探测器并行地移动到下一个位置。由于采用手动操作，因此工作效率低，而且由于不同的操作员的经验不同往往使得最后所得到的感兴趣区域的图像不精确。
美国US6944265，与美国US7177455类似，所公开的重叠是定义在图片平面上，即第一图片和第二图片有重叠，因此，其也存在与美国US7177455类似的会使最终得到的感兴趣区域的图像不精确的缺陷。
发明内容
本发明解决的主要问题是提供一种能够获得精确图像的图像获取方法及装置和X光机。
为了解决上述问题，本发明图像获取方法，用于通过基于X光的拍片机对病人的感兴趣区域进行成像，所述拍片机包括相对布置的球管和探测器，所述图像获取方法包括：
确定步骤：确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
计算步骤：根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄步骤：将所述球管和探测器移动到各个所述位置并在所述位置对所述感兴趣区域进行拍摄以得到若干图片；
拼接步骤：将所述若干图片拼接成所述感兴趣区域的图像。
其中，所述计算步骤进一步包括，
根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
其中，所述拼接步骤进一步包括：
把所述图片中无用信息剪切掉；
根据所述感兴趣区域在相邻两幅图片中的重叠大小确定相邻图像之间作配准时需要的搜索范围；
根据所述搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；
依据所述相对位置将相邻图像的对应像素作图像融合；
将融合后的图像进行垂直化均衡。
另外，所述重叠部分大小优选为5cm至7cm。
相应地，本发明图像获取装置，用于通过基于X光的拍片机对病人的感兴趣区域进行成像，所述拍片机包括相对布置的球管和探测器，所述图像获取装置包括：
确定单元，用于确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
计算单元，用于根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄单元，用于将所述球管和探测器移动到各个所述位置并在所述位置控制所述球管对所述感兴趣区域进行拍摄以得到若干图片；
拼接单元，用于将所述若干图片拼接成所述感兴趣区域的图像。
其中，所述计算单元进一步包括，
第一单元，用于根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
第二单元，用于根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
第三单元，用于根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
第四单元，用于根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
另外，所述拼接单元进一步包括：
剪切单元，用于把所述图片中无用信息剪切掉；
搜索范围确定单元，用于根据所述感兴趣区域在相邻两幅图片中的重叠大小确定相邻图像之间作配准时需要的搜索范围；
相对位置确定单元，用于根据所述搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；
融合单元，用于依据所述相对位置将相邻图像的对应像素作图像融合；
垂直化均衡单元，用于将融合后的图像进行垂直化均衡。
进一步地，所述重叠部分大小优选为5cm至7cm。
在本发明的另一方面，还提供了一种X光拍片机，所述X光拍片机包括相对布置的球管和探测器，其还包括图像获取装置，所述图像获取装置包括：
确定单元，用于确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
计算单元，用于根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄单元，用于将所述球管和探测器移动到各个所述位置并在所述位置对所述感兴趣区域进行拍摄以得到若干图片；
拼接单元，用于将所述若干图片拼接成所述感兴趣区域的图像。
其中，所述计算单元进一步包括，
第一单元，用于根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
第二单元，用于根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
第三单元，用于根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
第四单元，用于根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
与现有技术相比，本发明图像获取方法及装置和X光机的有益效果为：
首先，由于通过感兴趣区域将在相邻两幅图片中的重叠大小来计算得出所需拍摄的图片数目、球管和探测器移动的位置等，使得所拍摄出的各个相邻图像之间一定在感兴趣区域平面上有重叠，从而保证了诊断效果和图像拼接质量；
其次，由于本发明对于X光球管控制装置没有电动转动的要求，从而可起到降低成本；
最后，由于本发明通过确定起始位置和结束位置的方式然后自动确定曝光位置、X光视域、曝光次数等方式可以增大工作效率，节省操作者时间。
附图说明
下面通过结合附图所进行的下列描述，对于本领域技术人员来说将能更透彻地理解本发明，在附图中相同的附图标记指代相同的元素，其中：
图1图示了一种X光机的示意图；
图2图示了采用美国US7177455的方法获取感兴趣区域的图像的一种情况的示意图；
图3是本发明图像获取方法的流程图；
图4A是用球管转动方式确定感兴趣区域起始位置和结束位置的一个实施例的示意图；
图4B是用球管平行移动方式确定感兴趣区域起始位置和结束位置的一个实施例的示意图；
图4C是采用本发明技术方案得到的曝光位置与图片对应关系的示意图；
图5是图3中计算步骤的流程图；
图6是图3中的拼接步骤的流程图；
图7是确定病人的感兴趣区域起始位置的示意图；
图8是确定病人的结束位置的示意图；
图9是经过计算后球管以及探测器移动到第一个曝光位置的示意图；
图10是球管以及探测器移动到第二个曝光位置的示意图；
图11是球管以及探测器移动到第三个曝光位置的示意图；
图12是采用本发明的技术方案确保在感兴趣区域的固定重叠部分的示意图；
图13是采用本发明所获得的图片以及将所述图片拼接后的图像的图示；
图14是本发明图像获取装置的示意图。
具体实施方式
下面通过示例实施例来描述本发明的各特征及优点等。
再如图1所示，该X光机主要包括X光球管1，X光限束器2，患者固定器3和探测器4。下面以此X光机为基础来介绍本发明的技术方案。
图3图示了本发明图像获取方法的流程图。本发明图像获取方法用于通过基于X光的拍片机对病人的感兴趣区域进行成像，所述拍片机包括相对布置的球管1(见图1)和探测器4(见图1)，所述探测器4用于接收球管1发射的X光并生成图像。如图3所示，所述图像获取方法包括：
1)确定步骤：确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
2)计算步骤：根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
3)拍摄步骤：将所述球管1和探测器4移动到各个所述位置并在所述位置对所述感兴趣区域进行拍摄以得到若干图片；
4)拼接步骤：将所述若干图片拼接成所述感兴趣区域的图像。
从上述可知，本发明图像获取方法首先确定感兴趣区域起始位置、结束位置以及将在相邻两幅图片中感兴趣区域的重叠部分大小，确定感兴趣区域起始位置和结束位置可以有多种方式，如通过球管转角的方式或者球管平行移动的方式，如图4A和4B所示，其是确定感兴趣区域起始位置和结束位置的示例实施例的示意图。图4A是转角确定方式和图4B是平行移动确定方式。转角确定方式可以通过转动球管1来确定感兴趣区域平面13上感兴趣区域起始位置15和结束位置16；平行移动确定方式可以通过平行移动球管来确定感兴趣区域起始位置15和结束位置16。然后再根据所述感兴趣区域起始位置15、结束位置16以及所述重叠部分17大小计算得出所需拍摄的图片数目、视域在球管移动的方向的分量以及与每个图片相对应的球管1和探测器4的位置。确定了拍摄每一幅图片球管1和探测器4所应在的位置之后，将球管1和探测器4移动到所确定的每个位置以对感兴趣区域进行拍摄，即在每一个位置拍摄一幅图片，所拍摄的图片数目等于计算得出的所需拍摄的图片数目。如图4C所示，球管1沿球管平行移动平面14移动，探测器4沿探测器入射面12移动，其示出了曝光位置(球管和探测器位置)与所要拍摄的图片的对应关系。最后，把所拍摄的这些图片拼接在一起得到感兴趣区域的图像。
由于本发明图像获取方法的技术方案，能够保证感兴趣区域在两幅相邻图片中有重叠，而不是如美国专利US7177455那样，只是确保第一图片和第二图片有重叠，由于第一图片和第二图片有重叠并不能确保图片中的感兴趣区域有重叠，因此，采用本发明图片获取方法获得的图像更精确。
对于所述重叠部分17大小可以为5厘米至7厘米，也可以为其他的值。这是经过一系列的实验发现这一重叠大小可以使得曝光次数和图像质量达到了最好的平衡。
如图5所示，所述计算步骤进一步包括：
步骤20)根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
步骤21)根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
步骤22)根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出视域在球管移动的方向的分量；
步骤23)根据所述病人覆盖范围、所述视域在球管移动的方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
再如图6所示，所述拼接步骤4)进一步包括：
40)把所述图片中无用信息剪切掉；
41)根据所述感兴趣区域在相邻两幅图片中的重叠大小确定相邻图像之间作配准时需要的搜索范围；
42)根据所述搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；
43)依据所述相对位置将相邻图像的对应像素作图像融合；
44)将融合后的图像进行垂直化均衡。
下面以要获取病人脊柱图像为例来描述本发明图像获取方法的技术方案。如图7至12所示，球管在球管移动平面14上移动，探测器在探测器入射面12上移动，首先确定感兴趣区域起始位置topMarkedHt为1800mm；结束位置botMarkedHt为1250mm；感兴趣区域将在相邻两幅图片中的重叠大小overlap_anat为50mm；假设系统FOV(视域)最大值Hfov_prefer为250mm；
感兴趣平面13上病人覆盖范围为：
covAnatPlane＝topMarkedHt-botMarkedH＝550mm；
暂定每次X光球管及探测器移动的距离为：
DFS_tmp＝Hfov_prefer-Hfov_prefer  *(detAnatSep/acqSID)-overlap_anat，其中，detAnatSep表示从探测器入射面12到感兴趣区域平面13的距离(为常数)，acqSID表示球管焦点到探测器入射面12的垂直距离(为常数)；
球管和探测器的移动次数为：
N＝ceil((covAnatPlane-Hfov_prefer*(acqSID-detAnatSep)/acqSID)/DFS_tmp)，其中，ceil()函数表示实数向上取整；
最终的视域在球管移动方向的分量为：
VertColl＝(covAnatPlane+overlap_anat*N)/(N+(acqSID-detAnatSep)/acqSID-N*detAnatSep/acqSID)；
最终的球管以及探测器的移动距离为：
DFS＝(covAnatPlane-VertColl*(acqSID-detAnatSep)/acqSID)/N；
在探测器入射面12上的重叠大小为：
overlap＝VertColl-DFS；
最终的曝光次数为N+1；
与每幅图片相对应的球管和探测器的位置：
i从1至N+1；
location(i)＝topMarkedHt-DFS*(i-1)-(1/2)*
((acqSID-COI-detBarrierSep)/acqSID)*VertColl；
就本例而言，因为感兴趣区域起始位置topMarkedHt为1800mm；结束位置botMarkedHt为1250mm；感兴趣区域将在相邻两幅图片中的重叠大小overlap_anat为50mm；假设系统希望的视域大小在球管移动方向的分量为250mm；通过上面的计算公式，我们首先得到，球管以及探测器将移动的次数为2，然后最终的视域在球管移动方向的分量为243.75mm，接着能够得出最终的球管以及探测器的移动距离为166.67mm，最后得到最终的曝光次数为3，以及球管和探测器在每次曝光的位置分别为1691.67mm、1525mm、1358.33mm。然后，将球管和探测器移动到1691.67mm、1525mm、1358.33mm进行拍摄，如图13所示，左侧是所拍摄的三幅图片。
得到三幅图片之后，把这些图片中无用信息剪切掉，例如，把限束器2(见图1)限制范围之外的图像剪切掉；
根据所述探测器入射面上的重叠大小overlap确定相邻图像之间作配准时需要的搜索范围，大约为70mm；下面根据该搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；依据所述相对位置将相邻图像的对应像素作图像融合；然后再将融合后的图像进行垂直化均衡，从而将这三幅图片拼接成一幅图像，如图13的右侧所示。其中，所述垂直化均衡指给各个图片作规范化，使得各个图片的亮度、对比度等均衡化。
为了使所得到的图像进一步精确，还可以对图像进行进一步处理，如组织化均衡、多分辨率处理、对比度拉伸等图像增强方法。
对于拼接步骤还可以采用其他方式，例如mAs SCALING(毫安秒缩放)方法、BLENDING(混合)方法等等。
对于球管1移动方向而言，可以为水平移动、垂直移动或者成一定角度移动。
如图14所示，本发明图像获取装置包括：
确定单元100，用于确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
计算单元110，用于根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄单元120，用于将所述球管和探测器移动到各个所述位置并在所述位置控制所述球管对所述感兴趣区域进行拍摄以得到若干图片；
拼接单元130，用于将所述若干图片拼接成所述感兴趣区域的图像。
其中，所述计算单元110进一步包括，
第一单元，用于根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
第二单元，用于根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
第三单元，用于根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
第四单元，用于根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
另外，所述拼接单元130进一步包括：
剪切单元，用于把所述图片中无用信息剪切掉；
搜索范围确定单元，用于根据所述感兴趣区域平面上的重叠大小确定相邻图像之间作配准时需要的搜索范围；
相对位置确定单元，用于根据所述搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；
融合单元，用于依据所述相对位置将相邻图像的对应像素作图像融合；
垂直化均衡单元，用于将融合后的图像进行垂直化均衡。对于所述重叠部分大小而言，其可以为5厘米至7厘米或者其他值，优选为5厘米。
本发明还公开了一种X光拍片机，所述X光拍片机包括相对布置的球管和探测器，其中，所述X光拍片机还包括图像获取装置，所述图像获取装置包括：
确定单元100，用于确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
计算单元110，用于根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄单元120，用于将所述球管和探测器移动到各个所述位置并在所述位置控制所述球管对所述感兴趣区域进行拍摄以得到若干图片；
拼接单元130，用于将所述若干图片拼接成所述感兴趣区域的图像。
其中，所述计算单元110进一步包括，
第一单元，用于根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
第二单元，用于根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
第三单元，用于根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
第四单元，用于根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
另外，所述拼接单元130进一步包括：
剪切单元，用于把所述图片中无用信息剪切掉；
搜索范围确定单元，用于根据所述感兴趣区域平面上的重叠大小确定相邻图像之间作配准时需要的搜索范围；
相对位置确定单元，用于根据所述搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；
融合单元，用于依据所述相对位置将相邻图像的对应像素作图像融合；
垂直化均衡单元，用于将融合后的图像进行垂直化均衡。
综上所述，首先，由于通过感兴趣区域将在相邻两幅图片中的重叠大小来计算得出所需拍摄的图片数目、球管和探测器移动的位置等，使得所拍摄出的各个相邻图像之间一定在感兴趣区域平面上有重叠，从而保证了诊断效果和图像拼接质量；
其次，由于本发明对于X光球管控制装置没有电动转动的要求，从而可起到降低成本；
最后，由于本发明通过确定起始位置和结束位置的方式然后自动确定曝光位置、X光视域、曝光次数等方式可以增大工作效率，节省操作者时间。
虽然上述通过本发明的具体实施例对本发明的各特征进行了描述，但是对于本领域技术人员来说，在不脱离本发明的精神和范围的情况下可以对本发明进行各种修改和等价变换。所有这些修改和等价变换都意在被包括在随附的权利要求所限定的范围之内。

Claims (10)

1.一种图像获取方法，用于通过基于X光的拍片机对病人的感兴趣区域进行成像，所述拍片机包括相对布置的球管和探测器，其特征在于，所述图像获取方法包括：
确定步骤：确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
计算步骤：根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄步骤：将所述球管和探测器移动到各个所述位置并在所述位置对所述感兴趣区域进行拍摄以得到若干图片；
拼接步骤：将所述若干图片拼接成所述感兴趣区域的图像。
2.如权利要求1所述的图像获取方法，其特征在于，所述计算步骤进一步包括，
根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
3.如权利要求2所述的图像获取方法，其特征在于，所述拼接步骤进一步包括：
把所述图片中无用信息剪切掉；
根据所述感兴趣区域在相邻两幅图片中的重叠大小确定相邻图像之间作配准时需要的搜索范围；
根据所述搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；
依据所述相对位置将相邻图像的对应像素作图像融合；
将融合后的图像进行垂直化均衡。
4.如权利要求1至3所述的图像获取方法，其特征在于，所述重叠部分大小优选为5cm至7cm。
5.一种图像获取装置，用于通过基于X光的拍片机对病人的感兴趣区域进行成像，所述拍片机包括相对布置的球管和探测器，其特征在于，所述图像获取装置包括：
确定单元，用于确定感兴趣区域起始位置、结束位置以及感兴趣致区域将在相邻两幅图片中的重叠部分大小；
计算单元，用于根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄单元，用于将所述球管和探测器移动到各个所述位置并在所述位置控制所述球管对所述感兴趣区域进行拍摄以得到若干图片；
拼接单元，用于将所述若干图片拼接成所述感兴趣区域的图像。
6.如权利要求5所述的图像获取装置，其特征在于，所述计算单元进一步包括，
第一单元，用于根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
第二单元，用于根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
第三单元，用于根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
第四单元，用于根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
7.如权利要求6所述的图像获取装置，其特征在于，所述拼接单元进一步包括：
剪切单元，用于把所述图片中无用信息剪切掉；
搜索范围确定单元，用于根据所述感兴趣区域在相邻两幅图片中的重叠大小确定相邻图像之间作配准时需要的搜索范围；
相对位置确定单元，用于根据所述搜索范围计算相邻图像之间的相似度，从而确定相邻图像之间的匹配的相对位置；
融合单元，用于依据所述相对位置将相邻图像的对应像素作图像融合；
垂直化均衡单元，用于将融合后的图像进行垂直化均衡。
8.如权利要求5至7所述的图像获取装置，其特征在于，所述重叠部分大小优选为5cm至7cm。
9.一种X光拍片机，所述X光拍片机包括相对布置的球管和探测器，其特征在于，所述X光拍片机还包括图像获取装置，所述图像获取装置包括：
确定单元，用于确定感兴趣区域起始位置、结束位置以及感兴趣区域将在相邻两幅图片中的重叠部分大小；
计算单元，用于根据所述感兴趣区域起始位置、结束位置以及所述重叠部分大小计算得出所需拍摄的图片数目、视域在球管移动方向的分量以及与每个图片相对应的球管和探测器的位置；
拍摄单元，用于将所述球管和探测器移动到各个所述位置并在所述位置对所述感兴趣区域进行拍摄以得到若干图片；
拼接单元，用于将所述若干图片拼接成所述感兴趣区域的图像。
10.如权利要求9所述的X光拍片机，其特征在于，所述计算单元进一步包括，
第一单元，用于根据所述感兴趣区域起始位置和结束位置计算出在感兴趣区域平面上病人覆盖范围；
第二单元，用于根据所述病人覆盖范围、所述探测器入射面到感兴趣区域平面的距离、焦点到所述探测器入射面的距离和所述重叠部分大小计算出所需拍摄的图片数目；
第三单元，用于根据所需拍摄的图片数目、所述探测器入射面到感兴趣区域平面的距离、所述焦点到所述探测器入射面的距离和所述重叠部分大小计算得出所述视域在球管移动方向的分量；
第四单元，用于根据所述病人覆盖范围、所述视域在球管移动方向的分量、所述探测器入射面到感兴趣区域平面的距离和所述图片数目计算得出与每个图片相对应的球管和探测器的位置。
CN200910141210.1A 2009-05-12 2009-05-12 图像获取方法及装置和x光拍片机 Active CN101884544B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN200910141210.1A CN101884544B (zh) 2009-05-12 2009-05-12 图像获取方法及装置和x光拍片机
EP10161544.1A EP2250965B1 (en) 2009-05-12 2010-04-30 Image acquisition method, device and radiography system
US12/777,414 US20100290707A1 (en) 2009-05-12 2010-05-11 Image acquisition method, device and radiography system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910141210.1A CN101884544B (zh) 2009-05-12 2009-05-12 图像获取方法及装置和x光拍片机

Publications (2)

Publication Number Publication Date
CN101884544A true CN101884544A (zh) 2010-11-17
CN101884544B CN101884544B (zh) 2014-12-10

Family

ID=42395007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910141210.1A Active CN101884544B (zh) 2009-05-12 2009-05-12 图像获取方法及装置和x光拍片机

Country Status (3)

Country Link
US (1) US20100290707A1 (zh)
EP (1) EP2250965B1 (zh)
CN (1) CN101884544B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102551742A (zh) * 2010-12-15 2012-07-11 深圳迈瑞生物医疗电子股份有限公司 一种放射影像拼接装置和方法
CN103429157A (zh) * 2011-03-07 2013-12-04 爱克发医疗保健公司 放射照相成像方法和设备
CN103505229A (zh) * 2012-06-20 2014-01-15 三星电子株式会社 X射线成像设备及其控制方法
CN103873753A (zh) * 2012-12-13 2014-06-18 联想(北京)有限公司 一种拍摄方法及电子设备
CN104414660A (zh) * 2013-08-29 2015-03-18 深圳市蓝韵实业有限公司 一种dr图像获取拼接方法及系统
CN106324975A (zh) * 2016-09-09 2017-01-11 上海波城医疗科技有限公司 双x光机系统及其控制方法
CN107106097A (zh) * 2014-08-28 2017-08-29 三星电子株式会社 医疗成像装置及其操作方法
CN108652653A (zh) * 2018-05-28 2018-10-16 上海联影医疗科技有限公司 拍摄医学影像的方法、装置、医学影像系统及存储介质
CN110301934A (zh) * 2019-08-14 2019-10-08 晓智科技（成都）有限公司 基于关键点检测的待拍摄部位光野区域调节系统及方法
CN111603190A (zh) * 2019-02-25 2020-09-01 西门子医疗有限公司 借助可移动的医学x射线设备记录检查对象的全景数据组
CN113489863A (zh) * 2021-05-25 2021-10-08 深圳技术大学 一种移动拍摄系统和方法

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102106740B (zh) * 2011-03-11 2013-04-17 河海大学 X射线复式断层扫描成像系统及方法
US10098598B2 (en) 2013-06-13 2018-10-16 Samsung Electronics Co., Ltd. X-ray imaging apparatus and method for controlling the same
FI125206B (fi) * 2013-11-29 2015-07-15 Planmed Oy Anatomian osavolyymien asemointi
JP6271382B2 (ja) * 2014-09-18 2018-01-31 富士フイルム株式会社 放射線画像撮影システム、放射線画像撮影装置、制御装置、及び合成放射線画像生成方法
JP2017018160A (ja) * 2015-07-07 2017-01-26 株式会社島津製作所 Ｘ線撮影装置
CN106264587B (zh) 2016-07-25 2020-05-19 东软医疗系统股份有限公司 多序列扫描方法及装置
US10448907B2 (en) 2018-01-16 2019-10-22 Shimadzu Corporation X-ray imaging apparatus
WO2019228359A1 (en) * 2018-05-28 2019-12-05 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for taking x-ray images
US10835196B2 (en) 2019-01-24 2020-11-17 General Electric Company Method and systems for camera-aided x-ray imaging
JP7415571B2 (ja) * 2020-01-09 2024-01-17 コニカミノルタ株式会社 撮影制御装置、長尺撮影システム及びプログラム

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000049572A1 (en) * 1999-02-18 2000-08-24 Koninklijke Philips Electronics N.V. Image processing method, system and apparatus for forming an overview image of an elongated scene
US20030031290A1 (en) * 2001-08-09 2003-02-13 Naoki Sugihara X-ray computed tomographic imaging apparatus
EP1484016A1 (en) * 2003-06-03 2004-12-08 GE Medical Systems Global Technology Company LLC Method and apparatus for acquisition of a composite image with a digital detector
US20050232397A1 (en) * 2004-03-30 2005-10-20 Siemens Aktiengesellschaft Medical imaging device
US20080152088A1 (en) * 2006-12-20 2008-06-26 Xiaohui Wang Long length imaging using digital radiography

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69425416T2 (de) * 1993-11-26 2001-03-15 Koninkl Philips Electronics Nv Verfahren zur Bildzusammensetzung und Bilderzeugungsvorrichtung zur Durchführung dieses Verfahrens
US5844242A (en) * 1996-01-26 1998-12-01 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Digital mammography with a mosaic of CCD arrays
US6273606B1 (en) * 1997-12-01 2001-08-14 Agfa-Gevaert Method and assembly for recording a radiation image of an elongate body
US6694047B1 (en) * 1999-07-15 2004-02-17 General Electric Company Method and apparatus for automated image quality evaluation of X-ray systems using any of multiple phantoms
US6944265B2 (en) 2002-11-25 2005-09-13 Ge Medical Systems Global Technology Company, Llc Image pasting using geometry measurement and a flat-panel detector
US20060098897A1 (en) * 2004-11-10 2006-05-11 Agfa-Gevaert Method of superimposing images
FR2904750B1 (fr) * 2006-08-03 2008-10-17 Gen Electric Procede de reconstruction tridimensionnelle d'une enveloppe exterieure d'un corps d'un appareil a rayons x
JP5027714B2 (ja) * 2008-03-31 2012-09-19 富士フイルム株式会社 Ｘ線画像撮影装置及び方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000049572A1 (en) * 1999-02-18 2000-08-24 Koninklijke Philips Electronics N.V. Image processing method, system and apparatus for forming an overview image of an elongated scene
US20030031290A1 (en) * 2001-08-09 2003-02-13 Naoki Sugihara X-ray computed tomographic imaging apparatus
EP1484016A1 (en) * 2003-06-03 2004-12-08 GE Medical Systems Global Technology Company LLC Method and apparatus for acquisition of a composite image with a digital detector
US20050232397A1 (en) * 2004-03-30 2005-10-20 Siemens Aktiengesellschaft Medical imaging device
US20080152088A1 (en) * 2006-12-20 2008-06-26 Xiaohui Wang Long length imaging using digital radiography

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102551742B (zh) * 2010-12-15 2015-09-23 深圳迈瑞生物医疗电子股份有限公司 一种放射影像拼接装置和方法
CN102551742A (zh) * 2010-12-15 2012-07-11 深圳迈瑞生物医疗电子股份有限公司 一种放射影像拼接装置和方法
CN103429157A (zh) * 2011-03-07 2013-12-04 爱克发医疗保健公司 放射照相成像方法和设备
CN103429157B (zh) * 2011-03-07 2017-01-18 爱克发医疗保健公司 放射照相成像方法和设备
US10368827B2 (en) 2012-06-20 2019-08-06 Samsung Electronics Co., Ltd. X-ray imaging apparatus and control method thereof
CN103505229A (zh) * 2012-06-20 2014-01-15 三星电子株式会社 X射线成像设备及其控制方法
US9532763B2 (en) 2012-06-20 2017-01-03 Samsung Electronics Co., Ltd. X-ray imaging apparatus and control method thereof
US11202614B2 (en) 2012-06-20 2021-12-21 Samsung Electronics Co., Ltd. X-ray imaging apparatus and control method thereof
CN108542407A (zh) * 2012-06-20 2018-09-18 三星电子株式会社 X射线成像设备及其控制方法
CN103873753A (zh) * 2012-12-13 2014-06-18 联想(北京)有限公司 一种拍摄方法及电子设备
CN104414660A (zh) * 2013-08-29 2015-03-18 深圳市蓝韵实业有限公司 一种dr图像获取拼接方法及系统
CN107106097A (zh) * 2014-08-28 2017-08-29 三星电子株式会社 医疗成像装置及其操作方法
CN106324975B (zh) * 2016-09-09 2017-12-22 上海涛影医疗科技有限公司 双x光机系统及其控制方法
CN106324975A (zh) * 2016-09-09 2017-01-11 上海波城医疗科技有限公司 双x光机系统及其控制方法
CN108652653A (zh) * 2018-05-28 2018-10-16 上海联影医疗科技有限公司 拍摄医学影像的方法、装置、医学影像系统及存储介质
CN111603190A (zh) * 2019-02-25 2020-09-01 西门子医疗有限公司 借助可移动的医学x射线设备记录检查对象的全景数据组
CN111603190B (zh) * 2019-02-25 2023-10-27 西门子医疗有限公司 借助可移动的医学x射线设备记录检查对象的全景数据组
CN110301934A (zh) * 2019-08-14 2019-10-08 晓智科技（成都）有限公司 基于关键点检测的待拍摄部位光野区域调节系统及方法
CN110301934B (zh) * 2019-08-14 2022-11-29 晓智未来(成都)科技有限公司 基于关键点检测的待拍摄部位光野区域调节系统及方法
CN113489863A (zh) * 2021-05-25 2021-10-08 深圳技术大学 一种移动拍摄系统和方法
CN113489863B (zh) * 2021-05-25 2023-03-07 深圳技术大学 一种移动拍摄系统和方法

Also Published As

Publication number Publication date
EP2250965A1 (en) 2010-11-17
CN101884544B (zh) 2014-12-10
EP2250965B1 (en) 2020-03-25
US20100290707A1 (en) 2010-11-18

Similar Documents

Publication Publication Date Title
CN101884544B (zh) 图像获取方法及装置和x光拍片机
CN108982500B (zh) 一种宫颈液基细胞学智能辅助阅片方法和系统
CN104287756A (zh) X射线图像获取方法及装置
CN105184784B (zh) 基于运动信息的单目相机获取深度信息的方法
CN103767722A (zh) Ct或pet-ct系统及其进行扫描的定位方法
US7035450B1 (en) Method for locating an element of interest contained in a three-dimensional object, in particular during a stereotactic investigation in the X-ray examination of the breast
Piccinini et al. Automated image mosaics by non‐automated light microscopes: the MicroMos software tool
CN102053357A (zh) 用于具有增强景深的成像的系统和方法
CN106618615B (zh) 一种x光摄像定位系统、定位标识板组件和x光机系统
CN110507338B (zh) 定位方法、装置、设备及数字化x线摄影系统
JP2006141904A (ja) Ｘ線撮影装置
US20120027170A1 (en) Radiological image capturing and displaying method and apparatus
JP2002063563A (ja) 画像処理方法および画像処理システム
CN114788254A (zh) 辅助对焦方法、装置及系统
WO2018190092A1 (ja) Ｘ線検査装置
CN106023078B (zh) 一种dr影像的拼接方法
JP5149984B2 (ja) 撮像装置
KR101133503B1 (ko) 광학 및 ｘ선 융합 ｃｔ 장치와 그 데이터의 재구성 방법
JP5564385B2 (ja) 放射線画像撮影装置、放射線画像撮影方法およびプログラム
CN112241953A (zh) 基于多聚焦图像融合和hdr算法的样本图像融合方法及装置
JP2010081996A (ja) サブトラクション画像処理装置及び処理方法
Drazic et al. Optimal design and critical analysis of a high-resolution video plenoptic demonstrator
EP2285282B1 (en) Spatial orientation method of an immobilized subject
CN113362386B (zh) 确定x射线摄影系统的照射区域的方法、装置和存储介质
Al-Durgham et al. Rigorous accuracy assessment for 3D reconstruction using time-series Dual Fluoroscopy (DF) image pairs

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant