Blame | Last modification | View Log | Download
set table 'gauss1++.txt'
splot 'gauss1++.avs' binary format='%4uint8' record=180x180 skip=8 using 2
unset table
set xlabel 'pixels'
set ylabel 'pixels'
set zlabel 'intensity' rotate by 90;
set hidden3d
unset key
set terminal png
set output 'gauss1++_3d.png'
splot 'gauss1++.txt' with lines
unset output
unset key
set samples 5000
set xlabel 'pixels'
set ylabel 'intensity'
f(x)=a*exp(-(x-x0)**2/(2*w**2))
a=100
x0=90
w=30
fit f(x) 'gauss1++.txt' every ::90::100 using 2:3 via a,x0,w
set output 'gauss1++_fit1.png'
plot f(x) w lines linewidth 3, 'gauss1++.txt' every ::90::100 using 2:3
unset output
unset key
f(x)=a*exp(-(x-x0)**2/(2*w**2))
a=200
x0=80
w=20
fit f(x) 'gauss1++.txt' every :::90::100 using 1:3 via a,x0,w
set output 'gauss1++_fit2.png'
plot f(x) w lines linewidth 3, 'gauss1++.txt' every :::90::100 using 1:3
unset output