User Tools

Site Tools


PolyLX - the MATLAB™ toolbox for quantitative analysis of microstructures


The textural analysis is a powerful, but underused tool of petrostructural analysis. Except acquirement of common statistical parameters, this technique can significantly improve understanding of processes of grain nucleation and grain growth, can bring insights on the role of surface energies or quantify duration of metamorphic and magmatic cooling events as long as appropriate thermodynamical data for studied mineral exist. This technique also allows systematic evaluation of degree of preferred orientations of grain boundaries in conjunction with their frequencies. This may help to better understand the mobility of grain boundaries and precipitations or removal of different mineral phases.

object_model.jpgWe introduce a new open platform, object-oriented MATLAB® toolbox PolyLX providing several core routines for data exchange, visualization and analysis of microstructural data, which can be run on any platform supported by MATLAB®. Detailed descriptions of toolbox routines and methods of implementation of new techniques are given in PolyLX Reference Manual.

The quantitative textural analysis rely on detailed and precise description of grain sizes, grain shapes, grain boundaries as well as preferred orientations of grain and grain boundaries. The manual or semi-automatic digitalizing of scanned or digital photographs could be done in any software which is capable to export ESRI Shapefiles. The desktop GIS packages which allow create topologically correct geometries are best candidates. The grains are generated in form of adjacent polygons and each grain is labelled with a unique ID and a phase name. The unique ID is used to keep relations integrity with imported results of analysis, when visualization capabilities of GIS packages are in demand. Once all grains are digitized, PolyLX toolbox can build a set of polylines representing individual grain boundaries. Individual boundaries span from neighbouring triple points and are labelled by a unique ID and by ID’s and phase names of adjacent grains. For those having access to ESRI ArcView 3.x, the available PolyLX extension could be used to generate boundary files and perform some error-checking tasks.

The grain and/or grain boundary geometries together can be directly imported from ESRI ArcView shapefiles. In MATLAB® environment grain and grain boundaries are represented as objects, which hold all important properties, such as ID, Phase name, Area, Perimeter, Axial ration, Length, Width etc. All these values can be obtained on-the-fly and statistically analysed. Along with wide spectrum of mathematical and statistical methods of data exploration available in MATLAB™, PolyLX incorporates several geologically specific techniques:

Grain and grain boundary shape

shape_props.jpgShape is extremely difficult property to measure, or even to define in a precise manner. Perhaps this is why there are so many proposed shape measures, none of which has been proved as entirely satisfactory. A shape measure should possess several desirable properties. Obviously, objects with different shapes should yield different measures, and similarly shapes should yield similar values regardless of the size or orientation of the object. Unfortunately, a shape measure possessing these properties may be a chimera; it has been prove that no single measure can be unique to only one shape. Therefore, there is a wide spectrum of single value measuring methods available in PolyLX toolbox.

Grain and grain boundary preferred orientation

To obtain data of preferred orientation, several techniques are implemented in PolyLX toolbox. The most general one is method of analysis matrix of inertia. This method can be applied on individual grains or boundaries as well as on a set of grains or grain boundaries. In latter case, the result is weighted by size of objects, which is welcome in specific tasks and differs from the results obtained from orientation analysis based on histograms (rose diagrams) or Fisher distribution. Another group of routines using approach of direction dependent projection of grain or grain boundaries (PAROR, SURFOR and PARIS) and advertised by Pannozo (1983) are fully implemented.

Spatial statistics

One of the most important aspects of quantitative texture analysis is description of spatial characteristics of grains or grain boundaries. PolyLX contains several routines dealing with spatial distribution of grains or grain boundaries (grain density method, nearest neighbour analysis, spatial pattern index) or evaluating deviation from random distribution (contact frequency and contact area methods). Crystal size distributions: PolyLX offers implementation of method to construct the CSD plots using technique described by Peterson (1996).

Strain analysis

rfphi.jpgSeveral techniques to estimate strain are available. Classical ones as Rf/φ, centre-to-centre method and Harvey & Fergusson (1981) as well as some of their recent modifications like DTNNM (Mulchrone, 2000) or area weighted Rf/φ are implemented.

Want to know more ?

You can download my short presentation of PolyLX toolbox here.

You can look on short descriptions of routines or check the recent filelist linked to manual pages.

Download PolyLX

Click here to download latest R3.1 version of the PolyLX Toolbox.

