// This Macros are made for PETIS Image Analysis.
// Created by Naoki Kawachi.
// Version 2.1.1 08/May/2012

//////////Macor in Pull Down Menu//////////
macro "---PETIS Format Import Utilities---" {}
macro "PETIS_3 Text Format			[f1]" {importPETIS3TextImage();}
macro "ALL PETIS_3 Text Format		[f2]" {importAllPTETIS3TextImages();}
macro "PETIS_4 Analyze Format		[f3]" {importPETIS4AnaImage();}
macro "Gamma Camera Text Format		[f10]" {importGammaCTextImage();}
macro "-" {}
macro "---PETIS Images Analysis Utilities---" {}
macro "Set PETIS Image Frames		[f4]" {setPETIS_ImageFrame();}
macro "PETIS Image Multiplier 		[f5]" {PETISimageMultiplier();}
macro "Re-Stack PETIS Image 		[f6]" {reStackPETIS_Images();}
macro "-" {}
macro "---PETIS Numerical Analysis Utilities---" {}
macro "ROIs Submit					[f7]" {roiSubmit();}
macro "Get PETIS TAC 				[f8]" {getPETIS_TAC();}
macro "TACs Analysis				[f9]" {tacAnalysis();}
macro "-" {}
macro "---Presantation Utilities---" {}
macro "  PETIS Color Bar			" {setColorBar();}
macro "  PETIS Time Bar 			" {setTimeBar();}
macro "  PETIS Lap Time Stamp		" {setLapTimeStamp();}
macro "  PETIS Montage 				" {makePETIS_Montage();}
macro "-" {}
macro "---Look up Tables---" {}
macro "  RGB Color" {
    getLut(reds, greens, blues);
	G256_RGB_LUT(reds,greens,blues);
	setLut(reds, greens, blues);
}
macro "  RGBw Color" {
    getLut(reds, greens, blues);
	G256_RGB_LUT(reds,greens,blues);
	reds[255] = 255; greens[255] = 255; blues[255] = 255;	//Write White
	setLut(reds, greens, blues);
}
macro "  NIH Color" {
    getLut(reds, greens, blues);
	G256_NIH_LUT(reds,greens,blues);
	setLut(reds, greens, blues);
}
macro "  NIHw Color" {
    getLut(reds, greens, blues);
	G256_NIH_LUT(reds,greens,blues);
	reds[255] = 255; greens[255] = 255; blues[255] = 255;	//Write White
	setLut(reds, greens, blues);
}
macro "  FIRE Color" {run("Fire");}
macro "  Gray Color" {run("Grays");}
macro "  Invert LUT" {
    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);
}
//////////END of Pull Down Menu//////////


