// This Macros are made for PETIS Image Analysis.
// Created by Naoki Kawachi.
// Version1.0.3(b) 17/Jan/2007


macro "Import ALL PETIS Images of text format" {
// This Macro Make a stack image from the All PETIS Text Images in a Directory.
// Select the directory and check the image property.

	// Set Target Text Images Directory.
	dir=getDirectory("");
	dirlist = getFileList(dir);
	FILE_list = lengthOf(dirlist);

	j=0;
	for (i=0; i<FILE_list; i++){
		if ((substring(dirlist[i],0,1)=="I"||substring(dirlist[i],0,1)=="H"||substring(dirlist[i],0,1)=="G") && (substring(dirlist[i],lengthOf(dirlist[i])-1,lengthOf(dirlist[i]))=="/")){
			j=j+1;
			}
	}
	PETIS_list=j;
	PETIS_Images=newArray(PETIS_list);

	j=0;
	for (i=0; i<FILE_list; i++){
		if ((substring(dirlist[i],0,1)=="I"||substring(dirlist[i],0,1)=="H"||substring(dirlist[i],0,1)=="G") && (substring(dirlist[i],lengthOf(dirlist[i])-1,lengthOf(dirlist[i]))=="/")){
			//print("i="+i+", j="+j);
			PETIS_Images[j]=dirlist[i];
			j=j+1;
		}
	}
	
	if (PETIS_list==0){
		Dialog.create("Read Error");
		Dialog.addMessage("No PETIS Image Directory");
		Dialog.show();
		exit();
	}
	
	
	signal = newArray(PETIS_list);	
	Dialog.create("Open Text Images as one Stacks");
	Dialog.addMessage("Select Import PETIS_Image from TEXT format directories");
	for (i=0; i<PETIS_list; i++){
		Dialog.addCheckbox(PETIS_Images[i], true);
	}
	Dialog.addMessage(" ");
	Dialog.addMessage("Check Your PC Memory!!");
	Dialog.addCheckbox("Close the Windows of Created Images", true);
	Dialog.show();
	for (i=0; i<PETIS_list; i++){
		signal[i]=Dialog.getCheckbox();
	}

	check_close=Dialog.getCheckbox();
	
	dirs=newArray(PETIS_list);
	for (i=0; i<PETIS_list; i++){
		if (signal[i]==1){
			dirs[i]=dir+PETIS_Images[i];
			list=getFileList(dirs[i]);
			fileNameHead=substring(list[0],0,lengthOf(list[0])-5);
			headerFileName=dirs[i]+fileNameHead+0+".txt";
			period=getTagPETIS(headerFileName,"Measurement period");
			period = 0 + substring(period, 0, lengthOf(period)-4); // convert to number
			name=getTagPETIS(headerFileName,"File name");
			name=substring(name,1,lengthOf(name));
			comment = getTagPETIS(headerFileName,"Comment");
			comment ="PETIS Comment:\t"+comment;
			PETIS_date = getTagPETIS(headerFileName,"Measurement date");
			PETIS_start = getTagPETIS(headerFileName,"Measurement start time");
			PETIS_date = "PETIS Measurement Date:\t"+PETIS_date+PETIS_start;
			PETIS_RI = getTagPETIS(headerFileName,"Name of RI");
			PETIS_RI ="PETIS Positron RI:\t"+PETIS_RI;
			blank_name = getTagPETIS(headerFileName,"Blank Data Name");
			blank_region_a = getTagPETIS(headerFileName,"Blank Data Region (First)");
			blank_region_b = getTagPETIS(headerFileName,"Blank Data Region (End)");
			PETIS_Blank ="PETIS Blank:\t"+blank_name+" ("+blank_region_a+" -"+blank_region_b+" )";
			
			depth=lengthOf(list)-1;
			fileName=dirs[i]+fileNameHead+1+".txt";
			lines = File.openAsString(fileName);
			line = split(lines,"\n");
			height = line.length;
			counts = split(line[0]," ");
			width = counts.length;
			start = 1 + 0;
			end = depth;
			frame = end - start +1;
			
			// Make New Image for New Stack Image
			tmpname="stacking_"+name;
			newImage(tmpname, "32-bit Black", width, height, frame);
			run("Properties...", "width="+width+" height="+height+" channels=1 depth=1 frames="+frame+" unit=mm pixel_width=1.10000 pixel_height=1.10000 voxel_depth=1.00000 frame="+period+" origin=0,0");
			PETIS_header = "Info:"+comment+"\n"+PETIS_date+"\n"+PETIS_RI+"\n"+PETIS_Blank;
			setMetadata(PETIS_header);
			
			// Set Each Pixel Counts from ASCII Files
			for (h=start; h<=end; h++){ //Read File (Z)
				fileName=dirs[i]+fileNameHead+h+".txt";
				lines = File.openAsString(fileName);
				line = split(lines,"\n"); // Read Lines (Y)
				setSlice(h-start+1);
				for (j=0; j<line.length; j++){
					counts = split(line[j], " "); // Read Lines (X)
					for (k=0; k<counts.length; k++){
						setPixel(k,j,counts[k]);
					}
				}
			}
			rename(name);
			
			// Set New Image Contrast and Color.
			setslice = floor(frame/2)+1;
			setSlice(setslice);
			run("Grays");
			run("Brightness/Contrast...");
			run("Enhance Contrast", "saturated=0.5");
			setSlice(1);
			selectWindow("B&C");
			run("Close");
			
			selectWindow(name);
			print(dir+name+".tif");
			save(dir+name+".tif");
			if (check_close==1){close();}
		}			
	}
}

macro "Import PETIS Text Images 	[f1]" {
// This Macro Make a stack image from the PETIS text images.
// Select the directory and check the image property.
	//1. Set Target Text Images Directory.
	dir=getDirectory("");
	list = getFileList(dir);
	fileNameHead=substring(list[0],0,lengthOf(list[0])-5);
	
	//2. Read Image info from PETIS Image header file of "xx_0.txt".
	headerFileName=dir+fileNameHead+0+".txt";
	period=getTagPETIS(headerFileName,"Measurement period");
	period = 0 + substring(period, 0, lengthOf(period)-4); // convert to number
	name=getTagPETIS(headerFileName,"File name");
	name=substring(name,1,lengthOf(name));
	comment = getTagPETIS(headerFileName,"Comment");
	comment ="PETIS Comment:\t"+comment;
	PETIS_date = getTagPETIS(headerFileName,"Measurement date");
	PETIS_start = getTagPETIS(headerFileName,"Measurement start time");
	PETIS_date = "PETIS Measurement Date:\t"+PETIS_date+PETIS_start;
	PETIS_RI = getTagPETIS(headerFileName,"Name of RI");
	PETIS_RI ="PETIS Positron RI:\t"+PETIS_RI;
	blank_name = getTagPETIS(headerFileName,"Blank Data Name");
	blank_region_a = getTagPETIS(headerFileName,"Blank Data Region (First)");
	blank_region_b = getTagPETIS(headerFileName,"Blank Data Region (End)");
	PETIS_Blank ="PETIS Blank:\t"+blank_name+" ("+blank_region_a+" -"+blank_region_b+" )";
	
	//3. Read Image width, Heights and File Name from PETIS Image "xx_1.txt".
	depth=lengthOf(list)-1;
	fileName=dir+fileNameHead+1+".txt";
	lines = File.openAsString(fileName);
	line = split(lines,"\n");
	height = line.length;
	counts = split(line[0]," ");
	width = counts.length;
	
	//4. Reset Image info.
	Dialog.create("Open Text Images as one Stack");
	Dialog.addString("Image Name", name);
	Dialog.addNumber("width", width);
	Dialog.addNumber("height", height);
	Dialog.addNumber("frame start", 1);
	Dialog.addNumber("frame end", depth);
	Dialog.addNumber("frame period (sec)", period);
	Dialog.show();
	name = Dialog.getString();
	width = Dialog.getNumber();
	height = Dialog.getNumber();
	start =  Dialog.getNumber();
	end = Dialog.getNumber();
	period =  Dialog.getNumber();
	frame = end - start +1;
	
	//5. Make New Image for New Stack Image
	tmpname="stacking_"+name;
	newImage(tmpname, "32-bit Black", width, height, frame);
	run("Properties...", "width="+width+" height="+height+" channels=1 depth=1 frames="+frame+" unit=mm pixel_width=1.10000 pixel_height=1.10000 voxel_depth=1.00000 frame="+period+" origin=0,0");
	PETIS_header = "Info:"+comment+"\n"+PETIS_date+"\n"+PETIS_RI+"\n"+PETIS_Blank;
	setMetadata(PETIS_header);
	
	//6. Set Each Pixel Counts from ASCII Files
	for (i=start; i<=end; i++){ //Read File (Z)
		fileName=dir+fileNameHead+i+".txt";
		lines = File.openAsString(fileName);
		line = split(lines,"\n"); // Read Lines (Y)
		setSlice(i-start+1);
		for (j=0; j<line.length; j++){
			counts = split(line[j], " "); // Read Lines (X)
			for (k=0; k<counts.length; k++){
				//print("x="+k+" y="+j+" z="+i);
				setPixel(k,j,counts[k]);
			}
		}
	}
	rename(name);
	
	//7. Set New Image Contrast and Color.
	setslice = floor(frame/2)+1;
	setSlice(setslice);
	run("Grays");
	run("Brightness/Contrast...");
	run("Enhance Contrast", "saturated=0.5");
	setSlice(1);
	selectWindow("B&C");
	run("Close");

}//EOMacro


