I’m very sorry to say that I do not have any sample code for this question. The reason being, I’ve been looking for a way to do this and I’ve got no ideas. I’m not looking for specific code help, more like some general guidance on where to start.
I have lat/lon for people living around California. Based on the geographic groupings of this data (county, city, zip, etc) I’m supposed to figure out a density score. Basically, how many points per square mile, or something along those lines. I have been looking for ways to do this, and unfortunately I haven’t found anything that seems right. I have lat/lon/geographic boundary columns in my data, so I can group by the various geographic types, but I’d like to rely on the lat/lon for the density scoring.
Again, I’m sorry that I don’t have any specific code to share. Any suggestions on packages or tools for this problem would be greatly appreciated!
I’m comfortable with both R and Python, but my colleagues would prefer python driven solutions.
Below is some sample data using zip codes and lon/lat (in this order).
Zip Longitude Latitude95223 -120.045063 38.46730895223 -120.040889 38.46543695223 -120.072499 38.45420295223 -120.049251 38.46205895223 -120.041697 38.46219495223 -120.045757 38.47063796120 -119.959615 38.70396596120 -119.937276 38.74133796120 -119.9382 38.73934496120 -119.901794 38.77658496120 -119.936094 38.74186596120 -119.957587 38.70753396120 -119.93456 38.7419495646 -120.072087 38.68706195646 -120.066752 38.68409795646 -120.069591 38.68419395646 -120.071754 38.69973895646 -120.066111 38.68516495646 -120.067082 38.68388195646 -120.070923 38.69604995646 -120.068004 38.68361595646 -120.07161 38.69930995646 -120.07385 38.69071995646 -120.066131 38.68501995646 -120.071263 38.686228Advertisement
Answer
From the Pyviz Examples page for census mapping here: Example code for mapping population density: https://examples.pyviz.org/census/census.html
Creating your boundaries: The answer to this question covers this quite well, giving a list of boundary detection algorithms: https://gis.stackexchange.com/questions/5426/finding-boundary-co-ordinates-from-given-set-of-point-co-ordinates
Finding the area of your polygon(AreaOfZipcode): How to calculate the area of a polygon on the earth’s surface using python?
Calculate Zipcode Populations with an accumulator algorithm.
Then: PopDensity = PopOfZipcode/AreaOfZipcode
Define your high density/mid-density/low-density boundaries and then assign your zipcodes to each bucket.