divide_line_from_stack_jru_v1.java

/*******************************************************************************
 * Copyright (c) 2012 Jay Unruh, Stowers Institute for Medical Research.
 * All rights reserved. This program and the accompanying materials
 * are made available under the terms of the GNU Public License v2.0
 * which accompanies this distribution, and is available at
 * http://www.gnu.org/licenses/old-licenses/gpl-2.0.html
 ******************************************************************************/
import ij.*;
import ij.process.*;
import ij.gui.*;
import java.awt.*;
import ij.plugin.*;
import jguis.*;

public class divide_line_from_stack_jru_v1 implements PlugIn {

	public void run(String arg) {
		int[] wList = WindowManager.getIDList();
		String[] titles = new String[wList.length];
		for(int i=0;i<wList.length;i++){
			ImagePlus imp = WindowManager.getImage(wList[i]);
			if(imp!=null){
				titles[i]=imp.getTitle();
			} else {
				titles[i]="NA";
			}
		}
		GenericDialog gd = new GenericDialog("Options");
		gd.addChoice("Plot Window",titles,titles[0]);
		gd.addChoice("Stack",titles,titles[0]);
		gd.addCheckbox("Multiply",false);
		gd.showDialog();
		if(gd.wasCanceled()){return;}
		int index1 = gd.getNextChoiceIndex();
		int index2 = gd.getNextChoiceIndex();
		ImagePlus imp = WindowManager.getImage(wList[index1]);
		ImagePlus imp2 = WindowManager.getImage(wList[index2]);
		boolean mult=gd.getNextBoolean();

		ImageWindow iw=imp.getWindow();
		float[] yvals=((float[][])jutils.runPW4VoidMethod(iw,"getYValues"))[0];
		int length=yvals.length;
		ImageStack stack=imp2.getStack();
		int width=imp2.getWidth(); int height=imp2.getHeight();
		ImageStack retstack=new ImageStack(width,height);

		for(int i=0;i<length;i++){
			float[] pixels=((float[])stack.getPixels(i+1)).clone();
			for(int j=0;j<width*height;j++){
				if(mult) pixels[j]*=yvals[i];
				else pixels[j]/=yvals[i];
			}
			retstack.addSlice("",pixels);
		}
		String name="Divided";
		if(mult) name="Multiplied";
		(new ImagePlus(name,retstack)).show();
	}

}