//////////Macor in Tool Bar//////////
//1st Button: Blank
macro 'Unused Tool -' {}
//2nd Button: PETIS Import
var aCmds = newMenu("PETIS Data Import Menu Tool", newArray("PETIS_3 Text Format", "ALL PETIS_3 Text Format", "-", "PETIS_4 Analyze Format", "-", "Save as Tiff"));
macro "PETIS Data Import Menu Tool - C000T0809P C000T5809E C000T8809T C000Td809I C000Tf809S C030T4f09I C030T6f09/ C030Taf09O" {
	cmd = getArgument();
	if (cmd=="PETIS_3 Text Format") {importPETIS3TextImage();}
	if (cmd=="ALL PETIS_3 Text Format") {importAllPTETIS3TextImages();}
	if (cmd=="PETIS_4 Analyze Format") {importPETIS4AnaImage();}
	if (cmd=="Save as Tiff") {saveAs("Tiff", "");}
}
//3rd Button: PETIS Image Analysis
var bCmds = newMenu("PETIS Images Analysis Menu Tool", newArray("Set PETIS Image Frames","PETIS Image Multiplier","-","Average PETIS Image","Integrated PETIS Images","-","Re-Stack PETIS Image"));
macro "PETIS Images Analysis Menu Tool - C000T0809I C000T2809M C000T8809A C000Tb809G C000Tf809E C030T2f09A C030T7f09N C030Tdf09A" {
	cmd = getArgument();
	if (cmd=="Set PETIS Image Frames") {setPETIS_ImageFrame();}
	if (cmd=="PETIS Image Multiplier") {PETISimageMultiplier();}
	if (cmd=="Average PETIS Image") {
		n=nSlices;
		run("Z Project...", "start=1 stop="+n+" projection=[Average Intensity]");
	}
	if (cmd=="Integrated PETIS Images") {
		n=nSlices;
		run("Z Project...", "start=1 stop="+n+" projection=[Sum Slices]");
	}
	if (cmd=="Re-Stack PETIS Image") {reStackPETIS_Images();}
}
//4th Button: ROI Submit
macro "ROIs Submit Action Tool -  C000 T0a10R  C000 T6a10O  C000 Tda10I  C000 Tfa10s" {roiSubmit();}
//5th Button: TAC Analysis
macro "TACs Analysis Action Tool -  C000 T0a10T  C000 T5a10A  C000 Tba10C  C000 Tfa10s" {tacAnalysis();}
//6th Button: Presentation Utilities
var cCmds = newMenu("Presentation Utilities Menu Tool", newArray("PETIS Color Bar","PETIS Time Bar","PETIS Lap Time Stamp","PETIS Montage"));
macro "Presentation Utilities Menu Tool - C005F1155C40cF6155Ca0bFa155Cc05F1655Cf40F6655Cf80Fa655Cfb0F1a55Cff1F6a55CffeFaa55" {
	cmd = getArgument();
	if (cmd=="PETIS Color Bar") {setColorBar();}
	if (cmd=="PETIS Time Bar") {setTimeBar();}
	if (cmd=="PETIS Lap Time Stamp") {setLapTimeStamp();}
	if (cmd=="PETIS Montage") {makePETIS_Montage();}
}
//7th Button: Color LUT
var dCmds = newMenu("Color of Look up Tables Menu Tool", newArray("Brightness and Contrast","-","RGB Color","RGBw Color","-","NIH Color","NIHw Color","-","FIRE Color","Gray Color","-","Invert LUT"));
macro "Color of Look up Tables Menu Tool - Cf00 T0510L  C0f0 T7b10U  C50f Tef10T Cfff" {
	cmd = getArgument();
	if (cmd=="Brightness and Contrast") {run("Brightness/Contrast...");}
	if (cmd=="RGB Color") {getLut(reds, greens, blues);G256_RGB_LUT(reds,greens,blues);setLut(reds, greens, blues);}
	if (cmd=="RGBw Color") {
		getLut(reds, greens, blues);
		G256_RGB_LUT(reds,greens,blues);
		reds[255] = 255; greens[255] = 255; blues[255] = 255;	//Write White
		setLut(reds, greens, blues);
	}
	if (cmd=="NIH Color") {getLut(reds, greens, blues);G256_NIH_LUT(reds,greens,blues);setLut(reds, greens, blues);}
	if (cmd=="NIHw Color") {
		getLut(reds, greens, blues);
		G256_NIH_LUT(reds,greens,blues);
		reds[255] = 255; greens[255] = 255; blues[255] = 255;	//Write White
		setLut(reds, greens, blues);
	}
	if (cmd=="FIRE Color") {run("Fire");}
	if (cmd=="Gray Color") {run("Grays");}
	if (cmd=="Invert LUT") {
		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);
	}
}
//////////END of Tool Bar Menu//////////


/////////////////MACRO FUNCTIONS///////////////////////
// (1)
// This Function Make a stack image from the PETIS3 text images.
// Select the directory and check the image property.
function importPETIS3TextImage() {
//1. Set Target Text Images Directory.
	dir=getDirectory("Choose Text Images Data Folder");
	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");
}


// (2)
// This Function Make a stack image from the All PETIS Text Images in a Directory.
// Select the directory and check the image property.
function importAllPTETIS3TextImages() {
// Set Target Text Images Directory.
	dir=getDirectory("Choose Text Image Folders Directory");
	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]))=="/")){
			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();}
		}			
	}
}