macro "Set PETIS Image Frames		[f2]" {

	name=getTitle();
	width=getWidth();
	height=getHeight();
	frame=nSlices;
	
	Dialog.create("PETIS Frame Slecter");
	Dialog.addNumber("Select Start Frame: ",1);
	Dialog.addNumber("Select End Frame: ",frame);
	Dialog.show();
	start = Dialog.getNumber();
	end = Dialog.getNumber();
	
	
	selectWindow(name);
	
	for (i=frame;i>end;i--){
	setSlice(i);
	run("Delete Slice");
	}
	
	for (i=start-1;i>0;i--){
	setSlice(i);
	run("Delete Slice");
	}

}





macro "Re-Stack PETIS Image 		[f3]" {
// This Macro Make a New Stack Image.
// Choose a Target Image and Restacking.
	////1. Select Target Stack Image//////
	if (nImages>1){
		title = getTitle();
		Dialog.create("Select Target Image");
		Dialog.addMessage("Target Images is "+title);
		Dialog.show();
		selectWindow(title);
	}
	else if (nImages==1){
		title = getTitle();
		selectWindow(title);
	}
	else if (nImages==0){
		open("");
	}
	
	////2. Read Image Informations//////
	name=getTitle();
	name_head=split(name, ".");
	name_head=name_head[0];
	width=getWidth();
	height=getHeight();
	frame=nSlices;
	step = floor(frame/10);
	
	////3. Set New Stack Image Parameters//////
	selectWindow(name);
	Dialog.create("Set New Frame Info");
	Dialog.addNumber("Enter Start Frame No.", 1);
	Dialog.addNumber("How Many Steps?",step);
	Dialog.addString("Restack Image File Name", name_head+"_restack");
	Dialog.show();
	start = Dialog.getNumber();
	step = Dialog.getNumber();
	new_name = Dialog.getString();
	new_frame = floor((frame-start+1) / step);
	tag = "Frame interval";
	period = getTag(tag);
	period = split(period," ");
	period = (0+period[0])*step;
	
	tag="PETIS Comment";
	comment="PETIS Comment:"+getTag(tag);
	tag="PETIS Measurement Date";
	PETIS_date="PETIS Measurement Date:"+getTag(tag);
	tag="PETIS Positron RI";
	PETIS_RI="PETIS Positron RI:"+getTag(tag);
	tag="PETIS Blank";
	PETIS_Blank="PETIS Blank:"+getTag(tag);
	
	////4. Make New Stack Images //////
	newImage(new_name,"32-bit",width,height,new_frame);
	run("Properties...", "width="+width+" height="+height+" channels=1 depth=1 frames="+new_frame+" unit=mm pixel_width=1.10000 pixel_height=1.10000 voxel_depth=1.00000 frame="+period+" origin=0,0");
	PETIS_header = "Info:"+comment+"\n"+PETIS_date+"\n"+PETIS_RI+"\n"+PETIS_Blank;
	setMetadata(PETIS_header);
	
	for (l=1;l<=new_frame;l++){
		selectWindow(name);
		pixel_count=newArray(width*height);
		for (k=start;k<=start+step-1;k++){
			setSlice(k+step*(l-1));
			//print("Slice="+k+step*(l-1));
			for (j=0;j<height;j++){
				for (i=0;i<width;i++){
					//print(i+j*width+","+"("+i+","+j+")");
					pixel_count[i+j*width]=pixel_count[i+j*width]+getPixel(i,j)/step;
					}
				}
		}
	
		selectWindow(new_name);
		setSlice(l);
		for (j=0;j<height;j++){
			for (i=0;i<width;i++){
				setPixel(i,j,pixel_count[i+j*width]);
				}
		}	
	}
	//5. Set New Image Contrast and Color.
	setslice = floor(new_frame/2)+1;
	setSlice(setslice);
	run("Grays");
	run("Brightness/Contrast...");
	run("Enhance Contrast", "saturated=0.5");
	setSlice(1);
	selectWindow("B&C");
	run("Close");	
		
		
}//End of Macro


macro "Get PETIS TAC 				[f4]" {
// This Macro Generates a Time Acitvity Curve from PETIS Image(stack).
// After a ROI selected, start up this Macro.
	//1. Get Header Information from PETIS Image
	name=getTitle();
	dir=getDirectory("image");
	if (dir == "") dir = "This Image on the ImageJ Memory";
	tag="Frame interval";
	str=getTag(tag);
	period=split(str," ");
	period=0+period[0];
	if (period==0) period = 1;
	tag="PETIS Comment";
	comment=getTag(tag);
	tag="PETIS Measurement Date";
	mdate=getTag(tag);
	tag="PETIS Positron RI";
	positronRI=getTag(tag);
	n = nSlices;
	roiManager("Add");
	
	//2. Get Image Spec
	area = newArray(n);
	mean = newArray(n);
	min = newArray(n);
	max = newArray(n);
	time = newArray(n);
	info = getImageInfo();
	info=split(info,"\n");
	roi_shape=info[lengthOf(info)-5];
	roi_x=info[lengthOf(info)-4];
	roi_y=info[lengthOf(info)-3];
	roi_width=info[lengthOf(info)-2];
	roi_height=info[lengthOf(info)-1];
	a=split(roi_x," ");
	roi_x=kakko_hiki(a[lengthOf(a)-1]);
	a=split(roi_y," ");
	roi_y=kakko_hiki(a[lengthOf(a)-1]);	
	a=split(roi_width," ");
	roi_width=kakko_hiki(a[lengthOf(a)-1]);
	a=split(roi_height," ");
	roi_height=kakko_hiki(a[lengthOf(a)-1]);
	
	//3. Generate Time Activity Curve
	run("Clear Results");
	run("Plot Z-axis Profile");
	close();
	for (i=0; i<n; i++){
		mean[i] = getResult("Mean",i);
		area[i] = getResult("Area",i);
		time[i]=i*period;
	}
	
	//4. Set TAC Graph
	meanMin=getMin(mean);
	meanMax=getMax(mean);
	Plot.create("Time Activity Curve ("+name+")", "Time", "counts(mean)", time, mean);
	Plot.setLimits(time[0], time[n-1], 0, meanMax);
	Plot.show();
	Dialog.create("Time Activity Time Save Dialog");
	Dialog.addChoice("OK?",newArray("Save TAC", "Retry"),"Retry");
	Dialog.addCheckbox("Hold ROIs in ROI Manager", true)
	Dialog.show();
	ans = Dialog.getChoice();
	roi_save=Dialog.getCheckbox();
	if (roi_save==0) roiManager("reset");
	if(ans == "Retry"){
		selectWindow("Time Activity Curve ("+name+")");
		close();
		selectWindow("Results");
		run("Close");
		exit();
	}
	
	//5. SAVE TAC Results as a .txt file (cvs format)
	ijpath=getDirectory("startup");
	roipath=ijpath+"roi_files";
	exist=File.exists(roipath);
	if (exist==0){
		Dialog.create("Title");
		Dialog.addMessage("ROI directory will be created: "+roipath);
		Dialog.show();
		File.makeDirectory(roipath);
	}
	selectWindow(name);	
	roiSaveFileName=getTime1();
	roilist_ln=roiManager("count");
	//roiManager("Select",roilist_ln-1);
	roiManager("Save",roipath+"/"+roiSaveFileName+".zip");
	//roiManager("reset");
	
	f = File.open("");
	print(f,"##\tImage File Name = "+name);
	print(f,"##\tImage File Directory = "+dir);
	print(f,"##\tSaved ROI File = "+roipath+"/"+roiSaveFileName+".zip");
	print(f,"##\tSaved ROI Shape = "+roi_shape);
	print(f,"##\t\tX= "+roi_x+" Y= "+roi_y+" Width= "+roi_width+" Height= "+roi_height+" Area= "+area[0]);
	print(f,"##\tComment ="+comment);
	print(f,"##\tMeasurement date ="+mdate);
	print(f,"##\tPositron RI ="+positronRI);
	print(f,"##\tframe\tstart_time(sec)\tend_time(sec)\tmean_count\ttotal_count");
	for (i=0; i<n; i++){
		print(f,"\t"+i+1+" \t"+time[i]+" \t"+time[i]+period+" \t"+mean[i]+" \t"+mean[i]*area[i]);
	}
	File.close(f);

	//6. Macro End
	selectWindow("Time Activity Curve ("+name+")");
	close();
	selectWindow("Results");
	run("Close");
	exit();

}//End of Macro

macro "PETIS Image Multiplier 		[f5]"{
// This Macro Multiple PETIS Image (each Frames).
	name=getTitle();
	width=getWidth();
	height=getHeight();
	frame=nSlices;
	
	n=nImages;
	image_name=newArray(n);
	for (i=0; i<n; i++){
		image_name[i]=getTitle();
		run("Put Behind [tab]");
	}
	
	Dialog.create("PETIS Image Multiplier");
	Dialog.addChoice("Target Image: ", image_name, image_name[0]);
	Dialog.addNumber("Select Start Frame: ",1);
	Dialog.addNumber("Select End Frame: ",frame);
	Dialog.addMessage(" ");
	Dialog.addNumber("Multiple Number: ",1.0);
	Dialog.show();
	name=Dialog.getChoice();
	start = Dialog.getNumber();
	end = Dialog.getNumber();
	multi = Dialog.getNumber();
	selectWindow(name);
	
	for (k=start;k<=end;k++){
		setSlice(k);
		pixel_count=newArray(width*height);
		for (j=0;j<height;j++){
			for (i=0;i<width;i++){
				pixel_count[i+j*width]=getPixel(i,j)*multi;
				setPixel(i,j,pixel_count[i+j*width]);
			}
		}
	}
	setSlice(1);
}

