ex3

Author: kgjenkins

exercise1.md

@@ -66,7 +66,7 @@ Once we're satisfied with the output, we are ready to save it as a GeoJSON file.
 
 13. Right-click the "Refactored" layer > Make Permanent...
   - Set Format = GeoJSON
-  - Always click the '...' button to specify where you want the file saved!
+  - Always click the `...` button to specify where you want the file saved!
   - Under Layer Options, set RFC7946 = YES (this will force it into WGS84, and set decimal precision at 7 digits)
 
 Now you should have an index map saved as a GeoJSON file!  If you open it in a text editor, you'll see that it starts something like this:

exercise2.md

@@ -45,7 +45,11 @@ Use the "Refactor fields" tool to rename or delete several of the fields:
 - Always click the `...` button to specify where you want the file saved!
 - Under Layer Options, set RFC7946 = YES (this will force it into WGS84, and set decimal precision at 7 digits)
 
-Now we have all the points saved to a single GeoJSON file.  What if we wanted to save each county separately, based on the "county" column?
+Now we have all the points saved to a single GeoJSON file.
+
+## 6. Save a GeoJSON for each county
+
+What if we wanted to save each county separately, based on the "county" column?
 
 - Open the "Split Vector Layer" tool
 - Input layer = "Points from table"

exercise3.md

@@ -1,7 +1,60 @@
 # Exercise 3: Create a grid index map from scratch
-- Create grid in QGIS
-- formula for grid labels
-  - label expression like `char(45-bottom+65)||(left+81)`
-  - once the formula looks good, calculate new 'label' column using the expression
-- edit form for names - can hide columns (layer properties > attributes form > drag and drop designer)
 
+In this scenario, we will create an index map from scratch, using a grid of 15-minute quads for the state of Colorado.  We'll use the Stamen Toner basemap for reference, because the state boundaries are easier to see.
+
+- Start a new project in QGIS and add the Stamen Toner basemap
+
+## 1. Create the grid
+
+A 15-minute quad is 15/60 or 1/4 of a degree, so we need to create a rectangular grid where each cell is 0.25 degrees wide and tall.  To cover all of Colorado, so we'll need to figure out the lat/lon bounds of the state.
+
+- Zoom in to Colorado (shift-drag will zoom into the rectangle you select)
+- Click the Project CRS ![EPSG:3857](https://kgjenkins.github.io/openindexmaps-workshop/image/project-crs.png) and change it to EPSG:4326 "WGS 84"
+- Open the "Create Grid" processing tool
+- Grid type = Rectangle (Polygon)
+- Grid extent -- click `...` > Select extent on canvas, then drag a rectangle covering Colorado
+- Fine tune the coordinates by rounding to the nearest whole number ("-109, -102, 37, 41")
+- Set the horizontal and vertical spacing to 0.25 degrees
+- Click "Run"
+
+## 2. Label the cells
+
+Let's suppose we have a Colorado map series that uses these quads, and each map is number A1, A2, A3 across the top row, and A1, B1, C1 down the left column, like this:
+
+![grid labels](https://kgjenkins.github.io/openindexmaps-workshop/image/ex3-grid-labels.png)
+
+Although we could just add a new column and type in the value for every cell, that would be very tedious and prone to error.  With some experimentation, we can figure out a formula for calculating the code for each cell.
+
+Use the Identify tool ![identify tool](https://kgjenkins.github.io/openindexmaps-workshop/image/identify-tool.png) to look at the attributes for one of the grid cells.  Notice that it includes values for the left, top, right, and bottom coordinates of the cell.
+
+- Open the Layer Styling Panel, and select the "Grid" layer
+- Click the "ABC" button on the left of the panel to switch to label styling
+- Change "No labels" to "Single labels"
+- Set Value = `top`
+
+QGIS dynamic label expressions can help us to figure out the formula in an iterative manner.  The complete formula is below, but here is the sequence used to figure it out.  After typing to change the `top` expression, just press "tab" or click elsewhere and the labels will update on the map.
+
+- `top*4` (to get whole numbers)
+- `165-top*4` (so that the top row is 1, and the numbers increase southwards)
+- `char(64 + 165-top*4)` (this converts the numbers to letters -- char(65) is A, char(66) is B, etc.)
+- `char(64 + 165-top*4) || (left*4)` (this appends numbers based on the longitude)
+- `char(64 + 165-top*4) || (437 + left*4)` (so that the left column is 1, and increases to the east)
+
+Once the expression is generating the correct codes, we can copy it and use it to create a new column.
+
+- Copy the value expression for the labels
+- Open the Field Calculator ![field calculator button](https://kgjenkins.github.io/openindexmaps-workshop/image/field-calculator.png)
+- Output field name = "code"
+- Output field type = "Text (string)"
+- Paste the expression into the large area on the left of the dialog
+- Check the output preview beneath
+- Click "OK" -- the adds the column, but now we are in edit mode
+- Click the pencil icon to toggle editing off (save when prompted)
+
+Now the code is just another property that could be used when we run "Refactor Fields".
+
+When in edit mode, other fields could be added and edited while viewing the attribute table.
+
+----
+
+Next: [Exercise 4](exercise4)