1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
|
#!/bin/bash
# remove whitespace and braces,
# filter out panasonic and powershot etc because they usually blow up the plot ranges
# also only choose nikon d800 and canon 5d mark 2 for comparison
grep '"name"' data/noiseprofiles.json | grep -v '"skip"[ \t]*:[ \t]*true' | tr -d "{}()[]\"" | tr ":" "," | sed 's/\s*,\s*/,/g' | tr -s " " "_" | awk -F, "{print \$2 \",\" \$4 \",\" \$7 \",\" \$11;}" | sed 's/_*iso_*[0-9]*//' | \
grep -v "PowerShot" | \
grep -v "Panasonic" | \
grep -v "DYNAX" | \
grep -v "NEX-C3" | \
grep -v "pentax_k-x" | \
grep -v "D5100" | \
awk -F, "{if (\$1 == \"Nikon_D800\" || \$1 == \"Canon_EOS_5D_Mark_II\") { print \$0; }}" \
> trim.txt
# get all:
filter="cat"
# only canon:
# filter="grep Canon"
# only canon mark 2+3
# filter="grep Mark"
# only Nikon
# filter="grep NIKON"
# get list of cameras (use exif model only):
cams=$(cat trim.txt | $filter | awk -F, "{ print \$1; }" | sort | uniq | tr "\n" " ")
# these won't work any longer due to changed input format:
# cams=$(cat trim.txt | awk -F, "{ print \$1; }" | sed 's/_iso_.*$//' | sort | uniq | tr "/" "_" | tr "\n" " ") # use first field to get commented double measurements as separate data points
# filtered:
# cams=$(cat trim.txt | $filter | awk -F, "{ print \$1; }" | sed 's/_iso_.*$//' | sort | uniq | tr "/" "_" | tr "\n" " ")
# is actually num cams + 1, because the first column is iso
num_cams=$(( $(echo $cams | wc -w) + 1))
# get sorted list of iso values:
isos=$(cat trim.txt | $filter | awk -F, "{ print \$2; }" | sort -g | uniq)
# isos="200 400 800 1600"
# we want output files to look like:
# iso cam1 cam2 cam3 cam4
# iso1 ..
# iso2 ...
#
# so we first output the header:
echo "iso ${cams}" > data.txt
echo "iso ${cams}" > data2.txt
for iso in $isos
do
echo -n "$iso " >> data.txt
echo -n "$iso " >> data2.txt
echo "collecting iso $iso .."
for cam in $cams
do
# echo "looking for cam $cam"
# collect green poissonian value for this camera and iso ($6)
a=$(cat trim.txt | awk -F, "{if ((\$1 == \"$cam\") && (\$2 == $iso)) { print \$3; } }" | head -n 1 | tr -d "\n")
# same for gaussian one
b=$(cat trim.txt | awk -F, "{if ((\$1 == \"$cam\") && (\$2 == $iso)) { print \$4; } }" | head -n 1 | tr -d "\n")
if [ "$a" = "" ]
then
# echo "no value found for $cam iso $iso"
a="?"
fi
echo -n "$a " >> data.txt
if [ "$b" = "" ]
then
# echo "no value found for $cam iso $iso"
b="?"
fi
echo -n "$b " >> data2.txt
done
echo "" >> data.txt
echo "" >> data2.txt
done
# now for some plotting pleasure:
gnuplot << EOF
set term pdf fontscale 0.5 size 10, 10
set output 'poissonian.pdf'
set mxtics 10
set grid mxtics xtics ytics
set logscale xy
# set yrange [0:1e-5]
set datafile missing "?"
set key autotitle columnhead
set xlabel 'iso speed'
set ylabel 'photon noise'
plot for [i=2:${num_cams}] './data.txt' u 1:(column(i)) w lp title column(i)
set output 'gaussian.pdf'
unset logscale
set logscale x
set yrange [-1e-6:1e-5]
set ylabel 'sensor noise'
plot for [i=2:${num_cams}] './data2.txt' u 1:(column(i)) w lp title column(i)
EOF
rm -f trim.txt data.txt data2.txt
|
