DigitalMicrographでangファイルを読み込む

###angファイル
　EBSD (Electron BackScatter Diffraction) もしくは EBSP (Electron BackScatter Pattern) の測定データをテキストデータ化したファイルフォーマット。オイラー角、xy座標その他の測定値が記録されている。「その他」が具体的に何かは装置依存？ヘッダー中には特に説明がないので自力で解読する必要あり。外部プログラムでいろいろ解析する場合は、測定点のxy座標を正方格子状のデータに変換した方が扱いやすいと思う。

データ例
Materials Science and Engineering: A, 755 (2019), 295.
Site-specific quasi in situ investigation of primary static recrystallization in a low carbon steel
Martin Diehl, Lukas Kertsch, Konstantina Traka, Dirk Helm, Dierk Raabe
https://doi.org/10.1016/j.msea.2019.02.032
文献中のデータ
Data for: Site-Specific Quasi In Situ Investigation of Primary Static Recrystallization in a Low Carbon Steel
DOI: 10.17632/7dn4g8rgzc.1

上記のデータセットの一部で読み込みを試してみたが、ファイルによってはバグがあるかも。上記データの測定点座標は三角格子状。

import-ang.s

//read .ang file
image ReadAng(string filename)
{
	number fileID=OpenFileForReading(filename)
	object fstream=NewStreamFromFileReference(fileID,1)
	//read header
	TagGroup headerTg=NewTagList()
	number pos //stream position
	string text
	string CommentBegin="#"
	number CommentCheck=1
	while(CommentCheck)
	{
		pos=fstream.StreamGetPos()
		fstream.StreamReadTextLine(0,text)
		if(text.left(1)==CommentBegin)
		{
			headerTg.TagGroupInsertTagAsString(infinity(),text)
		}
		else
		{
			CommentCheck=0
		}
	}
	//read data
	fstream.StreamSetPos(0,pos)
	number bLinesAreRows=1
	number bSizeByCount=1
	number kREAL4_DATA=2
	object imgSizeObj=Alloc( "ImageData_ImageDataSize" )
	image ang:=ImageImportTextData(filename.PathExtractBaseName(0),fstream,kREAL4_DATA,imgSizeObj,bLinesAreRows,bSizeByCount)
	
	//set header Tag
	TagGroup angTg=ang.ImageGetTagGroup()
	angTg.TagGroupSetTagAsTagGroup(filename.PathExtractFileName(0),headerTg)
	return ang
}

//find string
number findkey(string input,string key)
{
	number n=0
	if(input.len()>=key.len())
	{
		if(input.left(key.len())==key)
		{
			n=1
		}
	}
	return n
}
//get string
string getkey(string input,string key)
{
	result(input)
	string output
	output=input.right(input.len()-key.len())
	return output
}

//shaping ang image
image SetAng(image ang)
{
	number angx,angy
	ang.GetSize(angx,angy)
	string name
	ang.GetName(name)
	//extract image information
	TagGroup angTg=ang.ImageGetTagGroup()
	TagGroup headerTg
	angTg.TagGroupGetTagAsTagGroup(name+".ang",headerTg)
	string grid
	number xstep,ystep,col_odd,col_even,row
	string text
	string comment
	string output
	for(number i=0;i<headerTg.TagGroupCountTags();i++)
	{
		headerTg.TagGroupGetTagAsString("["+i+"]",text)
		//grid
		comment="# GRID: "
		if(findkey(text,comment))
		{
			grid=getkey(text,comment)
		}
		//step
		comment="# XSTEP: "
		if(findkey(text,comment))
		{
			xstep=getkey(text,comment).val()
		}
		comment="# YSTEP: "
		if(findkey(text,comment))
		{
			ystep=getkey(text,comment).val()
		}
		//size
		comment="# NCOLS_ODD: "
		if(findkey(text,comment))
		{
			col_odd=getkey(text,comment).val()
		}
		comment="# NCOLS_EVEN: "
		if(findkey(text,comment))
		{
			col_even=getkey(text,comment).val()
		}
		comment="# NROWS: "
		if(findkey(text,comment))
		{
			row=getkey(text,comment).val()
		}
	}
	if(grid=="HexGrid\n")
	{
		result("Warning: Grid is NOT square.\n")
		result("X-position is different from the original file.\n")
		result("And Y-scale of image on screen is enlarged by a factor of 2/sqrt(3).\n")
	}
	//set image
	number x,y
	x=min(col_odd,col_even)
	y=row
	image angcube:=RealImage(name,4,x,y,angx)
	//odd col
	for(number i=0;i<y/2;i++)
	{
		angcube.slice3(0,2*i,0,0,x,1,1,1,1,2,angx,1)=ang.slice3(0,i*(col_odd+col_even),0,1,x,1,2,1,1,0,angx,1)
	}
	//even col
	for(number i=1;i<y/2;i++)
	{
		angcube.slice3(0,2*i-1,0,0,x,1,1,1,1,2,angx,1)=ang.slice3(0,i*(col_odd+col_even)+col_odd,0,1,x,1,2,1,1,0,angx,1)
	}
	//set header Tag
	TagGroup angcubeTg=angcube.ImageGetTagGroup()
	angcubeTg.TagGroupSetTagAsTagGroup(name,headerTg)
	//set step size
	angcube.ImageSetDimensionCalibration(0,0,xstep,"μm",0)
	angcube.ImageSetDimensionCalibration(1,0,ystep,"μm",0)
	
	return angcube
}

//get file path of ang file
string filename
if(!OpenDialog(Null, "Select ang file", "*.ang", filename)) exit(0)
result("opening the next file\n")
result(filename+"\n")

image ang:=ReadAng(filename)
image angcube:=SetAng(ang)
angcube.showimage()