Magnetic resonance image correction methods, apparatus, computer equipment, and readable storage media

By acquiring the sensitivity and positional relationship of the receiving coil of the magnetic resonance device, a correction image is generated to correct the brightness of the initial magnetic resonance image, which solves the problem of image brightness non-uniformity under high field and improves image quality.

CN115685033BActive Publication Date: 2026-06-30SHANGHAI UNITED IMAGING HEALTHCARE

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
SHANGHAI UNITED IMAGING HEALTHCARE
Filing Date
2021-07-30
Publication Date
2026-06-30

AI Technical Summary

Technical Problem

In high-field magnetic resonance imaging, the non-uniformity of the radio frequency receiving field sensitivity leads to non-uniform image brightness, which affects diagnosis and image post-processing analysis. Existing correction methods are ineffective or lose contrast under high field conditions.

Method used

By acquiring the receiving field sensitivity of the receiving coil in the magnetic resonance device, the positional relationship between the coil and the object being detected is determined. Based on the positional relationship and sensitivity, a correction image is generated, and the brightness of the initial magnetic resonance image is corrected.

Benefits of technology

It improves the uniformity of image brightness distribution, solves the problem of uneven image brightness in high-field conditions caused by traditional methods, and improves the accuracy of diagnosis.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN115685033B_ABST
    Figure CN115685033B_ABST
Patent Text Reader

Abstract

This application relates to a method, apparatus, computer device, and readable storage medium for magnetic resonance image correction. The method includes: acquiring the receiving field sensitivity of a receiving coil in a magnetic resonance device; determining the positional relationship between the receiving coil in the magnetic resonance device and the region of interest (ROI) of the target object; determining the receiving signal sensitivity of each point in the imaging region where the ROI is located based on the positional relationship and the receiving field sensitivity; generating a corrected image of the corresponding imaging region based on the receiving signal sensitivity of each point; and correcting the initial magnetic resonance image of the ROI using the corrected image to obtain the target magnetic resonance image. This method can improve the uniformity of brightness distribution in high-field magnetic resonance images.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This application relates to the field of medical image processing technology, and in particular to a magnetic resonance image correction method, apparatus, computer device, and readable storage medium. Background Technology

[0002] Magnetic Resonance Imaging (MRI) is an imaging technique that reconstructs images by utilizing the signals generated by the resonance of hydrogen nuclei within a main magnetic field. It primarily uses a radio frequency (RF) transmitting coil to generate RF pulses. These RF pulses excite atomic nuclei with non-zero spin placed in the main magnetic field. After the RF pulse stops, the nuclei relax. During this relaxation process, a receiving coil acquires the signal. The acquired signal undergoes analog-to-digital conversion, filtering, Fourier transform, and other reconstruction processes to form a magnetic resonance image reflecting the tissue structure of the examined object. With the continuous improvement of medical technology, MRI is increasingly used in clinical diagnosis, offering significant advantages for disease diagnosis: it can directly produce tomographic images in transverse, sagittal, coronal, and various oblique planes; it does not produce artifacts as seen in CT scans; angiography does not require contrast agents; and it involves no ionizing radiation, having no adverse effects on the body, among many other advantages.

[0003] In recent years, clinical needs have driven the development of magnetic resonance imaging (MRI) technology towards higher signal-to-noise ratios, higher resolutions and contrasts, and faster imaging speeds. Major manufacturers have increased the MRI capabilities of their equipment from 3T (Tesla) to 5T or even higher. However, in high magnetic field environments, the image brightness information received by the MRI system using local coils contains information about the expected non-uniform distribution of the receiving field sensitivity, contrast information weighted by different tissue structures, and non-uniformity of the emission field intensity. These three types of brightness information are intertwined, making it difficult to measure the receiving field sensitivity distribution independently. Therefore, calibrating the receiving field sensitivity separately in high fields presents certain challenges.

[0004] Given the non-uniformity in the sensitivity distribution of the radio frequency (RF) receiving field in a magnetic resonance imaging (MRI) system, which can lead to uneven brightness in MRI images and affect doctors' interpretation and diagnosis, as well as post-processing analysis such as image segmentation, calibration of the RF receiving field sensitivity is essential. However, as the magnetic resonance field strength increases, the non-uniformity of the RF receiving field sensitivity becomes more severe, thus increasing the difficulty of uniformity calibration.

[0005] One existing method for correcting image intensity variations caused by radio frequency (RF) receiver sensitivity involves: using a pre-scanned image of a volume coil with the same field of view (FOV) as a reference; obtaining a receiver sensitivity distribution map of the receiving coil to be corrected; and using the reference map to perform uniformity correction on the receiver sensitivity distribution map. More specifically, each pixel in the receiver sensitivity distribution map can be divided by its corresponding pixel in the reference map, thus obtaining the ratio of the receiver sensitivity distribution of the receiving coil based on the same transmission field information and contrast information in the two scans. However, this method is based on the assumption of a uniform receiver sensitivity distribution when the volume coil has receiving functionality. It shows good correction results under low field strength conditions. However, under high field strength conditions, due to the increased RF frequency and shorter wavelength, the dielectric effect with the human body is enhanced, leading to a decrease in the uniformity of the RF field of the volume coil itself. The image generated by the volume coil cannot meet the uniformity requirements of the reference map.

[0006] Another existing method for correcting image intensity changes caused by RF receiver sensitivity is based on extracting coil sensitivity distribution from the image itself, such as the B-spline offset field correction algorithm, the Discrete Cosine Transform (DCT) coefficient offset field correction algorithm, and the polynomial offset field correction algorithm. These methods cannot effectively distinguish between contrast and brightness uniformity; while correcting brightness, there is a certain loss in contrast. Summary of the Invention

[0007] Therefore, it is necessary to provide a magnetic resonance image correction method, apparatus, computer equipment, and readable storage medium that can improve the uniformity of image brightness distribution, in order to address the above-mentioned technical problems.

[0008] A magnetic resonance image correction method, the method comprising:

