QGIS » QGIS Application

#!/usr/bin/python

from qgis.core import * 

from PyQt4.QtGui import * 

from PyQt4.QtCore import * 

from pprint import pprint 

import sys, os 

if __name__ == '__main__': 

  layerType = { QgsMapLayer.VectorLayer: "Vector", QgsMapLayer.RasterLayer: "Raster"}  

  #oswgeo4w = os.environ.get("OSGEO4W_ROOT"); 

  #print oswgeo4w 

  QgsApplication.setPrefixPath("/usr", True); 

  QgsApplication.initQgis() 

  print 'QGis initialized' 

  a = QApplication( sys.argv ) 

  # load a project 

  p = QgsProject.instance() 

  qgsInfo = QFileInfo( "wms.qgs" ) 

  isLoaded = p.read( qgsInfo ) 

  print 'IsLoaded '+str(isLoaded) 

  print 'Title: ' + p.title() 

  # Open MapRegistry 

  mlr = QgsMapLayerRegistry.instance() 

  print 

  print 'Layers count: ' + str(mlr.count()) 

  layers = mlr.mapLayers(); 

  # Display some layer information 

  #pprint(layers) 

  for layerId in layers: 

    layer = layers[layerId] 

    print 'Layer name: ' + layer.name() 

    print 'Layer abstract: ' + layer.abstract() 

    print 'Layer type: ' + layerType[layer.type()] 

  # Add a layer 

  print 

  print "Add Layer" 

  tifInfo = QFileInfo( "img.tif" ) 

  if (tifInfo.exists()):

    print "File exists!"

    rasterLayer = QgsRasterLayer(tifInfo.absoluteFilePath(), "david" ) 

    crs = QgsCoordinateReferenceSystem(); 

    rasterLayer.setCrs(crs) 

    mlr.addMapLayer(rasterLayer, True) 

    print 

    print 'Layers count: ' + str(mlr.count()) 

    layers = mlr.mapLayers(); 

    for layerId in layers: 

      layer = layers[layerId] 

      print 'Layer name: ' + layer.name() 

      print 'Layer abstract: ' + layer.abstract() 

      print 'Layer type: ' + layerType[layer.type()]

  # Save project 

  print 

  print "Save project" 

  p.write() 

  QgsApplication.exit() 

  print "Done."