macro "--Presantation Utilities--" {}

macro "  PETIS Color Bar 			[f6]"{
// This Macro Add Color Bar at Leftside of Stack Image.
	name=getTitle();
	width=getWidth();
	height=getHeight();
	frame=nSlices;
	getMinAndMax(min,max);
	
	Dialog.create("Color Bar");
	Dialog.addChoice("Select",newArray("Add Color Bar", "Reset Color Bar", "Clear Color Bar"),"Add Color Bar");
	Dialog.show();
	ans = Dialog.getChoice();
	if(ans == "Add Color Bar"){
		tag = "Display range";
		str=getTag(tag);
		range=split(str," ");
		lower=0+range[0];
		upper=0+range[2];
		//print("("+lower+","+upper+")");
		new_width=width+10;
		run("Canvas Size...", "width="+new_width+" height="+height+" position=Center-Right zero");
		for (k=1;k<=frame;k++){
			setSlice(k);
			pixel_count=newArray(height);
			for (j=0;j<height;j++){
				pixel_count[j]=upper-(upper-lower)*j/(height-1);
				for (i=1;i<10;i++){
					setPixel(i,j,pixel_count[j]);
				}
			}
		}
	}
	else if (ans == "Reset Color Bar"){
		tag = "Display range";
		str=getTag(tag);
		range=split(str," ");
		lower=0+range[0];
		upper=0+range[2];
		for (k=1;k<=frame;k++){
			setSlice(k);
			pixel_count=newArray(height);
			for (j=0;j<height;j++){
				pixel_count[j]=upper-(upper-lower)*j/(height-1);
				//print(pixel_count[j]);
				for (i=1;i<10;i++){
					setPixel(i,j,pixel_count[j]);
				}
			}
		}
	}
	else{
		width=getWidth();
		new_width = width-10;
		run("Canvas Size...", "width="+new_width+" height="+height+" position=Center-Right zero");
	}
}//End of Macro



macro "  PETIS Time Bar 			[f7]"{
// This Macro Add Time Bar at Bottom of Stack Image.
	name=getTitle();
	width=getWidth();
	height=getHeight();
	frame=nSlices;
	getMinAndMax(min,max);
	new_height=height+10;
	
	Dialog.create("Time Bar");
	Dialog.addChoice("Select Time Bar Color",newArray("Gray Scale with White", "256RGB with Red", "256RGBw with White","256RGB with LightBlue","Clear Time Bar"),"Gray Scale with White");
	Dialog.addChoice("How Many Partitions",newArray("1","2", "3", "4","6","8"),"4");
	Dialog.show();
	ans = Dialog.getChoice();
	part = Dialog.getChoice();
	if (ans == "Clear Time Bar") {
		run("Canvas Size...", "width="+width+" height="+height-10+" position=Top-Center zero");
		exit();
	}
	if (ans == "256RGB with LightBlue") max=max*0.4;
	
	start=0; end=width;
	half=round(width/2);
	oneThird=round(width/3); twoThird=round(2*width/3);
	firstQ=round(width/4); thirdQ=round(3*width/4);
	one6=round(width/6);five6=round(5*width/6);
	one8=round(width/8);three8=round(3*width/8);five8=round(5*width/8);seven8=round(7*width/8);
	
	run("Canvas Size...", "width="+width+" height="+new_height+" position=Top-Center zero");
	
	for (k=1;k<=frame;k++){
		setSlice(k);
		pixel_count=newArray(width);
		for (i=0;i<width;i++){
			pixel_count[i]=max;
			setPixel(i,height+1,pixel_count[i]);
		}
		for (j=height+2;j<=new_height;j++){
			pixel_count=newArray(width);
			if (part == "1"){
			for (i=0;i<width;i++){
				if( i==start ) pixel_count[i]=max;
				if( i==end-1 ) pixel_count[i]=max;
				setPixel(i,j,pixel_count[i]);
				}
			}
			if (part == "2"){
			for (i=0;i<width;i++){
				if( i==start ) pixel_count[i]=max;
				if( i==half ) pixel_count[i]=max;
				if( i==end-1 ) pixel_count[i]=max;
				setPixel(i,j,pixel_count[i]);
				}
			}
			if (part == "3"){
			for (i=0;i<width;i++){
				if( i==start ) pixel_count[i]=max;
				if( i==oneThird ) pixel_count[i]=max;
				if( i==twoThird ) pixel_count[i]=max;
				if( i==end-1 ) pixel_count[i]=max;
				setPixel(i,j,pixel_count[i]);
				}
			}
			if (part == "4"){
			for (i=0;i<width;i++){
				if( i==start ) pixel_count[i]=max;
				if( i==firstQ ) pixel_count[i]=max;
				if( i==half ) pixel_count[i]=max;
				if( i==thirdQ ) pixel_count[i]=max;
				if( i==end-1 ) pixel_count[i]=max;
				setPixel(i,j,pixel_count[i]);
				}
			}
			if (part == "6"){
			for (i=0;i<width;i++){
				if( i==start ) pixel_count[i]=max;
				if( i==one6 ) pixel_count[i]=max;
				if( i==oneThird ) pixel_count[i]=max;
				if( i==half ) pixel_count[i]=max;
				if( i==twoThird ) pixel_count[i]=max;
				if( i==five6 ) pixel_count[i]=max;
				if( i==end-1 ) pixel_count[i]=max;
				setPixel(i,j,pixel_count[i]);
				}
			}
			if (part == "8"){
			for (i=0;i<width;i++){
				if( i==start ) pixel_count[i]=max;
				if( i==one8 ) pixel_count[i]=max;
				if( i==firstQ ) pixel_count[i]=max;
				if( i==three8 ) pixel_count[i]=max;
				if( i==half ) pixel_count[i]=max;
				if( i==five8 ) pixel_count[i]=max;
				if( i==thirdQ ) pixel_count[i]=max;
				if( i==seven8) pixel_count[i]=max;
				if( i==end-1 ) pixel_count[i]=max;
				setPixel(i,j,pixel_count[i]);
				}
			}
		}
		pixel_count=newArray(width);
		for (j=height+3;j<height+8;j++){
			for (i=0;i<width;i++){
				if(i<k/frame*width) pixel_count[i]=max;
				setPixel(i,j,pixel_count[i]);
			}
		}
	}
	if (ans == "Gray Scale with White") run("Grays");
	if (ans == "256RGB with Red") G256_RGBb_LUT();
	if (ans == "256RGB with LightBlue") G256_RGBb_LUT();
	if (ans == "256RGBw with White") G256_RGBw_LUT();
}//End of Macro

macro "  PETIS Lap Time Stamp 		[f8]"{
// This Macro Add Laptime on the Stack Image.
	name=getTitle();
	width=getWidth();
	height=getHeight();
	frame=nSlices;
	tag="Frame interval";
	str=getTag(tag);
	period=split(str," ");
	period=0+period[0];
	
	getMinAndMax(min,max);
	new_height=height+10;
	
	Dialog.create("Lap Time Stamp");
	Dialog.addChoice("Lap Time Location:",newArray("Top Left", "Top Right", "Bottom Left","Bottom Right","Location Set"),"Top Left");
	Dialog.addChoice("Select Time Unit:",newArray("hour", "min", "sec"),"min");
	Dialog.show();
	location = Dialog.getChoice();
	unit = Dialog.getChoice();
	if (location == "Top Left") {lo_x= width*0.1; lo_y= height*0.1;}
	else if (location == "Top Right") {lo_x = width-width*0.5; lo_y= height*0.1;}
	else if (location == "Bottom Left") {lo_x = width*0.1; lo_y= height-height*0.1;}
	else if (location == "Bottom Right") {lo_x = width-width*0.5; lo_y= height-height*0.1;}
	else {
		Dialog.create("Set Lap Time Indication");
		Dialog.addNumber("x-locate", 0);
		Dialog.addNumber("x-locate", 0);
		Dialog.show();
		lo_x=Dialog.getNumber();
		lo_y=Dialog.getNumber();
	}
	
	setColor(255, 255, 255);
	for (k=1;k<=frame;k++){
		setSlice(k);
		if (unit == "sec") {
			lap_time=k*period;
			drawString(lap_time+"(sec)", lo_x, lo_y);
		}
		if (unit == "min") {
			lap_time=k*period/60;
			lap_time=d2s(lap_time, 1);
			drawString(lap_time+"(min)", lo_x, lo_y);
		}
		if (unit == "hour") {
			lap_time=k*period/60/60;
			lap_time=d2s(lap_time, 1);
			drawString(lap_time+"(hour)", lo_x, lo_y);
		}
	}
}//End of Macro




