<div dir="ltr">(Just now noticed that by default replies *do not* go to the mailing list.  Oops.  In any case, Carlos insists that I try the point sampling plugin, which I'm doing now.  Given that I'm dealing with 5 million points, this could take a while.)<br>
<br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Mark Coletti</b> <span dir="ltr"><<a href="mailto:mcoletti@gmail.com">mcoletti@gmail.com</a>></span><br>Date: Tue, Aug 12, 2014 at 4:07 PM<br>
Subject: Re: [Qgis-user] Question regarding point to grid analysis<br>To: Carlos Cerdán <<a href="mailto:sig.upagu@gmail.com">sig.upagu@gmail.com</a>><br><br><br><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">
<div class="">On Tue, Aug 12, 2014 at 3:31 PM, Carlos Cerdán <span dir="ltr"><<a href="mailto:sig.upagu@gmail.com" target="_blank">sig.upagu@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi Mark:<br><br></div>I think that "point sampling tool" plugin is the tool what you need, With it, you take the corresponding raster value into the points attribute table, so in a new column you can calculate difference between point value and sampled value and make a statistic analysis of that difference.</div>

</div></blockquote><div><br></div></div><div>There are a number of problems with this.  First, you say "take the corresponding raster value into the points attribute table," which implies going from raster to vector when I would be doing the exact inverse.  I.e., I <i>*start with*</i> vector data to then derive raster data.  Second, "sampling" has the connotation of taking a subset of values from a larger set when I want to use <i>*all the values*</i> and not a subset.</div>

<div><br></div><div>I'll make my intent more explicit:</div><div><br></div><div><ol><li>convert point layer P_1 to raster layer R_1 using P_1's attribute M_1<br></li><li>convert point layer P_2 to raster layer R_2 using P_2's attribute M_2<br>

</li><li>perform measures of similarity between R_1 and R_2<br></li></ol></div><div>Moreover these constraints must be satisfied:</div><div><br></div><div><ul><li>R_1 and R_2 have identical minimum bounding rectangles; i.e., they cover the same area<br>

</li><li>R_1 and R_2 have the same cell sizes<br></li></ul></div><div>The first challenge is noting that there exist many ways of converting a point layer to a corresponding raster.  E.g., heatmaps, MMQGIS gridify, and several Processing "Vector -> Raster" algorithms.  The second challenge is that once I decide on a means of creating a grid layer from vector that those constraints are satisfied.</div>

<div><br></div><div>I know qgis is capable of doing this task, I just need some pointers on how to proceed.</div><div> </div><div>Cheers!</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Mark</div></font></span></div>
<span class="HOEnZb"><font color="#888888">-- <br><a href="mailto:mcoletti@gmail.com" target="_blank">mcoletti@gmail.com</a><br>
<div><br></div>
</font></span></div></div>
</div><br><br clear="all"><div><br></div>-- <br><a href="mailto:mcoletti@gmail.com" target="_blank">mcoletti@gmail.com</a><br><div><br></div>
</div>