Bone fractures are caused by diseases or accidents and are a widespread problem throughout the globe. Worldwide, 1.6 millions of hip fractures occur every year and are expected to rise to 6.3 millions in 2050. The current gold standard to assess fracture risk is the Dual-energy X-ray Absorptiometry (DXA), which provides a projected image of the bone from which areal bone mineral density is extracted. Ultrasound techniques have been proposed as non invasive alternatives. Recently, estimates of cortical thickness and porosity, obtained by Bi-Directional Axial Transmission (BDAT) in a pilot clinical study, have been shown to be associated with non-traumatic fractures in post menopausal women. Cortical parameters were derived from the comparison between experimental and theoretical guided modes. This model-based inverse approach failed for the patients associated with poor guided mode information. Moreover, even if the fracture discrimination ability was found similar to DXA, it remained moderate. The goal of this paper is to explore if these two limitations could be overcome by using automatic classification tools, independent of any waveguide model. To this end, a dynamic machine learning approach based on a Support Vector Machine (SVM) has been used to classify ultrasonic guided wave spectrum images measured by BDAT on post menopausal women with or without non-traumatic fractures. This approach has then been improved using parameters tuned by Bat Algorithm Optimization (BOA). The applied methodology focused on the extraction of texture features through a gray level co-occurrence matrix, structural comparison and histograms. The results accuracy was assessed using a confusion matrix. The effectiveness of the learning approach reached an accuracy of 92.31%. INDEX TERMS cortical bone, osteoporosis, fracture discrimination, quantitative ultrasound, guided waves, support vector machine, bat algorithm, guided wave spectrum image.