Spatial Interpolation of Regional Rainfall Data of Thailand Using Different Methods

Main Article Content

Usawan Upakhot
Chatchai Tantasirin


Distributed areal rainfall is important data for hydrological modeling. Spatial interpolation is used to predict the values of a continuous variable at unmeasured points. This study compared the distribution and prediction error in the spatial interpolation maps of rainfall derived from different methods. Inverse Distance Weighted (IDW) using a power (p) equal to 2 and 3, and number of neighbor point (n) equal to 15, 30, and 50 points point and Kriging Ordinary (KrigOrd) with Circular, Exponential, Gaussian, Linear, and Spherical model were used and their predictions compared. Cohen’s kappa coefficient and mean absolute error were used to compare distribution and error in spatial interpolation, respectively. The results indicated that the distribution and parameters obtained from different spatial interpolation methods were in good agreement as indicated by k values (0.9202 to 0.9999) but results of IDW predicted a gradual change, while that from KrigOrd had localized features and abrupt change. The lowest spatial interpolation error from measurements in the north, northeast, central, east and south (east coast) parts was obtained for IDW with p of 2 with differing number of neighbors equal to 15, 15, 50, 50, and 15 points respectively, while that for the south (west coast) was obtained using the KrigOrd with a linear setup. Mean absolute error (MAE) for each method was 133.88 132.49 114.83 184.01 181.54, and 142.52 millimeters or 11.30, 10.25, 10.58, 11.76, 10.90, and 5.44 percent of the mean annual rainfall, respectively.

Article Details

How to Cite
Upakhot, U., & Tantasirin, C. . (2023). Spatial Interpolation of Regional Rainfall Data of Thailand Using Different Methods. Thai Journal of Forestry, 42(2), 63–76. Retrieved from
Original Articles


Agdas, M.G., Yenen, Z. 2023. Determining land use/land cover (LULC) changes using remote sensing method in Lüleburgaz and LULC change's impacts on SDGs. European Journal of Sustainable Development, 12(1): 1-24. doi: 10.14207/ejsd.2023.v12n1p1.

Bailey, G., Li, Y., McKinney, N., Yoder, D., Wright, W. Herrero, H. 2022. Comparison of ground point giltering algorithms for high-density point clouds collected by terrestrial LiDAR. Remote Sensing, 14(4776): 1-20.

Burrough, P.A., McDonanell, R.A. 1998. Principles of Geographical Information Systems. Oxford University Press, New York, USA.

Chomwinyan, W. 2017. Geography of Thailand. Udon Thani Rajabhat University. Udon Thani, Thailand. (in Thai)

Donald, S. 1968. A two-dimensional interpolation function for irregularly-spaced data. In : Proceedings of the 1968 23rd ACM National Conference. New York, USA, pp. 517–524.

Eiamprasert, S. 2011. Comparison of Spatial Daily Rainfall Amount Interpolation Methods by Using Geographic Information System at Chao-Phra-Ya Flood Plai. M.S. Thesis, National Institute of Development Administration, Bangkok, Thailand. (in Thai)

Grubbs, F.E. 1969. Procedures for Detecting Outlying Observations in Samples. Technometrics, (11)1: 1-21. doi: :10.1080/00401706.1969.10490657.

Kumari, M., Singh, C.K., Bakimchandra, O., Basistha, A. 2016. DEM-based delineation for improving geostatistical interpolation of rainfall in mountainous region of Central Himalayas, India. Theoretical and Applied Climatology, 130(1-2): 51-58. doi:10.1007/s00704-016-1866-y.

Maddala, G.S. 1992. Introduction to Econometrics, 2nd ed. Macmillan. New York, USA.

Matheron, G. 1963. Principles of geostatistics. Economic Geology, 58(8): 1246–1266. doi: gsecongeo.58.8.1246.

McHugh, L.M. 2012. Interrater reliability: the Kappa statistic. Biochem Med (Zagreb)., 22 March 2022.

Meechamnan, T. 2017. An Areal Rainfall Interpolation derived from Ground Measurement Rainfall and GSMaP Satellite with Topography Data in Southern Thailand. M.Eng. Thesis, Kasetsart University (Kamphaeng Saen Campus), Nakhon Pathom, Thailand. (in Thai)

Piazza, A.D., Conti, F.L., Noto, L.V., Viola, F. Loggia, G.L. 2011. Comparative analysis of different techniques for spatial interpolation of rainfall data to create a serially complete monthly time series of precipitation for Sicily, Italy. International Journal of Applied Earth Observation and Geoinformation, 13(3): 396-408. doi:10.1016/j.jag.2011.01.005.

Rattanaburee, P. 2014. Comparison of the Spatial Rainfall in Upper Ping Basin and Tha Chin Basin by Spatial Interpolation Techniques. M.Eng. Thesis, Kasetsart University (Kamphaeng Saen Campus), Nakhon Pathom, Thailand. (in Thai)

Thai Meteorological Department. (n.d.). Climate of Thailand., 19 March 2023. (in Thai)

Yang, X., Xie, X., Liu, D.L., Ji, F., Wang, L. 2015. Spatial interpolation of daily rainfall data for local climate impact assessment over greater Sydney region. Advances in Meteorology, 2015: 563629. doi:10.1155/2015/563629.

Yodkongdee, E. 2017. A Comparative Study on Areal Rainfall Estimation by Spatial Interpolation Techniques in Bang Pakong Basin and Mae Klong Basin. M.Eng. Thesis, Kasetsart University (Kamphaeng Saen Campus), Nakhon Pathom, Thailand. (in Thai)