The present invention relates to a method for computing external crystal shapes from X-Ray Diffraction Data (XRD) of a substance. Each diffraction peak arises from a set of crystal planes and the peak width is related to the thickness of the crystal in a direction perpendicular to these set of planes. The crystal shape is actually given by the mathematical envelope of the pairs of planes corresponding to each diffraction peak.