Calculating Line Lengths and Statistics (QGIS3)

Overview of the task

Given a polyline layer of railroads in North America, we will determine the total length of railroads in the United States.

Get the data

Natural Earth has a public domain railroads dataset.

Download the North America supplement zip file from the portal.

For convenience, you may directly download a copy of the dataset from the link below:

ne_10m_railroads_north_america..zip

Data Source [NATURALEARTH]

Procedure

1. Locate the downloaded ne_10m_railroads_north_america.zip file in the Browser panel and expand it. Drag the ne_10m_railroads_north_america.shp file to the canvas.

2. You will see a new layer ne_10m_railroads_north_america loaded in the Layers panel. You will see that the layer has lines representing railroads for all of North America. Now, let’s calculate the lengths of each line feature. Go to Processing ‣ Toolbox.

3. Search for and locate the Vector geometry ‣ Add geometry attributes algorithm. Double-click to launch it.

4. In the Add Geometry Attributes dialog, select ne_10m_railroads_north_america as the Input layer. The input layer’s Coordinate Reference System (CRS) is EPSG:4326 WGS84. This is a Geographic CRS with Latitude and Longitude as coordinates, WGS84 as ellipsoid and degrees as units. Because latitude and longitude don’t have a standard length, you can’t measure distances or areas accurately using planar geometry functions. Fortunately, QGIS provides a better way to compute distances using ellipsoidal geometry, which is the most accurate method for layers spanning large areas such as this. Pick Ellipsoidal as the Calculate using option. Click Run. Once the process finishes, click Close.

5. You will see a new layer Added geom info loaded in the Layers panel. This is a copy of the input layer with a new column added for distance. Right-click the Added geom info layer and select Open Attribute Table.

6. In the Attribute Table, you will see a new column called distance. This contains the length of each line feature in meters. Also note that the sov_a3 attribute which contains the country code for each feature. Close the Attribute Table window.

7. Now that we have lengths of individual railroad line segments, we can add them up to find the total length of railroads. But as the problem statement demands we need total railroad length in the United States, we must use only the segments contained within USA. We can use the country code value in the sov_a3 column to filter the layer. Right-click the Added geom info layer and select Filter.

8. In the Query Builder dialog, enter the following expression and click OK.

"sov_a3" = 'USA'

9. You will see a Filter icon appear next to the Added geom info layer in the Layers panel indicating that a filter is applied to the layer. You can also visually confirm that the layer now contains line segments only for United States. Now we are ready to calculate the sum. Click the Show statistical summary button on the Attributes Toolbar.

10. A new Statistics panel will open. Select Added geom info layer and length column.

11. You will see various statistics displayed in the panel. The unit of the statistics is the same as the units of length column – meters. Let’s change the computation to use kilometers instead. Click the Expression icon next to the fields drop-down menu in the Statistics panel.

12. Enter the following expression in the Expression Dialog that converts the length to kilometers.

length / 1000

13. The Sum value displayed is the total length of railroads in USA.