macro "  Make Montage 			[f9]" {
// This Macro Make a montage from a PETIS stack image.
// Select the directory and check the image property.
	////1. Select Target Stack Image//////
	if (nImages>1){
		title = getTitle();
		Dialog.create("Select Target Image");
		Dialog.addMessage("Target Images is "+title);
		Dialog.show();
		selectWindow(title);
	}
	else if (nImages==1){
		title = getTitle();
		selectWindow(title);
	}
	else if (nImages==0){
		open("");
	}
	
	////2. Read Image Informations//////
	depth=nSlices;
	name = getTitle();
	if(lengthOf(name)>4) title = substring(name,0,lengthOf(name)-4);
	else title = name;
	tag="PETIS Comment";
	comment="PETIS Comment:"+getTag(tag);
	tag="PETIS Measurement Date";
	PETIS_date="PETIS Measurement Date:"+getTag(tag);
	tag="PETIS Positron RI";
	PETIS_RI="PETIS Positron RI:"+getTag(tag);
	tag="PETIS Blank";
	PETIS_Blank="PETIS Blank:"+getTag(tag);
	tag="Frame interval";
	str=getTag(tag);
	PETIS_Period="PETIS Stack Image Frame Interval:"+str;
	
	
	////3. Make Montage Image//////
	rows = floor(depth / 5);
	Dialog.create("Make Montage with PETIS Images");
	Dialog.addMessage("This Stack Image has "+depth+" frames.");
	Dialog.addNumber("columns", 5);
	Dialog.addNumber("rows", rows);
	Dialog.addNumber("first frame", 1);
	Dialog.addNumber("last frame", depth);
	Dialog.addNumber("border", 2);
	Dialog.show();
	columns = Dialog.getNumber();
	rows = Dialog.getNumber();
	first = Dialog.getNumber();
	last = Dialog.getNumber();
	boader = Dialog.getNumber();
	run("Make Montage...", "columns="+columns+" rows="+rows+" scale=1 first="+first+" last="+last+" increment=1 border="+boader+"");
	PETIS_header = "Info:"+comment+"\n"+PETIS_date+"\n"+PETIS_RI+"\n"+PETIS_Blank+"\n"+PETIS_Period;
	setMetadata(PETIS_header);
	
}//EOMacro


macro "--Look up Tables for Image Color--" {}

macro "  256G RGB" {
	reds = newArray(256); 
	greens = newArray(256); 
	blues = newArray(256);
	G256_RGB_LUT(reds,greens,blues);
	setLut(reds, greens, blues);
}

macro "  256G RGBw" {
	reds = newArray(256); 
	greens = newArray(256); 
	blues = newArray(256);
	G256_RGB_LUT(reds,greens,blues);
	reds[255] = 255; greens[255] = 255; blues[255] = 255;
	setLut(reds, greens, blues);
}

macro "  256G NIH" {
	reds = newArray(256); 
	greens = newArray(256); 
	blues = newArray(256);
	G256_NIH_LUT(reds,greens,blues);
	setLut(reds, greens, blues);
}

macro "  256G NIHw" {
	reds = newArray(256); 
	greens = newArray(256); 
	blues = newArray(256);
	G256_NIH_LUT(reds,greens,blues);
	reds[255] = 255; greens[255] = 255; blues[255] = 255;
	setLut(reds, greens, blues);
}

macro "  256G Red" {
    reds = newArray(256); greens = newArray(256); blues = newArray(256);
    for (i=0; i<256; i++) reds[i] = i;
    setLut(reds, greens, blues);
}

macro "  256G Green" {
    reds = newArray(256); greens = newArray(256); blues = newArray(256);
    for (i=0; i<256; i++) greens[i] = i;
    setLut(reds, greens, blues);
}

macro "  256G Blue" {
    reds = newArray(256); greens = newArray(256); blues = newArray(256);
    for (i=0; i<256; i++) blues[i] = i;
    setLut(reds, greens, blues);
}

macro "  Invert LUT" {
    requires("1.30j");
    getLut(reds, greens, blues);
    for (i=0; i<reds.length; i++) {reds[i] = 255-reds[i]; greens[i] = 255-greens[i]; blues[i] = 255-blues[i];}
    setLut(reds, greens, blues);
}

macro ' ' {}
macro 'Open text files (Old ver.)' {
   dir = getDirectory("Choose a Directory");
   list = getFileList(dir);
   n=getNumber("1, 2 or 3 ?", 1);
   path = dir + list[n];
   path = substring(path, 0, lengthOf(path)-5);
//   print(path);

   n=getNumber("Number of Images", 360);
   for (i=1; i<=10; i++){
   run("Text Image... ", "open=" + path+i+".txt");
  // print(i);
   }

   run("Convert Images to Stack");
   run("Set Slice...", "slice=10");

   for (i=11; i<=n; i++){
   run("Text Image... ", "open=" + path+i+".txt");
   run("Select All");
   run("Copy");
   run("Close");
   run("Add Slice");
   run("Paste");

   }
}


/////////////////MACRO FUNCTIONS///////////////////////

