text: #!/usr/bin/env python
  # -*- coding: utf-8 -*-
  ###############################################################################
  # $Id$
  #
  # Project:  GDAL/OGR Test Suite
  # Purpose:  KML Driver testing.
  # Author:   Matuesz Loskot <mateusz@loskot.net>
  #
  ###############################################################################
  # Copyright (c) 2007, Matuesz Loskot <mateusz@loskot.net>
  # Copyright (c) 2008-2014, Even Rouault <even dot rouault at mines-paris dot org>
  #
  # Permission is hereby granted, free of charge, to any person obtaining a
  # copy of this software and associated documentation files (the "Software"),
  # to deal in the Software without restriction, including without limitation
  # the rights to use, copy, modify, merge, publish, distribute, sublicense,
  # and/or sell copies of the Software, and to permit persons to whom the
  # Software is furnished to do so, subject to the following conditions:
  #
  # The above copyright notice and this permission notice shall be included
  # in all copies or substantial portions of the Software.
  #
  # THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
  # OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  # FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
  # THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  # LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
  # FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
  # DEALINGS IN THE SOFTWARE.
  ###############################################################################
  
  import os
  import sys
  
  sys.path.append('../pymod')
  
  import gdaltest
  import ogrtest
  from osgeo import ogr
  from osgeo import osr
  from osgeo import gdal
  
  ###############################################################################
  # Test basic open operation for KML datastore.
  #
  
  
  def ogr_kml_datastore():
  
      ogrtest.kml_ds = None
      ogrtest.have_read_kml = 0
      ogrtest.kml_drv = None
      ogrtest.libkml_drv = None
  
      ogrtest.kml_drv = ogr.GetDriverByName('KML')
      if ogrtest.kml_drv is None:
          return 'skip'
  
      ogrtest.libkml_drv = ogr.GetDriverByName('LIBKML')
      # Unregister LIBKML driver if present as it's behaviour is not identical
      # to old KML driver
      if ogrtest.libkml_drv is not None:
          print('Unregister LIBKML driver')
          ogrtest.libkml_drv.Deregister()
  
      ogrtest.kml_ds = ogr.Open('data/samples.kml')
      if ogrtest.kml_ds is not None:
          ogrtest.have_read_kml = 1
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds.GetLayerCount() != 6:
          gdaltest.post_reason('wrong number of layers')
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading attributes for first layer (point).
  #
  
  
  def ogr_kml_attributes_1():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is None:
          gdaltest.post_reason('kml_ds is none')
          return 'fail'
  
      lyr = ogrtest.kml_ds.GetLayerByName('Placemarks')
      feat = lyr.GetNextFeature()
  
      if feat.GetField('Name') != 'Simple placemark':
          gdaltest.post_reason('Wrong name field value')
          return 'fail'
  
      if feat.GetField('description')[:23] != 'Attached to the ground.':
          gdaltest.post_reason('Wrong description field value')
          print('got: ', feat.GetField('description')[:23])
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      if feat.GetField('Name') != 'Floating placemark':
          gdaltest.post_reason('Wrong name field value')
          return 'fail'
  
      if feat.GetField('description')[:25] != 'Floats a defined distance':
          gdaltest.post_reason('Wrong description field value')
          print('got: ', feat.GetField('description')[:25])
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      if feat.GetField('Name') != 'Extruded placemark':
          gdaltest.post_reason('Wrong name field value')
          return 'fail'
  
      if feat.GetField('description') != 'Tethered to the ground by a customizable \"tail\"':
          gdaltest.post_reason('Wrong description field value')
          print('got: ', feat.GetField('description'))
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading attributes for another layer (point).
  #
  
  
  def ogr_kml_attributes_2():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is None:
          gdaltest.post_reason('kml_ds is none')
          return 'fail'
  
      lyr = ogrtest.kml_ds.GetLayerByName('Highlighted Icon')
      feat = lyr.GetNextFeature()
  
      if feat.GetField('Name') != 'Roll over this icon':
          gdaltest.post_reason('Wrong name field value')
          return 'fail'
  
      if feat.GetField('description') != '':
          gdaltest.post_reason('Wrong description field value')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is not None:
          gdaltest.post_reason('unexpected feature found.')
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading attributes for another layer (linestring).
  #
  
  
  def ogr_kml_attributes_3():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is None:
          gdaltest.post_reason('kml_ds is none')
          return 'fail'
  
      lyr = ogrtest.kml_ds.GetLayerByName('Paths')
      feat = lyr.GetNextFeature()
  
      if feat.GetField('Name') != 'Tessellated':
          gdaltest.post_reason('Wrong name field value')
          return 'fail'
  
      if feat.GetField('description') != 'If the <tessellate> tag has a value of 1, the line will contour to the underlying terrain':
          gdaltest.post_reason('Wrong description field value')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      if feat.GetField('Name') != 'Untessellated':
          gdaltest.post_reason('Wrong name field value')
          return 'fail'
  
      if feat.GetField('description') != 'If the <tessellate> tag has a value of 0, the line follow a simple straight-line path from point to point':
          gdaltest.post_reason('Wrong description field value')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading attributes for another layer (polygon).
  #
  
  
  def ogr_kml_attributes_4():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is None:
          gdaltest.post_reason('kml_ds is none')
          return 'fail'
  
      lyr = ogrtest.kml_ds.GetLayerByName('Google Campus')
      feat = lyr.GetNextFeature()
  
      i = 40
      while feat is not None:
          name = 'Building %d' % i
          if feat.GetField('Name') != name:
              gdaltest.post_reason('Wrong name field value')
              print('Got: "%s"' % feat.GetField('name'))
              return 'fail'
  
          if feat.GetField('description') != '':
              gdaltest.post_reason('Wrong description field value')
              return 'fail'
  
          i = i + 1
          feat = lyr.GetNextFeature()
  
      return 'success'
  
  ###############################################################################
  # Test reading of KML point geometry
  #
  
  
  def ogr_kml_point_read():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is None:
          gdaltest.post_reason('kml_ds is none')
          return 'fail'
  
      lyr = ogrtest.kml_ds.GetLayerByName('Placemarks')
      lyr.ResetReading()
      feat = lyr.GetNextFeature()
  
      wkt = 'POINT(-122.0822035425683 37.42228990140251)'
  
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      wkt = 'POINT(-122.084075 37.4220033612141 50)'
  
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      wkt = 'POINT(-122.0857667006183 37.42156927867553 50)'
  
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading of KML linestring geometry
  #
  
  
  def ogr_kml_linestring_read():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is None:
          gdaltest.post_reason('kml_ds is none')
          return 'fail'
  
      lyr = ogrtest.kml_ds.GetLayerByName('Paths')
      lyr.ResetReading()
      feat = lyr.GetNextFeature()
  
      wkt = 'LINESTRING (-112.081423783034495 36.106778704771372 0, -112.087026775269294 36.0905099328766 0)'
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      wkt = 'LINESTRING (-112.080622229594994 36.106734600079953 0,-112.085242575314993 36.090495986124218 0)'
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      wkt = 'LINESTRING (-112.265654928602004 36.094476726025462 2357,-112.266038452823807 36.093426088386707 2357,-112.266813901345301 36.092510587768807 2357,-112.267782683444494 36.091898273579957 2357,-112.268855751095202 36.091313794118697 2357,-112.269481071721899 36.090367720752099 2357,-112.269526855561097 36.089321714872852 2357,-112.269014456727604 36.088509160604723 2357,-112.268152881533894 36.087538135979557 2357,-112.2670588176031 36.086826852625677 2357,-112.265737458732104 36.086463123013033 2357)'
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading of KML polygon geometry
  #
  
  
  def ogr_kml_polygon_read():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is None:
          gdaltest.post_reason('kml_ds is none')
          return 'fail'
  
      lyr = ogrtest.kml_ds.GetLayerByName('Google Campus')
      lyr.ResetReading()
      feat = lyr.GetNextFeature()
  
      wkt = 'POLYGON ((-122.084893845961204 37.422571240447859 17,-122.084958097919795 37.422119226268563 17,-122.084746957304702 37.42207183952619 17,-122.084572538096197 37.422090067296757 17,-122.084595488672306 37.422159327008949 17,-122.0838521118269 37.422272785643713 17,-122.083792243334997 37.422035391120843 17,-122.0835076656616 37.422090069571063 17,-122.083470946415204 37.422009873951609 17,-122.083122108574798 37.422104649494599 17,-122.082924737457205 37.422265039903863 17,-122.082933916938501 37.422312428430942 17,-122.083383735973698 37.422250460876178 17,-122.083360785424802 37.422341592287452 17,-122.083420455164202 37.42237075460644 17,-122.083659133885007 37.422512920110009 17,-122.083975843895203 37.422658730937812 17,-122.084237474333094 37.422651439725207 17,-122.0845036949503 37.422651438643499 17,-122.0848020460801 37.422611339163147 17,-122.084788275051494 37.422563950551208 17,-122.084893845961204 37.422571240447859 17))'
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      wkt = 'POLYGON ((-122.085741277148301 37.422270331552568 17,-122.085816976848093 37.422314088323461 17,-122.085852582875006 37.422303374697442 17,-122.085879994563896 37.422256861387893 17,-122.085886010140896 37.422231107613797 17,-122.085806915728796 37.422202501738553 17,-122.085837954265301 37.42214027058678 17,-122.085673264051906 37.422086902144081 17,-122.085602292640701 37.42214885429042 17,-122.085590277843593 37.422128290487002 17,-122.085584167223701 37.422081719672462 17,-122.085485206574106 37.42210455874995 17,-122.085506726435199 37.422142679498243 17,-122.085443071291493 37.422127838461719 17,-122.085099071490404 37.42251282407603 17,-122.085676981863202 37.422818153236513 17,-122.086016227378295 37.422449188587223 17,-122.085726032700407 37.422292396042529 17,-122.085741277148301 37.422270331552568 17))'
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      wkt = 'POLYGON ((-122.085786228724203 37.421362088869692 25,-122.085731299060299 37.421369359894811 25,-122.085731299291794 37.421409349109027 25,-122.085607707367899 37.421383901665649 25,-122.085580242651602 37.42137299550869 25,-122.085218622197104 37.421372995043157 25,-122.085227776563897 37.421616565082651 25,-122.085259818934702 37.421605658944031 25,-122.085259818549901 37.421682001560001 25,-122.085236931147804 37.421700178603459 25,-122.085264395782801 37.421761979825753 25,-122.085323903274599 37.421761980139067 25,-122.085355945432397 37.421852864451999 25,-122.085410875246296 37.421889218237339 25,-122.085479537935697 37.42189285337048 25,-122.085543622981902 37.421889217975462 25,-122.085626017804202 37.421860134999257 25,-122.085937287963006 37.421860134536047 25,-122.085942871866607 37.42160898590042 25,-122.085965546986102 37.421579927591438 25,-122.085864046234093 37.421471150029568 25,-122.0858548911215 37.421405713261841 25,-122.085809116276806 37.4214057134039 25,-122.085786228724203 37.421362088869692 25))'
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat is None:
          gdaltest.post_reason('expected feature not found.')
          return 'fail'
  
      wkt = 'POLYGON ((-122.084437112828397 37.421772530030907 19,-122.084511885574599 37.421911115428962 19,-122.0850470999805 37.421787551215353 19,-122.085071991339106 37.421436630231611 19,-122.084916406231997 37.421372378221157 19,-122.084219386816699 37.421372378016258 19,-122.084219386589993 37.421476171614962 19,-122.083808641999099 37.4214613409357 19,-122.083789972856394 37.421313064107963 19,-122.083279653469802 37.421293288405927 19,-122.083260981920702 37.421392139442979 19,-122.082937362173695 37.421372363998763 19,-122.082906242566693 37.421515697788713 19,-122.082850226966499 37.421762825764652 19,-122.082943578863507 37.421767769696352 19,-122.083217411188002 37.421792485526858 19,-122.0835970430103 37.421748007445601 19,-122.083945555677104 37.421693642376027 19,-122.084007789463698 37.421762838158529 19,-122.084113587521003 37.421748011043917 19,-122.084076247378405 37.421713412923751 19,-122.084144704773905 37.421678815345693 19,-122.084144704222993 37.421817206601972 19,-122.084250333307395 37.421817070044597 19,-122.084437112828397 37.421772530030907 19))'
      if ogrtest.check_feature_geometry(feat, wkt):
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Write test
  
  
  def ogr_kml_write_1():
  
      if ogrtest.kml_drv is None:
          return 'skip'
  
      srs = osr.SpatialReference()
      srs.SetWellKnownGeogCS('WGS72')
      ds = ogr.GetDriverByName('KML').CreateDataSource('tmp/kml.kml')
      lyr = ds.CreateLayer('test_wgs72', srs=srs)
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('POINT (2 49)'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
      if dst_feat.GetGeometryRef().ExportToWkt() != 'POINT (2 49)':
          print(dst_feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('CreateFeature changed the geometry.')
          return 'fail'
  
      lyr = ds.CreateLayer('test_wgs84')
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetField('name', 'my_name')
      dst_feat.SetField('description', 'my_description')
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('POINT (2 49)'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('POINT (2 49 1)'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('LINESTRING (0 1,2 3)'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('POLYGON ((0 1,2 3,4 5,0 1),(0 1,2 3,4 5,0 1))'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('MULTIPOINT (2 49,2 49)'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('MULTILINESTRING ((0 1,2 3),(0 1,2 3))'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('MULTIPOLYGON (((0 1,2 3,4 5,0 1),(0 1,2 3,4 5,0 1)),((0 1,2 3,4 5,0 1),(0 1,2 3,4 5,0 1)))'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      dst_feat = ogr.Feature(lyr.GetLayerDefn())
      dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('GEOMETRYCOLLECTION (POINT (2 49 1),LINESTRING (0 1 0,2 3 0))'))
      if lyr.CreateFeature(dst_feat) != 0:
          gdaltest.post_reason('CreateFeature failed.')
          return 'fail'
  
      ds = None
  
      return 'success'
  
  
  ###############################################################################
  # Check previous test
  
  def ogr_kml_check_write_1():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      content = open('tmp/kml.kml').read()
      if content.find('Schema') >= 0:
          gdaltest.post_reason('Did not expect Schema tags.')
          return 'fail'
  
      ds = ogr.Open('tmp/kml.kml')
      lyr = ds.GetLayerByName('test_wgs84')
      if lyr.GetFeatureCount() != 8:
          gdaltest.post_reason('Bad feature count.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetField('name') != 'my_name':
          print(feat.GetField('name'))
          gdaltest.post_reason('Unexpected name.')
          return 'fail'
      if feat.GetField('description') != 'my_description':
          print(feat.GetField('description'))
          gdaltest.post_reason('Unexpected description.')
          return 'fail'
      if feat.GetGeometryRef().ExportToWkt() != 'POINT (2 49)':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetGeometryRef().ExportToWkt() != 'POINT (2 49 1)':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetGeometryRef().ExportToWkt() != 'LINESTRING (0 1,2 3)':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetGeometryRef().ExportToWkt() != 'POLYGON ((0 1,2 3,4 5,0 1),(0 1,2 3,4 5,0 1))':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetGeometryRef().ExportToWkt() != 'MULTIPOINT (2 49,2 49)':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetGeometryRef().ExportToWkt() != 'MULTILINESTRING ((0 1,2 3),(0 1,2 3))':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetGeometryRef().ExportToWkt() != 'MULTIPOLYGON (((0 1,2 3,4 5,0 1),(0 1,2 3,4 5,0 1)),((0 1,2 3,4 5,0 1),(0 1,2 3,4 5,0 1)))':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      feat = lyr.GetNextFeature()
      if feat.GetGeometryRef().ExportToWkt() != 'GEOMETRYCOLLECTION (POINT (2 49 1),LINESTRING (0 1 0,2 3 0))':
          print(feat.GetGeometryRef().ExportToWkt())
          gdaltest.post_reason('Unexpected geometry.')
          return 'fail'
  
      return 'success'
  
  
  ###############################################################################
  # Test reading attributes with XML content in them
  #
  def ogr_kml_xml_attributes():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/description_with_xml.kml')
  
      lyr = ds.GetLayer(0)
      feat = lyr.GetNextFeature()
  
      if feat.GetField('description') != 'Description<br></br><i attr="val">Interesting</i><br></br>':
          gdaltest.post_reason('Wrong description field value')
          print('got: ', feat.GetField('description'))
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading all geometry types (#3558)
  
  
  def ogr_kml_read_geometries():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/geometries.kml')
  
      lyr = ds.GetLayer(0)
      feat = lyr.GetNextFeature()
      while feat is not None:
          feat = lyr.GetNextFeature()
  
      return 'success'
  
  ###############################################################################
  # Run test_ogrsf
  
  
  def ogr_kml_test_ogrsf():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      import test_cli_utilities
      if test_cli_utilities.get_test_ogrsf_path() is None:
          return 'skip'
  
      ret = gdaltest.runexternal(test_cli_utilities.get_test_ogrsf_path() + ' --config OGR_SKIP LIBKML -ro data/samples.kml')
  
      if ret.find("using driver `KML'") == -1 or ret.find('INFO') == -1 or ret.find('ERROR') != -1:
          print(ret)
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test fix for #2772
  
  
  def ogr_kml_interleaved_writing():
  
      ds = ogr.GetDriverByName('KML').CreateDataSource('/vsimem/ogr_kml_interleaved_writing.kml')
      lyr1 = ds.CreateLayer("lyr1")
      ds.CreateLayer("lyr2")
      feat = ogr.Feature(lyr1.GetLayerDefn())
      with gdaltest.error_handler():
          ret = lyr1.CreateFeature(feat)
      ds = None
  
      gdal.Unlink('/vsimem/ogr_kml_interleaved_writing.kml')
  
      # CreateFeature() should fail
      if ret == 0:
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading KML with only Placemark
  
  
  def ogr_kml_read_placemark():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/placemark.kml')
      lyr = ds.GetLayer(0)
      feat = lyr.GetNextFeature()
      if feat is None:
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading KML without any layer
  
  
  def ogr_kml_read_empty():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/empty.kml')
      if ds.GetLayerCount() != 0:
          gdaltest.post_reason('failed')
          print(ds.GetLayerCount())
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading KML with empty layers
  
  
  def ogr_kml_read_emptylayers():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/emptylayers.kml')
      if ds.GetLayerCount() != 2:
          gdaltest.post_reason('failed')
          print(ds.GetLayerCount())
          return 'fail'
  
      if ds.GetLayer(0).GetFeatureCount() != 0:
          gdaltest.post_reason('failed')
          print(ds.GetLayer(0).GetFeatureCount())
          return 'fail'
  
      if ds.GetLayer(1).GetFeatureCount() != 0:
          gdaltest.post_reason('failed')
          print(ds.GetLayer(1).GetFeatureCount())
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  
  
  def compare_output(content, expected_content):
      content_lines = content.strip().split('\n')
      expected_lines = expected_content.strip().split('\n')
  
      if len(content_lines) != len(expected_lines):
          gdaltest.post_reason('fail')
          print(content)
          return 'fail'
      for i, content_line in enumerate(content_lines):
          if content_line.strip() != expected_lines[i].strip():
              gdaltest.post_reason('fail')
              print(content)
              return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test that we can write a schema
  
  
  def ogr_kml_write_schema():
  
      ds = ogr.GetDriverByName('KML').CreateDataSource('/vsimem/ogr_kml_write_schema.kml')
      lyr = ds.CreateLayer("lyr")
      lyr.CreateField(ogr.FieldDefn('strfield', ogr.OFTString))
      lyr.CreateField(ogr.FieldDefn('intfield', ogr.OFTInteger))
      lyr.CreateField(ogr.FieldDefn('realfield', ogr.OFTReal))
      feat = ogr.Feature(lyr.GetLayerDefn())
      feat.SetField('strfield', 'strfield_val')
      feat.SetField('intfield', '1')
      feat.SetField('realfield', '2.34')
      lyr.CreateFeature(feat)
      ds = None
  
      f = gdal.VSIFOpenL('/vsimem/ogr_kml_write_schema.kml', 'rb')
      content = gdal.VSIFReadL(1, 1000, f).decode('ascii')
      gdal.VSIFCloseL(f)
  
      gdal.Unlink('/vsimem/ogr_kml_write_schema.kml')
  
      expected_content = """<?xml version="1.0" encoding="utf-8" ?>
  <kml xmlns="http://www.opengis.net/kml/2.2">
  <Document id="root_doc">
  <Schema name="lyr" id="lyr">
      <SimpleField name="strfield" type="string"></SimpleField>
      <SimpleField name="intfield" type="int"></SimpleField>
      <SimpleField name="realfield" type="float"></SimpleField>
  </Schema>
  <Folder><name>lyr</name>
    <Placemark>
      <ExtendedData><SchemaData schemaUrl="#lyr">
          <SimpleData name="strfield">strfield_val</SimpleData>
          <SimpleData name="intfield">1</SimpleData>
          <SimpleData name="realfield">2.34</SimpleData>
      </SchemaData></ExtendedData>
    </Placemark>
  </Folder>
  </Document></kml>"""
  
      return compare_output(content, expected_content)
  
  ###############################################################################
  #
  
  
  def ogr_kml_empty_layer():
  
      ds = ogr.GetDriverByName('KML').CreateDataSource('/vsimem/ogr_kml_empty_layer.kml')
      ds.CreateLayer("empty")
      ds = None
  
      f = gdal.VSIFOpenL('/vsimem/ogr_kml_empty_layer.kml', 'rb')
      content = gdal.VSIFReadL(1, 1000, f).decode('ascii')
      gdal.VSIFCloseL(f)
  
      gdal.Unlink('/vsimem/ogr_kml_empty_layer.kml')
  
      expected_content = """<?xml version="1.0" encoding="utf-8" ?>
  <kml xmlns="http://www.opengis.net/kml/2.2">
  <Document id="root_doc">
  <Folder><name>empty</name>
  </Folder>
  </Document></kml>"""
  
      return compare_output(content, expected_content)
  
  ###############################################################################
  # Empty layer followed by regular layer
  
  
  def ogr_kml_two_layers():
  
      ds = ogr.GetDriverByName('KML').CreateDataSource('/vsimem/ogr_kml_two_layers.kml')
      ds.CreateLayer("empty")
      lyr = ds.CreateLayer("lyr")
      lyr.CreateField(ogr.FieldDefn("foo", ogr.OFTString))
      feat = ogr.Feature(lyr.GetLayerDefn())
      feat.SetField("foo", "bar")
      lyr.CreateFeature(feat)
      ds = None
  
      f = gdal.VSIFOpenL('/vsimem/ogr_kml_two_layers.kml', 'rb')
      content = gdal.VSIFReadL(1, 1000, f).decode('ascii')
      gdal.VSIFCloseL(f)
  
      gdal.Unlink('/vsimem/ogr_kml_two_layers.kml')
  
      # FIXME: the schema for lyr should be written before the first Folter for XML compliance
      expected_content = """<?xml version="1.0" encoding="utf-8" ?>
  <kml xmlns="http://www.opengis.net/kml/2.2">
  <Document id="root_doc">
  <Folder><name>empty</name>
  </Folder>
  <Folder><name>lyr</name>
    <Placemark>
      <ExtendedData><SchemaData schemaUrl="#lyr">
          <SimpleData name="foo">bar</SimpleData>
      </SchemaData></ExtendedData>
    </Placemark>
  </Folder>
  <Schema name="lyr" id="lyr">
      <SimpleField name="foo" type="string"></SimpleField>
  </Schema>
  </Document></kml>"""
  
      return compare_output(content, expected_content)
  
  ###############################################################################
  #  Cleanup
  
  
  def ogr_kml_cleanup():
  
      os.remove('tmp/kml.kml')
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      if ogrtest.kml_ds is not None:
          ogrtest.kml_ds = None
  
      # Re-register LIBKML driver if necessary
      if ogrtest.libkml_drv is not None:
          print('Re-register LIBKML driver')
          ogrtest.libkml_drv.Register()
  
      return 'success'
  
  ###############################################################################
  # Test reading KML with folder with empty subfolder and placemark
  
  
  def ogr_kml_read_folder_with_subfolder_placemark():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/folder_with_subfolder_placemark.kml')
      if ds.GetLayerCount() != 1:
          gdaltest.post_reason('failed')
          print(ds.GetLayerCount())
          return 'fail'
  
      if ds.GetLayer(0).GetFeatureCount() != 0:
          gdaltest.post_reason('failed')
          print(ds.GetLayer(0).GetFeatureCount())
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading invalid KML (#6878)
  
  
  def ogr_kml_read_truncated():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      with gdaltest.error_handler():
          ds = ogr.Open('data/truncated.kml')
      if ds is not None:
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test fix for https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=1591
  
  
  def ogr_kml_read_weird_empty_folders():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/weird_empty_folders.kml')
      if ds.GetLayerCount() != 1:
          gdaltest.post_reason('failed')
          print(ds.GetLayerCount())
          return 'fail'
  
      if ds.GetLayer(0).GetFeatureCount() != 0:
          gdaltest.post_reason('failed')
          print(ds.GetLayer(0).GetFeatureCount())
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test fix for https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=1683
  
  
  def ogr_kml_read_junk_content_after_valid_doc():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      with gdaltest.error_handler():
          ds = ogr.Open('data/junk_content_after_valid_doc.kml')
      if ds is not None:
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading KML with kml: prefix
  
  
  def ogr_kml_read_placemark_with_kml_prefix():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/placemark_with_kml_prefix.kml')
      lyr = ds.GetLayer(0)
      feat = lyr.GetNextFeature()
      if feat is None:
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading KML with duplicated folder name
  
  
  def ogr_kml_read_duplicate_folder_name():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/duplicate_folder_name.kml')
      lyr = ds.GetLayer(0)
      if lyr.GetName() != 'layer':
          gdaltest.post_reason('failure')
          print(lyr.GetName())
          return 'fail'
      lyr = ds.GetLayer(1)
      if lyr.GetName() != 'layer (#2)':
          gdaltest.post_reason('failure')
          print(lyr.GetName())
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Test reading KML with a placemark in root document, and a subfolder (#7221)
  
  
  def ogr_kml_read_placemark_in_root_and_subfolder():
  
      if not ogrtest.have_read_kml:
          return 'skip'
  
      ds = ogr.Open('data/placemark_in_root_and_subfolder.kml')
      lyr = ds.GetLayerByName('TopLevel')
      if lyr is None:
          gdaltest.post_reason('failure')
          return 'fail'
      if lyr.GetFeatureCount() != 1:
          gdaltest.post_reason('failure')
          return 'fail'
  
      lyr = ds.GetLayerByName('SubFolder1')
      if lyr is None:
          gdaltest.post_reason('failure')
          return 'fail'
      if lyr.GetFeatureCount() != 1:
          gdaltest.post_reason('failure')
          return 'fail'
  
      return 'success'
  
  ###############################################################################
  # Build tests runner
  
  
  gdaltest_list = [
      ogr_kml_datastore,
      ogr_kml_attributes_1,
      ogr_kml_attributes_2,
      ogr_kml_attributes_3,
      ogr_kml_attributes_4,
      ogr_kml_point_read,
      ogr_kml_linestring_read,
      ogr_kml_polygon_read,
      ogr_kml_write_1,
      ogr_kml_check_write_1,
      ogr_kml_xml_attributes,
      ogr_kml_read_geometries,
      ogr_kml_test_ogrsf,
      ogr_kml_interleaved_writing,
      ogr_kml_read_placemark,
      ogr_kml_read_empty,
      ogr_kml_read_emptylayers,
      ogr_kml_write_schema,
      ogr_kml_empty_layer,
      ogr_kml_two_layers,
      ogr_kml_read_folder_with_subfolder_placemark,
      ogr_kml_read_truncated,
      ogr_kml_read_weird_empty_folders,
      ogr_kml_read_junk_content_after_valid_doc,
      ogr_kml_read_placemark_with_kml_prefix,
      ogr_kml_read_duplicate_folder_name,
      ogr_kml_read_placemark_in_root_and_subfolder,
      ogr_kml_cleanup]
  
  if __name__ == '__main__':
  
      gdaltest.setup_run('ogr_kml')
  
      gdaltest.run_tests(gdaltest_list)
  
      gdaltest.summarize()
  