你能计算出Gnuplot中轮廓的面积吗?

时间:2015-01-27 15:19:01

标签: gnuplot contour area

我已经使用gnuplot几个星期了。我有23个变量的大型数据文件,但我选择了特定的x-y坐标数据和荧光强度数据供我分析。

我想做的事情是我的荧光颗粒的等高线图。我应该补充一点,这个轮廓图是随着时间的推移,所以会有几个点几乎重叠,但实际上这是相同的粒子。我想在这些斑点周围绘制轮廓,根据强度绘制颜色代码,并在图形上显示轮廓区域。

我已经为我的轮廓图实现了除了其中一个目标之外的所有目标。我无法为gnuplot设计一种方法来计算和显示轮廓内的区域。如果可以的话,我会估计粒子的面积。我认识到我的目标可能超出了gnuplot的能力,但是如果有解决方案那么它会非常整洁。

这是我的轮廓图的脚本,正如我所说,我需要的所有内容都是轮廓内的区域。

坐标以纳米为单位,数据集上的每个点都是分子的中心。我采取了一小部分坐标,因为有太多的数据,否则就无法区分(有超过80 000个数据点)。我还设置了一个强度阈值,因为我只想要相对明亮的荧光粒子(用set cntrparam level增量8000,5000,100000完成)。 $ 23和$ 24分别是x和y坐标。 12美元是强度。

#Contour plot of Fluorescent Particle Location with Intensity
#Gnuplot script file for plotting data in file "1002 all.txt"
reset
set dgrid3d 100,1000,1
set pm3d
set isosample 30
set xlabel 'x (nm)'
set ylabel 'y (nm)'
set contour base
set cntrparam levels incremental 8000,5000,100000
unset key
unset surface
set view map
set xrange[20000:22000]
set yrange[7000:10000]
splot "1002 all.txt" using ($23<22000 && $23>20000 ?$23 : 1/0):$24<10000 && $24>7000 ?$24 : 1/0):12 with lines
set terminal push               


set terminal png                
set output "1002_all_fluorophores_section_contour.png" # set the output filename
set terminal png size 1280,760
replot    
set output                               

0 个答案:

没有答案