[0009] To obtain the receiving field sensitivity of the receiving coil in a magnetic resonance imaging (MRI) device;

[0010] Determine the positional relationship between the receiving coil and the region of interest of the object being detected in the magnetic resonance imaging (MRI) device;

[0011] Based on the positional relationship and the receiving field sensitivity, the receiving signal sensitivity of each point in the imaging region where the region of interest is located is determined.

[0012] Based on the received signal sensitivity at each location point, a corrected image of the corresponding imaging area is generated;

[0013] By correcting the image, the initial magnetic resonance image of the region of interest is corrected to obtain the target magnetic resonance image.

[0014] In one embodiment, determining the positional relationship between the receiving coil in the magnetic resonance apparatus and the region of interest of the object being detected includes:

[0015] Determine the positional relationships between each single-coil loop in the receiving coil and the region of interest of the detected object;

[0016] Based on the positional relationship and receiving field sensitivity, the receiving signal sensitivity of each point in the imaging region where the region of interest is located is determined, including:

[0017] Based on the positional relationships and receiving field sensitivity, the receiving signal sensitivity of each single-coil circuit at each position point in the imaging region is determined.

[0018] In one embodiment, determining the positional relationship between each single-coil loop in the receiving coil and the region of interest of the detected object includes:

[0019] The positioning coil is determined from the receiving coil, and the positional distance between the positioning coil and the region of interest of the detected object is determined;

[0020] Obtain the relative positions of the coils between the individual coil loops in the receiving coil;

[0021] Based on the relative positions and distances of the coils, the positional relationship between each single coil loop and the region of interest of the detection object is determined.

[0022] In one embodiment, a corrected image for the corresponding imaging region is generated based on the received signal sensitivity at each location point, including:

[0023] Obtain the signal sensitivity combining coefficient of each single coil loop at each location point;

[0024] Based on the combined sensitivity coefficients of each received signal and the received signal sensitivity of each single coil loop at each location point, the received field sensitivity of each location point in the imaging region is obtained.

[0025] Based on the receiving field sensitivity of each location point in the imaging region, a corrected image of the corresponding imaging region is generated.

[0026] In one embodiment, a corrected image for the corresponding imaging region is generated based on the receiving field sensitivity of each location point in the imaging region, including:

[0027] The receiver field sensitivity at each location point in the imaging region is smoothed to obtain the smoothed receiver field sensitivity.

[0028] Based on the sensitivity of each receiving field after smoothing, a corresponding corrected image is generated.

[0029] In one embodiment, the initial magnetic resonance image of the region of interest is corrected by image correction to obtain the target magnetic resonance image, including:

[0030] Align the coordinate positions of the corrected image and the initial magnetic resonance image of the region of interest to obtain the coordinate-aligned corrected image and the initial magnetic resonance image;

[0031] The initial magnetic resonance image, after coordinate alignment, is corrected based on the corrected image after coordinate alignment to generate the target magnetic resonance image.

[0032] A magnetic resonance image correction device, the device comprising:

[0033] The receiving field sensitivity acquisition module is used to acquire the receiving field sensitivity of the receiving coil in the magnetic resonance device;

[0034] The positional relationship determination module is used to determine the positional relationship between the receiving coil and the region of interest of the object being detected in the magnetic resonance imaging device.

[0035] The received signal sensitivity determination module is used to determine the received signal sensitivity of each point in the imaging region where the region of interest is located, based on the positional relationship and the received field sensitivity.

[0036] The corrected image generation module is used to generate a corrected image of the corresponding imaging area based on the received signal sensitivity of each location point;

[0037] The correction module is used to correct the initial magnetic resonance image of the region of interest by using the correction image to obtain the target magnetic resonance image.

[0038] In one embodiment, the positional relationship determination module is used to determine the positional relationship between each single coil loop in the receiving coil and the region of interest of the detected object;

[0039] The received signal sensitivity determination module is used to determine the received signal sensitivity of each single coil circuit at each position point in the imaging area based on the positional relationship and the received field sensitivity.

[0040] A computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the steps of the method described in any of the above embodiments.

[0041] A computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method described in any of the above embodiments.

[0042] The aforementioned magnetic resonance image correction method, apparatus, computer device, and readable storage medium acquire the receiving field sensitivity of the receiving coil in the magnetic resonance device and determine the positional relationship between the receiving coil and the region of interest (ROI) of the detected object. Then, based on the positional relationship and the receiving field sensitivity, the receiving signal sensitivity of each point in the imaging region where the ROI is located is determined. Further, based on the receiving signal sensitivity of each point, a corrected image of the corresponding imaging region is generated. This corrected image is then used to correct the initial magnetic resonance image of the ROI, resulting in the target magnetic resonance image. Therefore, the receiving signal sensitivity of the receiving coil in the imaging region can be determined based on the receiving field sensitivity of the receiving coil and the positional relationship between the receiving coil and the ROI. Based on the determined receiving signal sensitivity, a corresponding corrected image is generated and used for brightness correction of the initial magnetic resonance image, solving the problem of uneven brightness distribution in high-field magnetic resonance images in traditional methods. Attached Figure Description

[0043] Figure 1 This is an application scenario diagram of a magnetic resonance image correction method in one embodiment;

[0044] Figure 2 This is a flowchart illustrating a magnetic resonance image correction method in one embodiment;

[0045] Figure 3 This is a schematic diagram illustrating the positional relationship between the receiving coil and the imaging area in one embodiment;

[0046] Figure 4 This is a schematic diagram illustrating the positional relationship between the receiving coil and the imaging area in another embodiment;

[0047] Figure 5 This is a schematic diagram illustrating the positional relationship between the receiving coil and the imaging area in yet another embodiment;

[0048] Figure 6 This is a structural block diagram of a magnetic resonance image correction device in one embodiment;

