We analysed receiver functions from teleseismic events recorded at 11 broadband seismometers in the western part of Java Island, Indonesia. The stations are mostly located at three main geological environment including Northwest Java Basin, Bogor Zone, and Southern Mountains Arc. A total of about 341 receiver functions were computed using iterative time domain deconvolution. We derived shear-wave velocity structure and crustal $Vp/Vs$ ratio by inverting stacked radial receiver functions using non-linear neighbourhood algorithm. Inversion results show sediment thickness varies between 1 and 2 km thick in Western Java. Our inversion shows that crustal thickness in this region varies between 25 and 32 km. Average crustal $Vp/Vs$ ratio is estimated to be about 1.69–1.78. We hope the study may provide useful information for velocity model and crustal thickness for Indonesia region.