Macroseismic intensity, a useful measure of earthquake effects, is still applied in a wide range of seismological applications like seismic hazard assessments, attenuation relationships, etc. Isoseismals represent the spatial distribution of macroseismic intensities and their shapes depend on source properties, lithosphere structures, tectonic line orientations, site geology and topography. The applications ask for both the higher number of isoseismal maps and their standardization and homogenization. The point kriging gridding method for an automatic computer drawing of isoseismal maps was delivered. Smoothing rates and numerical parameters used in the kriging algorithm were tested on macroseismic data of Greek earthquakes representing different tectonic and geomorphological regimes. The optimum kriging default option was defined. Its application for four Greek earthquakes is presented and discussed from viewpoint of a broad use in recent macroseismology.