[0049] Figure 7 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0050] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0051] The magnetic resonance image correction method provided in this application can be applied to, for example... Figure 1In the application environment shown, terminal 102 communicates with server 104 via a network. The user can input commands through terminal 102, which are transmitted to server 104, enabling server 104 to perform subsequent processing operations. Server 104 can obtain the receiving field sensitivity of the receiving coil in the magnetic resonance imaging (MRI) device based on the received commands, and determine the positional relationship between the receiving coil and the region of interest (ROI) of the detected object. Then, based on the positional relationship and the receiving field sensitivity, server 104 can determine the receiving signal sensitivity of each point in the imaging region where the RIO is located, and generate a corrected image for the corresponding imaging region based on the receiving signal sensitivity of each point. Furthermore, server 104 can use the corrected image to correct the initial MRI image corresponding to the imaging region to obtain the target MRI image. The correction of the initial MRI image mainly involves correcting the brightness of the initial MRI image. Terminal 102 can be, but is not limited to, various personal computers, laptops, smartphones, tablets, and portable wearable devices. Server 104 can be implemented using a standalone server or a server cluster consisting of multiple servers.

[0052] In one embodiment, such as Figure 2 As shown, a magnetic resonance image correction method is provided, which can be applied to... Figure 1 Taking the server in the example, the following steps are included:

[0053] Step S202: Obtain the receiving field sensitivity of the receiving coil in the magnetic resonance device.

[0054] In this context, the receiving coil refers to the array of receiving coils used in a magnetic resonance imaging (MRI) device to receive magnetic resonance signals. In an MRI device, the receiving coils can include multiple single-coil loops (receiving loops formed by individual receiving coil units) arranged in rows or columns and mechanically connected, such as... Figure 3 As shown in 301 and 302, both receiving coil 301 and receiving coil 302 include multiple arrayed single coil loops, each single coil loop corresponds to a receiving coil unit, and multiple receiving coil units form a receiving array.

[0055] In this embodiment, the server can obtain the receiving field sensitivity of the receiving coil based on the equipment parameters of the magnetic resonance device and perform subsequent processing.

[0056] In this embodiment, the receiving field sensitivity of each receiving coil located in the same row or column can be the same.

[0057] In this embodiment, the receiving field sensitivity of the receiving coil in the magnetic resonance device, obtained by the server, can be described by function expansion parameters, for example, by a harmonic function or a higher-order polynomial function. Taking harmonic function expansion as an example, the receiving field sensitivity of the receiving coil... Expanded using harmonic functions, as shown in formula (1):

[0058]

[0059] in, This represents the magnetic resonance signal received by the receiving coil; γ, θ, and Let x, y, and z represent polar coordinates (spherical coordinates). Specifically, for a point (x, y, z) in the Cartesian coordinate system XYZ, its corresponding polar coordinates are: γ is the distance from a point (x, y, z) in space to the origin, 0 < γ < +∞; It is the angle formed by the half-plane passing through the Z-axis and a point in space, and the coordinate plane ZOX. θ is the angle between the origin and a point in space and the positive z-axis, 0 < θ < π; and The coefficients are spherical harmonic expansion coefficients; N takes a sufficiently large value (positive integer) to ensure... The accuracy of the harmonic function expansion is considered while also taking into account the computational complexity; n and m are positive integers, and m≤n.

[0060] Step S204: Determine the positional relationship between the receiving coil in the magnetic resonance device and the region of interest of the object being detected.

[0061] The region of interest (ROI) of the object being detected is typically located at the center of the field of view (FOV) of the magnetic resonance imaging (MRI) scanner. The center of the FOV is generally set at the center of the main magnetic field generated by the MRI scanner's magnets. Since the ROI is located within the imaging region, the MRI scanner's excitation scan sequence can acquire the magnetic resonance signal corresponding to the ROI when the object is within the imaging region. Performing a Fourier transform on this magnetic resonance signal yields the initial MRI image of the ROI.

[0062] Continue to refer to Figure 3 In magnetic resonance imaging (MRI) equipment, when testing an object, the detection area of ​​the object is located at... Figure 3 In the imaging region 303, radio frequency pulses are emitted by the magnetic resonance device to excite the detection area, and the corresponding magnetic resonance signal is acquired by the receiving coil.

[0063] Furthermore, the server can generate detection images of the corresponding detection sites based on the magnetic resonance signals collected by the receiving coil, using a certain image reconstruction algorithm, including but not limited to detection images of the upper limbs, lower limbs, head, and chest cavity of the detection object.

[0064] In this embodiment, before the hospital bed enters the magnet aperture, the server uses a camera to photograph the coil, obtaining the horizontal and vertical positions of the coil. The server further integrates the travel distance information of the hospital bed (hospital implementation position information) and the mechanical dimension information of the equipment (distance from the hospital bed to the center of the magnet) to determine the relative positional relationship between the receiving coil in the magnetic resonance imaging equipment and the region of interest of the object being detected, i.e., to determine... Figure 3 The positional relationship between receiving coil 301 and receiving coil 302 and the region of interest is determined, and subsequent processing is performed.

[0065] The server can also take pictures of the coil through a camera after the hospital bed enters the magnet aperture to obtain the horizontal and vertical positions of the coil. Furthermore, based on the camera installation position information (the relative position between the camera and the magnetic resonance device is a fixed value), the server can determine the relative positional relationship between the receiving coil in the magnetic resonance device and the region of interest of the object being detected.

[0066] During the positioning process, the server can also determine the relative positional relationship between the receiving coil in the magnetic resonance equipment and the region of interest of the detected object independently, or jointly with the camera, based on the equipment parameters of the magnetic resonance equipment (the positioning information of the coil itself, such as the fixed installation position information of the head coil, spine coil, ankle coil, etc.).

[0067] Step S206: Based on the positional relationship and the receiving field sensitivity, determine the receiving signal sensitivity of each location point in the imaging region where the region of interest is located.

