基于cellpose分割结果统计细胞膜荧光强度

1
// 测量给出细胞膜大致的厚度（同scale单位，此处是微米）
2
width = 3;
3
img = getTitle();
4
wks = getDirectory("image");
5
n = roiManager("count");
6

7

8
print("id,area,mean,min,max,std");
9
for (i = 0; i < n; i++) {
10
  roiManager("select", i);
11
  getMembraneIntensity(width, i);
12
  selectWindow(""+i);
13
  fname = img+"_p_"+i+".tif";
14
  fpath = wks + "/" + fname;
15
  save(fpath);
16
  close();
17
}
18
// 保存Log窗口中的结果为csv文件并关闭
19
fname2 = img+"_results.csv";
20
fpath2 = wks + "/" + fname2;
21
selectWindow("Log");
22
saveAs("text", fpath2);
23
run("Close");
24

25

26
function getMembraneIntensity(width, idx) {
27
  // function description
28
  run("Enlarge...", "enlarge=-"+width);
29
  run("Make Band...", "band="+width+0.5);
30
  run("Duplicate...", "title="+idx);
31
  setBackgroundColor(0, 0, 0);
32
  run("Clear Outside");
33
  run("Add Selection...");
34
  setAutoThreshold();
35
  run("Create Selection");
36
  getStatistics(area, mean, min, max, std, histogram);
37
  print(""+idx+","+area+","+mean+","+min+","+max+","+std);
38
}