Initial commit

First upload of /compton-fastfit/ code onto public Github

.gitattributes

@@ -0,0 +1,2 @@
+# Auto detect text files and perform LF normalization
+* text=auto

.gitignore

@@ -0,0 +1,2 @@
+*.csv
+

Perceptually uniform colormaps/Colormaps (5)/Colormaps/README.txt

@@ -0,0 +1,9 @@
+DISCLAIMER:
+
+The author of this Matlab File Exchange submission has not created any of the colormaps.
+
+Parula has been taken from Matlab >=2014b
+
+py_A-D has been taken from the discussion of the new colormaps for matlplotlib library in python.
+
+Read more about colormaps on python: https://bids.github.io/colormap/

Perceptually uniform colormaps/Colormaps (5)/Colormaps/demo1.m

@@ -0,0 +1,112 @@
+%% Colormaps that are pretty awesome: DEMO1
+% this is an exampe of 4 colormaps that are being considered as the default
+% colormap in python's matplotlib lybrary.
+%
+% All of them look quite good and they dont have any official name, so at
+% the moment they are A,B,C,D.
+%
+% colormaps from https://github.com/bids/colormap
+%
+% Ander Biguri
+%% Clear workspace and get screen data
+clear;
+clc
+close all;
+
+screen=get(0,'ScreenSize') ;
+w0=screen(1);
+h0=screen(2);
+w =screen(3);
+h =screen(4);
+%% Generate sample data
+load flujet
+X=X.';
+
+%% Plot original with jet and parula
+% Parula
+h2=figure('name','Sample data with "parula" colormap');
+set(h2,'position',[w0,h0,w/3,h/2])
+imagesc(X)
+
+if verLessThan('matlab', '8.4')
+   % if parula is the "future"
+   colormap(fake_parula());
+else
+   % if parula is already in Matlab
+   colormap('parula');
+end
+axis image
+xlabel('Parula Colormap')
+set(gca,'xtick',[],'ytick',[]) % This is axis off without offing the labels
+
+% Jet
+h1=figure('name','Sample data with "jet" colormap');
+set(h1,'position',[w0,h/2,w/3,h/2])
+imagesc(X)
+colormap('jet')
+xlabel('jet Colormap')
+set(gca,'xtick',[],'ytick',[]) % This is axis off without offing the labels
+axis image
+
+
+%% Load new colormaps
+
+m=100;
+cm_magma=magma(m);
+cm_inferno=inferno(m);
+cm_plasma=plasma(m);
+cm_viridis=viridis(m);
+
+
+%% Plot new colormaps
+h3=figure('name','Super-cool new colormaps that you can easily use');
+set(h3,'position',[w/3,h0,2*w/3,h])
+if verLessThan('matlab', '8.4')
+% If you are using old Matlab figure engine do  it this way
+% (some very lousy colormap problems before)
+    subplot(2,2,1,'Position',[0.05 0.55 0.4 0.4])
+    subimage(uint8(X/max(X(:))*255),cm_magma)
+    xlabel('MAGMA')
+    set(gca,'xtick',[],'ytick',[]) 
+
+
+    subplot(2,2,2,'Position',[0.55 0.55 0.4 0.4])
+    subimage(uint8(X/max(X(:))*255),cm_inferno)
+    xlabel('INFERNO')
+    set(gca,'xtick',[],'ytick',[]) 
+
+
+    subplot(2,2,3,'Position',[0.05 0.05 0.4 0.4])
+    subimage(uint8(X/max(X(:))*255),cm_plasma)
+    xlabel('PLASMA')
+    set(gca,'xtick',[],'ytick',[])
+
+    subplot(2,2,4,'Position',[0.55 0.05 0.4 0.4])
+    subimage(uint8(X/max(X(:))*255),cm_viridis)
+    xlabel('VIRIDIS')
+    set(gca,'xtick',[],'ytick',[])
+else
+    sp1=subplot(2,2,1,'Position',[0.05 0.55 0.4 0.4]);
+    imagesc(X)
+    colormap(sp1,cm_magma)
+    xlabel('MAGMA')
+    set(gca,'xtick',[],'ytick',[])
+
+    sp2=subplot(2,2,2,'Position',[0.55 0.55 0.4 0.4]);
+    imagesc(X)
+    colormap(sp2,cm_inferno)
+    xlabel('INFERNO')
+    set(gca,'xtick',[],'ytick',[])
+
+    sp3=subplot(2,2,3,'Position',[0.05 0.05 0.4 0.4]);
+    imagesc(X)
+    colormap(sp3,cm_plasma)
+    xlabel('PLASMA')
+    set(gca,'xtick',[],'ytick',[])
+
+    sp4=subplot(2,2,4,'Position',[0.55 0.05 0.4 0.4]);
+    imagesc(X)
+    colormap(sp4,cm_viridis)
+    xlabel('VIRIDIS')
+    set(gca,'xtick',[],'ytick',[])
+end

Perceptually uniform colormaps/Colormaps (5)/Colormaps/demo2.m