// (3)
// This Function Make a stack image from the PETIS4 analyze format images.
// Select the directory and check the image property.
function importPETIS4AnaImage() {
//1. Set Target Analyze Format Images Directory.
	path=File.openDialog("Choose Analyze Format Data");
	dirName=File.directory;
	fileName=File.getName(path);
	expTmpName=split(fileName,".");
	if(expTmpName[1]=="inf"){
		expName=split(fileName,"P.");
	} else if(expTmpName[1]=="hdr" || expTmpName[1]=="img" ){
		expName=split(expTmpName[0],"_");
	} else {
		exit("Please select a analyze format image file");
	}
	anaImgFile=dirName+expName[0]+"_ANA.img";
	anaHeadFile=dirName+expName[0]+"_ANA.hdr";
	headerFile =dirName+expName[0]+"P.inf";
//2. Read Image
	open(anaImgFile);
//3. Set New Image Header and Name
	comment = "PETIS Comment:\t"+getTagPETIS(headerFile,"COMMENT");
	period = getTagPETIS(headerFile,"FRAME TIME(0)");
	start_time = getTagPETIS(headerFile,"TIME");
	exp_date = getTagPETIS(headerFile,"DATE");
	PETIS_date ="PETIS Measurement Date:\t"+exp_date+" "+start_time;
	PETIS_RI = "PETIS Positron RI:\t"+getTagPETIS(headerFile,"NUCLIDE")+" ("+getTagPETIS(headerFile,"DECAY TIME")+")";
	PETIS_Head_Distance="PETIS Head Distance:\t"+getTagPETIS(headerFile,"DET DISTANCE");
	PETIS_frames=getTagPETIS(headerFile,"TOTAL FRAMES");
	PETIS_header ="Info:"+comment+"\n"+PETIS_date+"\n"+PETIS_RI+"\n"+PETIS_Head_Distance;
	setMetadata(PETIS_header);
	run("Properties...", "channels=1 slices=1 frames="+PETIS_frames+" unit=mm pixel_width=1.1000 pixel_height=1.1000 voxel_depth=1.0000 frame=["+period+" sec] origin=0,0");
	rename(expName[0]);
}


// (4)
// This Function Make a stack image from the Espec Inc. text format images.
function importGammaCTextImage() {
	path=File.openDialog("Choose Text Image File (APOLLOMEC)");
	fileName=File.getName(path);
	imageName=split(fileName,".");

	lines = File.openAsString(path);
	line = split(lines,"\n");
	coment_line = split(line[0],", ");
	frame_line = split(line[1],", ");
	dim_z = frame_line[1];
	dim_line = split(line[2]," ");
	dim_x = dim_line[3];
	dim_y = dim_line[5];
	newImage(imageName[0], "32-bit Black", dim_x, dim_y, dim_z);
	for (j=5; j<5+dim_y; j++){				// Read Lines (Y)
			counts = split(line[j],"\t");	// Read Lines (X)
			for (k=1; k<1+dim_x; k++){
			setPixel(k-1,j-5,counts[k]);	// Set Values (X,Y)
			}
	}
	run("Grays");
	if(dim_z==1){
	getStatistics(area, mean, min, max, std, histogram);
	} else {
	Stack.getStatistics(voxelCount, mean, min, max, stdDev);
	}
	setMinAndMax(0.0, max);
}


// (5)
function setPETIS_ImageFrame() {
	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");
	}
}


// (6)
// This Function Multiple PETIS Image (each Frames).
function 	PETISimageMultiplier() {
	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);
}


