DigitalMicrographでプラゴミのIRスペクトルを調べる

//read .csv as image
image readCSV(string path,number datatype,number header)//datatype=2 for real4
{
    number fileID=OpenFileForReading(path)
    object fstream=NewStreamFromFileReference(fileID,1)
    number n=0
    string text
    while(n<header)
    {
        fstream.StreamReadTextLine(0,text)
        n=n+1
    }
    number bLinesAreRows=1
    number bSizeByCount=1
    object imgSizeObj=Alloc("ImageData_ImageDataSize")
    image img:=ImageImportTextData(path,fstream,datatype,imgSizeObj,bLinesAreRows,bSizeByCount)
    return img
}

//add tag to an image
void addTg(TagGroup Tg,image img)
{
    TagGroup imgTg=img.ImageGetTagGroup()
    TagGroupCopyTagsFrom(imgTg,Tg)
}

//get path of target directory
string filename
string targetname="ATR.csv" //"ATR.CSV" for 2019 data (uppercase character)
if(!OpenDialog(Null, "Select "+targetname+" file", "*"+targetname, filename)) exit(0)
string dirpath=GetApplicationDirectory(0,0)
result("exploring the below directory\n")
result(dirpath+"\n")

//get file list in the target diectory
TagGroup fileList=GetFilesInDirectory(dirpath,1)
//filelist.TagGroupOpenBrowserWindow("file list",0)

//get file list of "*ATR.csv"
TagGroup targetList=NewTagList()
number n
string tempname
for(n=0;n<fileList.TagGroupCountTags();n++)
{
    //check file name
    fileList.TagGroupGetTagAsString("["+n+"]:Name",tempname)
    if(tempname.len()>targetname.len())
    {
        if(tempname.right(targetname.len())==targetname)
        {
            targetList.TagGroupInsertTagAsString(infinity(),tempname)
        }
    }
}
result("number of ***"+targetname+" files: "+n+"\n")
//targetList.TagGroupOpenBrowserWindow(targetname,0)

//opening csv files if file format is the same as the first one
//opened csv image is stored into a tag group "openListAsImg".
string path //file path
number datatype=2 //datatype=2 for real4 (float)
number header=1 //number of header rows in csv
image firstImg //image of the first csv file in the targetList
TagGroup openList=NewTagList() //list of the name of the opened files
TagGroup openListAsImg=NewTagList() //image list of the opened files
number x,y,x0,y0
//first csv image
targetList.TagGroupGetIndexedTagAsString(0,path)
firstImg:=readCSV(dirpath+path,datatype,header)
firstImg.GetSize(x0,y0)
openListAsImg.TagGroupInsertTagAsArray(infinity(),firstImg)
openList.TagGroupInsertTagAsString(infinity(),path)

//opening csv files
for(n=1;n<targetList.TagGroupCountTags();n++)
{
    targetList.TagGroupGetIndexedTagAsString(n,path)
    //check size of the csv file
    readCSV(dirpath+path,datatype,header).GetSize(x,y)
    if((x==x0)&&(y==y0))
    {
        openListAsImg.TagGroupInsertTagAsArray(infinity(),readCSV(dirpath+path,datatype,header))
        openList.TagGroupInsertTagAsString(infinity(),path)
    }
}

//stack all opened csv data
number z0=openListAsImg.TagGroupCountTags()
image stackImg:=RealImage(targetname,4,x0,y0,z0)
for(n=0;n<z0;n++)
{
    openListAsImg.TagGroupGetIndexedTagAsArray(n,stackImg[0,0,n,x0,y0,n+1])
}

//set openList tag to "file list" tag
TagGroup openFileList=NewTagGroup()
openFileList.TagGroupSetTagAsTagGroup("file list",openList)

//image stackImg should have (x,y,z)=(col,row,number of opened data)
result("number of opened data files: "+z0+"\n")

//wavelength as an image
image wavelength=stackImg.slice2(0,0,0,1,y0,1,2,z0,1)
wavelength.SetName("wavelenghth")
addTg(openFileList,wavelength)
wavelength.showimage()
result("wavelength range: "+wavelength.min()+" - "+wavelength.max()+"\n")

//spetrum as an image
image spc=stackImg.slice2(1,0,0,1,y0,1,2,z0,1)
spc.SetName("spectrum")
addTg(openFileList,spc)
spc.showimage()

//spectrum as an spestrum imaging type data (3D)
image spccube=stackimg.slice3(1,0,0,2,z0,1,0,1,1,1,y0,1)
spccube.SetName("spectrum cube")
addTg(openFileList,spccube)
spccube.showimage()

//delete data which has a value larger than one
image in:=getFrontImage() //"spectrum cube" should be the front.
number x,y,z
in.Get3DSize(x,y,z)
string name
in.GetName(name)

TagGroup inTg=in.ImageGetTagGroup()
TagGroup originalFileList
inTg.TagGroupGetTagAsTagGroup("file list",originalFileList)
TagGroup cpTg=originalFileList.TagGroupClone() //copy

//store image into tag group
TagGroup imgTg=NewTagList()
image temp
number n
for(n=x-1;n>=0;n--)
{
    temp=in.slice1(n,0,0,2,z,1)
    if(temp.abs().max()>1)
    {
        cpTg.TagGroupDeleteTagWithIndex(n)
    }
    else
    {
        imgTg.TagGroupInsertTagAsArray(0,temp)
    }
}

number nn=cpTg.TagGroupCountTags()
image delImg:=RealImage(name+" del",4,nn,y,z)

for(n=0;n<nn;n++)
{
    imgTg.TagGroupGetIndexedTagAsArray(n,temp)
    delImg[n,0,0,n+1,y,z]=temp.slice3(0,0,0,1,1,1,2,1,1,0,z,1)
}

TagGroup resultTg=delImg.ImageGetTagGroup()
resultTg.TagGroupSetTagAsTagGroup("file list",cpTg)

delImg.showimage()

image spc:=GetFrontImage() //2D image can be loaded

image score,load,scree
MSA_PCA(spc,score,load,scree,0)

scree.setname("scree")
score.setname("score")
load.setname("load")

scree.showimage()
score.showimage()
load.showimage()