Well, I just realized that I was using the wrong module for what I wanted to do. v.neighbors does exactly what I want. However, the man page fro v.kernel is still quite spare, and if I made the mistake of thinking it would make a normal density map, then I think others may do so too? I'm still not entirely sure what v.kernel is actually calculating. Since it's listed as "Generate Gaussian Kernel Density Surface" in the "generate surfaces" portion of the raster menu in the GUI, and v.neighbors is listed as "Neighborhood Points" in the 'Neighborhood Analysis" portion of the raster menu in the GUI, it seems at a casual glance that v.kernel would be the one to make density maps, and not v.neighbors, even though it is, in reality, vice versa. And there is no info in the man page of v.kernel that suggests otherwise. I hope other's read this and don't make the same mistake I did! :)<br>
<br>Cheers,<br><br>Isaac<br><br><br><div class="gmail_quote">On Tue, Jan 26, 2010 at 12:18 PM, Isaac Ullah <span dir="ltr"><<a href="mailto:isaac.ullah@asu.edu">isaac.ullah@asu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi all,<br><br>I'm running an experiment to simulate the distribution of microartifacts (less than 1mm sized artifacts) on ancient housefloors. The simulations essentially produce vector point files where each point is the location of a single microartifact. These are arranged as clusters across a rectangular space (region extents), and there are several hundred thousand points in each vector point file. To simulate archaeological recovery of these microartifacts, I am using v.kernel to create density maps of the points at different resolutions. My question here concerns the output of v.kernal. Do the numbers in an output v.kernal density map represent density over some standard unchanging unit (ie. map unit), or is it density per the current resolution? The reason I'm asking is that the actual numbers in the output maps do not seem to be changing even though the resolution is different. For example, if I set the resolution to 10cm (0.1m), I get values ranging from 0 to over 2000. When I then change the resolution to 1m, I still get values ranging from 0 to 2000. The pattern of distribution between the maps is similar, so that seems to be working. This all makes me think that v.kernel is calculating the density of points <i>per some standard unit</i> rather than just counting the number of points in one cell given the current raster resolution (set in g.region). If the later was occurring, then I would expect an output map with raster resolution of 1m to have much higher counts per cell than one with resolution of only 10cm. The man page for v.kernel does not mention anything about this (and in fact pretty mentions very little at all). So I guess my question is: Can anyone confirm for me that v.kernel output will be density per map unit, or density per current resolution settings? I just need to know so that I can quantify my results properly. And while we are on it, can anyone also explain to me what the "standard dev" option in v.kernel actually does? Like I said before, the man page is lacking in these details...<br>
<br>Any help will be greatly appreciated!<br><br>Cheers,<br clear="all"><br>-- <br>Isaac I Ullah, M.A.<br><br>Archaeology PhD Candidate,<br>ASU School of Evolution and Social Change<br><br>Research Assistant,<br>Mediterranean Landscape Dynamics Project<br>
***************************************************<br><a href="mailto:isaac.ullah@asu.edu" target="_blank">isaac.ullah@asu.edu</a><br><a href="mailto:ullah@archaeologist.com" target="_blank">ullah@archaeologist.com</a><br>
<br><a href="http://www.public.asu.edu/%7Eiullah" target="_blank">http://www.public.asu.edu/~iullah</a><br>
***************************************************<br>
</blockquote></div><br><br clear="all"><br>-- <br>Isaac I Ullah, M.A.<br><br>Archaeology PhD Candidate,<br>ASU School of Evolution and Social Change<br><br>Research Assistant,<br>Mediterranean Landscape Dynamics Project<br>
***************************************************<br><a href="mailto:isaac.ullah@asu.edu">isaac.ullah@asu.edu</a><br><a href="mailto:ullah@archaeologist.com">ullah@archaeologist.com</a><br><br><a href="http://www.public.asu.edu/~iullah">http://www.public.asu.edu/~iullah</a><br>
***************************************************<br>