//Color Lookup Table Functions/////////////////////////////
function G256_RGB_LUT(reds,greens,blues) {
	for (i=0; i<52; i++) {reds[i] = 0; greens[i] = 0; blues[i] = 5*i;}
	for (i=52; i<103; i++) {reds[i] = 0; greens[i] = 5*(i-51); blues[i] = 255;}
	for (i=103; i<154; i++) {reds[i] = 0 ;greens[i] = 255; blues[i] = 255-5*(i-102);}
	for (i=154; i<205; i++) {reds[i] = 5*(i-153); greens[i] = 255; blues[i] = 0;}
	for (i=205; i<256; i++) {reds[i] = 255; greens[i] = 255-5*(i-204); blues[i] = 0;}
}
function G256_NIH_LUT(reds,greens,blues) {
	for (i= 0 ; i< 1 ;i++) { reds[i]= 0.0 ; greens[i]= 0.0 ; blues[i]= 0.0 ;}
	for (i= 1 ; i< 2 ;i++) { reds[i]= 5.1 ; greens[i]= 0.0 ; blues[i]= 5.1 ;}
	for (i= 2 ; i< 3 ;i++) { reds[i]= 9.945 ; greens[i]= 0.0 ; blues[i]= 9.945 ;}
	for (i= 3 ; i< 4 ;i++) { reds[i]= 15.045 ; greens[i]= 0.0 ; blues[i]= 15.045 ;}
	for (i= 4 ; i< 5 ;i++) { reds[i]= 20.91 ; greens[i]= 0.0 ; blues[i]= 19.89 ;}
	for (i= 5 ; i< 6 ;i++) { reds[i]= 26.01 ; greens[i]= 0.0 ; blues[i]= 24.99 ;}
	for (i= 6 ; i< 7 ;i++) { reds[i]= 31.11 ; greens[i]= 0.0 ; blues[i]= 30.09 ;}
	for (i= 7 ; i< 8 ;i++) { reds[i]= 35.955 ; greens[i]= 0.0 ; blues[i]= 34.935 ;}
	for (i= 8 ; i< 9 ;i++) { reds[i]= 42.075 ; greens[i]= 0.0 ; blues[i]= 41.055 ;}
	for (i= 9 ; i< 10 ;i++) { reds[i]= 46.92 ; greens[i]= 0.0 ; blues[i]= 45.9 ;}
	for (i= 10 ; i< 11 ;i++) { reds[i]= 52.02 ; greens[i]= 0.0 ; blues[i]= 51.0 ;}
	for (i= 11 ; i< 12 ;i++) { reds[i]= 57.12 ; greens[i]= 0.0 ; blues[i]= 56.1 ;}
	for (i= 12 ; i< 13 ;i++) { reds[i]= 62.985 ; greens[i]= 0.0 ; blues[i]= 60.945 ;}
	for (i= 13 ; i< 14 ;i++) { reds[i]= 68.085 ; greens[i]= 0.0 ; blues[i]= 66.045 ;}
	for (i= 14 ; i< 15 ;i++) { reds[i]= 72.93 ; greens[i]= 0.0 ; blues[i]= 70.89 ;}
	for (i= 15 ; i< 16 ;i++) { reds[i]= 78.03 ; greens[i]= 0.0 ; blues[i]= 75.99 ;}
	for (i= 16 ; i< 17 ;i++) { reds[i]= 83.895 ; greens[i]= 0.0 ; blues[i]= 82.11 ;}
	for (i= 17 ; i< 18 ;i++) { reds[i]= 88.995 ; greens[i]= 0.0 ; blues[i]= 86.955 ;}
	for (i= 18 ; i< 19 ;i++) { reds[i]= 94.095 ; greens[i]= 0.0 ; blues[i]= 92.055 ;}
	for (i= 19 ; i< 20 ;i++) { reds[i]= 99.96 ; greens[i]= 0.0 ; blues[i]= 96.9 ;}
	for (i= 20 ; i< 21 ;i++) { reds[i]= 96.9 ; greens[i]= 0.0 ; blues[i]= 102.0 ;}
	for (i= 21 ; i< 22 ;i++) { reds[i]= 95.115 ; greens[i]= 0.0 ; blues[i]= 107.1 ;}
	for (i= 22 ; i< 23 ;i++) { reds[i]= 92.055 ; greens[i]= 0.0 ; blues[i]= 111.945 ;}
	for (i= 23 ; i< 24 ;i++) { reds[i]= 90.015 ; greens[i]= 0.0 ; blues[i]= 117.045 ;}
	for (i= 24 ; i< 25 ;i++) { reds[i]= 86.955 ; greens[i]= 0.0 ; blues[i]= 122.91 ;}
	for (i= 25 ; i< 26 ;i++) { reds[i]= 84.915 ; greens[i]= 0.0 ; blues[i]= 128.01 ;}
	for (i= 26 ; i< 27 ;i++) { reds[i]= 82.11 ; greens[i]= 0.0 ; blues[i]= 133.11 ;}
	for (i= 27 ; i< 28 ;i++) { reds[i]= 80.07 ; greens[i]= 0.0 ; blues[i]= 137.955 ;}
	for (i= 28 ; i< 29 ;i++) { reds[i]= 77.01 ; greens[i]= 0.0 ; blues[i]= 143.055 ;}
	for (i= 29 ; i< 30 ;i++) { reds[i]= 74.97 ; greens[i]= 0.0 ; blues[i]= 147.9 ;}
	for (i= 30 ; i< 31 ;i++) { reds[i]= 71.91 ; greens[i]= 0.0 ; blues[i]= 153.0 ;}
	for (i= 31 ; i< 32 ;i++) { reds[i]= 70.125 ; greens[i]= 0.0 ; blues[i]= 158.1 ;}
	for (i= 32 ; i< 33 ;i++) { reds[i]= 67.065 ; greens[i]= 0.0 ; blues[i]= 163.965 ;}
	for (i= 33 ; i< 34 ;i++) { reds[i]= 65.025 ; greens[i]= 0.0 ; blues[i]= 169.065 ;}
	for (i= 34 ; i< 35 ;i++) { reds[i]= 61.965 ; greens[i]= 0.0 ; blues[i]= 173.91 ;}
	for (i= 35 ; i< 36 ;i++) { reds[i]= 59.925 ; greens[i]= 0.0 ; blues[i]= 179.01 ;}
	for (i= 36 ; i< 37 ;i++) { reds[i]= 57.12 ; greens[i]= 0.0 ; blues[i]= 184.11 ;}
	for (i= 37 ; i< 38 ;i++) { reds[i]= 55.08 ; greens[i]= 0.0 ; blues[i]= 188.955 ;}
	for (i= 38 ; i< 39 ;i++) { reds[i]= 52.02 ; greens[i]= 0.0 ; blues[i]= 194.055 ;}
	for (i= 39 ; i< 40 ;i++) { reds[i]= 49.98 ; greens[i]= 0.0 ; blues[i]= 199.92 ;}
	for (i= 40 ; i< 41 ;i++) { reds[i]= 46.92 ; greens[i]= 0.0 ; blues[i]= 200.94 ;}
	for (i= 41 ; i< 42 ;i++) { reds[i]= 44.88 ; greens[i]= 0.0 ; blues[i]= 201.96 ;}
	for (i= 42 ; i< 43 ;i++) { reds[i]= 42.075 ; greens[i]= 0.0 ; blues[i]= 204.0 ;}
	for (i= 43 ; i< 44 ;i++) { reds[i]= 40.035 ; greens[i]= 0.0 ; blues[i]= 205.02 ;}
	for (i= 44 ; i< 45 ;i++) { reds[i]= 36.975 ; greens[i]= 0.0 ; blues[i]= 206.04 ;}
	for (i= 45 ; i< 46 ;i++) { reds[i]= 34.935 ; greens[i]= 0.0 ; blues[i]= 208.08 ;}
	for (i= 46 ; i< 47 ;i++) { reds[i]= 31.875 ; greens[i]= 0.0 ; blues[i]= 209.1 ;}
	for (i= 47 ; i< 48 ;i++) { reds[i]= 30.09 ; greens[i]= 0.0 ; blues[i]= 210.885 ;}
	for (i= 48 ; i< 49 ;i++) { reds[i]= 27.03 ; greens[i]= 0.0 ; blues[i]= 211.905 ;}
	for (i= 49 ; i< 50 ;i++) { reds[i]= 24.99 ; greens[i]= 0.0 ; blues[i]= 212.925 ;}
	for (i= 50 ; i< 51 ;i++) { reds[i]= 21.93 ; greens[i]= 4.08 ; blues[i]= 214.965 ;}
	for (i= 51 ; i< 52 ;i++) { reds[i]= 19.89 ; greens[i]= 7.905 ; blues[i]= 215.985 ;}
	for (i= 52 ; i< 53 ;i++) { reds[i]= 17.085 ; greens[i]= 11.985 ; blues[i]= 217.005 ;}
	for (i= 53 ; i< 54 ;i++) { reds[i]= 15.045 ; greens[i]= 16.065 ; blues[i]= 219.045 ;}
	for (i= 54 ; i< 55 ;i++) { reds[i]= 11.985 ; greens[i]= 19.89 ; blues[i]= 220.065 ;}
	for (i= 55 ; i< 56 ;i++) { reds[i]= 9.945 ; greens[i]= 23.97 ; blues[i]= 222.105 ;}
	for (i= 56 ; i< 57 ;i++) { reds[i]= 6.885 ; greens[i]= 28.05 ; blues[i]= 223.125 ;}
	for (i= 57 ; i< 58 ;i++) { reds[i]= 5.1 ; greens[i]= 31.875 ; blues[i]= 223.89 ;}
	for (i= 58 ; i< 59 ;i++) { reds[i]= 2.04 ; greens[i]= 35.955 ; blues[i]= 225.93 ;}
	for (i= 59 ; i< 60 ;i++) { reds[i]= 0.0 ; greens[i]= 40.035 ; blues[i]= 226.95 ;}
	for (i= 60 ; i< 61 ;i++) { reds[i]= 0.0 ; greens[i]= 44.115 ; blues[i]= 227.97 ;}
	for (i= 61 ; i< 62 ;i++) { reds[i]= 0.0 ; greens[i]= 47.94 ; blues[i]= 230.01 ;}
	for (i= 62 ; i< 63 ;i++) { reds[i]= 0.0 ; greens[i]= 52.02 ; blues[i]= 231.03 ;}
	for (i= 63 ; i< 64 ;i++) { reds[i]= 0.0 ; greens[i]= 56.1 ; blues[i]= 233.07 ;}
	for (i= 64 ; i< 65 ;i++) { reds[i]= 0.0 ; greens[i]= 59.925 ; blues[i]= 234.09 ;}
	for (i= 65 ; i< 66 ;i++) { reds[i]= 0.0 ; greens[i]= 64.005 ; blues[i]= 235.11 ;}
	for (i= 66 ; i< 67 ;i++) { reds[i]= 0.0 ; greens[i]= 68.085 ; blues[i]= 236.895 ;}
	for (i= 67 ; i< 68 ;i++) { reds[i]= 0.0 ; greens[i]= 71.91 ; blues[i]= 237.915 ;}
	for (i= 68 ; i< 69 ;i++) { reds[i]= 0.0 ; greens[i]= 75.99 ; blues[i]= 238.935 ;}
	for (i= 69 ; i< 70 ;i++) { reds[i]= 0.0 ; greens[i]= 80.07 ; blues[i]= 240.975 ;}
	for (i= 70 ; i< 71 ;i++) { reds[i]= 0.0 ; greens[i]= 83.895 ; blues[i]= 241.995 ;}
	for (i= 71 ; i< 72 ;i++) { reds[i]= 0.0 ; greens[i]= 87.975 ; blues[i]= 244.035 ;}
	for (i= 72 ; i< 73 ;i++) { reds[i]= 0.0 ; greens[i]= 92.055 ; blues[i]= 245.055 ;}
	for (i= 73 ; i< 74 ;i++) { reds[i]= 0.0 ; greens[i]= 95.88 ; blues[i]= 246.075 ;}
	for (i= 74 ; i< 75 ;i++) { reds[i]= 0.0 ; greens[i]= 99.96 ; blues[i]= 248.115 ;}
	for (i= 75 ; i< 76 ;i++) { reds[i]= 0.0 ; greens[i]= 104.04 ; blues[i]= 248.88 ;}
	for (i= 76 ; i< 77 ;i++) { reds[i]= 0.0 ; greens[i]= 108.12 ; blues[i]= 249.9 ;}
	for (i= 77 ; i< 78 ;i++) { reds[i]= 0.0 ; greens[i]= 111.945 ; blues[i]= 251.94 ;}
	for (i= 78 ; i< 79 ;i++) { reds[i]= 0.0 ; greens[i]= 116.025 ; blues[i]= 252.96 ;}
	for (i= 79 ; i< 80 ;i++) { reds[i]= 0.0 ; greens[i]= 120.105 ; blues[i]= 255.0 ;}
	for (i= 80 ; i< 81 ;i++) { reds[i]= 0.0 ; greens[i]= 123.93 ; blues[i]= 250.92 ;}
	for (i= 81 ; i< 82 ;i++) { reds[i]= 0.0 ; greens[i]= 128.01 ; blues[i]= 248.115 ;}
	for (i= 82 ; i< 83 ;i++) { reds[i]= 0.0 ; greens[i]= 132.09 ; blues[i]= 245.055 ;}
	for (i= 83 ; i< 84 ;i++) { reds[i]= 0.0 ; greens[i]= 135.915 ; blues[i]= 241.995 ;}
	for (i= 84 ; i< 85 ;i++) { reds[i]= 0.0 ; greens[i]= 139.995 ; blues[i]= 238.935 ;}
	for (i= 85 ; i< 86 ;i++) { reds[i]= 0.0 ; greens[i]= 144.075 ; blues[i]= 235.11 ;}
	for (i= 86 ; i< 87 ;i++) { reds[i]= 0.0 ; greens[i]= 147.9 ; blues[i]= 232.05 ;}
	for (i= 87 ; i< 88 ;i++) { reds[i]= 0.0 ; greens[i]= 151.98 ; blues[i]= 228.99 ;}
	for (i= 88 ; i< 89 ;i++) { reds[i]= 0.0 ; greens[i]= 156.06 ; blues[i]= 225.93 ;}
	for (i= 89 ; i< 90 ;i++) { reds[i]= 0.0 ; greens[i]= 159.885 ; blues[i]= 223.125 ;}
	for (i= 90 ; i< 91 ;i++) { reds[i]= 0.0 ; greens[i]= 163.965 ; blues[i]= 219.045 ;}
	for (i= 91 ; i< 92 ;i++) { reds[i]= 0.0 ; greens[i]= 168.045 ; blues[i]= 215.985 ;}
	for (i= 92 ; i< 93 ;i++) { reds[i]= 0.0 ; greens[i]= 172.125 ; blues[i]= 212.925 ;}
	for (i= 93 ; i< 94 ;i++) { reds[i]= 0.0 ; greens[i]= 175.95 ; blues[i]= 210.12 ;}
	for (i= 94 ; i< 95 ;i++) { reds[i]= 0.0 ; greens[i]= 180.03 ; blues[i]= 207.06 ;}
	for (i= 95 ; i< 96 ;i++) { reds[i]= 0.0 ; greens[i]= 184.11 ; blues[i]= 204.0 ;}
	for (i= 96 ; i< 97 ;i++) { reds[i]= 0.0 ; greens[i]= 187.935 ; blues[i]= 199.92 ;}
	for (i= 97 ; i< 98 ;i++) { reds[i]= 0.0 ; greens[i]= 192.015 ; blues[i]= 197.115 ;}
	for (i= 98 ; i< 99 ;i++) { reds[i]= 0.0 ; greens[i]= 196.095 ; blues[i]= 194.055 ;}
	for (i= 99 ; i< 100 ;i++) { reds[i]= 0.0 ; greens[i]= 199.92 ; blues[i]= 190.995 ;}
	for (i= 100 ; i< 101 ;i++) { reds[i]= 0.0 ; greens[i]= 200.94 ; blues[i]= 187.935 ;}
	for (i= 101 ; i< 102 ;i++) { reds[i]= 0.0 ; greens[i]= 202.98 ; blues[i]= 184.11 ;}
	for (i= 102 ; i< 103 ;i++) { reds[i]= 0.0 ; greens[i]= 205.02 ; blues[i]= 181.05 ;}
	for (i= 103 ; i< 104 ;i++) { reds[i]= 0.0 ; greens[i]= 207.06 ; blues[i]= 177.99 ;}
	for (i= 104 ; i< 105 ;i++) { reds[i]= 0.0 ; greens[i]= 209.1 ; blues[i]= 174.93 ;}
	for (i= 105 ; i< 106 ;i++) { reds[i]= 0.0 ; greens[i]= 210.885 ; blues[i]= 172.125 ;}
	for (i= 106 ; i< 107 ;i++) { reds[i]= 0.0 ; greens[i]= 211.905 ; blues[i]= 168.045 ;}
	for (i= 107 ; i< 108 ;i++) { reds[i]= 0.0 ; greens[i]= 213.945 ; blues[i]= 164.985 ;}
	for (i= 108 ; i< 109 ;i++) { reds[i]= 0.0 ; greens[i]= 215.985 ; blues[i]= 161.925 ;}
	for (i= 109 ; i< 110 ;i++) { reds[i]= 0.0 ; greens[i]= 218.025 ; blues[i]= 159.12 ;}
	for (i= 110 ; i< 111 ;i++) { reds[i]= 0.0 ; greens[i]= 220.065 ; blues[i]= 156.06 ;}
	for (i= 111 ; i< 112 ;i++) { reds[i]= 0.0 ; greens[i]= 222.105 ; blues[i]= 153.0 ;}
	for (i= 112 ; i< 113 ;i++) { reds[i]= 0.0 ; greens[i]= 223.125 ; blues[i]= 148.92 ;}
	for (i= 113 ; i< 114 ;i++) { reds[i]= 0.0 ; greens[i]= 224.91 ; blues[i]= 146.115 ;}
	for (i= 114 ; i< 115 ;i++) { reds[i]= 0.0 ; greens[i]= 226.95 ; blues[i]= 143.055 ;}
	for (i= 115 ; i< 116 ;i++) { reds[i]= 0.0 ; greens[i]= 228.99 ; blues[i]= 139.995 ;}
	for (i= 116 ; i< 117 ;i++) { reds[i]= 0.0 ; greens[i]= 231.03 ; blues[i]= 136.935 ;}
	for (i= 117 ; i< 118 ;i++) { reds[i]= 0.0 ; greens[i]= 233.07 ; blues[i]= 133.11 ;}
	for (i= 118 ; i< 119 ;i++) { reds[i]= 0.0 ; greens[i]= 234.09 ; blues[i]= 130.05 ;}
	for (i= 119 ; i< 120 ;i++) { reds[i]= 0.0 ; greens[i]= 235.875 ; blues[i]= 126.99 ;}
	for (i= 120 ; i< 121 ;i++) { reds[i]= 0.0 ; greens[i]= 237.915 ; blues[i]= 123.93 ;}
	for (i= 121 ; i< 122 ;i++) { reds[i]= 0.0 ; greens[i]= 239.955 ; blues[i]= 121.125 ;}
	for (i= 122 ; i< 123 ;i++) { reds[i]= 0.0 ; greens[i]= 241.995 ; blues[i]= 117.045 ;}
	for (i= 123 ; i< 124 ;i++) { reds[i]= 0.0 ; greens[i]= 244.035 ; blues[i]= 113.985 ;}
	for (i= 124 ; i< 125 ;i++) { reds[i]= 0.0 ; greens[i]= 245.055 ; blues[i]= 110.925 ;}
	for (i= 125 ; i< 126 ;i++) { reds[i]= 0.0 ; greens[i]= 247.095 ; blues[i]= 108.12 ;}
	for (i= 126 ; i< 127 ;i++) { reds[i]= 0.0 ; greens[i]= 248.88 ; blues[i]= 105.06 ;}
	for (i= 127 ; i< 128 ;i++) { reds[i]= 0.0 ; greens[i]= 250.92 ; blues[i]= 102.0 ;}
	for (i= 128 ; i< 129 ;i++) { reds[i]= 0.0 ; greens[i]= 252.96 ; blues[i]= 97.92 ;}
	for (i= 129 ; i< 130 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 95.115 ;}
	for (i= 130 ; i< 131 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 92.055 ;}
	for (i= 131 ; i< 132 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 88.995 ;}
	for (i= 132 ; i< 133 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 85.935 ;}
	for (i= 133 ; i< 134 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 82.11 ;}
	for (i= 134 ; i< 135 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 79.05 ;}
	for (i= 135 ; i< 136 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 75.99 ;}
	for (i= 136 ; i< 137 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 72.93 ;}
	for (i= 137 ; i< 138 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 70.125 ;}
	for (i= 138 ; i< 139 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 66.045 ;}
	for (i= 139 ; i< 140 ;i++) { reds[i]= 0.0 ; greens[i]= 255.0 ; blues[i]= 62.985 ;}
	for (i= 140 ; i< 141 ;i++) { reds[i]= 6.12 ; greens[i]= 255.0 ; blues[i]= 59.925 ;}
	for (i= 141 ; i< 142 ;i++) { reds[i]= 11.985 ; greens[i]= 255.0 ; blues[i]= 57.12 ;}
	for (i= 142 ; i< 143 ;i++) { reds[i]= 19.125 ; greens[i]= 255.0 ; blues[i]= 54.06 ;}
	for (i= 143 ; i< 144 ;i++) { reds[i]= 24.99 ; greens[i]= 255.0 ; blues[i]= 51.0 ;}
	for (i= 144 ; i< 145 ;i++) { reds[i]= 31.11 ; greens[i]= 255.0 ; blues[i]= 46.92 ;}
	for (i= 145 ; i< 146 ;i++) { reds[i]= 37.995 ; greens[i]= 255.0 ; blues[i]= 44.115 ;}
	for (i= 146 ; i< 147 ;i++) { reds[i]= 44.115 ; greens[i]= 255.0 ; blues[i]= 41.055 ;}
	for (i= 147 ; i< 148 ;i++) { reds[i]= 51.0 ; greens[i]= 255.0 ; blues[i]= 37.995 ;}
	for (i= 148 ; i< 149 ;i++) { reds[i]= 57.12 ; greens[i]= 255.0 ; blues[i]= 34.935 ;}
	for (i= 149 ; i< 150 ;i++) { reds[i]= 62.985 ; greens[i]= 255.0 ; blues[i]= 31.11 ;}
	for (i= 150 ; i< 151 ;i++) { reds[i]= 70.125 ; greens[i]= 255.0 ; blues[i]= 28.05 ;}
	for (i= 151 ; i< 152 ;i++) { reds[i]= 75.99 ; greens[i]= 255.0 ; blues[i]= 24.99 ;}
	for (i= 152 ; i< 153 ;i++) { reds[i]= 82.11 ; greens[i]= 255.0 ; blues[i]= 21.93 ;}
	for (i= 153 ; i< 154 ;i++) { reds[i]= 88.995 ; greens[i]= 255.0 ; blues[i]= 19.125 ;}
	for (i= 154 ; i< 155 ;i++) { reds[i]= 95.115 ; greens[i]= 255.0 ; blues[i]= 15.045 ;}
	for (i= 155 ; i< 156 ;i++) { reds[i]= 102.0 ; greens[i]= 255.0 ; blues[i]= 11.985 ;}
	for (i= 156 ; i< 157 ;i++) { reds[i]= 108.12 ; greens[i]= 255.0 ; blues[i]= 8.925 ;}
	for (i= 157 ; i< 158 ;i++) { reds[i]= 113.985 ; greens[i]= 255.0 ; blues[i]= 6.12 ;}
	for (i= 158 ; i< 159 ;i++) { reds[i]= 121.125 ; greens[i]= 255.0 ; blues[i]= 3.06 ;}
	for (i= 159 ; i< 160 ;i++) { reds[i]= 126.99 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 160 ; i< 161 ;i++) { reds[i]= 133.11 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 161 ; i< 162 ;i++) { reds[i]= 139.995 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 162 ; i< 163 ;i++) { reds[i]= 146.115 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 163 ; i< 164 ;i++) { reds[i]= 153.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 164 ; i< 165 ;i++) { reds[i]= 159.12 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 165 ; i< 166 ;i++) { reds[i]= 164.985 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 166 ; i< 167 ;i++) { reds[i]= 172.125 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 167 ; i< 168 ;i++) { reds[i]= 177.99 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 168 ; i< 169 ;i++) { reds[i]= 184.11 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 169 ; i< 170 ;i++) { reds[i]= 190.995 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 170 ; i< 171 ;i++) { reds[i]= 197.115 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 171 ; i< 172 ;i++) { reds[i]= 204.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 172 ; i< 173 ;i++) { reds[i]= 210.12 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 173 ; i< 174 ;i++) { reds[i]= 215.985 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 174 ; i< 175 ;i++) { reds[i]= 223.125 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 175 ; i< 176 ;i++) { reds[i]= 228.99 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 176 ; i< 177 ;i++) { reds[i]= 235.11 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 177 ; i< 178 ;i++) { reds[i]= 241.995 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 178 ; i< 179 ;i++) { reds[i]= 248.115 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 179 ; i< 180 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 180 ; i< 181 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 181 ; i< 182 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 182 ; i< 183 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 183 ; i< 184 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 184 ; i< 185 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 185 ; i< 186 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 186 ; i< 187 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 187 ; i< 188 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 188 ; i< 189 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 189 ; i< 190 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 190 ; i< 191 ;i++) { reds[i]= 255.0 ; greens[i]= 248.88 ; blues[i]= 0.0 ;}
	for (i= 191 ; i< 192 ;i++) { reds[i]= 255.0 ; greens[i]= 255.0 ; blues[i]= 0.0 ;}
	for (i= 192 ; i< 193 ;i++) { reds[i]= 255.0 ; greens[i]= 250.92 ; blues[i]= 0.0 ;}
	for (i= 193 ; i< 194 ;i++) { reds[i]= 255.0 ; greens[i]= 247.095 ; blues[i]= 0.0 ;}
	for (i= 194 ; i< 195 ;i++) { reds[i]= 255.0 ; greens[i]= 243.015 ; blues[i]= 0.0 ;}
	for (i= 195 ; i< 196 ;i++) { reds[i]= 255.0 ; greens[i]= 238.935 ; blues[i]= 0.0 ;}
	for (i= 196 ; i< 197 ;i++) { reds[i]= 255.0 ; greens[i]= 235.11 ; blues[i]= 0.0 ;}
	for (i= 197 ; i< 198 ;i++) { reds[i]= 255.0 ; greens[i]= 231.03 ; blues[i]= 0.0 ;}
	for (i= 198 ; i< 199 ;i++) { reds[i]= 255.0 ; greens[i]= 226.95 ; blues[i]= 0.0 ;}
	for (i= 199 ; i< 200 ;i++) { reds[i]= 255.0 ; greens[i]= 223.125 ; blues[i]= 0.0 ;}
	for (i= 200 ; i< 201 ;i++) { reds[i]= 255.0 ; greens[i]= 219.045 ; blues[i]= 0.0 ;}
	for (i= 201 ; i< 202 ;i++) { reds[i]= 255.0 ; greens[i]= 214.965 ; blues[i]= 0.0 ;}
	for (i= 202 ; i< 203 ;i++) { reds[i]= 255.0 ; greens[i]= 210.885 ; blues[i]= 0.0 ;}
	for (i= 203 ; i< 204 ;i++) { reds[i]= 255.0 ; greens[i]= 207.06 ; blues[i]= 0.0 ;}
	for (i= 204 ; i< 205 ;i++) { reds[i]= 255.0 ; greens[i]= 202.98 ; blues[i]= 0.0 ;}
	for (i= 205 ; i< 206 ;i++) { reds[i]= 255.0 ; greens[i]= 198.9 ; blues[i]= 0.0 ;}
	for (i= 206 ; i< 207 ;i++) { reds[i]= 255.0 ; greens[i]= 195.075 ; blues[i]= 0.0 ;}
	for (i= 207 ; i< 208 ;i++) { reds[i]= 255.0 ; greens[i]= 190.995 ; blues[i]= 0.0 ;}
	for (i= 208 ; i< 209 ;i++) { reds[i]= 255.0 ; greens[i]= 186.915 ; blues[i]= 0.0 ;}
	for (i= 209 ; i< 210 ;i++) { reds[i]= 255.0 ; greens[i]= 183.09 ; blues[i]= 0.0 ;}
	for (i= 210 ; i< 211 ;i++) { reds[i]= 255.0 ; greens[i]= 179.01 ; blues[i]= 0.0 ;}
	for (i= 211 ; i< 212 ;i++) { reds[i]= 255.0 ; greens[i]= 174.93 ; blues[i]= 0.0 ;}
	for (i= 212 ; i< 213 ;i++) { reds[i]= 255.0 ; greens[i]= 171.105 ; blues[i]= 0.0 ;}
	for (i= 213 ; i< 214 ;i++) { reds[i]= 255.0 ; greens[i]= 167.025 ; blues[i]= 0.0 ;}
	for (i= 214 ; i< 215 ;i++) { reds[i]= 255.0 ; greens[i]= 162.945 ; blues[i]= 0.0 ;}
	for (i= 215 ; i< 216 ;i++) { reds[i]= 255.0 ; greens[i]= 159.12 ; blues[i]= 0.0 ;}
	for (i= 216 ; i< 217 ;i++) { reds[i]= 255.0 ; greens[i]= 155.04 ; blues[i]= 0.0 ;}
	for (i= 217 ; i< 218 ;i++) { reds[i]= 255.0 ; greens[i]= 150.96 ; blues[i]= 0.0 ;}
	for (i= 218 ; i< 219 ;i++) { reds[i]= 255.0 ; greens[i]= 146.88 ; blues[i]= 0.0 ;}
	for (i= 219 ; i< 220 ;i++) { reds[i]= 255.0 ; greens[i]= 143.055 ; blues[i]= 0.0 ;}
	for (i= 220 ; i< 221 ;i++) { reds[i]= 255.0 ; greens[i]= 138.975 ; blues[i]= 0.0 ;}
	for (i= 221 ; i< 222 ;i++) { reds[i]= 255.0 ; greens[i]= 134.895 ; blues[i]= 0.0 ;}
	for (i= 222 ; i< 223 ;i++) { reds[i]= 255.0 ; greens[i]= 131.07 ; blues[i]= 0.0 ;}
	for (i= 223 ; i< 224 ;i++) { reds[i]= 255.0 ; greens[i]= 126.99 ; blues[i]= 0.0 ;}
	for (i= 224 ; i< 225 ;i++) { reds[i]= 255.0 ; greens[i]= 122.91 ; blues[i]= 0.0 ;}
	for (i= 225 ; i< 226 ;i++) { reds[i]= 255.0 ; greens[i]= 119.085 ; blues[i]= 0.0 ;}
	for (i= 226 ; i< 227 ;i++) { reds[i]= 255.0 ; greens[i]= 115.005 ; blues[i]= 0.0 ;}
	for (i= 227 ; i< 228 ;i++) { reds[i]= 255.0 ; greens[i]= 110.925 ; blues[i]= 0.0 ;}
	for (i= 228 ; i< 229 ;i++) { reds[i]= 255.0 ; greens[i]= 107.1 ; blues[i]= 0.0 ;}
	for (i= 229 ; i< 230 ;i++) { reds[i]= 255.0 ; greens[i]= 103.02 ; blues[i]= 0.0 ;}
	for (i= 230 ; i< 231 ;i++) { reds[i]= 255.0 ; greens[i]= 98.94 ; blues[i]= 0.0 ;}
	for (i= 231 ; i< 232 ;i++) { reds[i]= 255.0 ; greens[i]= 95.115 ; blues[i]= 0.0 ;}
	for (i= 232 ; i< 233 ;i++) { reds[i]= 255.0 ; greens[i]= 91.035 ; blues[i]= 0.0 ;}
	for (i= 233 ; i< 234 ;i++) { reds[i]= 255.0 ; greens[i]= 86.955 ; blues[i]= 0.0 ;}
	for (i= 234 ; i< 235 ;i++) { reds[i]= 255.0 ; greens[i]= 82.875 ; blues[i]= 0.0 ;}
	for (i= 235 ; i< 236 ;i++) { reds[i]= 255.0 ; greens[i]= 79.05 ; blues[i]= 0.0 ;}
	for (i= 236 ; i< 237 ;i++) { reds[i]= 255.0 ; greens[i]= 74.97 ; blues[i]= 0.0 ;}
	for (i= 237 ; i< 238 ;i++) { reds[i]= 255.0 ; greens[i]= 70.89 ; blues[i]= 0.0 ;}
	for (i= 238 ; i< 239 ;i++) { reds[i]= 255.0 ; greens[i]= 67.065 ; blues[i]= 0.0 ;}
	for (i= 239 ; i< 240 ;i++) { reds[i]= 255.0 ; greens[i]= 62.985 ; blues[i]= 0.0 ;}
	for (i= 240 ; i< 241 ;i++) { reds[i]= 255.0 ; greens[i]= 58.905 ; blues[i]= 0.0 ;}
	for (i= 241 ; i< 242 ;i++) { reds[i]= 255.0 ; greens[i]= 55.08 ; blues[i]= 0.0 ;}
	for (i= 242 ; i< 243 ;i++) { reds[i]= 255.0 ; greens[i]= 51.0 ; blues[i]= 0.0 ;}
	for (i= 243 ; i< 244 ;i++) { reds[i]= 255.0 ; greens[i]= 46.92 ; blues[i]= 0.0 ;}
	for (i= 244 ; i< 245 ;i++) { reds[i]= 255.0 ; greens[i]= 43.095 ; blues[i]= 0.0 ;}
	for (i= 245 ; i< 246 ;i++) { reds[i]= 255.0 ; greens[i]= 39.015 ; blues[i]= 0.0 ;}
	for (i= 246 ; i< 247 ;i++) { reds[i]= 255.0 ; greens[i]= 34.935 ; blues[i]= 0.0 ;}
	for (i= 247 ; i< 248 ;i++) { reds[i]= 255.0 ; greens[i]= 31.11 ; blues[i]= 0.0 ;}
	for (i= 248 ; i< 249 ;i++) { reds[i]= 255.0 ; greens[i]= 27.03 ; blues[i]= 0.0 ;}
	for (i= 249 ; i< 250 ;i++) { reds[i]= 255.0 ; greens[i]= 22.95 ; blues[i]= 0.0 ;}
	for (i= 250 ; i< 251 ;i++) { reds[i]= 255.0 ; greens[i]= 19.125 ; blues[i]= 0.0 ;}
	for (i= 251 ; i< 252 ;i++) { reds[i]= 255.0 ; greens[i]= 15.045 ; blues[i]= 0.0 ;}
	for (i= 252 ; i< 253 ;i++) { reds[i]= 255.0 ; greens[i]= 10.965 ; blues[i]= 0.0 ;}
	for (i= 253 ; i< 254 ;i++) { reds[i]= 255.0 ; greens[i]= 6.885 ; blues[i]= 0.0 ;}
	for (i= 254 ; i< 255 ;i++) { reds[i]= 255.0 ; greens[i]= 3.06 ; blues[i]= 0.0 ;}
	for (i= 255 ; i< 256 ;i++) { reds[i]= 255.0 ; greens[i]= 0.0 ; blues[i]= 0.0 ;}
}

