A modulation format identification (MFI) scheme based on multi-dimensional amplitude features is proposed for elastic optical networks. According to the multi-dimensional amplitude features, incoming polarization division multiplexed (PDM) signals can be identified as QPSK, 8QAM, 16QAM, 32QAM, 64QAM and 128QAM signals using the k-nearest neighbors (KNNs) algorithm in the digital coherent receivers. The proposed scheme does not require any prior training or optical signal-to-noise ratio (OSNR) information. The performance of the proposed MFI scheme is verified based on numerical simulations with 28GBaud PDM-QPSK/-8QAM/-16QAM/-32QAM/-64QAM/-128QAM signals. The results show that the proposed scheme can achieve 100% of the correct MFI rate for all six modulation formats when the OSNR values are greater than their thresholds corresponding to the 20% forward error correction (FEC) related to a BER of 2.4 × 10−2. Meanwhile, the effects of residual chromatic dispersion, polarization mode dispersion and fiber nonlinearities on the proposed scheme are also explored. Finally, the computational complexity of the proposed scheme is analyzed, which is compared with relevant MFI schemes. The work indicates that the proposed technique could be regarded as a good candidate for identifying modulation formats up to 128QAM.