Image Reconstruction: PET data, backprojection, sinograms, list-mode, convolution

Поділитися
Вставка
  • Опубліковано 16 гру 2024

КОМЕНТАРІ • 23

  • @learn2444
    @learn2444 Рік тому +2

    Simple,clear and highly intuitive. Thanks for this contribution

  • @aljamela-is8tb
    @aljamela-is8tb 2 роки тому +1

    amazing brief description of the convolution

  • @AndrewChicco
    @AndrewChicco Рік тому +1

    Thank you for this well thought out & delivered series of lectures (as all of your lectures are). Could you perhaps consider adding in the (i of j) to the title of each, as it would assist in watching them in sequence.

    • @AndrewJReader
      @AndrewJReader  Рік тому +1

      Many thanks for your feedback.
      In the first instance, perhaps this playlist could help: ua-cam.com/play/PL557uxcMh3xx4XdfAf7eoO94f6AEr8Qun.html

    • @AndrewJReader
      @AndrewJReader  Рік тому +1

      More up to date playlist:
      ua-cam.com/play/PL557uxcMh3xy_lmu7kgfKM_nY31MDc8UM.html

  • @deepaphysics
    @deepaphysics Рік тому +1

    Hello Professor Andrew
    I am beginner in image reconstruction and was trying to obtain sinogram from the LOR data I got from simulations. I know the underlying concept but putting it into code is becoming difficult for me. I have the coordinates of coincident points. How can I generate sinogram using that ? I am facing difficulty in making a matrix where we can assign value of intensity for a given point. Kindly help me in this regard, I will be really grateful.

    • @AndrewJReader
      @AndrewJReader  Рік тому

      Hi, thanks for your message. I suggest you take a look at the newer version of this video - towards the end of the video I do give equations at the 36 minutes mark. Here is the link: ua-cam.com/video/8LTi0auyz2k/v-deo.htmlsi=BoGpmDftQAm3RYMF
      Hope that helps, let me know

    • @deepaphysics
      @deepaphysics Рік тому

      @@AndrewJReader Thank you for your reply professor. I understand the equation, will try to embed it. But kindly tell me that when plotting sinogram do we use simple matplotlib function ? How should I plot to get appropriate inverse radon transformation?

  • @morhafhaidar4949
    @morhafhaidar4949 Рік тому

    Thank you for the Video. I have a question: how can i apply that to cone beam ct. what is the training data and the ground truth in this case?

    • @AndrewJReader
      @AndrewJReader  Рік тому

      Thanks for the comment. In this video there is no machine learning, so no training data or ground truth labels needed. Cone beam CT would have a different ray geometry, and you would need a CT reconstruction method, so all a bit different from this video. You could look up deep learning methods for CT image reconstruction which might more closely match what you need.

  • @engels_diary
    @engels_diary 2 роки тому

    Great presentation!!

  • @anshuchoudhary5454
    @anshuchoudhary5454 Рік тому

    hello professor
    this video is extremely useful. I request you to kindly make one video on Fourier Rebinning Algorithm. I am facing some difficulty in understanding and implementing that. With your clear explanation, it could really help me a ton.

    • @AndrewJReader
      @AndrewJReader  Рік тому

      Hi, thanks for the feedback, really appreciated. For info, the latest version of the video is here: ua-cam.com/video/8LTi0auyz2k/v-deo.html
      For Fourier rebinning, I did indeed write my own version many years ago, but I do not use it anymore, either doing fully 3D or 2D only. Sorry to say that I will not plan to cover that topic at the moment. I guess you plan to do FORE+FBP?

    • @anshuchoudhary5454
      @anshuchoudhary5454 Рік тому

      @@AndrewJReader Oh alright professor. I was validating a research article using GEANT4 simulations where they used FORE +FBP for converting 3D sinogram to 2D and then getting the image. What method do you use to get resolution of PET scanner. Is it directly calculable from 3D sinograms.

  • @Dalat_Instruments
    @Dalat_Instruments Рік тому

    Hello, How to calculate S and phi angle? 5:00

    • @AndrewJReader
      @AndrewJReader  Рік тому +1

      Hi, you need to use equations based on the 2 detection coordinates of an event. So, something like phi = arctan ( difference in y positions / difference in x positions). Then s would be something like s = - x sin phi + y cos phi. Hope that helps?

    • @Dalat_Instruments
      @Dalat_Instruments Рік тому

      @@AndrewJReader I mean determine phi and S based on different detector pairs, from which to construct sinogram.

    • @AndrewJReader
      @AndrewJReader  Рік тому +1

      @@Dalat_Instruments Yes, that's also what I meant. You need to know the coordinates for each detector, (x1,y1), (x2,y2), that make up a detector pair. Given the coordinates of the 2 detectors, that allows the line joining the 2 detectors to be found (i.e. calculating s and phi).

    • @Dalat_Instruments
      @Dalat_Instruments Рік тому

      @@AndrewJReader detectors are arranged in a circle, why are the coordinates (x,y)?

    • @AndrewJReader
      @AndrewJReader  Рік тому +1

      @@Dalat_Instruments Hi, using (x,y) coordinates is more general for various scanners (I used to work with rotating planar detector PET systems). You can of course use a different coordinate system, but if you do, then the way of calculating phi and s will differ. I gave equations above for the general case of specifying (x,y) coordinates. To find s and phi just involves solving a geometric problem in the coordinate system that you choose.