Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

Hi I'm new to interpolations.

I'm looking for documentation to understand how inverse distance algorithm works and how to code it.

My intention is to code IDW algorithm in JavaScript to see if it is feasible in performance meaning.

Thanks !!!

share|improve this question

2 Answers 2

up vote 2 down vote accepted

Well, Wikipedia is the place to look at.

what you do is:

1) compute the distance between your location and the observations

2) define the weights as 1/distance^p

3) sum the weighted value of each observation

4) normalise the result by the sum of the weights

Note that IDW relies on one parameter(p), the power of the distance. The larger it is, the more infuence you give to the closest observation. the extreme value 0 and infinity yield a simple average and a thiessen polygonsation, respectively. The optimal parameter can be found based on cross-validation methods (but this add computational cost, of course).

Note also that you can limit the number of neighbours that are taken into account to remove the influence of far away points.

IDW is one of the fastest interpolation method, but not the one I prefer...

share|improve this answer
    
+1. However, the formula is the easy part. The hard parts are (1) efficiently selecting neighboring data (especially when implementing barriers) and (2) dealing with the exceptional cases such as when distances are zero or individual locations have multiple values. (Even many commercial packages, such as ArcGIS, don't get that correct or cannot even handle such situations.) –  whuber Mar 19 at 20:38
    
for 2, the way to deal with zero distances is to ignore other values if one distance is zero (infinite weight : IDW is an exact interpolator). Concerning 2 points at the same location, IDW will give them the same weight by design. For intelligent handling of this particular case, one should use kriging(the difference is that kriging will use the average value once, but IDW uses it twice) –  radouxju Mar 19 at 21:10
    
On the contrary, IDW, when properly implemented, also averages repeated values at common locations. This drops directly out of the formula: no special treatment is necessary to handle this case. (Which is why it is baffling that some software cannot deal with it.) Kriging, on the other hand, usually does not use the average value (except when the variogram is a pure nugget). That is why Kriging requires a preliminary step to aggregate clusters of very close data points. –  whuber Mar 20 at 14:42
    
I agree that a proper way to process coincident points is important, and averaging is then the safest solution. However, by construction, IDW only accounts for the distance while kriging also account for the relative position of the observations. Kriging indeed "naturally" compensate for the effect of clustering, assigning more weight to the isolated points than to the points from a cluster. IDW does not, so you could have markedly different results if you have 2 coincident points or two "nearly coincident" points. –  radouxju Mar 26 at 20:44
1  
for my last statement, what I tried to say is that you have different results if you take the average of the two coincident points (then considered as a single point if I understand your proposed implementation) compared with two nearly coincident points that are counted twice. I am not talking about discontinuties of the prediction near the points, but about the different global prediction that you could have depending on the decision to merge points within a given tolerance. With kriging this decision would have less consequences. –  radouxju Mar 26 at 21:31

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.