核医学影像中的数据处理课件

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,核医学影像中的数据处理,中国科学院高能物理研究所,北京市射线成像技术与装备工程技术研究中心,贠明凯,核医学影像中的数据处理中国科学院高能物理研究所,Modern Nuclear Medical Imaging,Acquire,Process,Apply,Scanners Computers,Users,Modern Nuclear Medical Imaging,Outline,Data organization,Correction methods,Rebinning,Image reconstruction,Image registration and fusion,DICOM and PACS,Outline Data organization,Outline,Data organization,Correction methods,Rebinning,Image reconstruction,Image registration and fusion,DICOM and PACS,Outline Data organization,Data organization,List mode,Histgram,Sinogram,Linogram,Data organizationList mode,SinogramPET,Sinogram,r,Projections and Sinogram,SinogramPETSinogram rProjec,Sinogram,PET,Sinogram,r,Projections and Sinogram,SinogramPETSinogramrProject,Sinogram,SPECT,SinogramSPECT,2D,VS.,3D,Septa between crystal rings,Lower sensitivity,Lower random,Lower scatter,2D reconstruction,No septa,Higher sensitivity,Higher random,Higher scatter,3D reconstruction or hybrid reconstruction,2D VS. 3DSepta between crysta,Outline,Data organization,Correction methods,Rebinning,Image reconstruction,Image registration and fusion,DICOM and PACS,Outline Data organization,Scatter Coincidence,Trues Coincidence,Random Coincidence,True Counts & Noise,Scatter CoincidenceTrues Coinc,Normalization,A,B,C,D,Attenuation,A,B,C,D,Scatter,A,B,C,D,Need to correct the data,NormalizationABCDAttenuationAB,Correction methods,random,“dead time”,normalization,scatter,attenuation,decay,Arc correction,Depth of interaction,Motion correction,Partial volume,Axial of rotation,Camera head tilt,Correction methodsrandomArc co,Random,Finite time window with,Energy window,Coincidence timing window,Activity,RandomFinite time window with,Random,Tail fitting,simplest,Small changes in tail, great changes in estimate,Estimation from singles rates,Measure the single count rate on each detector for a given time window,Subtracting from the prompts between detector pair,Singles rate is much larger than that of coincidence events,Single rates change in the same way over time,Random Tail fitting,Delayed coincidence channel estimation,One channel is delayed before being sent to coincidence processing,Subtracted form prompt coincidences,Advantage,Accurate,Same dead time environment as prompt channel,Disadvantage,Increased system dead time,Doubling of the statistical noise due to random,Delayed coincidence channel es,Dead time correction,Decaying source experiment is performed,Dead time correctionDecaying s,Dead time correction (con),Look up table,Uniform source,Known quantity,Short lived,Linear extrapolation from count rate for a given level of activity,Dead time correction (con)Look,Normalization,Causes of sensitivity variations,Summing of adjacent data elements,Detector efficiency variations,Geometric and solid angle effects,Rotational sampling,Time window alignment,Structural alignment,septa,Normalization Causes of sensit,Summing of adjacent data elements,Geometric and solid angle effects,Summing of adjacent data eleme,Rotational sampling,LOR at the edge are sampled less than LOR close to the center,Rotational samplingLOR at the,Crystal interface factors,Time window alignment factors,Crystal interface factorsTime,Normalization methods (con),Direct normalization,Simplest approach,Adequate statistical quality,Very uniform activity sources,Scatter in normalization should be substantially different from normal imaging,Normalization methods (con)Dir,Normalization methods (con),Component-based normalization,Normalization methods (con)Com,Scatter correction,LORs recorded outside object boundary can only be explained by scatter,The scatter distribution is very broad,Scattered coincidences fall within the photo-peak window mainly due to scattered once,Scatter correctionLORs recorde,Scatter correction,Energy spectra distribution of scattered 511KeV photons according to the number of times each photon scatters,Scatter correctionEnergy spect,Scatter correction,Empirical scatter corrections,Fitting the scatter tails,Direct measurement technique,Energy window techniques,Dual energy window methods,Multiple energy window methods,Convolution and de-convolution,Simulation-based scatter correction,Analytical simulation,Monte Carlo simulation,Scatter correctionEmpirical sc,Fitting the scatter tails,Simplest approach,Fit an analytical function to scatter tails,Second order polynomial or 1D Gaussian,Coincidences outside the object are entirely scatter events,Not always well approximated, particularly in thorax,Fitting the scatter tailsSimpl,Direct measurement technique,Only applicable to PET with retractable septa,Steps,Make a measurement of the same object with and without septa,Scaling septa extended projections for different efficiency,Subtract from projections of polar angle 0,Estimate the oblique scatter by interpolation of the direct plane scatter,Direct measurement techniqueOn,Dual energy window methods,Dual energy window methods,Dual energy window methods,Dual energy window methods,Multiple energy window methods,Multiple energy window methods,Scatter Correction,Scatter Correction,Analytical simulation,Analytical simulation,Scatter Correction,A,B,Single Scatter - Model based correction,Calculate the contribution for an arbitrary scatter point using the Klein-Nishina equation,Before,Scatter,correction,After,Scatter,correction,Scatter CorrectionABSingle Sca,Attenuation correction,Attenuation in the body is equal to that of source lying along the same LOR,Attenuation correctionAttenuat,Zaidi H, Hasegawa B. J Nucl Med 2003; 44:291-315.,SPECT,PET,Zaidi H, Hasegawa B. J Nucl Me,Attenuation correction (con),Measured attenuation correction,Coincidence transmission data,Long-lived positron emitter,Normally more than one rod source are used,Sinogram windowing is applied provide location of rod,Impractical in 3D,Singles transmission data,Shielded point transmission source,Separate blank scan is needed,Significant scatter and broad beam,Attenuation correction (con)Me,Measured attenuation correction,Coincidence measurement using rod source,Transmission measurement using point source,Measured attenuation correctio,CT scan,Advantage,High statistical quality,High spatial resolution,Significant reduction in scan time,Disadvantage,Faster CT, slower PET,Smaller FOV of CT,Difficulty in registration,values do not scale linearly,CT scanAdvantage,Attenuation correction for PET,Types of transmission images,Coincident photon Ge-68/Ga-68,(511 keV),high noise,15-30 min scan time,low bias,low contrast,Single photon Cs-137,(662 keV),lower noise,5-10 min scan time,some bias,lower contrast,X-ray,(30-140kVp),no noise,1 min scan time,potential for bias,high contrast,Attenuation correction for PET,Other attenuation correction methods,Calculated attenuation,Regular geometric outline,Constant tissue,Segmented attenuation,Segment transmission image according to tissue type,Assigning known attenuation coefficients,Forward projection,Other attenuation correction m,attenuation correction,attenuation correction,Attenuation/Scatter correction,University of Pennsylvania PET Center,No AC or Scatter Corr,AC and Scatter Corr,Philips,Allegro,Attenuation/Scatter correction,Arc correction,Different sampling distance at different radial position,Equal sampling distance is required in analytical method,Interpolation method,Nearest interpolation,Linear interpolation,B-spline interpolation (negative values!),Arc correctionDifferent sampli,DOIdepth of interaction,DOIdepth of interaction,DOIdepth of interaction,（,con,）,DOIdepth of interaction（con）,Dual Layer,Dual Layer,A Point Spread Function (PSF) describes the response of an imaging system to a point source or point object. A system that knows the response of a point source from everywhere in its field of view can use this information to recover the original shape and form of imaged objects.,PSFs are used in precision imaging instruments, such as microscopy, ophthalmology, and astronomy (e.g. the Hubble telescope) to make geometric corrections to the final image.,Point Spread Function (PSF),A Point Spread Function (PSF),Motion correction,Cardiac motion and respiration,Motion correctionCardiac motio,Motion correction,（,con,）,Gated frames,List mode,Motion correction（con）Gated fr,Respiratory motion is distributed throughout the whole body,Impact is rarely on detection, but often affects quantitation,Static wholebody,Single respiratory phase,(1 of 7, so noisier), 1 cc lesion on CT,Whole-body respiratory gated PET/CT: Patients,Respiratory motion is distribu,Partial volume effect,Characters,Object or structure being imaged only partially occupies the sensitive volume of scanner,Signal amplitude becomes diluted with signals from surrounding structures,The degree of underestimation of radioactivity concentration will depend not only on its size but also on the relative concentration in surrounding structures,Correction methods,Resolution recovery,Use of anatomical imaging data,Partial volume effectCharacter,A Point Spread Function (PSF) describes the response of an imaging system to a point source or point object. A system that knows the response of a point source from everywhere in its field of view can use this information to recover the original shape and form of imaged objects.,PSFs are used in precision imaging instruments, such as microscopy, ophthalmology, and astronomy (e.g. the Hubble telescope) to make geometric corrections to the final image.,Point Spread Function (PSF),A Point Spread Function (PSF),核医学影像中的数据处理课件,Partial volume effectMAP,Partial volume effectMAP,assumptions:,camera moves along,circular orbit,orbit is,reproducible,calibration method finds system geometry,assumptions:,problem 1: tilting detector,assumption: camera moves along circular orbit,problem 1: tilting detectorass,AORAxial of rotation,Offset of AOR,Rotation of AOR,Nutation of AOR,AORAxial of rotationOffset of,Camera head tilt,Heads need to be exactly parallel to axis of rotation,Correct alignment,Head tilt,Camera head tiltHeads need to,pinhole calibration,Dirk Bequ, Kathleen Vunckx,pinhole calibrationDirk Bequ,circular orbit,circular orbit + new model,extension 2: circular orbit +,arbitrary,small deviations,measurement,model,Michel Defrise, Chris Vanhove,circular orbitcircular orbit +,extension 2: circular orbit +,arbitrary,small deviations,old,new,translations,rotations,1mm,-3mm,1.5mm,-1.5mm,1.5mm,-1mm,1,o,-2,o,1.5,o,-1.5,o,3.5,o,-2.5,o,1mm,1.2,1.4,1.6,1.8,2mm,extension 2: circular orbit +,Outline,Data organization,Correction methods,Rebinning,Image reconstruction,Image registration and fusion,DICOM and PACS,Outline Data organization,Rebinning,Convert 3D data to 2D,RebinningConvert 3D data to 2D,SSRB and MSRB,SSRB-,Single-slice,rebinning,Detection: center slice,Simple,Fast,Resolution loss,MSRB-,Multi-slice rebinning,Distribute along all intermediate slices,De-blurring along z-axis,SSRB and MSRBSSRB- Single-slic,Fourier rebinning,Fourier rebinning,Outline,Data organization,Correction methods,Rebinning,Image reconstruction,Image registration and fusion,DICOM and PACS,Outline Data organization,Image reconstruction,Analytical,FBP,BPF,FDK,3D RP,Iterative,ART,MLEM,OSEM,OSLS,MAP,Image reconstructionAnalytical,Analytical algorithms,For example, FBP (Filtered Back-projection),Treat the unknown image as continuous,Point-by-point reconstruction,Regular grid points are commonly chosen,Treat projection process as line integral theoretically,Analytical algorithmsFor examp,解析重建,-FBP,FBP,解析重建-FBPFBP,back projection (BP) = summation of projections,back projection (BP) = summati,filtered back projection (FBP),filtered back projection (FBP),FDK,Feldkamp,、,Davis,、,Kress,FDKFeldkamp、Davis、Kress,FDK,FDK,3D RPRe-projection,3D RPRe-projection,Steps of 3D RP,Extract 2D sinograms,Reconstruct each with 2D FBP and stack to form 3D image,Forward project to calculate missing LORs,Extract 2D projection data of all oblique slices,Take 2D Fourier transform,Back project data through 3D image matrix,Repeat for all angles and oblique slices,Steps of 3D RPExtract 2D sinog,What is iterative reconstruction,Discrete measurements, discrete image,Optimization,What is iterative reconstructi,Attractions of iterative methods,Either consistent or inconsistent is OK,Complex geometry,Physical effects and detection processes can be modeled,Non-negativity,Great reducing streaking artifacts,Better contrast recovery,Attractions of iterative metho,Classification of iteration reconstruction methods,ART (algebraic reconstruction techniques),MART (multiplicative ART),AART (additive ART),SIRT (simultaneous iterative reconstruction),SMART (simultaneously MART),BI-ART (block iterative ART),BI-SMART (block iterative SMART),RBI-SMART (rescaled BI-SMART),Classification of iteration re,Statistical algorithms,MAP:,Maximize the conditional probability,P(image|data),MLEM:,Maximize the probability,P(data|image),Statistical algorithmsMAP:,Statistical algorithms,Gaussian assumption,P is projection column matrix, A is system matrix, F image column matrix, C is the covariance matrix of the data,Assumed all standard deviations are identical and equal to 1, idealized parallel projection, perfect resolution and no attenuation or other degrading affects,Statistical algorithmsGaussia,Statistical algorithms,Poisson assumption,Statistical algorithmsPoisson,实测数据,迭代重建,-MLEM&OSEM,正投影,比较,更新,重建,MLEM, OSEM, .,likelihood,iteration,实测数据迭代重建-MLEM&OSEM正投影比较更新重建MLE,Sinogram,r,Subset 1,Subset 2,Subset 3,Subset 4,1, 3 2 4,Subset order,SinogramrSubset 1Subset 2Subs,0,1,2,3,4,10,40,ordered subsets,1 iteration of 40 subsets,(2 proj per subset),012341040ordered subsets1 iter,System matrix,Scan geometry,Collimator/detector response,Attenuation,Scatter (object, collimator, scintillator),Duty cycle (dwell time at each angle),Detector efficiency,Dead-time losses,Positron range,Non-colinearity,Crystal penetration,System matrixScan geometry,Considerations of system matrix,Quantitative accuracy,Spatial accuracy,Computation time,Storage space,Model uncertainties,Artifacts due to over simpleifications,Considerations of system matri,System matrix tricks,Factorize,Symmetry,Sparseness,Approximation,Partial Monte Carlo,System matrix tricksFactorize,System matrix model,System matrix model,Reconstruction image of uniform source,Reconstruction image of unifor,FBP,VS.,OSEM,FBPanalytical,Pros,：,Single pass,Linear,Fast,Cons,：,Streak artifact,Poor resolution,Correction not built-in,OSEMiteration,Pros,：,Better resolution,Better contrast,Lower noise,Cons,：,Extensive time consuming,Memory consuming,Required user training,FBP VS. OSEMFBPanalyticalOSEM,FBP,VS.,OSEM,Phantom test (left),Clinical results (right),FBP VS. OSEMPhantom test (left,Outline,Data organization,Correction methods,Rebinning,Image reconstruction,Image registration and fusion,DICOM and PACS,Outline Data organization,Image Registration,PET,CT,PET/CT,Image RegistrationPETCTPET/CT,Voxel based image registration,Image Registration,Voxel based image registration,Image Registration,算法流程图,相似性测量一般用到的函数有,：,相同模态图像：残差（,sum of square difference,）,不同模态图像：互信息（,mutual information,）,一般用来做配准的优化算法有：,六参数或十二参数的优化一般使用,Powell,优化算法,多参数优化一般使用,LBFGS,（,limited-memory BroydenFletcherGoldfarbShanno,）优化算法（由牛顿算法演变而来）,Image Registration算法流程图相似性测量,Image Fusion,Alpha Blending based,Adaptive alpha blending,Alpha blending,Adaptive Alpha blending,Image FusionAlpha Blending bas,Outline,Data organization,Correction methods,Rebinning,Image reconstruction,Image registration and fusion,DICOM and PACS,Outline Data organization,DICOM and PACS,DICOM,Digital image and Communication in Medicine,Created by ACR (American College of Radiology) and NEMA (National Electrical Manufacturers Association) in 1985,Two components:,Communication protocols and image format standards,DICOM and PACSDICOM,DICOM image,VS.,General image,Structure:,DICOM contains header and image data sections;,Other image file such as BMP, JPG, TIFF, which contain also two sections;,Size:,The header size of DICOM is variable;,The header size of many general image is constant,Contents:,DICOM contains additional patients data such as basic information, study information and so on;,General image header describes basic image parameters, such as image size, compression type,DICOM image VS. General imageS,DICOM logical layers,Patient information:,Patients name, patients ID, patients birth date, hospital information system,Study information:,Study information such as dose, injection time and additional examination information such as study name study date;,Series information:,Series ID, manufacturer and institution name,Image information:,The size of pixel, image size, pixel value and how it is encoded,DICOM logical layersPatient in,PACS,PACS,The Picture Archival and Communication System,A system for storage of images and transferring images between computers in different facilities through networks,PACS is helpful to provide comparative studies among different image modalities,PACSPACS,PACS (con),PACS (con),Thanks,Thanks,Our Group is,growing up,Our Group is growing up,TOF-PET,TOF,信息,TOF-PETTOF信息,PET,原始数据,偶然符合校正,发射图像,衰变校正,死时间校正,弓形几何校正,衰减校正,散射校正,归一化校正,图像重建,数据重组,数据校正,PET原始数据偶然符合校正发射图像衰变校正死时间校正弓形几何,