官术网_书友最值得收藏!

Measuring the distance along a line sample

In this recipe, we'll measure the distance along a line with multiple vertices.

Getting ready

For this recipe, we'll use a line shapefile with two features. You can download the shapefile as a .ZIP file from https://geospatialpython.googlecode.com/svn/paths.zip

Unzip the shapefile into a directory named qgis_data/shapes within your root or home directory.

How to do it...

The steps for this recipe are fairly straightforward. We'll extract the geometry from the first line feature and pass it to the measurement object, as shown here:

  1. First, we must load the QGIS constants library:
    from qgis.core import QGis
    
  2. Load the line layer:
    lyr = QgsVectorLayer("/qgis_data/shapes/paths.shp", "Route", "ogr")
    
  3. Grab the features:
    fts = lyr.getFeatures()
    
  4. Get the first feature:
    route = fts.next()
    
  5. Create the measurement object instance:
    d = QgsDistanceArea()
    
  6. Then, we must configure the QgsDistanceArea object to use the ellipsoidal mode for accurate measurements in meters:
    d.setEllipsoidalMode(True)
    
  7. Pass the line's geometry to the measureLine method:
    m = d.measureLine(route.geometry().asPolyline())
    
  8. Convert the measurement output to miles:
    d.convertMeasurement(m, QGis.Meters, QGis.NauticalMiles, False)
    

Ensure that your output looks similar to the following:

(2314126.583384674, 7)

How it works...

The QgsDistanceArea object can perform any type of measurement, based on the method you call. When you convert the measurement from meters (represented by 0) to miles (identified by the number 7), you will get a tuple with the measurement in miles and the unit identifier. The QGIS API documentation shows the values for all the unit constants

(http://qgis.org/api/classQGis.html).

主站蜘蛛池模板: 平和县| 杭州市| 深圳市| 河曲县| 桃江县| 普安县| 蓬溪县| 新蔡县| 卓资县| 盐山县| 祁东县| 原阳县| 屏边| 屯昌县| 会同县| 新邵县| 台东县| 五华县| 板桥市| 宾川县| 宿迁市| 奇台县| 吉水县| 太保市| 平邑县| 黔西县| 龙胜| 陆丰市| 汉阴县| 洞头县| 长兴县| 大渡口区| 延川县| 突泉县| 大田县| 日照市| 永福县| 左贡县| 汕头市| 神农架林区| 嵩明县|