仅在ROI中获取像素,如何?

时间:2015-04-23 17:52:12

标签: python jython image-segmentation imagej roi

抱歉,我是Image Analysis游戏的新手。尝试编写脚本,到目前为止,我有一个返回通过阈值处理得出的ROI的函数。

PA.setRoiManager(roim)
paOpt = PA.CLEAR_WORKSHEET +PA.SHOW_NONE +PA.INCLUDE_HOLES +PA.EXCLUDE_EDGE_PARTICLES + PA.ADD_TO_MANAGER 
# measurement options: clear the results, show nothing, include holes, exclude particles on edges and add ROIs to manager
measOpt = PA.AREA + PA.MEAN + PA.MIN_MAX # measure the minimum, maximum, mean and area of particles          
pa = PA(paOpt, measOpt, rt, minSize, maxSize)
pa.analyze(imp)
roi_list = roim.getRoisAsArray()
print roi_list[0] 

这给控制台一个输出:Roi [Traced,x =,y =,width =,height =] 所以我相信我找到了投资回报率。 roim是RoiManager的一个实例。

现在问题是当我去查看ROI中的像素时。

def pixelStats(roiPatch, imp):
  '''
  get the pixel value distribution for the patch ROI 
  return the number of Lo pixels and the number of Ld pixels 
  roi_Patch returns handle to the patch ROI and imp is the 
  instance of current ImagePlus to be measured
  '''
  mask = roiPatch.getMask() # mask of roiPatch, do you need this?
  ip = imp.getProcessor() 
  ip.setRoi(roiPatch)
  cal = imp.getCalibration()
  # options = IS.MEAN | IS.STD_DEV | IS.MODE
  stats = IS.getStatistics(ip, M.MEAN | M.MODE | M.STD_DEV, cal) # what does getCalibration do??? stats for ROI only applied to imp, why?
  # see the following --> http://fiji.sc/Jython_Scripting#Obtain.2FView_histogram_and_measurements_from_an_image
  domThreshold = stats.mode - stats.STD_DEV 
  backGround = stats.MEAN - stats.STD_DEV 
  # define the threshold for Lo/Ld phase 
  pixels = ip.getPixels() 
  above = filter(lambda i: pixels[i] > backGround, xrange(len(pixels)))
  below = filter(lambda i: above[i] < domThreshold, xrange(len(above)))
  # calculate the length of the arrays containing pixels of Lo/Ld phase 
  return len(above), len(below), len(pixels)

查看返回的测量结果,我确信ip.getPixels()只返回原始图像处理器的像素而不是ROI。换句话说,我只是得到了原始&#34;矩形&#34;中的像素数。图片。

任何人都可以提供原因吗?或者建议一个更好的方法来做同样的事情。

许多人非常感谢!你可以从我的评论中看到我有多少菜鸟。

1 个答案:

答案 0 :(得分:-1)

我熟悉Java API但不熟悉python,因此我的答案尚未经过测试。

方法

ip.getPixels() 

正在做我期望的并返回所有像素。是否有投资回报率不会影响这一点。您的测量结果考虑了ROI,因此平均值和模式是根据您设置的ROI计算的。

要解决此问题,我会裁剪ip。应该有一个

ip.crop() 

方法会返回ip考虑了投资回报率,然后您可以调用

ip.getPixels()

仅获取ROI中的像素。