// This function returns the PETIS Header Value of the specified 
// tag  (e.g., "Measurement period") as a string. Returns "" if the tag is not found.
function getTagPETIS(headerfile,tag) {
	header = File.openAsString(headerfile);
	index1 = indexOf(header,tag);
	if (index1 == -1) return "";
	index1 = indexOf(header, ":", index1);
	if (index1 == -1) return "";
	index2 = indexOf(header, "\n", index1);
	value = substring(header, index1+1, index2);
	return value;
}


// This function returns the value of the specified 
// tag  (e.g., "Depth") as a string. Returns "" if the tag is not found.
function getTag(tag) {
	info = getImageInfo();
	index1 = indexOf(info, tag);
	if (index1==-1) return "";
	index1 = indexOf(info, ":", index1);
	if (index1==-1) return "";
	index2 = indexOf(info, "\n", index1);
	value = substring(info, index1+1, index2);
	return value;
}

// This function returns the value of the specified 
// tag  (e.g., "Depth") as a string. Returns "" if the tag is not found.
function getTag2(tag) {
	info = getImageInfo();
	index1 = indexOf(info, tag);
	if (index1==-1) return "";
	index1 = indexOf(info, "(", index1);
	if (index1==-1) return "";
	index2 = indexOf(info, ")", index1);
	value = substring(info, index1+1, index2);
	return value;
}