[0068] In this embodiment, after determining the receiving field sensitivity of the receiving coil in the magnetic resonance device and the positional relationship between the regions of interest, the server can calculate the receiving signal sensitivity of each position point in the imaging region where the region of interest is located based on the receiving field sensitivity and the positional relationship, so as to determine the receiving signal sensitivity of each position point in the imaging region, that is, to determine the receiving field sensitivity of each position point.

[0069] In this embodiment, the receiving field sensitivity of the i-th single-coil circuit is i can take any value from 1, 2, 3, ..., K, where K is the total number of single-coil loops contained in the receiving coil; after coupling, the receiving field sensitivity of the i-th coil loop is... and The relationship can be expressed as the following formula (2):

[0070]

[0071] Where K is the number of single-coil loops, and K is a positive integer. Let C be the coupling matrix between single-coil loops (the single-coil loops contained in the receiving coil); where C 11 C 22 C 33 C KK C represents the self-inductive coupling coefficient of a single-coil circuit. 12 -C 1K C 21 -C 2K ..., C K1 -C K(K-1) This represents the mutual inductance coupling coefficient of a single-coil circuit.

[0072] In this embodiment, the received signal sensitivity at any location point in the imaging area It can be expressed as the following formula (3):

[0073]

[0074] Step S208: Based on the received signal sensitivity of each location point, generate a corrected image of the corresponding imaging area.

[0075] In this context, a correction image refers to an image used to correct the brightness of an initial magnetic resonance image generated within a region of interest in the imaging area. Alternatively, the brightness of a magnetic resonance image generated from a detected object within the imaging area can be corrected using a correction image.

[0076] In this embodiment, after determining the received signal sensitivity of each location point in the imaging area, the server can generate a corresponding corrected image by performing a Fourier transform on the received field sensitivity of each location point.

[0077] Step S210: Correct the initial magnetic resonance image of the region of interest by correcting the image to obtain the target magnetic resonance image.

[0078] As mentioned earlier, the server can generate a detection image, i.e., an initial magnetic resonance image, based on the magnetic resonance signal received by the receiving coil.

[0079] Furthermore, the server can correct the data at each location point in the generated initial magnetic resonance image by correcting the data at each location point in the correction image, thereby obtaining the target magnetic resonance image after image brightness correction. For example, the server can quotient the two images, that is, divide the pixels of each point in the initial magnetic resonance image by the pixels of the corresponding point in the correction image, to obtain the target magnetic resonance image.

[0080] In the aforementioned magnetic resonance image correction method, the receiving field sensitivity of the receiving coil in the magnetic resonance device is obtained, and the positional relationship between the receiving coil and the region of interest (ROI) of the magnetic resonance device is determined. Then, based on the positional relationship and the receiving field sensitivity, the receiving signal sensitivity of each point in the imaging region where the ROI is located is determined. Furthermore, based on the receiving signal sensitivity of each point, a corrected image for the corresponding imaging region is generated. This corrected image is then used to correct the initial magnetic resonance image corresponding to the imaging region, resulting in the target magnetic resonance image. Therefore, the receiving signal sensitivity of the receiving coil in the imaging region can be determined based on the receiving field sensitivity of the receiving coil and the positional relationship between the receiving coil and the ROI. Based on this determined receiving signal sensitivity, a corresponding corrected image is generated and used for brightness correction of the initial magnetic resonance image, thus solving the problem of uneven image brightness distribution in traditional methods.

[0081] In one embodiment, determining the positional relationship between the receiving coil in the magnetic resonance device and the region of interest of the object being detected may include: determining the positional relationship between each single coil loop in the receiving coil and the region of interest of the object being detected.

[0082] As mentioned earlier, the receiving coil in a magnetic resonance imaging device may include multiple single-coil loops. When the server determines the positional relationship between the receiving coil and the region of interest of the detection object, it may determine the positional relationship between each single-coil loop and the region of interest.

[0083] In this embodiment, the server determines the received signal sensitivity of each location point in the imaging region where the region of interest is located based on the positional relationship and the receiving field sensitivity. This may include: determining the received signal sensitivity of each single coil circuit at each location point in the imaging region based on the positional relationship and the receiving field sensitivity.

[0084] In this embodiment, as Figure 4 As shown, depending on the different positional relationships between each single-coil circuit and the region of interest of the magnetic resonance device, the received signal sensitivity at each position point in the imaging region is different. The server can superimpose the received signal sensitivity of each single-coil circuit at each position point in the imaging region according to the positional relationships and the received field sensitivity to obtain the received signal sensitivity at each position point in the imaging region.

[0085] In one embodiment, determining the positional relationship between each single coil loop in the receiving coil and the region of interest of the object to be detected by the magnetic resonance device may include: determining a positioning coil from the receiving coil and determining the positional distance between the positioning coil and the region of interest of the object to be detected; obtaining the relative positions of the coils between each single coil loop in the receiving coil; and determining the positional relationship between each single coil loop and the region of interest of the object to be detected based on the relative positions of the coils and the positional distance.

[0086] The positioning coil refers to the coil used for reference positioning in the receiving coil.

[0087] In one embodiment, the server can obtain the initial position of each single coil loop (coil unit) in the receiving coil, specifically through the following steps:

[0088] (a) A pre-scanning sequence is activated to acquire the magnetic resonance signals collected by each single-coil loop in the body coil and receiving coil, respectively. The frequency encoding direction of the magnetic resonance signals is along the direction of the main magnetic field. It should be noted that in this embodiment, the magnetic resonance signals acquired by the body coil and single-coil loops are one-dimensional, and the direction of the main magnetic field is along the Z-axis direction along the subject's axial direction, where I... v (x) represents the intensity value of the magnetic resonance signal acquired by the body coil at position x, I i (x) represents the intensity value of the magnetic resonance signal acquired at position x by each single coil loop in the receiving coil, and i represents the number of the single coil loop.

