CN114459397B - CT image size measuring method - Google Patents

CT image size measuring method Download PDF

Info

Publication number
CN114459397B
CN114459397B CN202210195927.XA CN202210195927A CN114459397B CN 114459397 B CN114459397 B CN 114459397B CN 202210195927 A CN202210195927 A CN 202210195927A CN 114459397 B CN114459397 B CN 114459397B
Authority
CN
China
Prior art keywords
image
screen
dimensional
point
mapping
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210195927.XA
Other languages
Chinese (zh)
Other versions
CN114459397A (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.)
First Research Institute of Ministry of Public Security
Beijing Zhongdun Anmin Analysis Technology Co Ltd
Original Assignee
First Research Institute of Ministry of Public Security
Beijing Zhongdun Anmin Analysis Technology 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 First Research Institute of Ministry of Public Security, Beijing Zhongdun Anmin Analysis Technology Co Ltd filed Critical First Research Institute of Ministry of Public Security
Priority to CN202210195927.XA priority Critical patent/CN114459397B/en
Publication of CN114459397A publication Critical patent/CN114459397A/en
Application granted granted Critical
Publication of CN114459397B publication Critical patent/CN114459397B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/60Rotation of whole images or parts thereof
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Electromagnetism (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Processing Or Creating Images (AREA)

Abstract

The invention discloses a method for measuring CT image size, which comprises the following specific processes: s1, after the three-dimensional CT image is displayed on a screen, a user can input a scaling rotation signal of the CT image through a mouse to control scaling and rotation of the CT image; s2, acquiring the distance from the center of a three-dimensional data sampling point of the image to the center of an observation point in real time according to the current rotation angle and the zoom level of the image, calculating the actual pixel of the three-dimensional CT image mapped to the screen according to the geometrical similarity principle, and calculating and refreshing the real-time display value of the graduated scale on the screen according to the actual pixel of the three-dimensional CT image mapped to the screen; s3, the interpreter measures the size of the content in the image according to the real-time display numerical value of the scale on the screen. The invention can adaptively adjust the display content of the measuring ruler according to the scaling of the three-dimensional CT image, thereby realizing the function of measuring the actual size of the content of the measured package.

Description

CT image size measuring method
Technical Field
The invention relates to the technical field of image processing, in particular to a CT image size measuring method.
Background
With the continuous improvement of the requirements of airports at home and abroad on security inspection precision, CT equipment is gradually replacing traditional perspective equipment, and three-dimensional CT images are gradually replacing two-dimensional perspective images, so that the CT equipment becomes the main basis for security inspectors to judge images. For object size measurement in an image, CN104598228A discloses a method for measuring object size, comprising: outputting a measurement area on a terminal display screen and receiving a touch instruction input by a user; scaling the measurement area in response to the touch instruction; acquiring pixel coordinates of an intersection point of the scaled measurement area; and determining the measurement data of the scaled measurement region according to the pixel coordinates. The scheme is only suitable for scaling of common two-dimensional application, and a three-dimensional image display scheme under a specific environment is not provided. In addition, in the scheme, the scaling of the measurement area is realized by responding to the touch instruction of the user, so that the realization of the method is required to depend on terminal equipment with a touch function and has certain limitation.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a measuring method of CT image size, which can adaptively adjust the display content of a measuring ruler according to the scaling of a three-dimensional CT image, thereby realizing the function of measuring the actual size of the content of the measured package.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a CT image size measuring method comprises the following specific steps:
s1, after the three-dimensional CT image is displayed on a screen, a user can input a scaling rotation signal of the CT image through a mouse to control scaling and rotation of the CT image;
s2, in the zooming rotation process of the three-dimensional CT image, the distance from the center of the three-dimensional data sampling point of the image to the center of the observation point is continuously changed, so that the distance from the center of the three-dimensional data sampling point of the image to the center of the observation point is obtained in real time according to the current rotation angle and the zooming level of the image, then the actual pixel of the three-dimensional CT image mapped to the screen is calculated according to the geometrical similarity principle, and then the real-time display value of the graduated scale on the screen is calculated and refreshed according to the actual pixel of the three-dimensional CT image mapped to the screen;
s3, the interpreter measures the size of the content in the image according to the real-time display numerical value of the scale on the screen.
Further, the specific implementation process of step S2 is as follows:
s2.1, performing volume rendering on three-dimensional CT data: starting from each pixel point on the screen, emitting a light ray along the direction of the line of sight, sampling the light ray equidistantly along the direction of the light ray when the light ray passes through the three-dimensional CT image data, and calculating the color value of the pixel point on the screen corresponding to the light ray;
s2.2, performing perspective projection: defining the region which can be observed by human eyes as a viewing cone, wherein the range between a near plane and a far plane is a visual range, the near plane is taken as a screen, the eyes of observers are recorded at point e, the observation direction is g, the distance between the screen and the eyes is d, the distance from the center of the observation point to the center of a sampling point is z, the actual side length of the content in a three-dimensional CT image is a, the sampling point of edge data of the content projects towards the point e, the intersection point position of a projection line and the screen is the position of the sampling point on the screen, and the projection length of the three-dimensional CT image on the screen is 2y s The geometric expression is as follows:
s2.3, along with the enlargement of the three-dimensional CT image, the distance between the far plane and the near plane is reduced, the corresponding z is reduced, and the corresponding y s The display length of the graduated scale on the screen is kept unchanged, and the total mapping graduation corresponding to the graduated scale is reduced; conversely, as the three-dimensional CT image is reduced, the distance between the far plane and the near plane increases, the corresponding z increases, and the corresponding y s The display length of the graduated scale on the screen is still kept unchanged, and the total mapping graduation corresponding to the graduated scale is increased;
s2.4, when the sampling data point rotates 90 degrees clockwise around the center of the sampling point, the following steps are included:
the rotation process of the image is approximately expressed as y=asinx+b, and the following is obtained:
alpha represents the angle of rotation of the image; y with rotational scaling of the three-dimensional CT image s Continuously changing according to the above method;
s2.5, the mapping scale of the graduated scale in the initial state is L, and the mapping length of the three-dimensional CT image on the screen is 2y s The method comprises the steps of carrying out a first treatment on the surface of the With scaling of the CT image, the mapping length of the three-dimensional CT image on the screen becomes 2y s ' at this time, the scale corresponds to the mapping scale
S2.6, adaptively refreshing the display content of the graduated scale on the screen according to the mapping graduation L ', wherein the maximum value of the graduated scale display interval is the mapping graduation L'.
Further, the specific implementation process of step S2 is as follows:
s2.1, performing volume rendering on three-dimensional CT data: starting from each pixel point on the screen, emitting a light ray along the direction of the line of sight, sampling the light ray equidistantly along the direction of the light ray when the light ray passes through the three-dimensional CT image data, and calculating the color value of the pixel point on the screen corresponding to the light ray;
s2.2, performing perspective projection: defining the region which can be observed by human eyes as a viewing cone, wherein the range between a near plane and a far plane is a visible range, the near plane is taken as a screen, the eyes of observers are recorded at point e, the observation direction is g, the distance between the screen and the eyes is d, the distance from the center of the observation point to the center of a sampling point is z, the actual side length of a package is a, the sampling point of package edge data projects to the point e, the intersection point position of a projection line and the screen is the position of the sampling point on the screen, and the projection length of a three-dimensional CT image on the screen is 2y s The method comprises the steps of carrying out a first treatment on the surface of the The geometric expression is as follows:
s2.3, along with the enlargement of the three-dimensional CT image, the distance between the far plane and the near plane is reduced, the corresponding z is reduced, and the corresponding y s The display length of the graduated scale on the screen is kept unchanged, and the total mapping graduation corresponding to the graduated scale is reduced; conversely, as the three-dimensional CT image is reduced, the distance between the far plane and the near plane increases, the corresponding z increases, and the corresponding y s The display length of the graduated scale on the screen is still kept unchanged, and the total mapping graduation corresponding to the graduated scale is increased;
s2.4, when the sampling data point rotates 90 degrees clockwise around the center of the sampling point, the following steps are included:
the rotation process of the image is approximately expressed as y=asinx+b, and the following is obtained:
y with rotational scaling of the three-dimensional CT image s Continuously changing according to the above method;
s2.5, recording the mapping scale of the graduated scale in the initial state as L, and the mapping length of the three-dimensional CT image on the screen as 2y s The method comprises the steps of carrying out a first treatment on the surface of the With scaling of the CT image, the mapping length of the three-dimensional CT image on the screen becomes 2y s ' at this time, the scale corresponds to the mapping scale
S2.6, adaptively refreshing the display content of the graduated scale according to the mapping graduation L ', wherein the maximum value of the graduated scale display interval is the mapping graduation L'.
The present invention also provides a computer readable storage medium having stored therein a computer program which when executed by a processor implements the above method.
The invention also provides an electronic device comprising a processor and a memory for storing a computer program; the processor is configured to implement the above-described method when executing the computer program.
The invention has the beneficial effects that: the invention can adaptively adjust the display content of the measuring ruler according to the scaling of the three-dimensional CT image, thereby realizing the function of measuring the actual size of the content of the measured package.
Drawings
FIG. 1 is a schematic flow chart of a method according to an embodiment of the invention;
FIG. 2 is a schematic diagram of the principle of volume rendering in the method of the embodiment of the present invention;
FIG. 3 is a schematic diagram of perspective projection in a method according to an embodiment of the present invention;
fig. 4 is a schematic diagram for explaining the geometry of fig. 3.
Detailed Description
The present invention will be further described with reference to the accompanying drawings, and it should be noted that, while the present embodiment provides a detailed implementation and a specific operation process on the premise of the present technical solution, the protection scope of the present invention is not limited to the present embodiment.
The embodiment provides a method for measuring the size of a CT image, as shown in fig. 1, which specifically includes the following steps:
s1, after the three-dimensional CT image is displayed on a screen, a user can input a scaling rotation signal of the CT image through a mouse to control scaling and rotation of the CT image; specifically, scaling of the three-dimensional CT image is controlled by a mouse wheel, and rotation of the three-dimensional CT image is controlled by left mouse button.
S2, in the zooming rotation process of the three-dimensional CT image, the distance from the center of the three-dimensional data sampling point of the image to the center of the observation point is continuously changed, so that the distance from the center of the three-dimensional data sampling point of the image to the center of the observation point is obtained in real time according to the current rotation angle and the zooming level of the image, then the actual pixel of the three-dimensional CT image mapped to the screen is calculated according to the geometrical similarity principle, and then the real-time display value of the graduated scale on the screen is calculated and refreshed according to the actual pixel of the three-dimensional CT image mapped to the screen.
S3, the interpreter measures the size of the content in the image according to the real-time display numerical value of the scale on the screen.
Specifically, the specific implementation process of step S2 is as follows:
s2.1, performing volume rendering on the three-dimensional CT data. The imaging principle of volume rendering is shown in fig. 2, a light ray is emitted from each pixel point on the screen along the direction of the line of sight, when the light ray passes through the three-dimensional CT image data, the light ray is sampled at equal intervals along the direction of the light ray, and the color value of the pixel point on the screen corresponding to the light ray is calculated;
s2.2, performing perspective projection, as shown in fig. 3. Perspective projection is a method of converting three-dimensional coordinates into two-dimensional coordinates in computer three-dimensional imaging. Perspective projection is the projection mode closest to the human visual system, the region that can be observed by human eyes is defined as a viewing cone, and the range between the near plane and the far plane is the visual range. In this embodiment, a near plane is understood as a screen. FIG. 4 is a geometric illustration of the perspective projection principle shown in FIG. 3, wherein the eye of the observer is at point e, the observation direction is g, the distance between the observation plane (screen) and the eye is d, the distance from the center of the observation point to the center of the sampling point is z, the actual side length of the content in the three-dimensional CT image is a (i.e., the side length of the longest side in the content), the edge data sampling point of the content projects toward point e, the intersection point position of the projection line and the observation plane (screen) is the position of the sampling point on the screen, and the projection length of the three-dimensional CT image on the screen is 2y s The geometric expression is as follows:
s2.3, with an enlargement of the three-dimensional CT image, the distance between the far and near planes in FIG. 3 decreases, corresponding to a decrease in z in FIG. 4 (a), corresponding to y s The display length (screen pixels) of the graduated scale on the screen is kept unchanged, and the total mapping scale corresponding to the graduated scale is reduced; conversely, as the three-dimensional CT image is reduced, the distance between the far plane and the near plane in FIG. 3 increases, corresponding to z in FIG. 4 (a) increasing, corresponding to y s Reduced, while still maintaining the screenThe display length (screen pixels) of the upper graduated scale is unchanged, and the total mapping graduation corresponding to the graduated scale is increased.
S2.4, as shown in FIG. 4 (b), when the sampled data point is rotated 90 degrees clockwise about the center of the sampled point, there are
The rotation process of the image is approximately expressed as y=asinx+b, and the following is obtained:
thus, α is the rotation angle of the image; y with rotational scaling of the three-dimensional CT image s Continuously changing according to the above method;
s2.5, the mapping scale of the graduated scale in the initial state is L, and the mapping length of the three-dimensional CT image on the screen is 2y s The method comprises the steps of carrying out a first treatment on the surface of the With scaling of the CT image, the mapping length of the three-dimensional CT image on the screen becomes 2y s ' at this time, the scale corresponds to the mapping scale
S2.6, adaptively refreshing the display content of the graduated scale on the screen according to the mapping graduation L ', wherein the maximum value of the graduated scale display interval is the mapping graduation L'.
Various modifications and variations of the present invention will be apparent to those skilled in the art in light of the foregoing teachings and are intended to be included within the scope of the following claims.

Claims (3)

1. The CT image size measuring method is characterized by comprising the following specific steps:
s1, after the three-dimensional CT image is displayed on a screen, a user can input a scaling rotation signal of the CT image through a mouse to control scaling and rotation of the CT image;
s2, in the zooming rotation process of the three-dimensional CT image, the distance from the center of the three-dimensional data sampling point of the image to the center of the observation point is continuously changed, so that the distance from the center of the three-dimensional data sampling point of the image to the center of the observation point is obtained in real time according to the current rotation angle and the zooming level of the image, then the actual pixel of the three-dimensional CT image mapped to the screen is calculated according to the geometrical similarity principle, and then the real-time display value of the graduated scale on the screen is calculated and refreshed according to the actual pixel of the three-dimensional CT image mapped to the screen;
s3, measuring the size of the content in the image according to the real-time display numerical value of the scale on the screen by an interpreter;
the specific implementation process of the step S2 is as follows:
s2.1, performing volume rendering on three-dimensional CT data: starting from each pixel point on the screen, emitting a light ray along the direction of the line of sight, sampling the light ray equidistantly along the direction of the light ray when the light ray passes through the three-dimensional CT image data, and calculating the color value of the pixel point on the screen corresponding to the light ray;
s2.2, performing perspective projection: defining the region which can be observed by human eyes as a viewing cone, wherein the range between a near plane and a far plane is a visual range, the near plane is taken as a screen, the eyes of observers are recorded at point e, the observation direction is g, the distance between the screen and the eyes is d, the distance from the center of the observation point to the center of a sampling point is z, the actual side length of the content in a three-dimensional CT image is a, the sampling point of edge data of the content projects towards the point e, the intersection point position of a projection line and the screen is the position of the sampling point on the screen, and the projection length of the three-dimensional CT image on the screen is 2y s The geometric expression is as follows:
s2.3, along with the enlargement of the three-dimensional CT image, the distance between the far plane and the near plane is reduced, the corresponding z is reduced, and the corresponding y s The display length of the graduated scale on the screen is kept unchanged, and the total mapping graduation corresponding to the graduated scale is reduced; on the contrary, the method comprises the steps of,with the reduction of the three-dimensional CT image, the distance between the far plane and the near plane increases, the corresponding z increases, and the corresponding y s The display length of the graduated scale on the screen is still kept unchanged, and the total mapping graduation corresponding to the graduated scale is increased;
s2.4, when the sampling data point rotates 90 degrees clockwise around the center of the sampling point, the following steps are included:
the rotation process of the image is approximately expressed as y=asinx+b, and the following is obtained:
alpha represents the angle of rotation of the image; y with rotational scaling of the three-dimensional CT image s Continuously changing according to the above method;
s2.5, the mapping scale of the graduated scale in the initial state is L, and the mapping length of the three-dimensional CT image on the screen is 2y s The method comprises the steps of carrying out a first treatment on the surface of the With scaling of the CT image, the mapping length of the three-dimensional CT image on the screen becomes 2y s ' at this time, the scale corresponds to the mapping scale
S2.6, adaptively refreshing the display content of the graduated scale on the screen according to the mapping graduation L ', wherein the maximum value of the graduated scale display interval is the mapping graduation L'.
2. A computer readable storage medium, characterized in that the computer readable storage medium has stored therein a computer program which, when executed by a processor, implements the method of claim 1.
3. An electronic device comprising a processor and a memory, the memory for storing a computer program; the processor is configured to implement the method of claim 1 when executing the computer program.
CN202210195927.XA 2022-03-01 2022-03-01 CT image size measuring method Active CN114459397B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210195927.XA CN114459397B (en) 2022-03-01 2022-03-01 CT image size measuring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210195927.XA CN114459397B (en) 2022-03-01 2022-03-01 CT image size measuring method

Publications (2)

Publication Number Publication Date
CN114459397A CN114459397A (en) 2022-05-10
CN114459397B true CN114459397B (en) 2024-02-20

Family

ID=81416525

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210195927.XA Active CN114459397B (en) 2022-03-01 2022-03-01 CT image size measuring method

Country Status (1)

Country Link
CN (1) CN114459397B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5297043A (en) * 1987-11-27 1994-03-22 Picker International, Inc. Rapid display of rotated and translated three dimensional image representations
CN101802868A (en) * 2007-09-17 2010-08-11 皇家飞利浦电子股份有限公司 A caliper for measuring objects in an image
CN104794758A (en) * 2015-04-17 2015-07-22 青岛海信医疗设备股份有限公司 Three-dimensional image cutting method
CN111681277A (en) * 2020-06-11 2020-09-18 四川大学 Method and device for measuring geometric attributes in image, electronic equipment and storage medium
CN111986088A (en) * 2020-08-28 2020-11-24 网易(杭州)网络有限公司 Image processing method, image processing device, storage medium and terminal equipment
CN113960075A (en) * 2021-10-14 2022-01-21 科大讯飞(苏州)科技有限公司 Security check equipment, article size determination method, storage medium and equipment

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5297043A (en) * 1987-11-27 1994-03-22 Picker International, Inc. Rapid display of rotated and translated three dimensional image representations
CN101802868A (en) * 2007-09-17 2010-08-11 皇家飞利浦电子股份有限公司 A caliper for measuring objects in an image
CN104794758A (en) * 2015-04-17 2015-07-22 青岛海信医疗设备股份有限公司 Three-dimensional image cutting method
CN111681277A (en) * 2020-06-11 2020-09-18 四川大学 Method and device for measuring geometric attributes in image, electronic equipment and storage medium
CN111986088A (en) * 2020-08-28 2020-11-24 网易(杭州)网络有限公司 Image processing method, image processing device, storage medium and terminal equipment
CN113960075A (en) * 2021-10-14 2022-01-21 科大讯飞(苏州)科技有限公司 Security check equipment, article size determination method, storage medium and equipment

Also Published As

Publication number Publication date
CN114459397A (en) 2022-05-10

Similar Documents

Publication Publication Date Title
CN107564089B (en) Three-dimensional image processing method, device, storage medium and computer equipment
US11275248B2 (en) Head mounted display apparatus, virtual reality display system and driving method thereof
US11176637B2 (en) Foveated rendering using eye motion
CN109040736A (en) A kind of scaling method, device, equipment and the storage medium of eye space position
CN111970454B (en) Shot picture display method, device, equipment and storage medium
CN113643414B (en) Three-dimensional image generation method and device, electronic equipment and storage medium
CN109815913B (en) Visual enhancement perception system and method based on eye movement information
US7773829B1 (en) Image-centric rulers
US20150104084A1 (en) Scanning system and image display method
Moeslund et al. A natural interface to a virtual environment through computer vision-estimated pointing gestures
US11205309B2 (en) Augmented reality system and anchor display method thereof
AU2020301254B2 (en) Sticker generating method and apparatus, and medium and electronic device
CN112738496A (en) Image processing method, apparatus, system, and computer-readable medium
CN114459397B (en) CT image size measuring method
CN111311659A (en) Calibration method based on three-dimensional imaging of oblique plane mirror
CN112308768B (en) Data processing method, device, electronic equipment and storage medium
CN114972514B (en) SLAM positioning method, SLAM positioning device, electronic equipment and readable storage medium
CN115546417A (en) Three-dimensional reconstruction method, system, electronic device and computer-readable storage medium
CN110910482B (en) Method, system and readable storage medium for video data organization and scheduling
CN114972526A (en) Method and device for measuring angle of field, computer device and medium
CN115410242A (en) Sight estimation method and device
WO2023162504A1 (en) Information processing device, information processing method, and program
EP3929878B1 (en) Uncertainty display for a multi-dimensional mesh
JP7464933B2 (en) Display device and display system
PL411337A1 (en) Method for interaction with volumetric images by means of gestures and the system for the interaction with volumetric images by means of gestures

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