MATLAB scripts for Medipix2 analysis

The Medipix2 USB interface produces output files that are just ASCII arrays containing the no. of photons detected by each pixel. These can be loaded into MATLAB and analysed.

Spectrum from Medipix 3D

THLScanSpectrumImages.m - Script takes a series of images taken during one or more THL scans, and uses them to find the spectrum measured by the detector. Additional code in the script can find the peak from a monoenergetic source. If we know the calibration of the detector, we can plot the spectrum vs energy (rather than versus THL values) and calculate total no of hits and fraction of charge shared events.

NoiseAnalysis.m - This script takes a series of images taken during a THL scan without any beam or source present (in a similar fashion to the script above). It records the THL value where each pixel sees 1, 10, 100… noise counts. Then, it plots a graph showing the no. of pixels with 1, 10, 100… noise counts at each THL value. This can be useful for deciding which pixels to mask when doing analysis of spectra.

ThresholdDistributions.m - After performing a threshold equalisation, there is a “save distributions” option which allows you to save the histograms produced during the equalisation. This script loads the histogram of threshold positions after equalisation, and does a Gaussian fit. This can, for example, allow you to find the position of the noise centroid, giving you the “zero” level of the detector after equalisation.

MeanImage.m - If you use Pixelman to take multiple acquisitions, an index will be added to the filename of each image. This is a very simple function to find the mean of the images, by looping through them.

FlatFieldCoefficients.m - After taking a series of flat field images and finding their mean, this generates an matrix containing flat field coefficients. Flat field correction can then be applied in MATLAB by multiplying an image by this matrix. Note that any pixels in the flat-field image containing zero counts cannot be corrected – this function sets the corresponding FF coefficients to 1.

FindMTFFromEdge_testbeam.m - This script takes an image of an edge, and finds the edge spread function from it. It then differentiates this to find the line spread function, and takes the Fourier transform to get the Modulation Transfer Function. (Insert reference.) There are also a few bits of code to extract data directly from the line spread function, such as the full-width half maximum. Ideally, the image of the edge should have been flat-field corrected.

ClusterAnalysis.m - If you use a short acquisition time or a lower-intensity source, so that most pixels have 0 hits and very few have more than 1, it becomes possible to distinguish between hits produced by alphas (large round clusters), betas (smaller, irregular clusters due to scattering) and gammas (photoelectric absorption, hitting only 1 or 2 pixels). See the image below. This example script uses histogramming and image erosion to distinguish between the different clusters. Example ASCII image files are available for a mixed alpha, beta, gamma source, and the same source with alphas removed using a slide.

Mixed source image:

Mixed source image

-- DavidPennicard - 30 Jul 2008

Topic attachments
ISorted ascending Attachment History Action Size Date Who Comment
JPEGjpg 3D1_15keV_keVspectrum_stderror.jpg r1 manage 33.0 K 2008-07-30 - 15:17 DavidPennicard Example spectrum from 3D Medipix
JPEGjpg MixedSourceImage.jpg r1 manage 38.4 K 2008-07-30 - 17:53 DavidPennicard Image with mixed source showing clusters
Unknown file formatm ClusterAnalysis.m r1 manage 2.8 K 2008-07-30 - 17:31 DavidPennicard Cluster analysis script
Unknown file formatm FindMTFFromEdge_testbeam.m r1 manage 15.6 K 2008-07-30 - 17:31 DavidPennicard Modulation Transfer Function script
Unknown file formatm FlatFieldCoefficients.m r1 manage 1.9 K 2008-07-30 - 17:30 DavidPennicard Flat field script
Unknown file formatm MeanImage.m r1 manage 1.2 K 2008-07-30 - 17:30 DavidPennicard Mean image script
Unknown file formatm NoiseAnalysis.m r1 manage 4.7 K 2008-07-30 - 17:30 DavidPennicard Noise analysis script
Unknown file formatm THLScanSpectrumImages.m r1 manage 13.9 K 2008-07-30 - 17:30 DavidPennicard Spectrum-finding script
Unknown file formatm ThresholdDistributions.m r1 manage 4.2 K 2008-07-30 - 17:29 DavidPennicard Threshold distribution script
Texttxt Image_aby.txt r1 manage 128.2 K 2008-07-30 - 17:53 DavidPennicard Alpha, beta, gamma ASCII file
Texttxt Image_by.txt r1 manage 128.2 K 2008-07-30 - 17:54 DavidPennicard Beta, gamma ASCII file
Edit | Attach | Print version | History: r6 | r4 < r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r2 - 2008-07-30 - DavidPennicard
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback