home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #20 / NN_1992_20.iso / spool / comp / graphics / 9686 < prev    next >
Encoding:
Internet Message Format  |  1992-09-14  |  2.7 KB
  1. Xref: sparky comp.graphics:9686 comp.graphics.explorer:207 comp.graphics.visualization:1332
  2. Newsgroups: comp.graphics,comp.graphics.explorer,comp.graphics.visualization
  3. Path: sparky!uunet!munnari.oz.au!bunyip.cc.uq.oz.au!brolga!vthrc
  4. From: vthrc@brolga.cc.uq.oz.au (Danny Thomas)
  5. Subject: (Q) calculating spatial density
  6. Message-ID: <vthrc.716467347@brolga>
  7. Keywords: spatial density algorithm statistics
  8. Sender: news@bunyip.cc.uq.oz.au (USENET News System)
  9. Organization: Prentice Centre, University of Queensland
  10. Date: Mon, 14 Sep 1992 10:42:27 GMT
  11. Lines: 49
  12.  
  13.  
  14. I've got a complex 2D shape with a set of points in it and need to 
  15. generate a map of density of those points. I'm sure this topic must be 
  16. well-explored but one complication is that there are holes in the shape 
  17. too. 
  18.     The relatively small size of this data-set does however allow points 
  19. to be snapped to a fine pixel grid and areas calculated by pixel-
  20. counting. A simple approach would be to:
  21. 1. count the number of data points in a circular radius around each grid-
  22. point
  23. 2. calculate the area of holes in that radius by counting hole-pixels
  24. 3. subtract that area from that of the circle
  25. 4. divide the number of data-points by the non-hole area to give 
  26.    a measure of the density
  27.  
  28. but I would be interested to learn of "more analytical" algorithms that 
  29. don't depend on quantizing coordinates. It would also be interesting to 
  30. read about techniques that can handle arbitrary weighting functions. eg 
  31. if an inverse square weighting is applied to distances from a grid point 
  32. to each data-point, what do you do when that vector crosses over a hole? 
  33. Perhaps more accurately, how could you estimate a contribution from the 
  34. points that would exist if that hole wasn't there?
  35.  
  36.  
  37. The final step of converting a density grid into a contour map should be 
  38. solvable with other tools we have.
  39.  
  40. PS I've already received a pointer to the bivar module from avs.ncsc.org, 
  41. and apart from the fact it doesn't seem to account for holes (or edges), 
  42. we don't have AVS. We only have Macs & PCs but are planning on an SGI 
  43. Indigo XS24 soon, which means we are particularly interested in Explorer 
  44. (& Khoros?).
  45.  
  46. because of server problems I'd appreciate an email copy of any replies 
  47. posted to net-news. It might be worth reducing the number of groups any 
  48. reply is sent to.
  49.  
  50. cheers,
  51. X500: @c=AU [I'm about to update VTHRC entries, soon...]
  52.        @o=University of Queensland
  53.         @ou=Physiology and Pharmacology Department
  54.          @ou=Vision Touch and Hearing
  55.           @cn=Danny Thomas
  56.            @rfc822Mailbox=vthrc@cc.uq.oz.au
  57.                           ^^^^^^^^^^^^^^^^^  <- return address
  58.  
  59.  
  60. PS  & apologies for cross-posting, you'd think this a basic GIS function.
  61.  
  62.