I am using ArcGIS 10.2. I have a vector layer with points, which represent houses and a raster layer with values of 0 or 1, according to a certain threshold. In a buffer of 5 km around each of these houses I want to calculate the average number of raster points above a certain thresholds (value=1), weighted according to the squared distance from the house.
Without the weighting, this would be quite simple: I would itteratively make buffers around the house (to have buffers which do not dissolve) and then use the 'zonal statistics as table' command. Now, the weighting is complicating it.
I tried the following procedure (itterating house by house in a model, see below): first rasterize the houses (point to raster) then calculate the euclidian distance, then use raster calculator to take the square of the distance and to divide my raster with thresholds by the raster with squared distances. In the meantime I make buffers around the houses, rasterize these buffers and then finally run a zonal statistics as table command to get the statistics of the output of my second rasterization within the 5000 km. Finally I use the append command to make one single table for all the houses.
Now, I have two questions:
1) Is this a good approach or am I fundamentally wrong somewhere?
2) The analysis seems to work fine, but I get strange results for some houses (eg a maximum value of 5 for one house). Moreover, if I follow each of the previous steps without a model (manually) for some houses, I get different results. I don't understand why.
Could anyone help me?