在斐济(ImageJ),我打开了两张图片(Img1和Img2)。我想运行一个脚本,添加两个图像并将结果存储在Img1中。我想在一系列图像中这样做,所以我想尽量避免创建和关闭许多图像。
这可能吗?我尝试了下面的代码,但是当我在第一次Sum3and50.show()
通话后调用第二个AddSlice()
时崩溃了。基本上我很想能够到Sum3and50+=imp[Slice]
from __future__ import division
from ij import IJ
from ij import plugin
import time
def AddSlice(Stack,SumImg,Slice):
Stack.setSlice(Slice)
ic = plugin.ImageCalculator()
SliceImg = ic.run("Copy create", Stack, Stack)
SliceImg.show()
time.sleep(SLEEP_TIME)
SumImg=ic.run("Add RGB", SumImg, SliceImg)
return SumImg
SLEEP_TIME=1 #seconds
#imp = IJ.getImage()
imp = IJ.openImage("http://imagej.nih.gov/ij/images/flybrain.zip");
W,H,NCh,NSl,NFr = imp.getDimensions()
imp.show()
Sum3and50 = IJ.createImage("Sum3and50", "RGB black", W, H, 1)
Sum3and50.show()
time.sleep(SLEEP_TIME)
Sum3and50 = AddSlice(imp,Sum3and50,3)
Sum3and50.show()
time.sleep(SLEEP_TIME)
Sum3and50 = AddSlice(imp,Sum3and50,5)
Sum3and50.show()
为了避免窗口弹出,我倾向于避免使用插件并直接使用ImageProcessor
。这种取两个图像的每个像素对的总和覆盖第一个输入的函数看起来像:
def pixel_pair_sum(pro1, pro2):
for x in range(pro1.getWidth()):
for y in range(pro1.getHeight()):
v1 = pro1.get(x, y)
v2 = pro2.get(x, y)
pro1.set(x, y, v1 + v2)
pro1
和pro1
是ImageProcessor
s [1]。所以你需要在调用上面的函数之前从ImagePlus
中获取第一个:
...
sum3and50 = IJ.createImage("Sum3and50", "RGB black", W, H, 1)
p1 = sum3and50.getProcessor()
stk = imp.getStack()
p2 = stk.getProcessor(1) # get the processor for the first slice [2]
pixel_pair_sum(p1, p2) # add the pixel values of slice 1 to sum3and50
p2 = stk.getProcessor(2) # add another slice to sum3and50
pixel_pair_sum(p1, p2)
...
sum3and50.show()
供参考:https://imagej.nih.gov/ij/developer/api/ij/ImageStack.html [1] https://imagej.nih.gov/ij/developer/api/ij/process/ImageProcessor.html [2]也看到:https://imagej.nih.gov/ij/developer/api/ij/ImagePlus.html#setPositionWithoutUpdate-int-int-int-