Difference: MedipixMATLABScripts (4 vs. 5)

Revision 52009-04-08 - DavidPennicard

Line: 1 to 1
 
META TOPICPARENT name="MedipixUSBReadout"

MATLAB scripts for Medipix2 analysis

Line: 10 to 10
  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.
Changed:
<
<
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.
>
>
FullThresholdDistributions.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 before and after equalisation, and does Gaussian fits. 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 you take a series of flat field images and find their mean, this can be used to generate a 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.

Changed:
<
<
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.
>
>
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. For more info on the MTF, see Fujita 1992 http://dx.doi.org/10.1109/42.126908 and Moy 2000 http://dx.doi.org/10.1118/1.598859 . 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.
  NoisePowerSpectrum.m - As discussed on the previous page, to measure the frequency response of a detector you need to consider the signal-to-noise ratio, rather than the MTF alone. For example, some processes in a detector such as charge sharing will "blur" both the signal and the noise, whereas others like poor optics will blur the signal without reducing the noise. This script is used to calculate the noise power spectrum from a large number of flood images.
Line: 34 to 34
 
META FILEATTACHMENT attachment="ThresholdDistributions.m" attr="" comment="Threshold distribution script" date="1217438977" name="ThresholdDistributions.m" path="ThresholdDistributions.m" size="4309" stream="ThresholdDistributions.m" tmpFilename="/usr/tmp/CGItemp28978" user="DavidPennicard" version="1"
META FILEATTACHMENT attachment="THLScanSpectrumImages.m" attr="" comment="Spectrum-finding script" date="1217439014" name="THLScanSpectrumImages.m" path="THLScanSpectrumImages.m" size="14226" stream="THLScanSpectrumImages.m" tmpFilename="/usr/tmp/CGItemp31628" user="DavidPennicard" version="1"
META FILEATTACHMENT attachment="MeanImage.m" attr="" comment="Mean image script" date="1217439028" name="MeanImage.m" path="MeanImage.m" size="1245" stream="MeanImage.m" tmpFilename="/usr/tmp/CGItemp31727" user="DavidPennicard" version="1"
Changed:
<
<
META FILEATTACHMENT attachment="FlatFieldCoefficients.m" attr="" comment="Flat field script" date="1217439041" name="FlatFieldCoefficients.m" path="FlatFieldCoefficients.m" size="1904" stream="FlatFieldCoefficients.m" tmpFilename="/usr/tmp/CGItemp31680" user="DavidPennicard" version="1"
>
>
META FILEATTACHMENT attachment="FlatFieldCoefficients.m" attr="" comment="" date="1239210111" name="FlatFieldCoefficients.m" path="FlatFieldCoefficients.m" size="2074" stream="FlatFieldCoefficients.m" tmpFilename="/usr/tmp/CGItemp44684" user="DavidPennicard" version="2"
 
META FILEATTACHMENT attachment="NoiseAnalysis.m" attr="" comment="Noise analysis script" date="1217439053" name="NoiseAnalysis.m" path="NoiseAnalysis.m" size="4825" stream="NoiseAnalysis.m" tmpFilename="/usr/tmp/CGItemp31658" user="DavidPennicard" version="1"
META FILEATTACHMENT attachment="FindMTFFromEdge_testbeam.m" attr="" comment="" date="1227700078" name="FindMTFFromEdge_testbeam.m" path="FindMTFFromEdge_testbeam.m" size="17086" stream="FindMTFFromEdge_testbeam.m" tmpFilename="/usr/tmp/CGItemp29954" user="DavidPennicard" version="2"
META FILEATTACHMENT attachment="ClusterAnalysis.m" attr="" comment="Cluster analysis script" date="1217439082" name="ClusterAnalysis.m" path="ClusterAnalysis.m" size="2865" stream="ClusterAnalysis.m" tmpFilename="/usr/tmp/CGItemp31677" user="DavidPennicard" version="1"
Line: 42 to 42
 
META FILEATTACHMENT attachment="Image_aby.txt" attr="" comment="Alpha, beta, gamma ASCII file" date="1217440437" name="Image_aby.txt" path="Image_aby.txt" size="131328" stream="Image_aby.txt" tmpFilename="/usr/tmp/CGItemp31848" user="DavidPennicard" version="1"
META FILEATTACHMENT attachment="Image_by.txt" attr="" comment="Beta, gamma ASCII file" date="1217440460" name="Image_by.txt" path="Image_by.txt" size="131328" stream="Image_by.txt" tmpFilename="/usr/tmp/CGItemp31788" user="DavidPennicard" version="1"
META FILEATTACHMENT attachment="NoisePowerSpectrum.m" attr="" comment="" date="1227700184" name="NoisePowerSpectrum.m" path="NoisePowerSpectrum.m" size="4804" stream="NoisePowerSpectrum.m" tmpFilename="/usr/tmp/CGItemp35005" user="DavidPennicard" version="1"
Added:
>
>
META FILEATTACHMENT attachment="FullThresholdDistributions.m" attr="" comment="" date="1239209876" name="FullThresholdDistributions.m" path="FullThresholdDistributions.m" size="5037" stream="FullThresholdDistributions.m" tmpFilename="/usr/tmp/CGItemp41216" user="DavidPennicard" version="1"
 
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