sigprocexperts.com Git - sie_pcm/blob - ReadMe.md

   1
   2 # Subspace Identification Extension to the Phase Correlation Method (SIE-PCM) demo
   3
   4 This short bit of  Matlab code was the original implementation of the method described in
   5 ``A subspace identification extension to the phase correlation method,''
   6 IEEE Trans. Medical Imaging, 22(2):277-280, Feb. 2003.
   7 http://dx.doi.org/10.1109/TMI.2002.808359
   8
   9 Images from the original 5 test frame data are included as well.
  10
  11 # File List
  12
  13     run_example.m - start here
  14     siepcm.m      - a function to estimate the row and column translation between two images
  15
  16 The images are from 5 different MRI acquisition of a grapefruit, with
  17 the grapefruit at different locations within the field-of-view.  In
  18 the table below, the coordinates of one corner of the field of view
  19 are given.
  20
  21      Image File | Source file |  dim1 loc    | dim2 loc
  22     ------------|-------------|--------------|------------
  23      img1.png   | P64512.7    | y = 8.30 cm  | z= 2.00 cm
  24      img2.png   | P00000.7    | y = 8.15 cm  | z= 2.25 cm
  25      img3.png   | P00512.7    | y = 8.00 cm  | z= 2.50 cm
  26      img4.png   | P01024.7    | y = 7.85 cm  | z= 2.27 cm
  27      img5.png   | P01536.7    | y = 7.85 cm  | z= 2.75 cm
  28
  29 To determine the expected translational shift between the images, use the
  30 difference between the coordinates divided by the distance per pixel.
  31
  32                      dy = (y2 - y1) / n  where n=0.0625
  33
  34 The original data was acquired on a GE Signa Lx 1.5T MR Scanner using
  35 the following parameters:
  36  TR: 500  TE: 15   protocol: FSE  coil: HEAD
  37
  38 ----------------------------------------------------------------------
  39
  40 FOV length : Image Size :: 160mm : 256 pixels; => 0.0625 mm/pixel
  41
  42  compare images 1 and 2.  S: -2.5 mm, A: 1.5 mm
  43         expect z =  -4.0000,  y =  +2.4000
  44
  45  compare images 1 and 3.  S: -5 mm, A: 3 mm
  46         expect z =  -8.0000,  y =  +4.8000
  47
  48  compare images 1 and 4.  S: -2.7 mm, A: 4.5 mm
  49         expect z =  -4.3200,  y =  +7.2000
  50
  51  compare images 1 and 5.  S: -7.5 mm, A: 4.5 mm
  52         expect z = -12.0000,  y =  +7.2000
  53
  54  compare images 2 and 3.  S: -2.5 mm, A: 1.5 mm
  55         expect z =  -4.0000,  y =  +2.4000
  56
  57  compare images 2 and 4.  S: -0.2 mm, A: 3 mm
  58         expect z =  -0.3200,  y =  +4.8000
  59
  60  compare images 2 and 5.  S: -5 mm, A: 3 mm
  61         expect z =  -8.0000,  y =  +4.8000
  62
  63  compare images 3 and 4.  S: 2.3 mm, A: 1.5 mm
  64         expect z =  +3.6800,  y =  +2.4000
  65
  66  compare images 3 and 5.  S: -2.5 mm, A: 1.5 mm
  67         expect z =  -4.0000,  y =  +2.4000
  68
  69  compare images 4 and 5.  S: -4.8 mm, A: 0 mm
  70         expect z =  -7.6800,  y =  +0.0000
  71
  72 ----------------------------------------------------------------------
  73
  74 Copyright 2003, W. Scott Hoge (wsh032580 at proton dot me)
  75
  76 This images and text in this zip file are licensed under the Creative
  77 Commons Attribution License. To view a copy of this license, visit
  78 http://creativecommons.org/licenses/by/1.0/ or send a letter to Creative
  79 Commons, 559 Nathan Abbott Way, Stanford, California 94305, USA.
  80
  81 You are free:
  82
  83     * to copy, distribute, display, and perform the work
  84     * to make derivative works
  85     * to make commercial use of the work
  86
  87 Under the following conditions:
  88 by Attribution. You must give the original author credit.
  89
  90     * For any reuse or distribution, you must make clear to others the
  91       license terms of this work.
  92     * Any of these conditions can be waived if you get permission from the
  93       author.
  94
  95 Your fair use and other rights are in no way affected by the above.
  96
  97 This is a human-readable summary of the Legal Code
  98 (the full license: http://creativecommons.org/licenses/by/1.0/legalcode)
  99