@@ -0,0 +1,54 @@
+%% Colormaps that are pretty awesome: DEMO2
+%  In this demo we will show how to supereasily use the new colormaps
+%
+%
+% this is an exampe of 4 colormaps that are being considered as the default
+% colormap in python's matplotlib lybrary.
+%
+% All of them look quite good and they dont have any official name, so at
+% the moment they are A,B,C,D.
+%
+% colormaps from https://github.com/bids/colormap
+%
+% Ander Biguri
+%% Clear workspace and get screen data
+clear;
+clc
+close all;
+
+
+%% Generate sample data
+X=peaks(200);
+
+%% Load Colomaps
+
+jet=colormap('jet');
+parula=fake_parula();
+magma=magma();
+inferno=inferno();
+plasma=plasma();
+viridis=viridis();
+
+
+
+%% Chose colormap
+% Use 1 only, else it will  just use the last
+% CTRL+R -> comment line
+% CTRL+T -> Uncomment line
+
+% colormap(jet);
+% colormap(parula);
+% colormap(magma);
+% colormap(inferno);
+% colormap(plasma);
+colormap(viridis);
+
+%% Plot
+
+for ii=1:0.3:20
+    surf(cos(2*pi*ii/20)*X,'linestyle','none');
+    axis off
+    axis([0 200 0 200 -10 10])
+    set(gca, 'CLim', [-8, 8]);
+    drawnow
+end

Perceptually uniform colormaps/Colormaps (5)/Colormaps/fake_parula.m

@@ -0,0 +1,76 @@
+function cm_data=fake_parula(m)
+
+cm = [[0.2081, 0.1663, 0.5292], 
+    [0.2116238095, 0.1897809524, 0.5776761905], 
+ [0.212252381, 0.2137714286, 0.6269714286],
+ [0.2081, 0.2386, 0.6770857143], 
+ [0.1959047619, 0.2644571429, 0.7279], 
+ [0.1707285714, 0.2919380952, 0.779247619], 
+ [0.1252714286, 0.3242428571, 0.8302714286], 
+ [0.0591333333, 0.3598333333, 0.8683333333], 
+ [0.0116952381, 0.3875095238, 0.8819571429],
+ [0.0059571429, 0.4086142857, 0.8828428571], 
+ [0.0165142857, 0.4266, 0.8786333333],
+ [0.032852381, 0.4430428571, 0.8719571429],
+  [0.0498142857, 0.4585714286, 0.8640571429], 
+ [0.0629333333, 0.4736904762, 0.8554380952],
+ [0.0722666667, 0.4886666667,  0.8467], 
+ [0.0779428571, 0.5039857143, 0.8383714286], 
+ [0.079347619, 0.5200238095, 0.8311809524], 
+ [0.0749428571, 0.5375428571, 0.8262714286],
+  [0.0640571429, 0.5569857143, 0.8239571429], 
+ [0.0487714286, 0.5772238095, 0.8228285714],
+ [0.0343428571, 0.5965809524, 0.819852381], 
+  [0.0265, 0.6137, 0.8135], 
+  [0.0238904762, 0.6286619048, 0.8037619048],
+  [0.0230904762, 0.6417857143, 0.7912666667], 
+ [0.0227714286, 0.6534857143, 0.7767571429],
+ [0.0266619048, 0.6641952381, 0.7607190476], 
+  [0.0383714286, 0.6742714286, 0.743552381], 
+ [0.0589714286, 0.6837571429, 0.7253857143], 
+ [0.0843, 0.6928333333, 0.7061666667], 
+ [0.1132952381, 0.7015, 0.6858571429], 
+ [0.1452714286, 0.7097571429, 0.6646285714], 
+ [0.1801333333, 0.7176571429,  0.6424333333],
+  [0.2178285714, 0.7250428571, 0.6192619048], 
+ [0.2586428571, 0.7317142857, 0.5954285714], 
+ [0.3021714286, 0.7376047619, 0.5711857143],
+ [0.3481666667, 0.7424333333, 0.5472666667], 
+ [0.3952571429, 0.7459, 0.5244428571], 
+ [0.4420095238, 0.7480809524,  0.5033142857],
+ [0.4871238095, 0.7490619048, 0.4839761905], 
+ [0.5300285714, 0.7491142857, 0.4661142857], 
+ [0.5708571429, 0.7485190476, 0.4493904762], 
+  [0.609852381, 0.7473142857, 0.4336857143], 
+ [0.6473, 0.7456, 0.4188],
+ [0.6834190476, 0.7434761905, 0.4044333333], 
+ [0.7184095238, 0.7411333333, 0.3904761905], 
+ [0.7524857143, 0.7384, 0.3768142857], 
+ [0.7858428571, 0.7355666667,  0.3632714286], 
+ [0.8185047619, 0.7327333333, 0.3497904762], 
+ [0.8506571429, 0.7299, 0.3360285714],
+ [0.8824333333, 0.7274333333, 0.3217], 
+ [0.9139333333, 0.7257857143, 0.3062761905],
+ [0.9449571429, 0.7261142857,  0.2886428571],
+  [0.9738952381, 0.7313952381, 0.266647619], 
+ [0.9937714286, 0.7454571429, 0.240347619],
+ [0.9990428571, 0.7653142857,  0.2164142857], 
+  [0.9955333333, 0.7860571429, 0.196652381], 
+ [0.988, 0.8066, 0.1793666667],
+ [0.9788571429, 0.8271428571, 0.1633142857], 
+ [0.9697, 0.8481380952, 0.147452381], 
+ [0.9625857143, 0.8705142857, 0.1309], 
+ [0.9588714286, 0.8949, 0.1132428571], 
+ [0.9598238095, 0.9218333333,  0.0948380952],
+ [0.9661, 0.9514428571, 0.0755333333], 
+ [0.9763, 0.9831, 0.0538]];
+
+if nargin < 1
+    cm_data = cm;
+else
+    hsv=rgb2hsv(cm);
+    cm_data=interp1(linspace(0,1,size(cm,1)),hsv,linspace(0,1,m));
+    cm_data=hsv2rgb(cm_data);
+
+end
+end