Ok. If you plan to use Matlab, I suggest you to use Mtex. It can save your time. Any way, the ebsd data depend on your EBSD system. For TSL's ang file. The format is like this, and you can find (x,y) in 4,5th columns.
# The fields of each line in the body of the file are as follows:
# j1 F j2 x y IQ CI Phase ID Detector Intensity Fit
# where:
# j1,F,j2: Euler angles (in radians) in Bunge's notation for describing the lattice orientations and are given in radians.
# x,y: The horizontal and vertical coordinates of the points in the scan, in microns. The origin (0,0) is defined as the top-left corner of the scan.
# IQ: The image quality parameter that characterizes the contrast of the EBSP associated with each measurement point.
# CI: The confidence index that describes how confident the software is that it has correctly indexed the EBSP, i.e., confidence that the angles are correct.
# Phase ID: The material phase identifier. This field is 0 for single phase OIM scans or 1,2,3... for multi-phase scans.
# Detector Intensity: An integer describing the intensity from whichever detector was hooked up to the OIM system at the time of data collection, typically a forward scatter detector.
# Fit: The fit metric that describes how the indexing solution matches the bands detected by the Hough transform or manually by the user.
Thanks Xinfu but that is giving avg XY position but I need absolute XY position,as you can see in uploaded pic ...i tried in TSL software by goint to export>grain file>xy position but it didn't help.