[0089] (b) Using the magnetic resonance signal acquired by the bulk coil as a reference signal, a sensitivity curve for the magnetic resonance signal acquired by each single-coil circuit is fitted and obtained. The intensity value at each coordinate point of the sensitivity curve is the ratio of the intensity value of the magnetic resonance signal acquired by the single-coil circuit to the intensity value of the reference signal. In this embodiment, the intensity value of the sensitivity curve of the single-coil circuit corresponding to position x at x is I = I i (x) / I v (x). Using the magnetic resonance signal acquired by the volume coil as a reference to calculate the sensitivity characteristics of the coil unit can minimize the influence of the scanned object's tissue structure or contrast on the coil characteristics. The coil sensitivity curve reflects the ability of a single coil loop to receive signals in the direction of the main magnetic field, and the signal receiving ability is closely related to the position of the single coil loop. The closer to the center of the single coil loop, the stronger the signal receiving ability and the larger the intensity value of the sensitivity curve; while the farther away from the center of the single coil loop, the weaker the signal receiving ability and the smaller the intensity value of the sensitivity curve.

[0090] (c) Obtain the effective maximum value on the sensitivity curve. The position of the effective maximum value is the initial positioning information of the coil unit, wherein the number of single coil loops is greater than 1.

[0091] Furthermore, the server can obtain the relative zero-point position of the receiving coil. This relative zero-point position can be obtained based on the interrelationships between coil units within a local coil and the initial positioning information of several coil units. Specifically, the interrelationships between individual coil loops in the receiving coil are determined based on the subject's positioning and the inherent structure of the receiving coil. In this embodiment, the interrelationships are such that the arrangement direction and spacing of each individual coil loop are fixed. Within this fixed arrangement direction, based on the initial positioning information of each individual coil loop, the distances of each individual coil loop from the center of the magnet are compared, and the initial position of the individual coil loop with the smallest distance from the center of the magnet is taken as the relative zero-point position.

[0092] Furthermore, the server can obtain the calculated position of each coil unit based on the relative zero point position.

[0093] In this embodiment, reference Figure 5 The positioning coil can be any single coil loop in the receiving coil, and the server can determine it based on the equipment parameters of the magnetic resonance device.

[0094] Furthermore, the server can determine the relative positions of the individual coils based on the connection arrangement between the individual coil circuits in the receiving coil. Specifically, the connection arrangement can be uniform or non-uniform, and this application does not limit this.

[0095] In this embodiment, the server can determine the positional relationship between each individual receiving coil and the region of interest of the detected object based on the relative position of the coils between the positioning receiving coil and the region of interest of the detected object, as well as the positional distance between each individual coil loop in the receiving coil. For example, continuing to refer to... Figure 5 After determining the positional distance between the positioning coil 501 and the region of interest of the detected object, such as the center point of the imaging region 502, and after determining the relative positions of the coils between the individual coil loops in the receiving coil, the server can determine the positional relationships between each individual coil loop in the single coil loops 503 to 509 and the region of interest of the detected object, so as to obtain the positional relationships between each individual coil loop in the receiving coil and the region of interest of the detected object.

[0096] In one embodiment, generating a corrected image of the corresponding region of interest based on the received signal sensitivity at each location point may include: obtaining the received signal sensitivity combining coefficient of each single coil circuit at each location point; obtaining the received field sensitivity of each location point in the imaging region based on the received signal sensitivity combining coefficient and the received signal sensitivity of each single coil circuit at each location point; and generating a corrected image of the corresponding imaging region based on the received field sensitivity of each location point in the imaging region.

[0097] In this embodiment, the server can determine the received signal sensitivity combining coefficient of each single coil loop at each position point in the imaging region based on the relative positions of the coils between the single coil loops in the receiving coil, the positional relationships between each single coil loop and the region of interest of the detected object, and the position of the region of interest of the detected object in the imaging region. For example, the farther the single coil loop is from the position point in the imaging region, the smaller the received signal sensitivity combining coefficient is, and vice versa.

[0098] Furthermore, the server can obtain the receiving field sensitivity of each location point in the imaging area based on the combined coefficient of each received signal sensitivity and the received signal sensitivity of each single coil circuit at each location point, for example, by calculating it using the following formula (4).

[0099]

[0100] The above formula is illustrated using the calculation of the receiving field sensitivity at one location point as an example. P represents the received signal sensitivity of the single-coil loop numbered i at its current position in the imaging region. i This represents the receiver signal sensitivity combining coefficient of the single-coil loop numbered i at its current position in the imaging region. This represents the receiving field sensitivity at the current location point in the imaging region.

[0101] Furthermore, the server can perform a Fourier transform based on the determined receiving field sensitivity of each location point to generate a corrected image for the corresponding imaging area.

[0102] In one embodiment, generating a corrected image of the corresponding imaging region based on the receiving field sensitivity of each location point in the region of interest may include: smoothing the receiving field sensitivity of each location point in the imaging region to obtain smoothed receiving field sensitivities; and generating a corresponding corrected image based on the smoothed receiving field sensitivities.

[0103] In this embodiment, the server can filter the receiving field sensitivity of each location point in the imaging area to remove spikes or interfering locations in the receiving field sensitivity of the imaging area, so as to obtain the smoothed receiving field sensitivity.

[0104] Furthermore, the server can generate corresponding corrected images by performing Fourier transform on the sensitivity of each receiving field obtained after filtering.

[0105] It will be understood by those skilled in the art that this is merely an illustrative example, and other methods may be used in other embodiments, which are not limited in this application.

[0106] In the above embodiments, by smoothing the receiving field sensitivity of each location point in the imaging area, smoothed receiving field sensitivity is obtained. Then, based on the smoothed receiving field sensitivity, a corresponding correction image is generated. Thus, interference points in the receiving field sensitivity of the imaging area can be eliminated, the accuracy of the correction image generation can be improved, and the accuracy of subsequent brightness correction can be improved.

[0107] In one embodiment, the initial magnetic resonance image of the region of interest is corrected by correcting the image to obtain the target magnetic resonance image. This may include: aligning the coordinate positions of the corrected image and the initial magnetic resonance image corresponding to the region of interest within the imaging area to obtain the coordinate-aligned corrected image and the initial magnetic resonance image; and correcting the image brightness of the coordinate-aligned initial magnetic resonance image based on the coordinate-aligned corrected image to generate the target magnetic resonance image.