Publications where PolyLX has been used

  1. Ulrich S., Schulmann K. and Casey M., 2002. Microstructural evolution and rheological behaviour of marbles deformed at different crustal levels, Journal of Structural Geology, 24/5, 979-995. DOI: 10.1016/S0191-8141(01)00132-8
  2. Baratoux, L., Lexa, O., Cosgrove, J. W. & Schulmann, K., 2005. The quantitative link between fold geometry, mineral fabric and mechanical anisotropy: as exemplified by the deformation of amphibolites across a regional metamorphic gradient. Journal of Structural Geology, 27(4), 707-730. DOI: 10.1016/j.jsg.2005.01.001
  3. Baratoux, L., Schulmann, K., Ulrich, S. & Lexa, O., 2005. Contrasting microstructures and deformation mechanisms in metagabbro mylonites contemporaneously deformed under different temperatures (c. 650 and c. 750 °C). In: Deformation Mechanisms, Rheology and Tectonics: from Minerals to the Lithosphere (eds Gapais, D., Brun, J. P. & Cobbold, P. R.) Geological Society Special Publications, pp. 97-125, Geological Society of London, London, United Kingdom.
  4. Lexa, O., Štípská, P., Schulmann, K., Baratoux, L. & Kröner, A., 2005. Contrasting textural record of two distinct metamorphic events of similar P-T conditions and different durations. Journal of Metamorphic Geology, 23(8), 649-666. DOI: 10.1111/j.1525-1314.2005.00601.x
  5. Barraud, J., 2006. The use of watershed segmentation and GIS software for textural analysis of thin sections, Journal of Volcanology and Geothermal Research, 154, 1-2, 17-33. DOI: 10.1016/j.jvolgeores.2005.09.017
  6. Závada, P., Schulmann, K., Konopásek, J., Ulrich, S. & Lexa, O., 2007. Extreme ductility of feldspar aggregates—Melt-enhanced grain boundary sliding and creep failure: Rheological implications for felsic lower crust. Journal of Geophysical Research, B, Solid Earth and Planets, 112, B10210, doi:10.1029/2006JB004820.DOI: 10.1029/2006JB004820
  7. Jeřábek, P., Stunitz, H., Heilbronner, R., Lexa, O. & Schulmann, K., 2007. Microstructural-deformation record of an orogen-parallel extension in the Vepor Unit, West Carpathians. Journal of Structural Geology, 29(11), 1722-1743. DOI: 10.1016/j.jsg.2007.09.002
  8. Machek, M., Špaček, P., Ulrich, S. & Heidelbach, F., 2007. Origin and orientation of microporosity in eclogites of different microstructure studied by ultrasound and microfabric analysis, Engineering Geology, 89, 266-277. DOI: 10.1016/10.1016/j.enggeo.2006.11.001
  9. Hasalová, P., Schulmann, K., Lexa, O., Štípská, P., Hrouda, F., Ulrich, S., Haloda, J. & Týcová, P., 2008. Origin of migmatites by deformation enhanced melt infiltration of orthogneiss: a new model based on quantitative microstructural analysis. Journal of Metamorphic Geology, 26, 29-53. DOI: 10.1111/j.1525-1314.2007.00743.x
  10. Schulmann, K., Martelat, J-E., Ulrich, S., Lexa, O., Štípská, P. & Becker, J.K., 2008. Evolution of microstructure and melt topology in partially molten granitic mylonite: Implications for rheology of felsic middle crust, Journal of Geophysical Research, B, Solid Earth and Planets, 113, B10406, doi:10.1029/2007JB005508. DOI: 10.1029/2007JB005508
  11. Závada, P., Schulmann, K., Lexa, O., Hrouda F., Haloda, J. & Týcová, P., 2009. The mechanism of flow and fabric development in mechanically anisotropic trachyte lava, Journal of Structural Geology, 31, 11, 1295-1307. DOI:10.1016/j.jsg.2009.04.002
  12. Kratinová, Z., Ježek J., Schulmann K., Hrouda F., Shail R. K. & Lexa, O., 2010. Noncoaxial K-feldspar and AMS subfabrics in the Land's End granite, Cornwall: Evidence of magmatic fabric decoupling during late deformation and matrix crystallization, Journal of Geophysical Research, B, Solid Earth and Planets, 115, B09104, DOI:10.1029/2009JB006714
  13. Franěk, J., Schulmann, K., Lexa, O., Ulrich, S., Štípská, P., Haloda, J. & Týcová, P., 2011. Origin of felsic granulite microstructure by heterogeneous decomposition of alkali feldspar and extreme weakening of orogenic lower crust during the Variscan orogeny, Journal of Metamorphic Geology, 29, 1, 103-130. DOI:10.1111/j.1525-1314.2010.00911.x

Any comments ?

polylx/home.txt · Last modified: 2017/06/02 11:49 (external edit)