0,0 → 1,38 |
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 |
|