diff --git a/processing/pre/FromEBSD/Hex2Cub.cpp b/processing/pre/FromEBSD/Hex2Cub.cpp
new file mode 100644
index 000000000..e0ec3f581
--- /dev/null
+++ b/processing/pre/FromEBSD/Hex2Cub.cpp
@@ -0,0 +1,205 @@
+#include <iostream>
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <math.h>
+
+struct
+{
+	float phi1, Phi, phi2;
+	float ci, iq, fit, avgIQ;
+	int phase, ds;
+} data[3000][3000];
+
+int round2( double x )
+{
+	double fl, ce;
+
+	fl = floor( x );
+	ce = ceil( x );
+	if( fabs( fl-x ) < fabs( ce-x ) )
+		return fl;
+	else
+		return ce;
+}
+
+char lineBuffer[200];
+
+int ReadFileInfo(FILE *oimStream, int *xmax, int *ymax, double *stepSize )
+{
+	double x, y, dxmax;
+	long count;
+	
+	*stepSize = 0;
+	dxmax = 0;
+	*xmax = 0;
+
+	while( fgets( lineBuffer, 200, oimStream ) != NULL )
+	{
+		if( lineBuffer[0] == '#' )
+		{
+			if( strcmp( lineBuffer, "# GRID: SqrGrid" ) == 0 )
+			{
+				printf("\nThe file is already a square grid file.\nProgram terminated.");
+				return 0;
+			}
+			count = 0;
+			continue;
+		}
+		if( sscanf( lineBuffer, "%*lf %*lf %*lf %lf %lf %*lf %*lf %*i %*i %*lf %*lf", &x, &y ) != 2 )
+			return 0;
+		if( *stepSize == 0 && x != 0 )
+			*stepSize = x;
+		if( x > dxmax )
+		{
+			dxmax = x;
+			(*xmax)++;
+		}
+		count++;
+	}
+	(*xmax)++;
+	*ymax = (int)(count / *xmax );
+
+	return 1;
+}
+
+
+int main(int argc, char* argv[])
+{
+	int xx, yy, zz, xlimit, ylimit, zlimit, xlimitOut, xOut;
+	double stepSize;
+	char zFilename[50], zOutFilename[50], filename[50];
+	FILE *inStream, *outStream;
+	int zStartNumber, zEndNumber;
+
+	printf( "\nFilenames must have the format ""root_xxx.ang"""
+		"\nwith xxx indicating a 3-digit integer"
+		"\nEnter oim map filename-root, the start integer and the end integer number of the files: " );
+
+	scanf( "%s %i %i", filename, &zStartNumber, &zEndNumber );
+	zlimit = zEndNumber-zStartNumber+1;
+
+	//read the first data file and get all necessary start data
+	sprintf( zFilename, "%s_%03i.ang", filename, zStartNumber );
+	if( (inStream = fopen( zFilename, "r" )) == NULL )
+	{
+		printf( "\nCan't open %s", zFilename );
+		exit( 1 );
+	}
+	if( ReadFileInfo( inStream, &xlimit, &ylimit, &stepSize ) == 0 )
+	{
+		printf( "\nWrong file format in %s", filename );
+		exit( 1 );
+	}
+	fclose( inStream );
+	
+	for( zz=0; zz<zlimit; zz++ )
+	{
+		printf("\nReading");
+		sprintf( zFilename, "%s_%03i.ang", filename, zz+zStartNumber );
+		sprintf( zOutFilename, "%s_cub_%03i.ang", filename, zz+zStartNumber );
+		if( (inStream = fopen( zFilename, "r" )) != NULL )
+		{
+			outStream = fopen( zOutFilename, "w" );
+			//read file header
+			do
+			{
+				if( fscanf( inStream, "%[^\n]\n", lineBuffer ) == EOF )
+				{
+					printf( "\nEarly end of file encountered in ANG file" ); 
+					exit(1);
+				}
+				//write the file header
+				if( lineBuffer[0] == '#' )
+				{
+					if( strcmp( lineBuffer, "# GRID: HexGrid" ) == 0 )
+						fprintf( outStream, "# GRID: SqrGrid\n" );
+					else
+						fprintf( outStream, "%s\n", lineBuffer );
+				}
+
+			}
+			while( lineBuffer[0] == '#' );
+
+			for( yy=0; yy<ylimit; yy++)
+			{
+				printf(".");
+				for( xx=0; xx<xlimit; xx++)
+				{
+					//t1: pattern quality, iq: confidence index, avgIQ:  average Image Quality
+					if( sscanf( lineBuffer, "%f %f %f %*f %*f %f %f %i %i %f %f", 
+										&data[xx][yy].phi1, 
+										&data[xx][yy].Phi, 
+										&data[xx][yy].phi2, 
+										&data[xx][yy].iq, 
+										&data[xx][yy].ci, 
+										&data[xx][yy].phase, 
+										&data[xx][yy].ds, 
+										&data[xx][yy].fit, 
+										&data[xx][yy].avgIQ ) != 9 )
+					{
+						printf( "\nWrong file format in %s", filename );
+						exit( 1 );
+					}
+
+					//read the next line buffer if there is any.
+					//ylimit%2 only for hexagonal grid data (odd lines are shorter by 1 pixel)
+					if( yy%2 == 1 && xx == xlimit-1 )
+					{
+						data[xx][yy].phi1 = data[xx-1][yy].phi1; 
+						data[xx][yy].Phi = data[xx-1][yy].Phi; 
+						data[xx][yy].phi2 = data[xx-1][yy].phi2; 
+						data[xx][yy].iq = data[xx-1][yy].iq; 
+						data[xx][yy].ci = data[xx-1][yy].ci; 
+						data[xx][yy].phase = data[xx-1][yy].phase; 
+						data[xx][yy].ds = data[xx-1][yy].ds; 
+						data[xx][yy].fit = data[xx-1][yy].fit;
+						data[xx][yy].avgIQ = data[xx-1][yy].avgIQ;
+					}
+					else 
+					{
+						if( fscanf( inStream, "%[^\n]\n", lineBuffer ) == EOF )
+						{
+							printf( "\nEarly end of file encountered in ANG file" ); 
+							exit(1);
+						}
+					}
+				}//end for(x...
+			}//end for(y...
+			fclose( inStream );
+			printf("\nWriting");
+			//the step size in y-direction (=0.866*stepSizeX) 
+			//is the new step size for x and y
+			xlimitOut = round2( (double)xlimit/0.866); 
+			for( yy=0; yy<ylimit; yy++)
+			{
+				printf(".");
+				for( xx=0; xx<xlimitOut; xx++)
+				{
+					xOut = round2( (double)xx * 0.866 );
+					fprintf( outStream, "%f %f %f %f %f %f %f %i %i %f %f\n", 
+										data[xOut][yy].phi1, 
+										data[xOut][yy].Phi, 
+										data[xOut][yy].phi2, 
+										xx * stepSize * 0.866,
+										yy * stepSize * 0.866,
+										data[xOut][yy].iq, 
+										data[xOut][yy].ci, 
+										data[xOut][yy].phase, 
+										data[xOut][yy].ds, 
+										data[xOut][yy].fit, 
+										data[xOut][yy].avgIQ );
+				}//end for( xx...
+			}//end for( yy...
+			fclose( outStream );
+		}
+		else
+		{
+			printf( "\nExpected file %s does not exist", zFilename );
+			exit( 1 );
+		}
+	}
+
+	return 0;
+}
+
diff --git a/processing/pre/FromEBSD/SpectralMethodFromEBDS b/processing/pre/FromEBSD/SpectralMethodFromEBDS
new file mode 100644
index 000000000..8d68d8cb1
--- /dev/null
+++ b/processing/pre/FromEBSD/SpectralMethodFromEBDS
@@ -0,0 +1,42 @@
+#!/usr/bin/env python
+# -*- coding: UTF-8 no BOM -*-
+import sys, array, re, math
+file_in = open(sys.argv[1]+'.ang','r')
+file_out = open(sys.argv[1]+'_material.config_pt1','w')
+file_out.write('<homogenization>\n[SX]\ntype isostrain\nNgrains 1 \n\n<microstructure>\n')
+file_out2 = open(sys.argv[1]+'_material.config_pt2','w')
+file_out2.write('<texture>\n')
+data=array.array('f',[0.0,0.0,0.0])
+resY = 1
+resX = 0
+dimX = 0.0
+dimY = 0.0
+i=0
+for line in file_in:
+  if line[0]!='#':
+    i = i +1
+    data[0] = float(re.findall('\S*',str(line))[0])/2.0/math.pi*360.0
+    data[1] = float(re.findall('\S*',str(line))[2])/2.0/math.pi*360.0
+    data[2] = float(re.findall('\S*',str(line))[4])/2.0/math.pi*360.0
+    phasef = float(re.findall('\S*',str(line))[20])
+    if phasef>10000.0:
+      phase=1
+    else:
+      phase=2
+    if float(re.findall('\S*',str(line))[6]) < dimX:
+      dimX=0.0
+      resX=1
+      resY= resY + 1
+    else:
+      resX=resX+1
+      dimX = float(re.findall('\S*',str(line))[6])
+    file_out.write('[Grain'+str(i)+']\ncrystallite 1\n(constituent)'+\
+      ' phase '+str(phase)+' texture '+str(i)+' fraction 1.0\n')
+    file_out2.write('[Grain'+str(i)+']\n(gauss) phi1 '+str(data[0])+\
+      ' Phi '+str(data[1])+' phi2 '+str(data[2])+' scatter 0.0 fraction 1.0\n')
+dimY = float(re.findall('\S*',str(line))[8])
+dimZ = min(dimX/resX,dimY/resY)
+file_out3 = open(sys.argv[1]+'.geom','w')
+file_out3.write('resolution x '+str(resX)+' y '+str(resY)+' z 1 \ndimension a '+str(dimX)+' b '+str(dimY)+' c '+str(dimZ)+'\nhomogenization 1\n')
+for x in xrange(resX*resY):
+  file_out3.write(str(x+1)+'\n')
diff --git a/processing/pre/patchFromReconstructedBoundaries b/processing/pre/FromEBSD/patchFromReconstructedBoundaries
old mode 100755
new mode 100644
similarity index 100%
rename from processing/pre/patchFromReconstructedBoundaries
rename to processing/pre/FromEBSD/patchFromReconstructedBoundaries