Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Embodiment in the present invention, those of ordinary skill in the art are not having
Every other embodiment obtained under the premise of creative work is made, shall fall within the protection scope of the present invention.
In order to overcome the above problem of the prior art, the embodiment of the present invention provides a kind of assessment of ablative surgery therapeutic effect
Method, inventive concept are that the contraction situation actually occurred in art by analyzing tumour predicts the center of compression point of tumour,
According to the center of compression of tumour point, the tumor region in pre-operative image is mapped on postoperative image, is considered in mapping swollen
The actual change of tumor heat shrinkable, so as to more accurately reflect that ablation areas to the package situation of tumor region, obtains more quasi-
The assessment result of true ablative surgery therapeutic effect.
Fig. 1 is the flow diagram according to the appraisal procedure of ablative surgery therapeutic effect provided in an embodiment of the present invention.Such as
Shown in Fig. 1, a kind of appraisal procedure of ablative surgery therapeutic effect includes: step S101, respectively to the tumor area in pre-operative image
Ablation areas in domain and postoperative image is pre-processed so that in pretreated pre-operative image tumor region intensity profile
It is identical as the intensity profile of tumor region peripheral organs tissue regions, the gray scale of ablation areas point in pretreated postoperative image
Cloth is identical as the intensity profile of ablation areas peripheral organs tissue regions.
It should be noted that the appraisal procedure of ablative surgery therapeutic effect provided in an embodiment of the present invention is suitable for liver
The aftertreatment effect that the tumour of each diseased organ such as dirty, stomach, lung, pancreas, thyroid gland, mammary gland, enteron aisle implements ablative surgery is commented
Estimate, be not used to the implementation of ablative surgery, without after implementing ablative surgery to the diagnosing and treating of patient.
It is understood that carry out ablative surgery therapeutic effect assessment before, first obtain patient pre-operative image and
Postoperative image.Pre-operative image includes tumor region;Postoperative image includes ablation areas.Tumor region exists only in pre-operative image
In, be not present in postoperative image, and ablation areas is existed only in postoperative image, is not present in pre-operative image.
It should be noted that tumor region, refer to tumour shared region in the picture.
It should be noted that after ablative surgery, from inside to outside, three layers are sequentially formed: carbonized region, coagulation region and not disappearing
Tabetisol domain.Carbonized region and coagulation region are all shunk, thus ablation areas includes carbonized region and coagulation region.
It should be noted that due to the intensity profile of tumor region and ablation areas approximation, and with normal organ tissue area
The intensity profile in domain is significantly different, and will lead to tumor region and the edge of ablation areas in subsequent image registration process can tend to pair
Together, error hiding is generated, to have an impact to the elastic deformation field around tumour, so that the actual shrinkage situation finally predicted is not
Accurately.
Therefore, in order to avoid error hiding occurs in registration for tumor region and ablation areas, it can use image repair skill
Art pre-processes pre-operative image and postoperative image, specifically respectively to the tumor region and postoperative image in pre-operative image
In ablation areas pre-processed so that in pretreated pre-operative image the intensity profile of tumor region and tumor region week
The intensity profile for enclosing organ-tissue region is identical, the intensity profile and ablation areas of ablation areas in pretreated postoperative image
The intensity profile of peripheral organs tissue regions is identical, to will not lead to tumor region and zone of ablation in subsequent image registration process
The edge alignment in domain, the actual shrinkage situation of prediction are inaccurate.
The ash of the intensity profile of tumor region and tumor region peripheral organs tissue regions in pretreated pre-operative image
Degree distribution is identical, refers to the intensity profile of tumor region and tumor region peripheral organs tissue area in pretreated pre-operative image
The intensity profile in domain has the probability-distribution function of same form, and only the parameter of probability-distribution function has differences.
The ash of the intensity profile of ablation areas and ablation areas peripheral organs tissue regions in pretreated postoperative image
Degree distribution is identical, refers to the intensity profile of ablation areas and ablation areas peripheral organs tissue area in pretreated postoperative image
The intensity profile in domain has the probability-distribution function of same form, and only the parameter of probability-distribution function has differences.
Step S102, to the ablation in the tumor region and pretreated postoperative image in pretreated pre-operative image
Region carries out image registration, obtains elastic deformation field, and carry out vector field analysis to elastic deformation field, determines actual shrinkage center
Point;Wherein, elastic deformation field, for describing the contraction of tumour.
It should be noted that tumour will drive peripheral organs during heat shrinkable when carrying out microblogging ablative surgery
Therefore tissue, such as blood vessel, can derive tumour reality also to contract according to the elastic deformation of tumour peripheral organs tissue
Contraction process, which can be indicated by the elastic deformation field of tumour peripheral organs tissue.Tumour peripheral organs group
The elastic deformation field knitted, it is believed that identical as the elastic deformation field of tumour.
Specifically, after being pre-processed to pre-operative image and postoperative image, according to method for registering images to pretreated
The ablation areas in tumor region and pretreated postoperative image in pre-operative image is registrated, to be able to achieve to lesion
The estimation of organ internal motion obtains the elastic deformation field for describing the contraction of tumour.
Elastic deformation field reflects the elastic shape of the tumour peripheral organs tissue comprising diseased organ internal motion information
Become.Since peripheral organs tissue is located at around tumour, thus the elastic deformation of tumour peripheral organs tissue can also reflect tumour
It shrinks.
Elastic deformation field reflects the elastic movement inside diseased organ, due to wherein including tumour peripheral organs tissue
Motion information, therefore can be used to analysis actual shrinkage situation.It is shunk since tumour will drive surrounding tissue when heated,
So that elastic deformation field can be directed toward the one-point or multi-point of inside tumor jointly, the point being directed toward jointly is the center of compression of tumour
Point.Therefore, by finding actual shrinkage central point, it can analyze the contractile motion of tumour, and then state the contraction situation of tumour.
Since elastic deformation field is vector field, vector field analysis can be carried out to elastic deformation field, to determine actual shrinkage
Central point.Vector field analysis is carried out to elastic deformation field, a kind of suitable side can be selected from existing vector field analysis method
Method.
Step S103, according to actual shrinkage central point, the tumor region in pretreated pre-operative image is mapped to figure
As on the pretreated postoperative image after registration, and according to being mapped on the pretreated postoperative image after image registration
The ablation areas in pretreated postoperative image after tumor region and image registration obtains commenting for ablative surgery therapeutic effect
Estimate result.
It should be noted that the implementation process of microwave ablation operation is that ablation needle is inserted into tumour, produced by ablation needle
Raw microwave leads to tumour cell coagulation necrosis by strenuous exercise's frictional heat generation of hydrone in microwave concussion electric field.Cause
This, it is believed that the center of compression point of tumour is the position in art where ablation needle.According to the actual situation, ablation needle can be one
It is a or multiple, because the center of compression point of the tumour obtained is one or more.
After obtaining actual shrinkage central point by step S102, the conduct of actual shrinkage central point will be obtained by step S102
Position where ablation needle, so as to according to where ablation needle position and existing any actual shrinkage model, predict tumour
Contraction situation, by pretreated pre-operative image tumor region combine prediction obtain tumour contraction situation, mapping
Onto pretreated postoperative image, so as to obtain the ablation areas in the pretreated postoperative image after image registration,
Package situation to the tumor region being mapped on the pretreated postoperative image after image registration, and according to above-mentioned package feelings
Condition assesses the therapeutic effect of ablative surgery, obtains the assessment result of ablative surgery therapeutic effect.
After the embodiment of the present invention is by pre-processing pre-operative image and postoperative image, to pretreated pre-operative image
Image registration is carried out with pretreated postoperative image, obtains elastic deformation field, vector field analysis is carried out to elastic deformation field, really
Determine actual shrinkage central point, tumor region is mapped on postoperative image according to actual shrinkage central point, obtains ablative surgery and control
The assessment result of therapeutic effect can accurately reflect the contraction of tumour in art, so as to obtain more accurate ablative surgery therapeutic effect
Assessment result.
Content based on the above embodiment, carrying out pretreated specific steps to the tumor region in pre-operative image includes:
For each pixel in the tumor region in pre-operative image, the gray scale of pixel is replaced with to the ash of each pixel in the neighborhood of pixel
The weighted sum of degree.
Specifically, by the gray scale of each pixel in the tumor region in pre-operative image, with picture each in the neighborhood of the pixel
The weighted sum of the gray scale of element is replaced, and realizes the pretreatment to the tumor region in pre-operative image.
The neighborhood of each pixel, can be 4 neighborhoods or 8 neighborhoods, but not limited to this.
Pretreated specific formula is carried out to the tumor region in pre-operative image are as follows:
Wherein, Ω1Indicate tumor region;X is tumor region Ω1In pixel;I (x) indicates the ash of pixel x after pretreatment
Degree;Neighborhood Ω0Indicate the neighborhood of pixel x;The neighborhood Ω of y expression pixel x0In pixel;The gray scale of I (y) expression pixel y;W table
Show weight;α0For preset normal number, regulation coefficient is indicated.
For α0Value, lower value can be reduced noise, and higher value can retain image detail.
The gray scale of pixel each in tumor region is replaced with the ash of each pixel in the neighborhood of the pixel by the embodiment of the present invention
The weighted sum of degree is avoided that error hiding occurs in registration for tumor region and ablation areas, so as to obtain more accurate tumour
Situation is shunk, the assessment result of more accurate ablative surgery therapeutic effect can be obtained.
Content based on the above embodiment, carrying out pretreated specific steps to the ablation areas in postoperative image includes:
For each pixel in the ablation areas in postoperative image, the gray scale of pixel is replaced with to the ash of each pixel in the neighborhood of pixel
The weighted sum of degree.
Specifically, by the gray scale of each pixel in the ablation areas in postoperative image, with picture each in the neighborhood of the pixel
The weighted sum of the gray scale of element is replaced, and realizes the pretreatment to the tumor region in postoperative image.
The neighborhood of each pixel, can be 4 neighborhoods or 8 neighborhoods, but not limited to this.
Pretreated method can be carried out using to the tumor region in pre-operative image, to the ablation areas in postoperative image
It is pre-processed.Pretreated specific formula is carried out to the ablation areas in postoperative image are as follows:
Wherein, Ω2Indicate ablation areas;X is ablation areas Ω2In pixel;I (x) indicates the ash of pixel x after pretreatment
Degree;Neighborhood Ω0Indicate the neighborhood of pixel x;The neighborhood Ω of y expression pixel x0In pixel;The gray scale of I (y) expression pixel y;W table
Show weight;α0For preset normal number, regulation coefficient is indicated.
For α0Value, lower value can be reduced noise, and higher value can retain image detail.Pre-operative image and art
α in the specific formula of image preprocessing afterwards0It is identical.
The neighborhood used when it is understood that being pre-processed to the tumor region in pre-operative image, can with to art
The neighborhood used when the ablation areas in image is pre-processed afterwards is identical or different.For example, 8 neighborhoods can be all used, it can also
4 neighborhoods and 8 neighborhoods are respectively adopted.
The gray scale of pixel each in ablation areas is replaced with the ash of each pixel in the neighborhood of the pixel by the embodiment of the present invention
The weighted sum of degree is avoided that error hiding occurs in registration for tumor region and ablation areas, so as to obtain more accurate tumour
Situation is shunk, the assessment result of more accurate ablative surgery therapeutic effect can be obtained.
Content based on the above embodiment, to tumor region in pretreated pre-operative image and pretreated postoperative
Ablation areas in image carries out image registration, and the specific steps for obtaining elastic deformation field include: to pretreated preoperative figure
The ablation areas in tumor region and pretreated postoperative image as in carries out Rigid Registration, so that pretreated preoperative
Mutual information between pretreated postoperative image after image and Rigid Registration is maximum;It is right according to differomorphism deformation model
The ablation areas in pretreated postoperative image after tumor region and Rigid Registration in pretreated pre-operative image into
Row elastic registrating obtains elastic deformation field.
Specifically, to the zone of ablation in the tumor region and pretreated postoperative image in pretreated pre-operative image
Domain carries out image registration, is divided into two steps: Rigid Registration and elastic registrating.
Due to for assessing ablative surgery therapeutic effect pre-operative image and postoperative image pick up from different when and wheres,
There are the pose differences of large scale between image, can not only handle by elastic registrating.Thus need first predict two images between
Rigid deformation corrects the pose difference of large scale by Rigid Registration.
Pose refers to position and posture.
Rigid Registration can be guided using the grayscale information of entire image, and use mutual information (Mutual
Information, abbreviation MI) similarity degree between Lai Hengliang two images.The target of Rigid Registration is exactly to find an ideal
Rigid transformation matrixSo that the pretreated postoperative image I after rigid transformationR-postWith pretreated preoperative figure
As IpreBetween mutual information reach maximum value.Pretreated postoperative image can use IpostIt indicates.
Mutual information can use normalised mutual information (Normalized Mutual Information, abbreviation NMI), but
It is without being limited thereto.
Rigid Registration is carried out based on mutual information, pre-operative image is can be suitably used for and postoperative image is adopted by different medical images
The case where set method obtains.For example, carrying out rigidity based on mutual information when pre-operative image is CT images, postoperative image is MRI image
Registration can not be influenced by image collecting method is different, obtain better registration effect.
Rigid transformation matrixFor Euler's transformation (Euler Transform) matrix.Rigid transformation, by rotation and translation
Two parts form.Rigid transformation matrixIt can be indicated by following formula:
Wherein, α, beta, gamma are respectively indicated around x, y, the rotation angle of z-axis;tx, ty, tzIt respectively indicates along x, y, the translation of z-axis
Amount.
After Rigid Registration, pre-operative image and postoperative image have reached alignment spatially on the whole.But due to needing
The non-rigid motion inside diseased organ is analyzed, needs further to obtain corresponding motion information, i.e. bullet using elastic registrating
Property Deformation Field
When carrying out elastic registrating, it can be matched according to differomorphism deformation (diffeomorphic demons) model
Standard obtains elastic deformation field.The elastic deformation field that prediction acquisition is carried out according to differomorphism deformation model, with differomorphism
Characteristic.The characteristic makes elastic deformation field everywhere can be micro-, i.e. the variation of elastic deformation field can continuously be led, thus is grasping tumour week
On the basis of enclosing organ-tissue motion information, actual shrinkage center position can be derived.
For example, the deformation of log-domain differomorphism (log-domain diffeomorphic demons) model can be used
Elastic registrating is carried out, but not limited to this.
Pretreated postoperative image IpostThrough final Deformation FieldAfter deformation, with pretreated pre-operative image IpreHighly
Alignment.
Wherein, rigid transformation matrixReflect the pose rigidity deformation of patient's preoperative and postoperative, elastic deformation fieldInstead
The organ-tissue elastic deformation comprising diseased organ internal motion information is reflected.
Elastic registrating of the embodiment of the present invention by Rigid Registration based on mutual information and based on differomorphism deformation model,
Elastic deformation field is obtained, pose difference and tumour peripheral organs of the postoperative image relative to pre-operative image can be more accurately isolated
The elastic deformation of tissue can obtain more accurately for describing the elastic deformation field of the contraction of tumour, more acurrate so as to obtain
Actual shrinkage situation, the assessment result of more accurate ablative surgery therapeutic effect can be obtained.Further, based on mutual information
Rigid Registration can be suitably used for pre-operative image and postoperative image is the case where acquisition by different medical image acquisition methods, using model
It encloses wider.
Content based on the above embodiment carries out vector field analysis to elastic deformation field, determines actual shrinkage central point
Specific steps include: to decompose to elastic deformation field, obtain the potential function of the irrotational field and irrotational field in elastic deformation field;Root
According to the potential function of irrotational field, the central point that propagation is reduced in irrotational field is determined;Obtain the central point that propagation is respectively reduced in irrotational field
Shrinkage degree, according to shrinkage degree in irrotational field diminution propagate central point cluster, shrinkage degree is maximum
One kind reduces the central point propagated and is determined as actual shrinkage central point.
Specifically, vector field analysis can be carried out to elastic deformation field by following step, determines actual shrinkage center
Point.
After obtaining elastic deformation field, elastic deformation field is decomposed, elastic deformation field is decomposed into irrotational field, without end of a performance
With reconciliation field.
?It is middle that there are many key points, such as source (to expand the central point propagated, center of an
Expanding propagation), sink (reduces the central point propagated), rotational center (central point of rotation)
Deng.The contraction of tumour is reflected as the propagation of the diminution in Deformation Field, thus sink can indicate actual shrinkage central point.
Wherein,Indicate elastic deformation field;Indicate irrotational field;Indicate no end of a performance;Indicate reconciliation field.
Irrotational field, only comprising diverging and/or contraction;Without end of a performance, only comprising rotation;Reconciliation field, approximation translation.
Since need to obtain is center of compression point, center of compression point obviously only exist in irrotational field.
During decomposing elastic Deformation Field, the potential function D of irrotational field can also be obtained.
Wherein, G∞It indicates free space (free-space Green's function);▽ indicates divergence operator;Ω indicates pretreatment
Pre-operative image I afterwardspreSpace.
According to the potential function of irrotational field, the sink (reducing the central point propagated) in irrotational field can be determined.
Due toIn complicated movement, will lead to and generate many sink.Sink includes two classes: noise sink and correct
sink.Noise sink is the shrinkage degree very little as caused by noise;Correct sink is to shrink journey caused by actual shrinkage
It spends larger.
In order to eliminate the influence of noise, obtain correct actual shrinkage central point, propagation is reduced in determining irrotational field
After central point, the shrinkage degree of each central point for reducing and propagating is obtained.
It can reflect the shrinkage degree of each sink point by Jacobian.The expression of Jacobian
Are as follows:
Wherein,It respectively indicatesIn x, y, the component in three directions of z-axis.
JAC illustrates that volume remains unchanged for 1;The value of JAC is smaller, illustrates that volume contraction degree is bigger.
The expression that the coordinate of each sink point is substituted into Jacobian, can be obtained the receipts of each sink point
Contracting degree.
After the shrinkage degree for obtaining each sink point, carried out using shrinkage degree (JAC value) of the clustering algorithm to sink point
The maximum a kind of central point propagated that reduces of shrinkage degree is determined as correct sink, i.e. actual shrinkage central point by cluster.
And the lesser a kind of or multiclass of shrinkage degree reduces the central point propagated, then is noise sink.
The embodiment of the present invention is not specifically limited the clustering algorithm of use.
For example, can be clustered using K-means algorithm, the central point propagated will be reduced and be polymerized to two classes: shrinking journey
Biggish a kind of and lesser one kind of shrinkage degree is spent, the biggish one kind of shrinkage degree is correct sink, and shrinkage degree is lesser
One kind is noise sink.
The embodiment of the present invention reduces the central point of propagation by decompose to elastic deformation field in determining irrotational field, according to
The central point that shrinkage degree propagates the diminution in irrotational field clusters, will be in the maximum a kind of diminution propagation of shrinkage degree
Heart point is determined as actual shrinkage central point, can eliminate the interference of noise, more accurate actual shrinkage central point is obtained, so as to obtain
Obtain the assessment result of more accurate ablative surgery therapeutic effect.
Content based on the above embodiment determines the central point that propagation is reduced in irrotational field according to the potential function of irrotational field
Specific steps include: the maximum point for obtaining the potential function of irrotational field, the point in the corresponding irrotational field of maximum point determines
For the central point for reducing propagation in irrotational field.
Specifically, sink is the point in the corresponding irrotational field of maximum point of the potential function of irrotational field.
According to the potential function of irrotational field, the maximum point of the potential function of irrotational field can be determined, so as to according to irrotational field
The maximum point of potential function determine the central point that propagation is reduced in irrotational field.
Point in the corresponding irrotational field of the maximum point of the potential function of irrotational field is determined as irrotational field by the embodiment of the present invention
It is middle to reduce the central point propagated, the central point for accurately reducing and propagating can be obtained, so as to obtain in more accurate actual shrinkage
Heart point, the assessment result for obtaining more accurate ablative surgery therapeutic effect.
Content based on the above embodiment decomposes elastic deformation field, obtain elastic deformation field in irrotational field and
The specific steps of the potential function of irrotational field include: to be decomposed according to Huo Qi Helmholtz decomposition algorithm to elastic deformation field, are obtained
Take the potential function of the irrotational field and irrotational field in elastic deformation field.
Specifically, can be decomposed using Huo Qi Helmholtz (Hodge Helmholtz decomposition, referred to as
HHD) algorithm decomposes elastic deformation field.
For example, can be using Natural HHD (natural Huo Qi Helmholtz decomposes) algorithm or Efficient HHD (effect
Rate Huo Qi Helmholtz decomposes), elastic deformation field is decomposed.
The embodiment of the present invention decomposes elastic deformation field according to Huo Qi Helmholtz decomposition algorithm, can obtain more acurrate
Irrotational field and irrotational field potential function, so as to which the assessment result of more accurate ablative surgery therapeutic effect can be obtained.
Fig. 2 is the functional block diagram according to the assessment device of ablative surgery therapeutic effect provided in an embodiment of the present invention.It is based on
The content of above-described embodiment, as shown in Fig. 2, a kind of assessment device of ablative surgery therapeutic effect includes preprocessing module 201, figure
As registration module 202 and recruitment evaluation module 203, in which:
Preprocessing module 201, for respectively to the ablation areas in the tumor region and postoperative image in pre-operative image into
Row pretreatment, so that the intensity profile of tumor region and tumor region peripheral organs tissue regions in pretreated pre-operative image
Intensity profile it is identical, the intensity profile of ablation areas and ablation areas peripheral organs tissue area in pretreated postoperative image
The intensity profile in domain is identical;
Image registration module 202, for tumor region in pretreated pre-operative image and pretreated postoperative
Ablation areas in image carries out image registration, obtains elastic deformation field, and carry out vector field analysis to elastic deformation field, determines
Actual shrinkage central point;
Recruitment evaluation module 203 is used for according to actual shrinkage central point, by the tumor area in pretreated pre-operative image
On pretreated postoperative image after domain mapping to image registration, and it is swollen on pretreated postoperative image according to being mapped to
The ablation areas in pretreated postoperative image behind tumor region and image registration obtains the assessment of ablative surgery therapeutic effect
As a result;
Wherein, elastic deformation field, for describing the contraction of tumour.
Specifically, after the pre-operative image and postoperative image that obtain patient, preprocessing module 201 is respectively in pre-operative image
Ablation areas in tumor region and postoperative image is pre-processed so that in pretreated pre-operative image tumor region ash
Degree is distributed, in pretreated postoperative image ablation areas identical as the intensity profile of tumor region peripheral organs tissue regions
Intensity profile is identical as the intensity profile of ablation areas peripheral organs tissue regions, avoids tumor area in subsequent image registration process
Domain and the alignment of the edge of ablation areas, error hiding occurs in registration for tumor region and ablation areas.
Image registration module 202 is in the tumor region and pretreated postoperative image in pretreated pre-operative image
Ablation areas carry out image registration, obtain the elastic deformation field for describing the contraction of tumour;It is right after obtaining elastic deformation field
Elastic deformation field carries out vector field analysis, determines actual shrinkage central point.
Recruitment evaluation module 203 reflects the tumor region in pretreated pre-operative image according to actual shrinkage central point
On pretreated postoperative image after being mapped to image registration, the assessment result of ablative surgery therapeutic effect is obtained.
The assessment device of ablative surgery therapeutic effect provided in an embodiment of the present invention is provided for executing the embodiment of the present invention
Ablative surgery therapeutic effect appraisal procedure, each module that the assessment device of ablative surgery therapeutic effect includes realizes corresponding function
The specific method and process of energy are detailed in the embodiment of the appraisal procedure of above-mentioned ablative surgery therapeutic effect, and details are not described herein again.
Assessment of the assessment device of the ablative surgery therapeutic effect for the ablative surgery therapeutic effect of foregoing embodiments
Method.Therefore, the description and definition in the appraisal procedure of the ablative surgery therapeutic effect in foregoing embodiments, can be used for
The understanding of each execution module in the embodiment of the present invention.
After the embodiment of the present invention is by pre-processing pre-operative image and postoperative image, to pretreated pre-operative image
Image registration is carried out with pretreated postoperative image, obtains elastic deformation field, vector field analysis is carried out to elastic deformation field, really
Determine actual shrinkage central point, tumor region is mapped on postoperative image according to actual shrinkage central point, obtains ablative surgery and control
The assessment result of therapeutic effect can accurately reflect the contraction of tumour in art, so as to obtain more accurate ablative surgery therapeutic effect
Assessment result.
Fig. 3 is the structural block diagram according to electronic equipment provided in an embodiment of the present invention.Content based on the above embodiment, such as
Shown in Fig. 3, which may include: processor (processor) 301, memory (memory) 302 and bus 303;Its
In, processor 301 and memory 302 pass through bus 303 and complete mutual communication;Processor 301 is stored in for calling
In reservoir 302 and the computer program instructions that can be run on processor 301, to execute provided by above-mentioned each method embodiment
Method, for example, the ablation areas in the tumor region and postoperative image in pre-operative image is pre-processed respectively, so that
The intensity profile phase of the intensity profile of tumor region and tumor region peripheral organs tissue regions in pretreated pre-operative image
Together, in pretreated postoperative image the intensity profile of ablation areas and ablation areas peripheral organs tissue regions intensity profile
It is identical;Image is carried out to the ablation areas in the tumor region and pretreated postoperative image in pretreated pre-operative image
Registration obtains elastic deformation field, and carries out vector field analysis to elastic deformation field, determines actual shrinkage central point;According to tumour
Center of compression point, the pretreated postoperative figure after the tumor region in pretreated pre-operative image to be mapped to image registration
As upper, and according to the tumor region and the pre- place after image registration being mapped on the pretreated postoperative image after image registration
The ablation areas in postoperative image after reason obtains the assessment result of ablative surgery therapeutic effect.
Another embodiment of the present invention discloses a kind of computer program product, and computer program product is non-transient including being stored in
Computer program on computer readable storage medium, computer program include program instruction, when program instruction is held by computer
When row, computer is able to carry out method provided by above-mentioned each method embodiment, for example, respectively to swollen in pre-operative image
Ablation areas in tumor region and postoperative image is pre-processed so that in pretreated pre-operative image tumor region gray scale
It is distributed, in the pretreated postoperative image ash of ablation areas identical as the intensity profile of tumor region peripheral organs tissue regions
Degree distribution is identical as the intensity profile of ablation areas peripheral organs tissue regions;To the tumor area in pretreated pre-operative image
Ablation areas in domain and pretreated postoperative image carries out image registration, obtains elastic deformation field, and to elastic deformation field
Vector field analysis is carried out, determines actual shrinkage central point;It, will be in pretreated pre-operative image according to actual shrinkage central point
Tumor region is mapped on the pretreated postoperative image after image registration, and according to being mapped to the pretreatment after image registration
The ablation areas in the tumor region on postoperative image and the pretreated postoperative image after image registration afterwards obtains ablation
The assessment result of surgical result.
In addition, the logical order in above-mentioned memory 302 can be realized by way of SFU software functional unit and conduct
Independent product when selling or using, can store in a computer readable storage medium.Based on this understanding, originally
The technical solution of the inventive embodiments substantially part of the part that contributes to existing technology or the technical solution in other words
It can be embodied in the form of software products, which is stored in a storage medium, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes the present invention respectively
The all or part of the steps of a embodiment the method.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory
(ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic or disk
Etc. the various media that can store program code.
Another embodiment of the present invention provides a kind of non-transient computer readable storage medium, non-transient computer readable storages
Medium storing computer instruction, computer instruction makes computer execute method provided by above-mentioned each method embodiment, such as wraps
It includes: the ablation areas in the tumor region and postoperative image in pre-operative image being pre-processed respectively, so that pretreated
The intensity profile of tumor region is identical as the intensity profile of tumor region peripheral organs tissue regions in pre-operative image, after pretreatment
Postoperative image in ablation areas intensity profile it is identical as the intensity profile of ablation areas peripheral organs tissue regions;To pre- place
Ablation areas in the tumor region and pretreated postoperative image in pre-operative image after reason carries out image registration, obtains bullet
Property Deformation Field, and to elastic deformation field carry out vector field analysis, determine actual shrinkage central point;According to actual shrinkage central point,
On pretreated postoperative image after tumor region in pretreated pre-operative image to be mapped to image registration, and according to
Tumor region on pretreated postoperative image after being mapped to image registration and pretreated postoperative after image registration
Ablation areas in image obtains the assessment result of ablative surgery therapeutic effect.
The apparatus embodiments described above are merely exemplary, wherein unit can be as illustrated by the separation member
Or may not be and be physically separated, component shown as a unit may or may not be physical unit, i.e.,
It can be located in one place, or may be distributed over multiple network units.It can select according to the actual needs therein
Some or all of the modules achieves the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying creative labor
In the case where dynamic, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Such understanding, above-mentioned skill
Substantially the part that contributes to existing technology can be embodied in the form of software products art scheme in other words, the calculating
Machine software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used
So that a computer equipment (can be personal computer, server or the network equipment etc.) executes above-mentioned each implementation
The method of certain parts of example or embodiment.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features;
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.