[0108] In this embodiment, after obtaining the corrected image and the initial magnetic resonance image, the server can align the coordinate positions of the corrected image and the initial magnetic resonance image according to the positional relationship between each single coil loop and the region of interest of the detection object, so as to obtain the coordinate-aligned corrected image and the initial magnetic resonance image.

[0109] Furthermore, the server can correct the image brightness of the initial magnetic resonance image based on the aligned and corrected image to generate the target magnetic resonance image.

[0110] In the above embodiments, by aligning the coordinate positions of the corrected image and the initial magnetic resonance image, and then correcting the image brightness after alignment, the accuracy of image brightness correction can be further improved, thereby improving the accuracy of the generated target magnetic resonance image.

[0111] In one embodiment, obtaining the receiving field sensitivity of a single coil circuit in a magnetic resonance imaging (MRI) device may include: obtaining device parameters of the MRI device, including coil shape parameters and the number of coil turns of the single coil circuit; and determining the receiving field sensitivity of the single coil circuit in the MRI device based on the coil shape parameters and the number of coil turns.

[0112] Among them, equipment parameters refer to the parameters corresponding to the magnetic resonance equipment. For example, they may include the equipment's signal, type, coil type, coil shape parameters, number of single coil loops of the receiving coil, and number of coil turns of each single coil loop.

[0113] In this embodiment, the server can determine the receiving field sensitivity of a single coil circuit in the corresponding magnetic resonance imaging (MRI) device based on the coil shape parameters and the number of coil turns obtained from the device parameters. For example, the receiving field sensitivity of a single coil circuit can be simulated and calculated using Maxwell's equations to obtain the receiving field sensitivity.

[0114] It will be understood by those skilled in the art that this is merely an illustrative example, and that other methods may be used to determine the same properties in other embodiments. This application does not impose any limitations.

[0115] It should be understood that, although Figure 2 The steps in the flowchart are shown sequentially as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order in which these steps are executed, and they can be performed in other orders. Figure 2 At least some of the steps in the process may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least some of the sub-steps or stages of other steps.

[0116] In one embodiment, such as Figure 6 As shown, a magnetic resonance image correction device is provided, comprising: a receiving field sensitivity acquisition module 100, a positional relationship determination module 200, a receiving signal sensitivity determination module 300, a correction image generation module 400, and a correction module 500, wherein:

[0117] The receiving field sensitivity acquisition module 100 is used to acquire the receiving field sensitivity of the receiving coil in the magnetic resonance device.

[0118] The positional relationship determination module 200 is used to determine the positional relationship between the receiving coil and the region of interest of the detection object in the magnetic resonance device.

[0119] The received signal sensitivity determination module 300 is used to determine the received signal sensitivity of each position point in the imaging region where the region of interest is located, based on the positional relationship and the received field sensitivity.

[0120] The image correction generation module 400 is used to generate a correction image of the corresponding imaging area based on the received signal sensitivity of each location point.

[0121] The correction module 500 is used to correct the initial magnetic resonance image of the region of interest by using the correction image to obtain the target magnetic resonance image.

[0122] In one embodiment, the position relationship determination module 200 is used to determine the position relationship between each single coil loop in the receiving coil and the region of interest of the detection object.

[0123] In this embodiment, the receiving signal sensitivity determination module 300 determines the receiving signal sensitivity of each single coil circuit at each position point in the imaging area based on the positional relationship and the receiving field sensitivity.

[0124] In one embodiment, the position relationship determination module 200 may include:

[0125] The position distance determination submodule is used to determine the positioning coil from the receiving coil and to determine the position distance between the positioning coil and the region of interest of the detected object.

[0126] The coil relative position acquisition submodule is used to acquire the coil relative position between each single coil loop in the receiving coil.

[0127] The positional relationship determination submodule is used to determine the positional relationship between each single coil loop and the region of interest of the detection object based on the relative position and positional distance of the coils.

[0128] In one embodiment, the corrected image generation module 400 may include:

[0129] The received signal sensitivity combining coefficient acquisition submodule is used to obtain the received signal sensitivity combining coefficient of each single coil circuit at each location point.

[0130] The receiving field sensitivity determination submodule is used to obtain the receiving field sensitivity of each location point in the imaging area based on the combined sensitivity coefficient of each received signal and the received signal sensitivity of each single coil circuit at each location point.

[0131] The image correction generation submodule is used to generate a corrected image of the corresponding imaging area based on the receiving field sensitivity of each location point in the imaging area.

[0132] In one embodiment, the corrected image generation submodule may include:

[0133] The smoothing unit is used to smooth the receiving field sensitivity of each location point in the imaging area to obtain the smoothed receiving field sensitivity.

[0134] The correction image generation unit is used to generate a corresponding correction image based on the sensitivity of each receiving field after smoothing.

[0135] In one embodiment, the correction module 500 may include:

[0136] The alignment submodule is used to align the coordinate positions of the corrected image and the initial magnetic resonance image of the region of interest, resulting in the coordinate-aligned corrected image and the initial magnetic resonance image.

[0137] The correction submodule is used to correct the brightness of the initial magnetic resonance image after coordinate alignment based on the corrected image after coordinate alignment, and generate the target magnetic resonance image.

[0138] In one embodiment, the position relationship determination module 200 is used to determine the position relationship between each single coil loop in the receiving coil and the region of interest of the detection object.

[0139] In this embodiment, the receiving signal sensitivity determination module 300 is used to determine the receiving signal sensitivity of each single coil circuit at each position point in the imaging area based on the positional relationship and the receiving field sensitivity.

[0140] Specific limitations regarding the magnetic resonance image correction device can be found in the limitations of the magnetic resonance image correction method described above, and will not be repeated here. Each module in the aforementioned magnetic resonance image correction device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in hardware or independently of the processor in a computer device, or stored in software in the memory of a computer device, so that the processor can call and execute the corresponding operations of each module.