// This function returns the value of the frame number (in info, Slice: 1/xxxx)
// as a string. Returns "" if the frame number is not found.
function getFrame(tag) {
	info = getImageInfo();
	index1 = indexOf(info, tag);
	if (index1==-1) return "";
	index1 = indexOf(info, "/", index1);
	if (index1==-1) return "";
	index2 = indexOf(info, "\n", index1);
	value = substring(info, index1+1, index2);
	return value;
}

function getMin(array){
	n=lengthOf(array);
	minValue=array[0];
	for (i=0; i<n; i++){
		//print(array[i]);
		if (minValue > array[i]) {
			minValue=array[i];
		}
	}
	return minValue;
}

function getMax(array){
	n=lengthOf(array);
	maxValue=array[0];
	for (i=0; i<n; i++){
		//print(array[i]);
		if (maxValue < array[i]) {
			maxValue=array[i];
		}
	}
	return maxValue;
}

function getTime1(){
	getDateAndTime(year, month, dayOfWeek, dayOfMonth, hour, minute, second, msec);
	TimeString ="time= "+year;
	month = month + 1;
	if (month<10) {TimeString = TimeString+"0";}
	TimeString =TimeString+month;
	if (dayOfMonth<10) {TimeString = TimeString+"0";}
	TimeString = TimeString+dayOfMonth+"_";
	if (hour<10) {TimeString = TimeString+"0";}
	TimeString = TimeString+hour;
	if (minute<10) {TimeString = TimeString+"0";}
	TimeString = TimeString+minute;
	if (second<10) {TimeString = TimeString+"0";}
	TimeString = TimeString+second;
	TimeString=split(TimeString," ");
	time=TimeString[1];
	return time;
}

