mapping

Introduction to Mapping using QGIS

View the Project on GitHub dhsouthbend/mapping

«< Previous Next »>

Additional Layers and Ideas for Analysis

Now that we have visualizations of some of the types of populations within the vulnerable areas, we can add other layers and analyze the relationships between population and the data in these other layers. For this example, we will look at a layer of Hurricane Evacuation Centers in New York City. Are these centers capable of sheltering all of the population in the vulnerable areas?

Opening the Hurricane Evacuation Centers layer, we see the following:

Opening the Layer of Hurricane Evacuation Centers

We can remove the points that are outside of the study area by Clipping it with the “Clipped Census Blocks to Study Area” layer (Just like we clipped other layers before). However, when you try to do it, you might get this error:

Unmatching CRS error: Layers do not all use the same CRS. This can cause unexpected results. Do you want to continue?

When you get this error, know that whenever QGIS warns that there may be unexpected results, it means there WILL BE unexpected results. To avoid this you need to set the layer CRS to the project’s. To do it:

Note: Sometimes setting the CRS this way won’t work, then you have to “Save as…” and choose the proper projection from the Save dialog box.

Set Layer CRS button when right-clicking on the Layer

Now you can run the Clipping tool. Note that there are some centers that are within the study area, but are not in Manhattan, so it’s best to use the layer of the Blocks in the Study Area to clip the Hurricane Evacuation Centers. Once you do, rename this new layer “Evacuation Centers in Study Area”, and open the Attributes Table. You will notice that there are only eight centers for the whole population. Here are some ideas of types of anaylsis you could do:

About creating buffers

Something important that you need to know is that QGIS works in “map units”, determined by the project’s CRS (that’s why to avoid this, we used the MMQGIS plugin earlier). For geospatial operations to be done properly, all layers involved must be in the same CRS, and this CRS is ideally the same of the project. Some CRS use degrees as units, some use meters, or feet, or miles. If you are working with a CRS that uses degrees, such as EPSG:4269, a trick to do buffers properly is to save your dataset into a CRS that uses your desired unit (for example, EPSG:102718 is a proper projection of NYC that uses feet), then change your project’s CRS to the same one than the layer, and then do the buffer, and “Save as…” a layer in the original CRS of your project (in this case, EPSG:4269). Then revert your project to the original CRS and you will have the properly projected buffers.

Comparing layers

A final layer that we can overlay with our data is the one of Hurricane Inundation Zones - Worst Case. This is a layer prepared by USGS to determine what areas could potentially be affected in the case of Hurricanes (basically what we did, but with a more complex algorythm).

Overlay of Hurricane Inundation Zones - Worst Case over our Layer

Here are some ideas to use this layer:

Adding a Basemap

A useful tool for mapping, especially if you are starting from scratch or specific streets are important, is to use a basemap to add context to your points, lines, and polygons in space. Let’s learn how to add a basemap.

As you may notice, there are many basemap options, including those with bike lanes. If you want to keep your basemap, change the opacity of your census layer so that you can see it more clearly.

«< Previous Next »>