// (7)
// This Function Make a New Stack Image.
// Choose a Target Image and Restacking.
function reStackPETIS_Images() {
////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));
			for (j=0;j<height;j++){
				for (i=0;i<width;i++){
					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");
}


// (8)
// ROI Submit and Rename
function roiSubmit() {
	run("ROI Manager...");
	setKeyDown("alt");
	roiManager("add");
}


// (9)
// This Function Generates a Time Acitvity Curve from PETIS Image(stack).
// After a ROI selected, start up this Function.
function getPETIS_TAC() {
//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();
}


// (10)
// This Function Generates Time Acitvity Curves from PETIS Image(stack).
// All ROIs are selected to generate TACs.
function tacAnalysis() {
//<<Set Colors Array>>
	lineColor=newArray("black","red","blue","green","orange","darkGray","magenta","yellow","gray","cyan","pink","lightGray");
//<<Classify Image Windows>>
	windowList = getList("window.titles");
	if (windowList.length==0)
		print("No non-image windows are open");
	else {
		print("Non-image windows:");
		for (i=0; i<windowList.length; i++)
			print("   "+windowList[i]);
	}
	imageTitle=newArray(nImages);
	if (nImages==0)print("No images are open");
	else {
		print("Image windows:");
		setBatchMode(true);
		for (i=1; i<=nImages; i++) {
			selectImage(i);
			print("   "+getTitle);
			imageTitle[i-1]=getTitle;
		}
	setBatchMode(false);
	}
//<<Classify ROIs Name>>
	nRoi=roiManager("count");
	indexesRoi=newArray(nRoi);
	indexesRoiName=newArray(nRoi);
	for (i=0; i<nRoi; i++) {
		roiManager("select", i);
		indexesRoi[i]=i;
		indexesRoiName[i]=getInfo("roi.name");
	}
//<<Set Target Image and ROIs for Analysis>>
	//1. Selection Dialog
	Dialog.create("Select Target Image and ROIs for Analysis");
	Dialog.addMessage("Select Target Image");
	Dialog.addChoice("        ", imageTitle, imageTitle[0]);
	Dialog.addMessage("Select Submitted ROIs for Analysis");
	for (i=0; i<nRoi; i++) Dialog.addCheckbox(indexesRoiName[i], true);
	Dialog.show();
	//2. Set Target Image
	targetImage = Dialog.getChoice();
	//3. Reset ROIs for Analysis
	roiSelection=newArray(nRoi);
	for (i=0; i<nRoi; i++) roiSelection[i]=Dialog.getCheckbox();
	newNRoi=0;
	for (i=0; i<nRoi; i++) {if (roiSelection[i]==1) newNRoi=newNRoi+1;}
	indexesRoi=newArray(newNRoi);
	indexesRoiName=newArray(newNRoi);
	j=0;
	for (i=0; i<nRoi; i++) {
		if (roiSelection[i]==1){
			roiManager("select", i);
			indexesRoi[j]=i;
			indexesRoiName[j]=getInfo("roi.name");
			j=j+1;
		}
	}
//<<Read Target Image Info>>
	selectWindow(targetImage);
	name=getTitle();				//Image Name
	dir=getDirectory("image");		//Image Dir
	if (dir == "") dir = "This Image on the ImageJ Memory";
	tag="Frame interval";
	str=getTag(tag);
	period=split(str," ");
	period=0+period[0];				//Frame interval (sec)
	if (period==0) period = 1;
	tag="PETIS Comment";
	comment=getTag(tag);			//PETIS Comment
	tag="PETIS Measurement Date";
	mdate=getTag(tag);				//Measurement Date
	tag="PETIS Positron RI";
	positronRI=getTag(tag);			//Name of RI
	ns = nSlices;					//Number of Frames
	nAD = ns*newNRoi;				//Number of All Data
//<<ROIs Analysis with Multi Measure to the Target Image>>
	area=newArray(nAD);mean=newArray(nAD);
	intDen=newArray(nAD);
	time=newArray(ns);mtime=newArray(ns);
	for (i=0; i<ns; i++) {
		time[i]=(i+1)*period;
		mtime[i]=(i+1)*period/60;
	}
	selectWindow(targetImage);
	setBatchMode(true);
	roiManager("select", indexesRoi);
	roiManager("Multi Measure");
	run("Set Measurements...", "area mean integrated redirect=None decimal=9");
	for (j=0; j<newNRoi; j++){
		firstRow=getResult("Area("+indexesRoiName[j]+")",0);
		if(isNaN(firstRow)==true) {
			tmpRoiName=j+1;
			for (i=0; i<ns; i++){
				area[j*ns+i] = getResult("Area"+tmpRoiName+"",i);
				mean[j*ns+i] = getResult("Mean"+tmpRoiName+"",i);
				intDen[j*ns+i] = getResult("IntDen"+tmpRoiName+"",i);
		//print(i,"+",ns,"x",j,"=",i+j*ns,": IntDen"+indexesRoiName[j]+"(NaN)= ",intDen[i+j*ns]);
			}
		}
		else {
			for (i=0; i<ns; i++){
				area[j*ns+i] = getResult("Area("+indexesRoiName[j]+")",i);
				mean[j*ns+i] = getResult("Mean("+indexesRoiName[j]+")",i);
				intDen[j*ns+i] = getResult("IntDen("+indexesRoiName[j]+")",i);
		//print(i,"+",ns,"x",j,"=",i+j*ns,": IntDen("+indexesRoiName[j]+")= ",intDen[i+j*ns]);
			}
		}
	}
//<<Plot Time Activity Curves for Each ROIs>>
	setBatchMode(false);
	Array.getStatistics(intDen, min_0, max_0, mean_0, stdDev_0);
	run("Profile Plot Options...", "width=600 height=400 minimum=0 maximum=0 interpolate draw");
	Plot.create("ROIs Analysis Result: TACs of "+targetImage+" Image", "Time(min)", "Integrated Density (cps)", mtime, newArray(ns));
	Plot.setLimits(0,period*ns/60,0,max_0*1.1);
	for (j=0; j<newNRoi; j++){
		Y_Values=newArray(ns);
		for (i=0; i<ns; i++){
			Y_Values[i]=intDen[j*ns+i];
			}
		Plot.setColor(lineColor[j%12]);
		Plot.setLineWidth(1);
		Plot.add("line", mtime, Y_Values);
		Plot.setLineWidth(2);
		Plot.add("circles", mtime, Y_Values);
		Plot.addText("--------ROI "+(j+1)+" = "+indexesRoiName[j]+" ", 0.1, 0.1+0.05*j);
	}
	Plot.show();
//<<Save Data>>
	Dialog.create("Time Activity Time Save Dialog");
	Dialog.addChoice("This TAC_Graph Save as .xls file?",newArray("Save TAC", "No"),"Save TAC");
	Dialog.addNumber("Decimal places (0-9)",4);
	//Dialog.addCheckbox("Save Area Data", false);
	Dialog.addCheckbox("Save Mean Data", false);
	Dialog.addCheckbox("Save Integrated Density Data", true);
	Dialog.addCheckbox("Save with ROIs File", true);
	Dialog.addMessage("");
	Dialog.addCheckbox("Save TAC PNG File", true);
	Dialog.show();
	ans = Dialog.getChoice();
	dn = Dialog.getNumber();
	//area_save=Dialog.getCheckbox();
	mean_save=Dialog.getCheckbox();
	intDen_save=Dialog.getCheckbox();
	roi_save=Dialog.getCheckbox();
	tacpng_save=Dialog.getCheckbox();
	if (roi_save==0) print("ROIs_will not be saved.");
	if (ans == "No"){
		exit();
	}
	saveAs("Results","");
	saveFileDir=File.directory;
	saveFileName=File.name;
	saveFilePath= saveFileDir +File.name;
	saveFileNameOnly=split(saveFileName,".");
	setBatchMode(true);
	selectWindow("Results");
	run("Close");
//<<ROI Save>>
	roiManager("Select",indexesRoi);
	roiManager("Save",saveFileDir+"/"+saveFileNameOnly[0]+"_ROIs.zip");
	f = File.open(saveFilePath);
	print(f,"##\tImage File Name = "+name);
	print(f,"##\tImage File Directory = "+dir);
	print(f,"##\tComment ="+comment);
	print(f,"##\tMeasurement date ="+mdate);
	print(f,"##\tPositron RI ="+positronRI);
	print(f,"##\tSaved ROI File = "+saveFileDir+""+saveFileNameOnly[0]+"_ROIs.zip");
//<<Write Area Info>>
	line="";
	for (j=0; j<newNRoi; j++){
		line= line +"##\tArea Info\t"+j+".  ROI Area ("+indexesRoiName[j]+") =";
		line= line +"\t"+d2s(area[j*ns],dn);
		line= line +"\t(mm2)\n";
	}
//<<Write 1st Column>>
	print(f,line);
	line="Frame"; line= line +"\tTime (sec)";
	for (j=0; j<newNRoi; j++){
		//if (area_save==1) line= line +"\tArea("+indexesRoiName[j]+")";
		if (mean_save==1) line= line +"\tMean("+indexesRoiName[j]+")";
		if (intDen_save==1) line= line +"\tIntDen("+indexesRoiName[j]+")";
	}
	print(f,line);
//<<Write Selected Data>>
	for (i=0; i<ns; i++) {
		line=d2s((i+1),0);
		line=line+"\t"+d2s(time[i],0);
		for (j=0; j<newNRoi; j++){
			//if (area_save==1) line= line +"\t"+d2s(area[j*ns+i],dn);
			if (mean_save==1) line= line +"\t"+d2s(mean[j*ns+i],dn);
			if (intDen_save==1) line= line +"\t"+d2s(intDen[j*ns+i],dn);
		}
		print(f,line);
	}
	File.close(f);
	//<<TAC_PNG_Save>>
	if (tacpng_save==1) saveAs("PNG", saveFileDir+"/"+saveFileNameOnly[0]+"_TACs.png");
}


// (11)
// This Function Add Color Bar at Leftside of Stack Image.
function setColorBar() {
	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];
		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);
				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");
	}
}


// (12)
// This Function Add Time Bar at Bottom of Stack Image.
function setTimeBar() {
	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();
}


// (13)
// This Function Add Laptime on the Stack Image.
function 	setLapTimeStamp() {
	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);
		}
	}
}


//(14)
// This Function Make a montage from a PETIS stack image.
// Select the directory and check the image property.
function makePETIS_Montage() {
////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);
}





///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////


// 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++){
		if (minValue > array[i]) {
			minValue=array[i];
		}
	}
	return minValue;
}


function getMax(array){
	n=lengthOf(array);
	maxValue=array[0];
	for (i=0; i<n; 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;
}


////////////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 ;}
}