I have created a hexagonal grid shapefile using MMQGIS and want to use this for hexagonal binning.
I have another point shapefile which contains patient records for three different types of hospital activity (let's call them A,B and C) for two different years.
What would be the best way of quickly generating point counts (via Points in Polygon) for each subset of the activity shapefile so that I ended up with 9 additional columns in my hexagonal grid shapefile?
- Year 1, activity A
- Year 2, activity A
Difference between years 1 and 2 for activity A
Year 1, activity B
- Year 2, activity B
Difference between years 1 and 2 for activity B
Year 1, activity C
- Year 2, activity C
- Difference between years 1 and 2 for activity C
Is it possible to script this using pyqgis? I've started reading up on pyqgis but I'm a long way off being able to tell if this is possible.
If it is possible, are there any tutorials out there that would cover this sort of processing or at least point me in the right direction to start figuring it out?
I really want to learn scripting in QGIS and this feels like something which would lend itself to getting started.