GIS5100 - Module 1 - Crime Analysis
In this first module 1, we learn how to hot spot analysis using grid overlay, kernel density, and Local Moran's I. This first lab is definitely a step up from the previous classes, There were a lot of trials and errors for me to self-experiment since the directions were more vague and I had to relearn a lot of the tools we used. Belows are a quick summary of the steps I used which I hope prove to be useful to others as I personally struggled a lot with this first module.
I created a blank map called Burglaries and added data from the WashingtonDC dataset. I created a new table using “Select by Attribute” for ‘Burglary’ from the “OFFENSE” field and named it “Burglary”. I then used a spatial join to count the number of crimes per census tract. Next, I joined the generated “Burlaries_Tracts” w/ “dctract_data” using “TRACT” and “Tract Number”. Then I added a field to the “Burglaries_Tracts” named “BurlgaryRate”. Next, I calculated the burglary rate using the provided python expression. I added a graduated colors symbology which excluded the burglary rate of 125 or 1500 from the provided SQL expression. I then created a layout and added a legend, then exported the layout as a JPEG.
I created a new map called “Assault” and added the “CrimpeRep2018_DC” and “Boundary_DC” to the map. Next, I select by attributes where offense is equal to assult w/ dangerous weapon and created a new layer with the selection. I then open the Environment setting and set the extent to the Boundary_DC and mask to Boundary_DC. Next, I ran the Kernel Density tool with the provided inputs. I applied a Classify symbology to the generated assaults_kd. I excluded 0 from the layer’s data. Then, I set the upper values to the required mean plus means * X, and labeled it accordingly. I then created a AssaultRate layout and added a legend and exported the layout as a JPEG.
For the grid-based hot spot, I spatial joined the Chicago_Grid and the 2017 homicide which made a new homicide_grid layer. The homicide_grid layer now have a Join_count field which I use select by attribute to select any rows with higher than zero. I exported all the rows previously selected to a new layer call “homicide_count”. I then open the “homicide_count” and sorted the “Joint_count” from largest to smallest and then I selected the top 62 rows and then exported the 62 rows as “hom_top20per”. Next, I added a field called dissolve which use a short interger. I added the number “1” to the entire column. I then dissolved the layer using the provided settings to make “hom_top20per_diss” layer.
For the kernel density hot spot mapping, I used the provided inputs in the instruction and modified the symbology so that there are only two break values: 3 * mean and the maximum value. I also excluded zero from the symbology as well. I then reclassify the raster and used the provided inputs and exported “hom_kd_rc”. Next, I converted the rasted to a polygon and named the feature class “hom_kd_fc”. Next I selected by attribute for the gridcode value of ‘2’. I exported the selection as “hom_kd”.
For the local Moran’s I Hotspot Mapping, I spatial joined the census tract and the 2017 homicides and named the output feature class hom_morans_sj. I then added a new “Crime_Rate” field with a float data type. I use a python expression to do calculate the homicides per 1000 housing unit using the join_count and total_households. Next I ran the Cluster and Outlier Analysis spatial statistics tool and named the output hom_morans_run. After, I used SQL query to select the high-high cluster and output the selection to a new hom_morans_hh layer and dissolved it.



Comments
Post a Comment