Objective: Parkinson's disease (PD) is a common neurological disorder with variable clinical manifestations and magnetic resonance imaging (MRI) findings. We propose a handcrafted image classification model that can accurately (i) classify different PD stages, (ii) detect comorbid dementia, and (iii) discriminate PD-related motor symptoms. Methods: Selected image datasets from three PD studies were used to develop the classification model. Our proposed novel automated system was developed in four phases: (i) texture features are extracted from the non-fixed size patches. In the feature extraction phase, a pyramid histogram-oriented gradient (PHOG) image descriptor is used. (ii) In the feature selection phase, four feature selectors: neighborhood component analysis (NCA), Chi2, minimum redundancy maximum relevancy (mRMR), and ReliefF are used to generate four feature vectors. (iii) Two classifiers: k-nearest neighbor (kNN) and support vector machine (SVM) are used in the classifica