Computer-implemented methods and systems for image analysis of multiband images of geographic regions are described, including a method by one or more computer executing executable instructions stored in one or more non-transitory, tangible, computer readable media, the method comprising: receiving one or more multiband image of a geographic region, the one or more multiband image having pixels; generating a grey level co-occurrence matrix for the pixels in the one or more multiband image; generating a surface index for the one or more multiband image containing information indicative of a surface type represented by one or more of the pixels in the one or more multiband image; and classifying the pixels of the one or more multiband image into one of a group of predefined land cover classes, based on the surface index in combination with the grey level co-occurrence matrix.