function kakko_hiki(string) {
	string=split(string,"(");
	string=split(string[0],")");
	string=string[0];
	return string;
}


/*
macro "Import PETIS Text Images" {
// This Macro Make a stack image from the PETIS text images.
// Select the directory and check the image property.
//1. Set Target Text Image.
dir = getDirectory("");
list = getFileList(dir);
depth = lengthOf(list) - 1;
fileNameHead=substring(list[0],0,lengthOf(list[0])-5);

//2. Read Image info from PETIS Image header file of "xx_0.txt".
headerFileName=dir+fileNameHead+0+".txt";
tag =  "Measurement period";
period=getTagPETIS(headerFileName,tag);
//print("period="+period); 
period = 0 + substring(period, 0, lengthOf(period)-4); // convert to number
tag =  "File name";
name=getTagPETIS(headerFileName,tag);
name=substring(name,1,lengthOf(name));

//3. Read Image width, Heights and File Name from PETIS Image "xx_1.txt".
run("Text Image... ", "open="+dir+fileNameHead+1+".txt");
width=getWidth();
height=getHeight();
close();

//4. Reset Image info.
Dialog.create("Open Text Images as one Stack");
Dialog.addString("Image Name", name);
Dialog.addNumber("width", width);
Dialog.addNumber("height", height);
Dialog.addNumber("frame start", 1);
Dialog.addNumber("frame end", depth);
Dialog.addNumber("frame period (sec)", period);
Dialog.show();
name = Dialog.getString();
width = Dialog.getNumber();
height = Dialog.getNumber();
start =  Dialog.getNumber();
end = Dialog.getNumber();
period =  Dialog.getNumber();
frame = end - start +1;

//5. Make New Image for New Stack Image
tmpname="stacking_"+name;
newImage(tmpname, "32-bit Black", width, height, frame);
run("Properties...", "width="+width+" height="+height+" channels=1 depth=1 frames="+frame+" unit=mm pixel_width=1.10000 pixel_height=1.10000 voxel_depth=1.00000 frame="+period+" origin=0,0");

//6. Copy and Past
for (i=start; i<=end; i++){
	//print("open="+dir+fileNameHead+i+".txt");
	run("Text Image... ", "open="+dir+fileNameHead+i+".txt");
	run("Select All");
	run("Copy");
	close();
	selectWindow(tmpname);
	setSlice(i-start+1);
	run("Paste");
}
rename(name);

//7. Set New Image Contrast and Color.
setslice = floor(frame/2)+1;
setSlice(setslice);
run("Grays");
run("Brightness/Contrast...");
run("Enhance Contrast", "saturated=0.5");
setSlice(1);
selectWindow("B&C");
run("Close");
}//EOMacro
*/