Image temp, spectrum, spectrumimage, datacube

GetTwoLabeledImagesWithPrompt("Please select the rotated datacube and the spectrum.","Spectral Image","Datacube",datacube,"Spectrum",spectrum)

number height, width, depth
get3dsize(datacube, width, height, depth)

spectrumimage = slice2( datacube, 0, 0, 0, 0, width, 1, 1, height, 1 )

number top, left, bottom, right

//GetFrontImage(spectrum)
//GetSelection(spectrum,top,left,bottom,right)

for (number n=1006/*left*//*Change to first channel*/; n<=1030/*right*//*Change to last channel*/; n++)
{

number cuben = n

temp = slice2( datacube, 0, 0, cuben, 0, width, 1, 1, height, 1 )

spectrumimage = spectrumimage + temp
}

showimage(spectrumimage)