Tuesday, March 9, 2010

Lab #7 - Interpolation


Interpolation is a powerful procedure that can be used in ArcGIS to predict the value of cells at locations that do not have sample points. It is based on the idea of spatial autocorrelation which is premised upon Tobler's first law of geography: "Everything is related to everything else, but near things are more related than distant things..." Of the different interpolation methods, I chose to use inverse distance weighted (IDW) and kriging to create surface area maps of precipitation in Los Angeles county. The IDW method of interpolation estimates cell values by averaging the values of sample data points in the area of each processing cell. The closer the point is to the center of the cell being estimated, the more influence, or weight, it has in the averaging process. Kriging is similar to IDW in that it weights the surrounding measured values to derive a prediction for an unmeasured location. I used ordinary kriging in my maps, which assumes that the variation in z-values is free of any structural component (drift).

My first map using IDW shows that the most rainfall is in the central region of the map, while there is less precipitation in the northern and southern parts of Los Angeles County. Overall, this season's precipitation is not too different from normal precipitation patterns. I think that IDW is an appropriate interpolation technique for this data set because IDW requires a set of points that is dense enough to capture the extent of local surface variation. Since there are over sixty automatic rain gage points, I feel that there was enough to produce a quality surface map for analysis of LA County's precipitation.

The second map I created using kriging showed overall less variance in rainfall from this season compared to the season normal. It also emphasizes in both the normal season and this current season-to-date that the most rainfall is in the central part of the county, while the least is in the north (northeast especially) and in the south. Compared to the IDW map, the kriging maps does not show as much of the local variation in precipitation, causing me to favor the IDW procedure for interpolation since I find it yields more detailed results. Ultimately, while it is predicted to be an El Nino year (which means higher levels of precipitation), these maps show that for the 2009-2010 season there does not seem to be significantly higher levels of rainfall--in Los Angeles County at least.


Source: Los Angeles County Department of Public Works and ESRI