[0141] In one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as follows: Figure 7 As shown, the computer device includes a processor, memory, network interface, and database connected via a system bus. The processor provides computational and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores the operating system, computer programs, and the database. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage medium. The database stores data such as receiving field sensitivity, positional relationship, received signal sensitivity, corrected image, initial magnetic resonance image, and target magnetic resonance image. The network interface is used for communication with external terminals via a network connection. When executed by the processor, the computer program implements a magnetic resonance image correction method.

[0142] Those skilled in the art will understand that Figure 7 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0143] In one embodiment, a computer device is provided, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to perform the following steps: acquiring the receiving field sensitivity of a receiving coil in a magnetic resonance imaging device; determining the positional relationship between the receiving coil in the magnetic resonance imaging device and the region of interest of the detected object; determining the receiving signal sensitivity of each position point in the imaging region where the region of interest is located based on the positional relationship and the receiving field sensitivity; generating a corrected image of the corresponding imaging region based on the receiving signal sensitivity of each position point; and correcting the initial magnetic resonance image of the region of interest using the corrected image to obtain a target magnetic resonance image.

[0144] In one embodiment, when the processor executes a computer program to determine the positional relationship between the receiving coil in the magnetic resonance device and the region of interest of the object being detected, it may include: determining the positional relationship between each single coil loop in the receiving coil and the region of interest of the object being detected.

[0145] In this embodiment, when the processor executes the computer program, it determines the received signal sensitivity of each position point in the imaging region where the region of interest is located based on the positional relationship and the receiving field sensitivity. This may include: determining the received signal sensitivity of each single coil circuit at each position point in the imaging region based on the positional relationship and the receiving field sensitivity.

[0146] In one embodiment, when the processor executes a computer program to determine the positional relationships between each single coil loop in the receiving coil and the region of interest of the detected object, it may include: determining a positioning coil from the receiving coil and determining the positional distance between the positioning coil and the region of interest of the detected object; obtaining the relative positions of the coils between each single coil loop in the receiving coil; and determining the positional relationship between each single coil loop and the region of interest of the detected object based on the relative positions of the coils and the positional distance.

[0147] In one embodiment, when the processor executes a computer program, it generates a corrected image of the corresponding imaging region based on the received signal sensitivity at each location point. This may include: obtaining the received signal sensitivity merging coefficient of each single coil circuit at each location point; obtaining the received field sensitivity of each location point in the imaging region based on the received signal sensitivity merging coefficient and the received signal sensitivity of each single coil circuit at each location point; and generating a corrected image of the corresponding imaging region based on the received field sensitivity of each location point in the imaging region.

[0148] In one embodiment, when the processor executes a computer program, it generates a corrected image of the corresponding imaging region based on the receiving field sensitivity of each location point in the imaging region. This may include: smoothing the receiving field sensitivity of each location point in the imaging region to obtain smoothed receiving field sensitivities; and generating a corresponding corrected image based on the smoothed receiving field sensitivities.

[0149] In one embodiment, when the processor executes a computer program, it corrects an initial magnetic resonance image corresponding to the imaging region by correcting the image to obtain a target magnetic resonance image. This may include: aligning the coordinate positions of the corrected image and the initial magnetic resonance image of the region of interest to obtain a coordinate-aligned corrected image and an initial magnetic resonance image; and correcting the image brightness of the coordinate-aligned initial magnetic resonance image based on the coordinate-aligned corrected image to generate the target magnetic resonance image.

[0150] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon. When executed by a processor, the computer program performs the following steps: acquiring the receiving field sensitivity of a receiving coil in a magnetic resonance imaging device; determining the positional relationship between the receiving coil in the magnetic resonance imaging device and the region of interest of the object being detected; determining the receiving signal sensitivity of each point in the imaging region where the region of interest is located based on the positional relationship and the receiving field sensitivity; generating a corrected image based on the receiving signal sensitivity of each point; and correcting the initial magnetic resonance image of the region of interest using the corrected image to obtain a target magnetic resonance image.

[0151] In one embodiment, when the computer program is executed by the processor, it determines the positional relationship between the receiving coil in the magnetic resonance device and the region of interest of the object being detected, which may include: determining the positional relationship between each single coil loop in the receiving coil and the region of interest of the object being detected.

[0152] In this embodiment, when the computer program is executed by the processor, it determines the received signal sensitivity of each position point in the imaging region where the region of interest is located based on the positional relationship and the receiving field sensitivity. This may include: determining the received signal sensitivity of each single coil circuit at each position point in the imaging region based on the positional relationship and the receiving field sensitivity.

[0153] In one embodiment, when the computer program is executed by the processor, determining the positional relationships between each single coil loop in the receiving coil and the region of interest of the detection object may include: determining a positioning coil from the receiving coil and determining the positional distance between the positioning coil and the region of interest of the detection object; obtaining the relative positions of the coils between each single coil loop in the receiving coil; and determining the positional relationship between each single coil loop and the region of interest of the detection object based on the relative positions of the coils and the positional distance.

[0154] In one embodiment, when the computer program is executed by the processor, it generates a corrected image of the corresponding imaging region based on the received signal sensitivity at each location point. This may include: obtaining the received signal sensitivity combining coefficient of each single coil circuit at each location point; obtaining the received field sensitivity of each location point in the imaging region based on the received signal sensitivity combining coefficient and the received signal sensitivity of each single coil circuit at each location point; and generating a corrected image of the corresponding imaging region based on the received field sensitivity of each location point in the imaging region.

[0155] In one embodiment, when the computer program is executed by the processor, it generates a corrected image of the corresponding imaging region based on the receiving field sensitivity of each location point in the imaging region. This may include: smoothing the receiving field sensitivity of each location point in the imaging region to obtain smoothed receiving field sensitivities; and generating a corresponding corrected image based on the smoothed receiving field sensitivities.

