Abstract Three-dimensional (3D) ultrasonic imaging enables the viewing of internal features in a more accurate way than cross-sectional imaging. 3D imaging requires two-dimensional (2D) ultrasonic phased arrays to resolve all three spatial dimensions through beamforming along azimuthal and elevation angles. 2D phased arrays pose a manufacturing and control challenge due to the requirement of large number of elements to satisfy the Nyquist sampling limit. This problem can be alleviated through the use of sparse ultrasonic array designs. The optimization process is critical, as this dictates the efficiency of suppression of grating lobes. The aim of this work is to establish a framework for design of optimized sparse 2D phased arrays. This is enabled by exploiting the reconfigurability of Laser Induced Phased Arrays (LIPAs). LIPAs are synthetic arrays produced by scanning one laser for ultrasound generation and another for ultrasound detection, independently of each other. The reconfigurability and decoupled ultrasonic generation and detection of this systems enables easy and fast implementation of any arbitrary array layout. The framework consists of an analytical model for parametric evaluation of designs. The model performs a parametric sweep on the generation and detection element positions for each array design in order to identify the optimal parameters, based on mean and maximum side-lobe level. In the presented framework, four array designs are utilized for the optimization process: matrix, random, Vernier and a novel array design the Rotated Array (RA).