Fast GPU-based computation of spatial multigrid multiframe LMEM for PET.

Significant efforts were invested during the last decade to accelerate PET list-mode reconstructions, notably with GPU devices. However, the computation time per event is still relatively long, and the list-mode efficiency on the GPU is well below the histogram-mode efficiency. Since list-mode data...

Descripción completa

Detalles Bibliográficos
Publicado en:Medical & Biological Engineering & Computing Vol. 53; no. 9; pp. 791 - 804
Autores principales: Nassiri, Moulay, Carrier, Jean-François, Després, Philippe, Nassiri, Moulay Ali
Formato: diagnostic images equations & formulas tables/charts Journal Article
Publicado: Springer Nature Sep2015
Acceso en línea:Ver este registro en EBSCOhost
fields @attributes:
  recordID: 1
pdfLink:
plink: https://search.ebscohost.com/login.aspx?direct=true&db=ccm&AN=109466534&site=ehost-live
header:
  @attributes:
    shortDbName: ccm
    uiTerm: 109466534
    longDbName: CINAHL Complete
    uiTag: AN
  controlInfo:
    bkinfo:
    dissinfo:
    jinfo:
      jid:
        01400118
        PO0
      jtl: Medical & Biological Engineering & Computing
      issn: 01400118
      maglogo: N
    pubinfo:
      dt: Sep2015
      vid: 53
      iid: 9
      pid: 237
      pub: Springer Nature
      place: New York, New York
    artinfo:
      ui:
        109466534
        109466534
        NLM25850980
        109466534
        10.1007/s11517-015-1284-9
        NLM25850980
        109466534
      ppf: 791
      ppct: 13
      formats:
        fmt:
          @attributes:
            type: P
      tig:
        atl: Fast GPU-based computation of spatial multigrid multiframe LMEM for PET.
      aug:
        au:
          Nassiri, Moulay
          Carrier, Jean-François
          Després, Philippe
          Nassiri, Moulay Ali
          Carrier, Jean-François
          Després, Philippe
        affil: Département de radio-oncologie, Centre hospitalier de l'Université de Montréal (CHUM), Montréal Canada
      sug:
        subj:
          Computer Simulation
          Computer Graphics
          Tomography, Emission-Computed
          Algorithms
          Probability
          Time Factors
          Systems Analysis
          Reproducibility of Results
          Image Processing, Computer Assisted
          Sensitivity and Specificity
          Phantoms, Imaging
      ab: Significant efforts were invested during the last decade to accelerate PET list-mode reconstructions, notably with GPU devices. However, the computation time per event is still relatively long, and the list-mode efficiency on the GPU is well below the histogram-mode efficiency. Since list-mode data are not arranged in any regular pattern, costly accesses to the GPU global memory can hardly be optimized and geometrical symmetries cannot be used. To overcome obstacles that limit the acceleration of reconstruction from list-mode on the GPU, a multigrid and multiframe approach of an expectation-maximization algorithm was developed. The reconstruction process is started during data acquisition, and calculations are executed concurrently on the GPU and the CPU, while the system matrix is computed on-the-fly. A new convergence criterion also was introduced, which is computationally more efficient on the GPU. The implementation was tested on a Tesla C2050 GPU device for a Gemini GXL PET system geometry. The results show that the proposed algorithm (multigrid and multiframe list-mode expectation-maximization, MGMF-LMEM) converges to the same solution as the LMEM algorithm more than three times faster. The execution time of the MGMF-LMEM algorithm was 1.1 s per million of events on the Tesla C2050 hardware used, for a reconstructed space of 188 x 188 x 57 voxels of 2 x 2 x 3.15 mm3. For 17- and 22-mm simulated hot lesions, the MGMF-LMEM algorithm led on the first iteration to contrast recovery coefficients (CRC) of more than 75 % of the maximum CRC while achieving a minimum in the relative mean square error. Therefore, the MGMF-LMEM algorithm can be used as a one-pass method to perform real-time reconstructions for low-count acquisitions, as in list-mode gated studies. The computation time for one iteration and 60 millions of events was approximately 66 s.
      pubtype: Academic Journal
      doctype:
        diagnostic images
        equations & formulas
        tables/charts
        Journal Article
      ougenre: Article
    language: English
    refInfo:
    holdings:
      @attributes:
        islocal: N