[0156] In one embodiment, when the computer program is executed by the processor, it corrects an initial magnetic resonance image of the region of interest by correcting the correction image to obtain a target magnetic resonance image. This may include: aligning the coordinate positions of the correction image and the initial magnetic resonance image corresponding to the imaging region to obtain a coordinate-aligned correction image and an initial magnetic resonance image; and correcting the image brightness of the coordinate-aligned initial magnetic resonance image based on the coordinate-aligned correction image to generate the target magnetic resonance image.

[0157] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0158] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0159] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the invention patent. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this patent application should be determined by the appended claims.

Claims

1. A magnetic resonance image correction method, characterized in that, The method includes: Based on the equipment parameters of the magnetic resonance device, the receiving field sensitivity of the receiving coil in the magnetic resonance device is obtained; The positioning coil is determined from the receiving coil, and the positional distance between the positioning coil and the region of interest of the detected object is determined; Obtain the initial position of each single coil loop in the receiving coil, and obtain the relative zero point position based on the interrelationship between each coil loop and the initial position of the coil loop; The calculated position of each coil loop and the connection arrangement between each single coil loop are obtained based on the relative zero point position, so as to obtain the relative position of each single coil loop in the receiving coil. Based on the relative positions of the coils and the positional distances, the positional relationship between each single coil loop and the region of interest of the detection object is determined; Based on the positional relationships and the receiving field sensitivity, the receiving signal sensitivity of each single coil loop at each position point in the imaging region where the region of interest is located is determined. Based on the received signal sensitivity of each single coil circuit at each of the said locations and the combined coefficient of the received signal sensitivity of each single coil circuit at each of the said locations, the received field sensitivity of each location in the imaging region is obtained, so as to generate a corrected image corresponding to the imaging region. The initial magnetic resonance image of the region of interest is corrected using the corrected image to obtain the target magnetic resonance image.

2. The method according to claim 1, characterized in that, The step of obtaining the relative zero point position based on the interrelationship between each coil circuit and the initial position of the coil circuit includes: The interrelationships between the coil circuits in the receiving coil are determined based on the subject's positioning and the inherent structure of the receiving coil; the arrangement direction and spacing of the coil circuits are fixed in the interrelationships. In a fixed arrangement direction, based on the initial positioning information of each single coil loop, the distance of the single coil loop from the center of the magnet is compared, and the relative zero point position is determined based on the initial position of the single coil loop with the smallest distance from the center of the magnet.

3. The method according to claim 1, characterized in that, Obtaining the received signal sensitivity combining coefficients of each of the single-coil loops at each of the specified locations includes: Based on the relative positions of the coils between each single coil loop, the positional relationships between each single coil loop and the region of interest of the detected object, and the position of the region of interest of the detected object in the imaging region, the receiving signal sensitivity combining coefficient of each single coil loop at each position point is determined. The farther each single-coil circuit is from the location point in the imaging area, the smaller the signal sensitivity combining coefficient of each single-coil circuit at each location point.

4. The method according to claim 1, characterized in that, Based on the receiving field sensitivity of each location point in the imaging region, a corrected image corresponding to the imaging region is generated, including: The receiving field sensitivity of each location point in the imaging region is smoothed to obtain the smoothed receiving field sensitivity. Based on the sensitivity of each receiving field after smoothing, a corresponding corrected image is generated.

5. The method according to claim 1, characterized in that, The step of correcting the initial magnetic resonance image of the region of interest using the corrected image to obtain the target magnetic resonance image includes: Align the coordinate positions of the corrected image and the initial magnetic resonance image of the region of interest to obtain the coordinate-aligned corrected image and the initial magnetic resonance image; The initial magnetic resonance image after coordinate alignment is corrected for brightness based on the corrected image after coordinate alignment to generate the target magnetic resonance image.

6. A magnetic resonance image correction device, characterized in that, The device includes: The receiving field sensitivity acquisition module is used to acquire the receiving field sensitivity of the receiving coil in the magnetic resonance device based on the device parameters of the magnetic resonance device; A positional relationship determination module is used to determine a positioning coil from the receiving coil and to determine the positional distance between the positioning coil and the region of interest of the detection object; to obtain the initial position of each single coil loop in the receiving coil, and to obtain a relative zero-point position based on the mutual relationship between each coil loop and the initial position of the coil loop; to obtain the calculated position of each coil loop and the connection arrangement between each single coil loop based on the relative zero-point position, so as to obtain the relative coil position between each single coil loop in the receiving coil; and to determine the positional relationship between each single coil loop and the region of interest of the detection object based on the relative coil position and the positional distance. The received signal sensitivity determination module is used to determine the received signal sensitivity of each of the single coil loops at each position point in the imaging region where the region of interest is located, based on the positional relationship and the received field sensitivity. The corrected image generation module is used to obtain the receiving field sensitivity of each location point in the imaging region based on the receiving signal sensitivity of each single coil circuit at each location point and the receiving signal sensitivity combination coefficient of each single coil circuit at each location point, so as to generate a corrected image corresponding to the imaging region. The correction module is used to correct the initial magnetic resonance image of the region of interest using the correction image to obtain the target magnetic resonance image.

7. The apparatus according to claim 6, characterized in that, The device further includes a correction image generation submodule, which includes a smoothing processing unit and a correction image generation unit; The smoothing unit is used to smooth the receiving field sensitivity of each location point in the imaging area to obtain the smoothed receiving field sensitivity. The corrected image generation unit is used to generate a corresponding corrected image based on the sensitivity of each receiving field after smoothing.

8. The apparatus according to claim 6, characterized in that, The correction module includes an alignment submodule and a correction submodule; The alignment submodule is used to align the coordinate positions of the corrected image and the initial magnetic resonance image of the region of interest, so as to obtain the coordinate-aligned corrected image and the initial magnetic resonance image. The correction submodule is used to correct the image brightness of the initial magnetic resonance image after coordinate alignment based on the corrected image after coordinate alignment, and generate the target magnetic resonance image.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 5.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 5.