CN114863017A - Data processing method, equipment and device - Google Patents

Data processing method, equipment and device Download PDF

Info

Publication number
CN114863017A
CN114863017A CN202210457789.8A CN202210457789A CN114863017A CN 114863017 A CN114863017 A CN 114863017A CN 202210457789 A CN202210457789 A CN 202210457789A CN 114863017 A CN114863017 A CN 114863017A
Authority
CN
China
Prior art keywords
image
dicom
dicom image
image sequence
dimensional
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
CN202210457789.8A
Other languages
Chinese (zh)
Other versions
CN114863017B (en
Inventor
孙林
吴乙荣
吴海燕
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.)
Qingdao Hisense Medical Equipment Co Ltd
Original Assignee
Qingdao Hisense Medical Equipment Co Ltd
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 Qingdao Hisense Medical Equipment Co Ltd filed Critical Qingdao Hisense Medical Equipment Co Ltd
Priority to CN202210457789.8A priority Critical patent/CN114863017B/en
Publication of CN114863017A publication Critical patent/CN114863017A/en
Application granted granted Critical
Publication of CN114863017B publication Critical patent/CN114863017B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Architecture (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The present application relates to the field of data processing technologies, and in particular, to a data processing method, device and apparatus. The method comprises the steps of obtaining a plurality of medical digital imaging and communication DICOM image sequences in different stages needing three-dimensional reconstruction; determining unified pixel spacing and layer spacing corresponding to a plurality of DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; for any stage of DICOM image sequence, carrying out interpolation processing on image data in the DICOM image sequence according to position information of image reference points in the DICOM image sequence and unified pixel spacing and layer spacing to obtain a processed DICOM image sequence; and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.

Description

Data processing method, equipment and device
Technical Field
The present application relates to the field of data processing, and in particular, to a data processing method, device and apparatus.
Background
With the continuous development of hospital information platforms, the medical image cloud platform in the construction institute has very important significance for image data sharing among different departments. On the cloud platform, organs, lesions and blood vessels can be reconstructed in three dimensions according to a DICOM (Digital Imaging and Communications in Medicine) image sequence, and browsing and operation of a three-dimensional model are realized. Compared with a stand-alone workstation (client-server architecture), the cloud platform (browser-server) deploys large-scale algorithm calculation and artificial intelligence reasoning processes on a server side, and a client can access the server side only through the browser. The server side can be connected with a PACS (Picture Archiving and Communication Systems) system, and can process medical image data more conveniently and safely.
Compared with two-dimensional image information, the three-dimensional reconstruction model more intuitively reflects three-dimensional individual anatomical information, helps doctors to diagnose some diseases more accurately, such as some lung diseases and blood vessel diseases, and can help doctors to obtain higher accuracy rate clinically.
Disclosure of Invention
The embodiment of the application provides a data processing method, equipment and a device, which are used for unifying interlayer spacing and corresponding space size of DICOM image sequences in different stages.
In a first aspect, an embodiment of the present application provides a data processing method, including:
acquiring a plurality of medical digital imaging and communication DICOM image sequences in different stages needing three-dimensional reconstruction; wherein the one DICOM image sequence contains image data of a plurality of DICOM images in the same period;
determining unified pixel spacing and layer spacing corresponding to the plurality of DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
for any stage of DICOM image sequence, carrying out interpolation processing on image data in the DICOM image sequence according to position information of image reference points in the DICOM image sequence and the unified pixel spacing and layer spacing to obtain a processed DICOM image sequence;
and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
Optionally, the determining, according to the pixel interval and the layer interval corresponding to the DICOM image sequence in each period, the unified pixel interval and layer interval corresponding to the multiple DICOM image sequences in different periods specifically includes:
and selecting the minimum pixel interval and the layer interval from the pixel intervals and the layer intervals corresponding to the plurality of DICOM image sequences in different stages as unified pixel intervals and layer intervals corresponding to the plurality of DICOM image sequences in different stages.
Optionally, before performing interpolation processing on image data in the DICOM image sequence to obtain a processed DICOM image sequence, the method further includes:
aiming at a DICOM image sequence of any period, determining a corresponding diagonal point of an image reference point in a three-dimensional image space formed by the DICOM image sequence according to the image reference point corresponding to the DICOM image sequence, and the image resolution and the number of images corresponding to the DICOM image sequence;
determining target image reference points corresponding to the plurality of different stages of DICOM image sequences according to image reference points corresponding to the DICOM image sequences of each stage; determining target diagonal points corresponding to the plurality of DICOM image sequences in different stages according to the diagonal points corresponding to the image reference points;
determining the minimum three-dimensional image space corresponding to the DICOM image sequences in different phases according to the position information of the target image reference points and the position information of the target diagonal points;
the three-dimensional reconstruction according to the processed DICOM image sequence of each stage specifically includes:
and mapping the processed DICOM image sequences of each period to the minimum three-dimensional image space for three-dimensional reconstruction.
Optionally, the determining, according to the image reference point corresponding to the DICOM image sequence in each phase, the target image reference points corresponding to the multiple DICOM image sequences in different phases specifically includes:
acquiring three-dimensional coordinate information of an image reference point corresponding to the DICOM image sequence of each stage;
selecting target coordinate information with the minimum coordinate value in the coordinate information of the dimensionality of an image reference point corresponding to the DICOM image sequence of each period aiming at any dimensionality in the three-dimensional coordinate information;
and combining the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target image reference point.
Optionally, the determining, according to the diagonal point corresponding to each image reference point, a target diagonal point corresponding to the multiple DICOM image sequences in different periods includes:
acquiring three-dimensional coordinate information of diagonal points corresponding to each image reference point;
selecting target coordinate information with the maximum coordinate value in the coordinate information of the dimensionality of the diagonal point corresponding to each image reference point aiming at any dimensionality in the three-dimensional coordinate information;
and combining the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target diagonal point.
Optionally, the performing interpolation processing on the image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence and the unified pixel pitch and layer pitch to obtain a processed DICOM image sequence specifically includes:
determining the size information of a three-dimensional image space formed by the DICOM image sequence according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the interlayer spacing corresponding to the DICOM image sequence;
determining the target image resolution and the number of target images corresponding to the DICOM image sequence after interpolation processing according to the size information of a three-dimensional image space formed by the DICOM image sequence and the unified pixel spacing and layer spacing;
and carrying out interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence, the resolution of the target image and the number of the target images to obtain a processed DICOM image sequence.
Optionally, the size information of the three-dimensional image space composed of the DICOM image sequence includes a width of the three-dimensional image space, a height of the three-dimensional image space, and a depth of the three-dimensional image space; the image resolution comprises the number of pixels of an image in the width direction and the number of pixels of the image in the height direction;
the determining, according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the layer spacing corresponding to the DICOM image sequence, the size information of the three-dimensional image space formed by the DICOM image sequence specifically includes:
taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the width direction and the pixel pitch as the width of the three-dimensional image space; taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the height direction and the pixel pitch as the height of the three-dimensional image space; and taking the product of the number of images corresponding to the DICOM image sequence and the interlayer distance as the depth of the three-dimensional image space.
Optionally, the target image resolution includes a target pixel number of the image in the width direction and a target pixel number of the image in the height direction;
the determining, according to the size information of the three-dimensional image space composed of the DICOM image sequence and the unified pixel pitch and layer pitch, the resolution and number of the target images corresponding to the DICOM image sequence after interpolation processing includes:
taking the ratio of the width of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the width direction; taking the ratio of the height of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the height direction; and taking the ratio of the depth of the three-dimensional image space to the unified interlayer spacing as the number of target images corresponding to the DICOM image sequence after interpolation processing.
In a second aspect, an embodiment of the present application provides a data processing apparatus, which includes at least one processor and at least one memory; wherein the memory stores program code that, when executed by the processor, causes the processor to perform the following:
acquiring a plurality of medical digital imaging and communication DICOM image sequences in different stages needing three-dimensional reconstruction; wherein the one DICOM image sequence contains image data of a plurality of DICOM images in the same period;
determining unified pixel spacing and layer spacing corresponding to the plurality of DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
for any stage of DICOM image sequence, carrying out interpolation processing on image data in the DICOM image sequence according to position information of image reference points in the DICOM image sequence and the unified pixel spacing and layer spacing to obtain a processed DICOM image sequence;
and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
In a third aspect, an embodiment of the present application provides a data processing apparatus, including:
the system comprises an acquisition module, a processing module and a display module, wherein the acquisition module is used for acquiring a plurality of medical digital imaging and communication DICOM image sequences in different stages which need to be subjected to three-dimensional reconstruction; wherein the one DICOM image sequence contains image data of a plurality of DICOM images in the same period;
the determining module is used for determining unified pixel spacing and layer spacing corresponding to the DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
the processing module is used for carrying out interpolation processing on image data in the DICOM image sequence according to the position information of image reference points in the DICOM image sequence and the unified pixel spacing and interlayer spacing aiming at the DICOM image sequence in any period to obtain a processed DICOM image sequence; and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
In a fourth aspect, an embodiment of the present application provides a computer-readable and writable storage medium, which includes program code for causing an electronic device to perform the steps of the method of the first aspect when the program code runs on the electronic device.
In the embodiment of the application, interpolation processing is performed on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence and the unified pixel spacing and layer spacing to obtain a processed DICOM image sequence, and three-dimensional reconstruction is performed according to the processed DICOM image sequence of each stage. Based on the data processing method provided by the embodiment of the application, the pixel spacing and the layer spacing corresponding to a plurality of DICOM image sequences in different stages can be unified, and three-dimensional reconstruction is performed based on the DICOM image sequences with unified pixel spacing and layer spacing, so that after the three-dimensional model is generated, the problem that DICOM images of small organs cannot be intercepted by enclosing frames of large organs and the like to cause incapability of reasoning due to inconsistent layer spacing of two-stage images of the large organs and the small organs in the process of browsing and operating the three-dimensional model can be avoided, and the reliability of the three-dimensional model obtained after three-dimensional reconstruction is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments of the present invention will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic view of an application scenario of a data processing method according to an embodiment of the present application;
FIG. 2 is a schematic diagram of a CT two-dimensional image of a liver according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a three-dimensional reconstructed model of a liver according to an embodiment of the present application;
FIG. 4 is a flowchart of a data processing method according to an embodiment of the present application;
fig. 5 is a flowchart illustrating a method for determining a minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different phases according to an embodiment of the present disclosure;
fig. 6 is a schematic diagram of image reference points corresponding to a DICOM image sequence according to an embodiment of the present application;
fig. 7 is a schematic diagram of image reference points and target three-dimensional coordinate information corresponding to a plurality of different DICOM image sequences according to an embodiment of the present disclosure;
fig. 8 is a schematic diagram of three-dimensional coordinate information of diagonal points and target diagonal points corresponding to image reference points corresponding to a plurality of different DICOM image sequences according to an embodiment of the present application;
fig. 9 is a schematic diagram illustrating a method for determining a minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different phases according to position information of a target image reference point and position information of a target diagonal point in an embodiment of the present application;
fig. 10 is a flowchart illustrating an interpolation process for image data in a DICOM image sequence according to an embodiment of the present application;
fig. 11 is a schematic diagram illustrating a uniform pixel pitch of images corresponding to a venous phase DICOM image sequence according to an embodiment of the present application;
fig. 12 is a schematic diagram illustrating a uniform interlayer spacing of images corresponding to a venous DICOM image sequence according to an embodiment of the present application;
FIG. 13 is a flowchart illustrating an overall data processing method according to an embodiment of the present application;
fig. 14 is a schematic structural diagram of a data processing apparatus according to an embodiment of the present application;
fig. 15 is a schematic structural diagram of a data processing apparatus according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application clearer, the present application will be described in further detail with reference to the accompanying drawings, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
Some terms appearing in the present application are explained below:
1. in the embodiment of the present application, the term "and/or" describes an association relationship of associated objects, and means that there may be three relationships, for example, B and/or B, which may mean: b exists alone, B and B exist simultaneously, and B exists alone. The character "/" generally indicates that the former and latter associated objects are in an "or" relationship.
2. The term "DICOM" in the embodiments of the present application is an international standard for medical images and related information. It defines a medical image format that can be used for data exchange with a quality that meets clinical needs. DICOM is widely used in radiomedicine, cardiovascular imaging and radiodiagnosis (X-ray, nuclear magnetic resonance, ultrasound, etc.) and is increasingly used in ophthalmology, dentistry, and other medical fields. DICOM is one of the most widely deployed standards for medical information among the tens of thousands of in-use medical imaging devices. There are currently about billions of medical images that comply with the DICOM standard for clinical use.
3. The terms "arterial phase", "venous phase" and "equilibrium phase" in the embodiments of the present application generally refer to the observed vascular dissimilarity in a CT (Computed Tomography) scan or angiographic examination; wherein, the artery phase is the arterial vessel filling developing phase, and the arterial vessel is observed to be clearer; the venous phase is the venous vessel filling developing phase, and the portal vein is observed to be clearer; the equilibrium period is that all sets of blood vessels are full and developed, and the development of the whole blood vessel system is observed to be clearer.
The application scenario described in the embodiment of the present application is for more clearly illustrating the technical solution of the embodiment of the present application, and does not form a limitation on the technical solution provided in the embodiment of the present application, and it can be known by a person skilled in the art that with the occurrence of a new application scenario, the technical solution provided in the embodiment of the present application is also applicable to similar technical problems. In the description of the present application, the term "plurality" means two or more unless otherwise specified.
After three-dimensional reconstruction of organs, focuses and blood vessels, doctors need to perform multi-layer brush modification reasoning results on DICOM image sequences in different stages. The different phases of the DICOM image sequence are divided into a venous phase, an arterial phase and a balance phase; in order to improve efficiency in the three-dimensional reconstruction process, a large organ is usually inferred, then a DICOM image of a small organ is intercepted by using a bounding box range of the large organ, and corresponding inference is carried out, for example, the liver is inferred first, then DICOM data of hepatic artery or hepatic vein is intercepted, and corresponding inference is carried out.
However, due to the DICOM images screened by the user, there is a difference in DICOM image sequence periods for observing the corresponding large organ and the small organ. If the large organ and the small organ correspond to different two-stage DICOM image sequences, and the interlayer spacing of the two-stage DICOM image sequences is inconsistent, the DICOM image of the small organ cannot be intercepted by the enclosing frame of the large organ, so that the problem of incapability of reasoning is caused, and the problem cannot be solved on a single-machine edition and three-dimensional reconstruction cannot be carried out.
In order to solve the above problems, the present application provides a data processing method, which may perform interpolation processing on a plurality of different stages of DICOM image sequences, perform three-dimensional reconstruction using the processed DICOM image sequences as a data source, unify pixel pitches and layer pitches of the plurality of different stages of DICOM image sequences, and perform three-dimensional reconstruction according to the processed plurality of different stages of DICOM image sequences.
Fig. 1 is a schematic view of an application scenario of a data processing method according to an embodiment of the present application. The application scenario includes a plurality of medical imaging devices 101, a terminal device 102, and a server 103.
The medical imaging apparatus 101 includes, but is not limited to, a CT scanning apparatus in the medical field. The terminal device 102 includes, but is not limited to, various desktop computers, notebook computers, and tablet computers. The server 103 may be a server, a server cluster composed of several servers, or a cloud computing center. The server 103 may be an independent physical server, a server cluster or a distributed system formed by a plurality of physical servers, or a cloud server providing basic cloud computing services such as a cloud service, a cloud database, cloud computing, a cloud function, cloud storage, a network service, cloud communication, middleware service, a domain name service, a security service, a CDN, a big data and artificial intelligence platform, and the like.
The medical imaging equipment 101 shoots a plurality of different-stage DICOM image sequences, a user selects the plurality of different-stage DICOM image sequences needing three-dimensional reconstruction through the terminal equipment 102, and the server 103 acquires the plurality of different-stage medical digital imaging and communication DICOM image sequences needing three-dimensional reconstruction; wherein, a DICOM image sequence contains image data of a plurality of DICOM images in the same period; determining unified pixel spacing and layer spacing corresponding to a plurality of DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlamellar spacing is the distance between adjacent scanning layers when the DICOM image is obtained by scanning; for any stage of DICOM image sequence, carrying out interpolation processing on image data in the DICOM image sequence according to position information of image reference points in the DICOM image sequence and unified pixel spacing and layer spacing to obtain a processed DICOM image sequence; and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
The server 103 sends the three-dimensional model obtained after the three-dimensional reconstruction to the terminal device 102, and the terminal device 102 displays the three-dimensional model obtained through the three-dimensional reconstruction in a user interface so that a user can perform medical diagnosis according to the three-dimensional model.
The three-dimensional reconstruction of organs, focuses and blood vessels can be carried out on the cloud platform according to the DICOM image sequence, the browsing and the operation of a three-dimensional model are realized, and the three-dimensional reconstruction model more intuitively embodies three-dimensional individual anatomical information than two-dimensional image information. For example, after CT scanning, obtaining a two-dimensional CT image of the liver as shown in fig. 2 does not help a doctor to diagnose some diseases more accurately, such as some liver diseases, blood vessel diseases, etc. After three-dimensional reconstruction, a three-dimensional reconstruction model of the liver as shown in fig. 3 is obtained based on a plurality of two-dimensional images obtained by CT scanning, so that doctors can be helped to diagnose some diseases more accurately, and higher accuracy is helped to be obtained clinically.
As shown in fig. 4, a flowchart of a data processing method provided in the embodiment of the present application may specifically include the following steps:
s401, acquiring a plurality of medical digital imaging and communication DICOM image sequences in different stages which need to be subjected to three-dimensional reconstruction; wherein, a DICOM image sequence contains image data of a plurality of DICOM images in the same period;
step S402, determining unified pixel pitches and layer pitches corresponding to a plurality of DICOM image sequences in different stages according to the pixel pitches and layer pitches corresponding to the DICOM image sequences in each stage; wherein, the interlamellar spacing is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
step S403, aiming at the DICOM image sequence of any period, carrying out interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence and the unified pixel spacing and interlayer spacing to obtain a processed DICOM image sequence;
and S404, performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
In step 401, after acquiring a plurality of different stages of digital imaging and communications in medicine and communications DICOM image sequences that need to be three-dimensionally reconstructed, the embodiment of the present application may further acquire image information corresponding to the DICOM image sequence in each stage;
it should be noted that, in the embodiment of the present application, image information corresponding to multiple DICOM images in the same period is the same; the image information corresponding to the DICOM images in different stages may be the same or different.
The image information of the embodiment of the present application includes but is not limited to:
the image processing method comprises the following steps of (1) pixel spacing of DICOM images, resolution of the DICOM images, layer spacing between two adjacent DICOM images scanned in the same period, number of DICOM images in the same period, and position information of image reference points corresponding to the plurality of DICOM images in the same period;
the resolution of the DICOM image includes the number of pixels in the width direction of the DICOM image and the number of pixels in the height direction of the DICOM image.
In step 402, the embodiment of the present application may determine a unified pixel pitch and a layer pitch corresponding to a plurality of DICOM image sequences in different phases by the following methods:
and selecting the pixel interval and the layer interval corresponding to one stage DICOM image sequence from the pixel intervals and the layer intervals corresponding to a plurality of different stages of DICOM image sequences as unified pixel intervals and layer intervals corresponding to the plurality of different stages of DICOM image sequences.
In an alternative embodiment, the minimum pixel pitch and layer pitch are selected from the pixel pitches and layer pitches corresponding to the multiple different-stage DICOM image sequences as the unified pixel pitch and layer pitch corresponding to the multiple different-stage DICOM image sequences.
The pixel spacing and the interlayer spacing corresponding to a plurality of different-stage DICOM image sequences are different, the pixel spacing and the interlayer spacing corresponding to the image data of a plurality of DICOM images in the same stage in one DICOM image sequence are the same, and the smallest pixel spacing and the smallest interlayer spacing in the pixel spacing and the interlayer spacing corresponding to the plurality of different-stage DICOM image sequences are selected as the unified pixel spacing and the unified interlayer spacing corresponding to the plurality of different-stage DICOM image sequences.
For example, suppose the corresponding pixel spacing of the venous phase DICOM image sequence is 0.8mm, and the layer spacing is 1.0 mm; the pixel pitch corresponding to the artery DICOM image sequence is 1.0mm, and the interlayer pitch is 2.0 mm; the pixel spacing corresponding to the DICOM image sequence in the venous phase is 0.6mm, and the interlamellar spacing is 1.6 mm; and selecting the smallest pixel pitch and layer pitch in the pixel pitches and layer pitches corresponding to the plurality of DICOM image sequences in different stages, namely selecting the pixel pitch of 0.6mm and the layer pitch of 1.0mm as the unified pixel pitch and layer pitch corresponding to the plurality of DICOM image sequences in different stages.
In addition, in the related technology, when a multi-layer brush is used to modify a specific three-dimensional model generated by a multi-period DICOM image sequence with different space sizes, if the DICOM period used by a user when the three-dimensional model is generated is different from the corresponding space size of the DICOM period used when the multi-layer brush is used, the problem that the drawn content exceeds the space corresponding to the model easily occurs, and the canvas in the space coordinate is different, so that the multi-layer brush reports errors is caused.
Based on the above problem, the embodiment of the present application provides a method for determining a minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different phases, where the determined minimum three-dimensional image space is used as a space corresponding to a DICOM image sequence in each phase; when three-dimensional reconstruction is carried out, the DICOM image sequence of each period is mapped to the minimum three-dimensional image space for three-dimensional reconstruction.
Before step 403 is performed, the embodiment of the present application may determine the minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different phases by the following method:
as shown in fig. 5, a flowchart for determining a minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different phases provided in the embodiment of the present application may specifically include the following steps:
step 501, aiming at a DICOM image sequence of any period, determining a corresponding diagonal point of an image reference point in a three-dimensional image space formed by the DICOM image sequence according to the image reference point corresponding to the DICOM image sequence, and the image resolution and the image number corresponding to the DICOM image sequence;
step 502, determining target image reference points corresponding to a plurality of DICOM image sequences in different phases according to image reference points corresponding to the DICOM image sequences in each phase; determining target diagonal points corresponding to a plurality of DICOM image sequences in different phases according to the diagonal points corresponding to the image reference points;
step 503, determining the minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different phases according to the position information of the target image reference point and the position information of the target diagonal point.
In step 501, for a DICOM image sequence of any period, determining a corresponding diagonal point of an image reference point in a three-dimensional image space formed by the DICOM image sequence according to the image reference point corresponding to the DICOM image sequence, and the image resolution and the number of images corresponding to the DICOM image sequence;
as shown in fig. 6, a schematic diagram of an image reference point corresponding to a DICOM image sequence provided in the embodiment of the present application, for example, a point O in fig. 6 is an image reference point corresponding to a venous DICOM image sequence; the image reference point may be a point where coordinate values of three dimensions in each DICOM image of the vein-period DICOM image sequence are all minimum.
In implementation, the following formula can be used to determine the corresponding diagonal points of the image reference points of each DICOM image sequence:
diagnal=(x1+width,y1+height,z1+depth);
the digital is three-dimensional coordinate information of a diagonal point corresponding to an image reference point, (x1, y1, z1) is the three-dimensional coordinate information of the image reference point corresponding to the DICOM image sequence in the period, width is the width of the three-dimensional image space formed by the DICOM image sequence, height is the height of the three-dimensional image space formed by the DICOM image sequence, and depth is the depth of the three-dimensional image space formed by the DICOM image sequence.
In step 502, acquiring three-dimensional coordinate information of an image reference point corresponding to the DICOM image sequence of each period; selecting target coordinate information with the minimum coordinate value in the coordinate information of the dimensionality of the image reference point corresponding to the DICOM image sequence of each period aiming at any dimensionality in the three-dimensional coordinate information;
and forming the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target image reference point.
For example, assuming that three-dimensional coordinate information of image reference points corresponding to DICOM image sequences in a vein phase, an artery phase and a balance phase is obtained, three-dimensional coordinate information of an image reference point a corresponding to the DICOM image sequence in the vein phase is (5, 8, 11), three-dimensional coordinate information of an image reference point B corresponding to the DICOM image sequence in the artery phase is (3, 12, 4), three-dimensional coordinate information of an image reference point C corresponding to the DICOM image sequence in the balance phase is (6, 2, 7), and for any one dimension in the three-dimensional coordinate information, target coordinate information with the smallest coordinate value in the coordinate information of the image reference points corresponding to the DICOM image sequences in the vein phase, the artery phase and the balance phase is selected, and target coordinate information corresponding to each dimension is combined into target three-dimensional coordinate information corresponding to the target image reference point, target three-dimensional coordinate information corresponding to the target image reference point D is (3, 2,4). As shown in fig. 7, an embodiment of the present application provides a schematic diagram of image reference points and target three-dimensional coordinate information corresponding to a plurality of different DICOM image sequences.
Acquiring three-dimensional coordinate information of diagonal points corresponding to each image reference point; selecting target coordinate information with the maximum coordinate value in the coordinate information of the dimensionality of the diagonal point corresponding to each image reference point aiming at any dimensionality in the three-dimensional coordinate information; and forming target three-dimensional coordinate information corresponding to the target diagonal points by using the target coordinate information corresponding to each dimension.
For example, suppose that three-dimensional coordinate information of a diagonal point corresponding to an image reference point corresponding to a DICOM image sequence in a vein phase, an artery phase and a balance phase is obtained, three-dimensional coordinate information of a diagonal point a 'corresponding to an image reference point corresponding to the DICOM image sequence in the vein phase is a' (169, 281, 450), three-dimensional coordinate information of a diagonal point B 'corresponding to an image reference point corresponding to the DICOM image sequence in the artery phase is (180, 175, 554), three-dimensional coordinate information of a diagonal point C' corresponding to an image reference point corresponding to the DICOM image sequence in the balance phase is (376, 241, 372), and target coordinate information with the maximum coordinate value in coordinate information of a diagonal point corresponding to each image reference point is selected for any one dimension in the three-dimensional coordinate information; the target coordinate information corresponding to each dimension is composed into target three-dimensional coordinate information corresponding to a target diagonal point, and the target three-dimensional coordinate information corresponding to the target diagonal point is E (376, 281, 554), as shown in fig. 8, an embodiment of the present application provides a schematic diagram of diagonal points corresponding to image reference points corresponding to multiple DICOM image sequences in different periods and target diagonal point three-dimensional coordinate information.
In step 503, according to the position information of the target image reference point and the position information of the target diagonal point, a minimum three-dimensional image space corresponding to a plurality of different DICOM image sequences is determined, wherein the minimum three-dimensional image space is a minimum three-dimensional image space capable of surrounding the plurality of different DICOM image sequences.
Assuming that the target three-dimensional coordinate information corresponding to the target image reference point D is (3, 2, 4) and the target three-dimensional coordinate information corresponding to the target diagonal point is E (376, 281, 554), the minimum three-dimensional image space corresponding to a plurality of different DICOM image sequences is determined according to the position information of the target image reference point D and the position information of the target diagonal point E. As shown in fig. 9, an embodiment of the present application provides a schematic diagram for determining a minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different phases according to position information of a target image reference point and position information of a target diagonal point.
After determining the minimum three-dimensional image space which can surround a plurality of DICOM image sequences in different stages, carrying out interpolation processing on image data in the DICOM image sequences according to the position information of image reference points in the DICOM image sequences and the unified pixel spacing and interlayer spacing aiming at the DICOM image sequences in any stage.
As shown in fig. 10, a flowchart for performing interpolation processing on image data in a DICOM image sequence according to an embodiment of the present application may specifically include the following steps:
step 1001, determining size information of a three-dimensional image space formed by the DICOM image sequence according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the interlayer spacing corresponding to the DICOM image sequence;
step 1002, determining the resolution and the number of target images corresponding to the DICOM image sequence after interpolation processing according to the size information of a three-dimensional image space formed by the DICOM image sequence and the unified pixel spacing and interlayer spacing;
and 1003, performing interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence, the resolution of the target image and the number of the target images to obtain a processed DICOM image sequence.
When the interpolation processing is performed on the image data in the DICOM image sequence, the interpolation processing can be performed by adopting a cubic spline interpolation algorithm.
It should be noted that, in the following, in the process of interpolating the image data in the DICOM image sequence according to the unified pixel pitch and interlayer pitch, the minimum pixel pitch and interlayer pitch may be selected from the pixel pitches and interlayer pitches corresponding to the DICOM image sequences in multiple different phases, and are used as the unified pixel pitch and interlayer pitch corresponding to the DICOM image sequences in multiple different phases.
In step 1001, the size information of the three-dimensional image space composed of the DICOM image sequence includes a width of the three-dimensional image space, a height of the three-dimensional image space, and a depth of the three-dimensional image space; the image resolution comprises the number of pixels of the image in the width direction and the number of pixels of the image in the height direction;
the width of the three-dimensional image space is determined by the number of pixels and the pixel pitch of the images corresponding to the DICOM image sequence in the width direction, the number of pixels and the pixel pitch of the images corresponding to the DICOM image sequence in the height direction of the three-dimensional image space are determined by the number of images and the interlayer pitch corresponding to the DICOM image sequence, and the depth of the three-dimensional image space is determined by the number of images and the interlayer pitch corresponding to the DICOM image sequence. The corresponding pixel pitch in the width direction of the DICOM image sequence is different from the corresponding pixel pitch in the height direction of the DICOM image sequence.
Taking the product of the pixel number and the pixel pitch of the images corresponding to the DICOM image sequence in the width direction as the width of a three-dimensional image space; taking the product of the pixel number and the pixel pitch of the images corresponding to the DICOM image sequence in the height direction as the height of the three-dimensional image space; and taking the product of the number of images corresponding to the DICOM image sequence and the interlayer distance as the depth of the three-dimensional image space.
For example, the width of the three-dimensional image space is calculated according to the following formula:
width=n 1 *spacing-x;
wherein width is the width of the three-dimensional image space, n 1 The spacing-x is the corresponding pixel pitch of the DICOM image sequence in the width direction.
For example, the height of the three-dimensional image space is calculated according to the following formula:
height=n 2 *spacing-y;
where height is the height of the three-dimensional image space, n 2 The spacing-y is the corresponding pixel pitch of the DICOM image sequence in the height direction, and is the pixel number of the image corresponding to the DICOM image sequence in the height direction.
For example, the depth of the three-dimensional image space is calculated according to the following formula:
depth=slice*spacing-z;
the depth is the height of a three-dimensional image space, slice is the number of images corresponding to a DICOM image sequence, and spacing-z is the interlayer distance corresponding to the DICOM image sequence.
In step 1002, the target image resolution includes the target pixel number of the image in the width direction and the target pixel number of the image in the height direction;
taking the ratio of the width of the three-dimensional image space to the unified pixel pitch as the target pixel number of the image corresponding to the DICOM image sequence after interpolation in the width direction; taking the ratio of the height of the three-dimensional image space to the unified pixel pitch as the number of target pixels of the image corresponding to the DICOM image sequence after interpolation processing in the height direction; and taking the ratio of the depth of the three-dimensional image space to the unified interlayer spacing as the number of target images corresponding to the DICOM image sequence after interpolation processing.
For example: calculating the number of target pixels of the image corresponding to the DICOM image sequence after interpolation processing in the width direction according to the following formula:
width/spacing-X;
wherein, width is the width of the three-dimensional image space, and spacing-X is the corresponding unified pixel pitch of the DICOM image sequence in the width direction.
For example: calculating the number of target pixels of the images corresponding to the DICOM image sequence after interpolation in the height direction according to the following formula:
height/spacing-Y;
wherein, height is the height of the three-dimensional image space, and spacing-Y is the corresponding unified pixel pitch of the DICOM image sequence in the height direction.
For example: calculating the number of target images corresponding to the DICOM image sequence after interpolation processing according to the following formula:
depth/spacing-Z;
wherein depth is the height of the three-dimensional image space, and spacing-Z is the interlayer distance corresponding to the DICOM image sequence.
In step 1003, the image data in the DICOM image sequence is interpolated according to the position information of the image reference point in the DICOM image sequence, the resolution of the target image and the number of the target images, so as to obtain a processed DICOM image sequence.
For example, assuming that the pixel pitch of each image in the venous phase is 1.0mm, the unified pixel pitch is 0.5mm, which indicates that the pixel interval on each image in the original venous phase is 1.0mm, and after the pixel pitch is unified, there is one pixel every 0.5mm, using an interpolation method to select two pixels adjacent to a new pixel position corresponding to the interval of 0.5mm for interpolation to generate a gray value of the new pixel with the pixel pitch of 0.5mm, so as to obtain a venous phase image with the unified pixel pitch, as shown in fig. 11, the embodiment of the present application provides a schematic diagram of the unified pixel pitch of the images corresponding to the venous phase DICOM image sequence;
for example, assuming that the interlayer spacing of each image of the venous phase is 1.5mm, the unified interlayer spacing is 1.0mm, which indicates that there is one venous phase image every 1.5mm originally, and after the interlayer spacing is unified, there is one venous phase image every 1.0mm, using an interpolation method to select two adjacent images of the new venous phase image position corresponding to the interval of 1.0mm for interpolation, and generate a new venous phase image with the interlayer spacing of 1.0mm, thereby obtaining a venous phase image with the unified interlayer spacing, as shown in fig. 12, the embodiment of the present application provides a schematic diagram of the unified interlayer spacing of the images corresponding to the venous phase DICOM image sequence.
After interpolation processing is carried out on image data in the DICOM image sequence to obtain a processed DICOM image sequence, the processed DICOM image sequence of each period is mapped to the minimum three-dimensional image space for three-dimensional reconstruction.
Determining image reference points corresponding to the DICOM image sequences of each period to be mapped to three-dimensional coordinate information in a minimum three-dimensional image space according to image reference points corresponding to the DICOM image sequences of each period, target image reference points and unified pixel spacing and interlayer spacing;
the formula can be used: (X2, Y2, Z2) ((X1-xn)/spacing-X, (Y1-yn)/spacing-Y, (Z1-zn)/spacing-Z) determines that the image reference point corresponding to the DICOM image sequence in any period is mapped to the three-dimensional coordinate information in the minimum three-dimensional image space, wherein (X1, Y1, Z1) is the three-dimensional coordinate information of the image reference point corresponding to the DICOM image sequence in the period, (xn, yn, zn) is the three-dimensional coordinate information of the target image reference point corresponding to a plurality of different DICOM image sequences, and (X2, Y2, Z2) is the three-dimensional coordinate information of the image reference point corresponding to the DICOM image sequence in the period which is mapped to the minimum three-dimensional image space.
And mapping the processed DICOM image sequence of each stage to the minimum three-dimensional image space for three-dimensional reconstruction according to the three-dimensional coordinate information mapped to the minimum three-dimensional image space by the image reference point corresponding to the DICOM image sequence of each stage.
As shown in fig. 13, an overall flowchart of a data processing method according to an embodiment of the present application may specifically include the following steps:
step 1301, acquiring a plurality of medical digital imaging and communication DICOM image sequences in different stages needing three-dimensional reconstruction;
step 1302, selecting a minimum pixel pitch and a layer pitch from the pixel pitches and layer pitches corresponding to the plurality of different DICOM image sequences as unified pixel pitches and layer pitches corresponding to the plurality of different DICOM image sequences;
step 1303, determining corresponding diagonal points of the image reference points in a three-dimensional image space formed by the DICOM image sequence according to the image reference points corresponding to the DICOM image sequence, and the image resolution and the number of images corresponding to the DICOM image sequence aiming at the DICOM image sequence in any period;
step 1304, determining target image reference points corresponding to a plurality of DICOM image sequences in different stages according to image reference points corresponding to the DICOM image sequences in each stage; determining target diagonal points corresponding to a plurality of DICOM image sequences in different phases according to the diagonal points corresponding to the image reference points;
step 1305, determining a minimum three-dimensional image space corresponding to a plurality of DICOM image sequences in different stages according to the position information of the target image reference point and the position information of the target diagonal point;
step 1306, determining the size information of a three-dimensional image space formed by the DICOM image sequence according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the interlayer spacing corresponding to the DICOM image sequence;
step 1307, determining the resolution and the number of target images corresponding to the DICOM image sequence after interpolation processing according to the size information of the three-dimensional image space formed by the DICOM image sequence and the unified pixel spacing and interlayer spacing;
step 1308, performing interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence, the resolution of the target image and the number of the target images to obtain a processed DICOM image sequence;
and step 1309, mapping the processed DICOM image sequence of each period to the minimum three-dimensional image space for three-dimensional reconstruction.
Based on the same inventive concept, in some possible implementations, embodiments of the present application further provide a data processing device, which may include at least one processor and at least one memory. Wherein the memory stores program code which, when executed by the processor, causes the processor to perform the steps of the data processing method according to various exemplary embodiments of the present application described above in the present specification.
A data processing apparatus 1400 according to this embodiment of the present application is described below with reference to fig. 14. The data processing apparatus 1400 of fig. 14 is only an example, and should not bring any limitation to the function and the scope of use of the embodiments of the present application.
As shown in fig. 14, the data processing apparatus 1400 is embodied in the form of a general purpose computing device. The components of the data processing device 1400 may include, but are not limited to: the at least one processor 1401, the at least one memory 1402, and a bus 1403 connecting the various system components including the memory 1402 and the processor 1401.
Bus 1403 represents one or more of any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, a processor, or a local bus using any of a variety of bus architectures.
The memory 1402 may include readable media in the form of volatile memory, such as Random Access Memory (RAM)1421 or cache memory 1422, and may further include Read Only Memory (ROM) 1423.
Memory 1402 can also include a program/utility 1425 having a set (at least one) of program modules 1424, such program modules 1424 including, but not limited to: an operating system, one or more application programs, other program modules, and program data, each of which, or some combination thereof, may comprise an implementation of a network environment.
The data processing device 1400 may also communicate with one or more external devices 1404 (e.g., keyboard, pointing device, etc.), with one or more devices that enable a user to interact with the data processing device 1400, or with any device (e.g., router, modem, etc.) that enables the data processing device 1400 to communicate with one or more other computing devices. Such communication may occur via an input/output (I/O) interface 1405. Also, the data processing device 1400 may communicate with one or more networks (e.g., a Local Area Network (LAN), a Wide Area Network (WAN), or a public network such as the internet) through the network adapter 1406. As shown in fig. 14, the network adapter 1406 communicates with other modules for the data processing device 1400 through a bus 1403. It should be understood that although not shown in the figures, other hardware or software modules may be used in conjunction with the data processing device 1400, including but not limited to: microcode, device drivers, redundant processors, external disk drive arrays, RAID systems, tape drives, and data backup storage systems, among others.
Processor 1401 is specifically configured to perform the following processes:
acquiring a plurality of medical digital imaging and communication DICOM image sequences in different stages needing three-dimensional reconstruction; wherein the one DICOM image sequence contains image data of a plurality of DICOM images in the same period;
determining unified pixel spacing and layer spacing corresponding to the plurality of DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
for any stage of DICOM image sequence, carrying out interpolation processing on image data in the DICOM image sequence according to position information of image reference points in the DICOM image sequence and the unified pixel spacing and layer spacing to obtain a processed DICOM image sequence;
and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
Optionally, the processor 1401 is specifically configured to:
and selecting the minimum pixel interval and the layer interval from the pixel intervals and the layer intervals corresponding to the plurality of DICOM image sequences in different stages as unified pixel intervals and layer intervals corresponding to the plurality of DICOM image sequences in different stages.
Optionally, before the processor 1401 performs interpolation processing on image data in the DICOM image sequence to obtain a processed DICOM image sequence, the processor 1401 is further configured to:
aiming at a DICOM image sequence of any period, determining a corresponding diagonal point of an image reference point in a three-dimensional image space formed by the DICOM image sequence according to the image reference point corresponding to the DICOM image sequence, and the image resolution and the number of images corresponding to the DICOM image sequence;
determining target image reference points corresponding to the plurality of different stages of DICOM image sequences according to image reference points corresponding to the DICOM image sequences of each stage; determining target diagonal points corresponding to the plurality of DICOM image sequences in different stages according to the diagonal points corresponding to the image reference points;
determining the minimum three-dimensional image space corresponding to the plurality of different-stage DICOM image sequences according to the position information of the target image reference point and the position information of the target diagonal point;
optionally, the processor 1401 is specifically configured to:
and mapping the processed DICOM image sequences of each period to the minimum three-dimensional image space for three-dimensional reconstruction.
Optionally, the processor 1401 is specifically configured to:
acquiring three-dimensional coordinate information of an image reference point corresponding to the DICOM image sequence of each stage;
selecting target coordinate information with the minimum coordinate value in the coordinate information of the dimensionality of an image reference point corresponding to the DICOM image sequence of each period aiming at any dimensionality in the three-dimensional coordinate information;
and forming the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target image reference point.
Optionally, the processor 1401 is specifically configured to:
acquiring three-dimensional coordinate information of diagonal points corresponding to each image reference point;
selecting target coordinate information with the maximum coordinate value in the coordinate information of the dimensionality of the diagonal point corresponding to each image reference point aiming at any dimensionality in the three-dimensional coordinate information;
and forming the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target diagonal point.
Optionally, the processor 1401 is specifically configured to:
determining the size information of a three-dimensional image space formed by the DICOM image sequence according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the interlayer spacing corresponding to the DICOM image sequence;
determining the target image resolution and the number of target images corresponding to the DICOM image sequence after interpolation processing according to the size information of a three-dimensional image space formed by the DICOM image sequence and the unified pixel spacing and layer spacing;
and carrying out interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence, the resolution of the target image and the number of the target images to obtain a processed DICOM image sequence.
Optionally, the size information of the three-dimensional image space composed of the DICOM image sequence includes a width of the three-dimensional image space, a height of the three-dimensional image space, and a depth of the three-dimensional image space; the image resolution comprises the number of pixels of an image in the width direction and the number of pixels of the image in the height direction; the processor 1401 is specifically configured to:
taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the width direction and the pixel pitch as the width of the three-dimensional image space; taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the height direction and the pixel pitch as the height of the three-dimensional image space; and taking the product of the number of images corresponding to the DICOM image sequence and the interlayer distance as the depth of the three-dimensional image space.
Optionally, the target image resolution includes a target pixel number of the image in the width direction and a target pixel number of the image in the height direction; the processor 1401 is specifically configured to:
taking the ratio of the width of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the width direction; taking the ratio of the height of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the height direction; and taking the ratio of the depth of the three-dimensional image space to the unified interlayer spacing as the number of target images corresponding to the DICOM image sequence after interpolation processing.
As shown in fig. 15, an embodiment of the present application provides a data processing apparatus 1500, including:
the acquiring module 1501 is configured to acquire a plurality of different phases of digital imaging in medicine and communication DICOM image sequences that need to be three-dimensionally reconstructed; wherein the one DICOM image sequence contains image data of a plurality of DICOM images in the same period;
a determining module 1502, configured to determine unified pixel pitches and layer pitches corresponding to the multiple different stages of DICOM image sequences according to the pixel pitches and layer pitches corresponding to the DICOM image sequences in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
a processing module 1503, which performs interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence and the unified pixel pitch and interlayer pitch to obtain a processed DICOM image sequence; and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
Optionally, the determining module 1502 is specifically configured to:
according to the pixel spacing and the layer spacing corresponding to the DICOM image sequence of each stage, determining the unified pixel spacing and layer spacing corresponding to the DICOM image sequences of different stages, specifically comprising:
and selecting the minimum pixel pitch and the minimum interlayer pitch from the pixel pitches and the interlayer pitches corresponding to the DICOM image sequences in the different phases as the unified pixel pitch and the unified interlayer pitch corresponding to the DICOM image sequences in the different phases.
Optionally, before the processing module 1503 performs interpolation processing on the image data in the DICOM image sequence to obtain a processed DICOM image sequence, the determining module 1502 is further configured to:
aiming at a DICOM image sequence of any period, determining a corresponding diagonal point of an image reference point in a three-dimensional image space formed by the DICOM image sequence according to the image reference point corresponding to the DICOM image sequence, and the image resolution and the number of images corresponding to the DICOM image sequence;
determining target image reference points corresponding to the plurality of different stages of DICOM image sequences according to image reference points corresponding to the DICOM image sequences of each stage; determining target diagonal points corresponding to the plurality of DICOM image sequences in different phases according to the diagonal points corresponding to the image reference points;
determining the minimum three-dimensional image space corresponding to the plurality of different-stage DICOM image sequences according to the position information of the target image reference point and the position information of the target diagonal point;
optionally, the processing module 1503 is specifically configured to:
and mapping the processed DICOM image sequences of each period to the minimum three-dimensional image space for three-dimensional reconstruction.
Optionally, the determining module 1502 is specifically configured to:
acquiring three-dimensional coordinate information of an image reference point corresponding to the DICOM image sequence of each stage;
selecting target coordinate information with the minimum coordinate value in the coordinate information of the dimensionality of an image reference point corresponding to the DICOM image sequence of each period aiming at any dimensionality in the three-dimensional coordinate information;
and forming the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target image reference point.
Optionally, the determining module 1502 is specifically configured to:
acquiring three-dimensional coordinate information of diagonal points corresponding to each image reference point;
selecting target coordinate information with the maximum coordinate value in the coordinate information of the dimensionality of the diagonal point corresponding to each image reference point aiming at any dimensionality in the three-dimensional coordinate information;
and forming the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target diagonal point.
Optionally, the processing module 1503 is specifically configured to:
determining the size information of a three-dimensional image space formed by the DICOM image sequence according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the interlayer spacing corresponding to the DICOM image sequence;
determining the target image resolution and the number of target images corresponding to the DICOM image sequence after interpolation processing according to the size information of a three-dimensional image space formed by the DICOM image sequence and the unified pixel spacing and layer spacing;
and carrying out interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence, the resolution of the target image and the number of the target images to obtain a processed DICOM image sequence.
Optionally, the size information of the three-dimensional image space composed of the DICOM image sequence includes a width of the three-dimensional image space, a height of the three-dimensional image space, and a depth of the three-dimensional image space; the image resolution comprises the number of pixels of an image in the width direction and the number of pixels of the image in the height direction; the processing module 1503 is specifically configured to:
taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the width direction and the pixel pitch as the width of the three-dimensional image space; taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the height direction and the pixel pitch as the height of the three-dimensional image space; and taking the product of the number of images corresponding to the DICOM image sequence and the interlayer distance as the depth of the three-dimensional image space.
Optionally, the target image resolution includes a target pixel number of the image in the width direction and a target pixel number of the image in the height direction; the processing module 1503 is specifically configured to:
taking the ratio of the width of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the width direction; taking the ratio of the height of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the height direction; and taking the ratio of the depth of the three-dimensional image space to the unified interlayer spacing as the number of target images corresponding to the DICOM image sequence after interpolation processing.
In an exemplary embodiment, a computer-readable storage medium comprising instructions, such as a memory comprising instructions, executable by a processor to perform the above-described data processing method is also provided. Alternatively, the storage medium may be a non-transitory computer readable storage medium, which may be, for example, a ROM, a Random Access Memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, and the like.
In an exemplary embodiment, a computer program product is also provided, comprising a computer program which, when executed by a processor, implements any of the data processing methods as provided herein.
In an exemplary embodiment, aspects of a data processing method provided in the present application may also be implemented in the form of a program product including program code for causing a computer device to perform the steps in the data processing method according to various exemplary embodiments of the present application described above in this specification when the program product is run on the computer device.
The program product may employ any combination of one or more readable media. The readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the readable storage medium include: an electrical connection having one or more wires, a portable disk, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.
The program product for the data processing method of the embodiments of the present application may employ a portable compact disc read only memory (CD-ROM) and include program code, and may be executed on an electronic device. However, the program product of the present application is not limited thereto, and in this document, a readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
A readable signal medium may include a propagated data signal with readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A readable signal medium may also be any readable medium that is not a readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
Program code embodied on a readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
Program code for carrying out operations of the present application may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, C + + or the like and conventional procedural programming languages, such as the "for example" programming language or similar programming languages. The program code may execute entirely on the consumer electronic device, partly on the consumer electronic device, as a stand-alone software package, partly on the consumer electronic device and partly on a remote electronic device, or entirely on the remote electronic device or server. In the case of remote electronic devices, the remote electronic devices may be connected to the consumer electronic device through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external electronic device (for example, through the internet using an internet service provider).
It should be noted that although several units or sub-units of the apparatus are mentioned in the above detailed description, such division is merely exemplary and not mandatory. Indeed, the features and functions of two or more units described above may be embodied in one unit, according to embodiments of the application. Conversely, the features and functions of one unit described above may be further divided into embodiments by a plurality of units.
Further, while the operations of the methods of the present application are depicted in the drawings in a particular order, this does not require or imply that these operations must be performed in this particular order, or that all of the illustrated operations must be performed, to achieve desirable results. Additionally or alternatively, certain steps may be omitted, multiple steps combined into one step execution, and/or one step broken down into multiple step executions.

Claims (10)

1. A method of data processing, the method comprising:
acquiring a plurality of medical digital imaging and communication DICOM image sequences in different stages needing three-dimensional reconstruction; wherein the one DICOM image sequence contains image data of a plurality of DICOM images in the same period;
determining unified pixel spacing and layer spacing corresponding to the DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
for any stage of DICOM image sequence, carrying out interpolation processing on image data in the DICOM image sequence according to position information of image reference points in the DICOM image sequence and the unified pixel spacing and layer spacing to obtain a processed DICOM image sequence;
and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
2. The method according to claim 1, wherein the determining unified pixel pitches and layer pitches corresponding to the plurality of DICOM image sequences in different phases according to the pixel pitches and layer pitches corresponding to the DICOM image sequence in each phase specifically comprises:
and selecting the minimum pixel interval and the layer interval from the pixel intervals and the layer intervals corresponding to the plurality of DICOM image sequences in different stages as unified pixel intervals and layer intervals corresponding to the plurality of DICOM image sequences in different stages.
3. The method of claim 1, wherein before interpolating the image data in the DICOM image sequence to obtain a processed DICOM image sequence, the method further comprises:
aiming at a DICOM image sequence of any period, determining a corresponding diagonal point of an image reference point in a three-dimensional image space formed by the DICOM image sequence according to the image reference point corresponding to the DICOM image sequence, and the image resolution and the number of images corresponding to the DICOM image sequence;
determining target image reference points corresponding to the plurality of different stages of DICOM image sequences according to image reference points corresponding to the DICOM image sequences of each stage; determining target diagonal points corresponding to the plurality of DICOM image sequences in different stages according to the diagonal points corresponding to the image reference points;
determining the minimum three-dimensional image space corresponding to the plurality of different-stage DICOM image sequences according to the position information of the target image reference point and the position information of the target diagonal point;
the three-dimensional reconstruction according to the processed DICOM image sequence of each stage specifically includes:
and mapping the processed DICOM image sequences of each period to the minimum three-dimensional image space for three-dimensional reconstruction.
4. The method according to claim 3, wherein the determining the target image reference points corresponding to the plurality of different stages of DICOM image sequences according to the image reference points corresponding to each stage of DICOM image sequences specifically comprises:
acquiring three-dimensional coordinate information of an image reference point corresponding to the DICOM image sequence of each stage;
selecting target coordinate information with the minimum coordinate value in the coordinate information of the dimensionality of an image reference point corresponding to the DICOM image sequence of each period aiming at any dimensionality in the three-dimensional coordinate information;
and forming the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target image reference point.
5. The method according to claim 3, wherein the determining the target diagonal points corresponding to the plurality of different DICOM image sequences according to the diagonal points corresponding to each image reference point specifically comprises:
acquiring three-dimensional coordinate information of diagonal points corresponding to each image reference point;
selecting target coordinate information with the maximum coordinate value in the coordinate information of the dimensionality of the diagonal point corresponding to each image reference point aiming at any dimensionality in the three-dimensional coordinate information;
and forming the target coordinate information corresponding to each dimension into target three-dimensional coordinate information corresponding to the target diagonal point.
6. The method according to claim 3, wherein the interpolating the image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence and the unified pixel pitch and layer pitch to obtain a processed DICOM image sequence specifically includes:
determining the size information of a three-dimensional image space formed by the DICOM image sequence according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the interlayer spacing corresponding to the DICOM image sequence;
determining the target image resolution and the number of target images corresponding to the DICOM image sequence after interpolation processing according to the size information of a three-dimensional image space formed by the DICOM image sequence and the unified pixel spacing and layer spacing;
and carrying out interpolation processing on image data in the DICOM image sequence according to the position information of the image reference points in the DICOM image sequence, the resolution of the target image and the number of the target images to obtain a processed DICOM image sequence.
7. The method of claim 6, wherein the dimension information of the three-dimensional image space composed of the DICOM image sequence comprises a width of the three-dimensional image space, a height of the three-dimensional image space, and a depth of the three-dimensional image space; the image resolution comprises the number of pixels of an image in the width direction and the number of pixels of the image in the height direction;
the determining, according to the image resolution and the number of images corresponding to the DICOM image sequence and the pixel spacing and the layer spacing corresponding to the DICOM image sequence, the size information of the three-dimensional image space formed by the DICOM image sequence specifically includes:
taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the width direction and the pixel pitch as the width of the three-dimensional image space; taking the product of the number of pixels of the images corresponding to the DICOM image sequence in the height direction and the pixel pitch as the height of the three-dimensional image space; and taking the product of the number of images corresponding to the DICOM image sequence and the interlayer distance as the depth of the three-dimensional image space.
8. The method of claim 6, wherein the target image resolution comprises a target number of pixels of the image in a width direction and a target number of pixels of the image in a height direction;
the determining, according to the size information of the three-dimensional image space composed of the DICOM image sequence and the unified pixel pitch and layer pitch, the resolution and number of the target images corresponding to the DICOM image sequence after interpolation processing includes:
taking the ratio of the width of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the width direction; taking the ratio of the height of the three-dimensional image space to the unified pixel pitch as the number of target pixels of an image corresponding to the DICOM image sequence after interpolation processing in the height direction; and taking the ratio of the depth of the three-dimensional image space to the unified interlayer spacing as the number of target images corresponding to the DICOM image sequence after interpolation processing.
9. A data processing apparatus, characterized in that the apparatus comprises at least one processor, and at least one memory; wherein the memory stores program code that, when executed by the processor, causes the processor to perform the following:
acquiring a plurality of different stages of medical digital imaging and communication DICOM image sequences needing three-dimensional reconstruction; wherein the one DICOM image sequence contains image data of a plurality of DICOM images in the same period;
determining unified pixel spacing and layer spacing corresponding to the plurality of DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
for any stage of DICOM image sequence, carrying out interpolation processing on image data in the DICOM image sequence according to position information of image reference points in the DICOM image sequence and the unified pixel spacing and layer spacing to obtain a processed DICOM image sequence;
and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
10. A data processing apparatus, comprising:
the system comprises an acquisition module, a processing module and a display module, wherein the acquisition module is used for acquiring a plurality of different medical digital imaging and communication DICOM image sequences which need to be subjected to three-dimensional reconstruction; wherein the DICOM image sequence contains image data of a plurality of DICOM images in the same period;
the determining module is used for determining unified pixel spacing and layer spacing corresponding to the plurality of DICOM image sequences in different stages according to the pixel spacing and layer spacing corresponding to the DICOM image sequence in each stage; wherein, the interlayer distance is the distance between adjacent scanning layers when the DICOM image is obtained by scanning;
the processing module is used for carrying out interpolation processing on image data in the DICOM image sequence according to the position information of image reference points in the DICOM image sequence and the unified pixel spacing and interlayer spacing aiming at the DICOM image sequence in any period to obtain a processed DICOM image sequence; and performing three-dimensional reconstruction according to the processed DICOM image sequence of each period.
CN202210457789.8A 2022-04-27 2022-04-27 Data processing method, device and apparatus Active CN114863017B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210457789.8A CN114863017B (en) 2022-04-27 2022-04-27 Data processing method, device and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210457789.8A CN114863017B (en) 2022-04-27 2022-04-27 Data processing method, device and apparatus

Publications (2)

Publication Number Publication Date
CN114863017A true CN114863017A (en) 2022-08-05
CN114863017B CN114863017B (en) 2024-07-12

Family

ID=82633192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210457789.8A Active CN114863017B (en) 2022-04-27 2022-04-27 Data processing method, device and apparatus

Country Status (1)

Country Link
CN (1) CN114863017B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106652017A (en) * 2016-09-18 2017-05-10 青岛海信医疗设备股份有限公司 Method and apparatus for determining the completeness of DICOM image file in three-dimensional reconstruction
US20170278281A1 (en) * 2016-03-25 2017-09-28 L-3 Communications Security & Detection Systems, Inc. Systems and methods for reconstructing projection images from computed tomography volumes
CN108053433A (en) * 2017-11-28 2018-05-18 浙江工业大学 A kind of multi-modal arteria carotis MRI method for registering based on physical alignment and outline
CN110738597A (en) * 2018-07-19 2020-01-31 北京连心医疗科技有限公司 Size self-adaptive preprocessing method of multi-resolution medical image in neural network
JP2020101988A (en) * 2018-12-21 2020-07-02 大日本印刷株式会社 Three-dimensional reconstructed image display device, three-dimensional reconstructed image display method, program, and image generation method
WO2020173052A1 (en) * 2019-02-28 2020-09-03 未艾医疗技术(深圳)有限公司 Three-dimensional image measurement method, electronic device, storage medium, and program product

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170278281A1 (en) * 2016-03-25 2017-09-28 L-3 Communications Security & Detection Systems, Inc. Systems and methods for reconstructing projection images from computed tomography volumes
CN106652017A (en) * 2016-09-18 2017-05-10 青岛海信医疗设备股份有限公司 Method and apparatus for determining the completeness of DICOM image file in three-dimensional reconstruction
CN108053433A (en) * 2017-11-28 2018-05-18 浙江工业大学 A kind of multi-modal arteria carotis MRI method for registering based on physical alignment and outline
CN110738597A (en) * 2018-07-19 2020-01-31 北京连心医疗科技有限公司 Size self-adaptive preprocessing method of multi-resolution medical image in neural network
JP2020101988A (en) * 2018-12-21 2020-07-02 大日本印刷株式会社 Three-dimensional reconstructed image display device, three-dimensional reconstructed image display method, program, and image generation method
WO2020173052A1 (en) * 2019-02-28 2020-09-03 未艾医疗技术(深圳)有限公司 Three-dimensional image measurement method, electronic device, storage medium, and program product

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孟齐文;赵卫东;王志成;: "医学图像快速插值算法的设计与实现", 计算机仿真, no. 02, 15 February 2011 (2011-02-15) *

Also Published As

Publication number Publication date
CN114863017B (en) 2024-07-12

Similar Documents

Publication Publication Date Title
US20210210179A1 (en) Evolving contextual clinical data engine for medical information
US9754390B2 (en) Reconstruction of time-varying data
JP6275876B2 (en) An evolutionary contextual clinical data engine for medical data processing
CN111540025B (en) Predicting images for image processing
JP2019530488A (en) Computer-aided diagnostic system for medical images using deep convolutional neural networks
AU2017290600A1 (en) Systems and methods for anonymization of health data and transmition of health data for analysis across geographic regions
US10083511B2 (en) Angiographic roadmapping mask
US11037659B2 (en) Data-enriched electronic healthcare guidelines for analytics, visualization or clinical decision support
CA3131877A1 (en) A platform for evaluating medical information and method for using the same
US10679350B2 (en) Method and apparatus for adjusting a model of an anatomical structure
WO2015106339A1 (en) Method and system for generating pre-scaled images for a series of mammography images
US8639008B2 (en) Mobile architecture using cloud for data mining application
CN107170021B (en) Refined reconstruction of time-varying data
CN114863017B (en) Data processing method, device and apparatus
CN115861283A (en) Medical image analysis method, device, equipment and storage medium
US11170889B2 (en) Smooth image scrolling
US10977792B2 (en) Quantitative evaluation of time-varying data
CN113823385B (en) Method, device, equipment and medium for modifying DICOM image
US20230154067A1 (en) Output Validation of an Image Reconstruction Algorithm
US12100075B2 (en) Image reconstruction by modeling image formation as one or more neural networks
US20220215601A1 (en) Image Reconstruction by Modeling Image Formation as One or More Neural Networks
CN117541742A (en) Image processing method, device, computing equipment and storage medium
WO2024171153A1 (en) Method and system for determining a deformation of a blood vessel
Torres Bisbal Web-based visualisation of MRI cardiac images

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant