New

git-svn-id: https://svn.osgeo.org/gdal/trunk@7032 f0d54148-0727-0410-94bb-9a71ac55c965

pymod/gcps2vec.py

@@ -0,0 +1,146 @@
+#!/usr/bin/env python
+#******************************************************************************
+#  $Id$
+# 
+#  Project:  GDAL
+#  Purpose:  Convert GCPs to a point layer.
+#  Author:   Frank Warmerdam, [email protected]
+# 
+#******************************************************************************
+#  Copyright (c) 2005, Frank Warmerdam <[email protected]>
+# 
+#  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.
+#******************************************************************************
+# 
+# $Log$
+# Revision 1.1  2005/02/07 14:36:12  fwarmerdam
+# New
+#
+#
+
+import gdal
+import ogr
+import sys
+
+def Usage():
+    print 'Usage: gcp2vec.py [-of <ogr_drivername>] <raster_file> <vector_file>'
+    sys.exit(1)
+
+# =============================================================================
+# 	Mainline
+# =============================================================================
+
+format = 'GML'
+in_file = None
+out_file = None
+
+gdal.AllRegister()
+argv = gdal.GeneralCmdLineProcessor( sys.argv )
+if argv is None:
+    sys.exit( 0 )
+
+# Parse command line arguments.
+i = 1
+while i < len(argv):
+    arg = argv[i]
+
+    if arg == '-of':
+        format = argv[i]
+
+    elif in_file is None:
+        in_file = argv[i]
+
+    elif out_file is None:
+        out_file = argv[i]
+
+    else:
+        Usage()
+
+    i = i + 1
+
+if out_file is None:
+    Usage()
+
+# ----------------------------------------------------------------------------
+# Open input file, and fetch GCPs. 
+# ----------------------------------------------------------------------------
+ds = gdal.Open( in_file )
+if ds is None:
+    print 'Unable to open ', filename
+    sys.exit(1)
+
+gcp_srs = ds.GetGCPProjection()
+gcps = ds.GetGCPs()
+
+ds = None
+
+if gcps is None or len(gcps) == 0:
+    print 'No GCPs on file %s!' % in_file
+    sys.exit(1)
+
+# ----------------------------------------------------------------------------
+# Create output file, and layer.
+# ----------------------------------------------------------------------------
+
+drv = ogr.GetDriverByName( format )
+if drv is None:
+    print 'No driver named %s available.' % format
+    sys.exit(1)
+
+ds = drv.CreateDataSource( out_file )
+
+srs = None
+
+layer = ds.CreateLayer( 'gcps', srs, geom_type = ogr.wkbPoint25D )
+
+fd = ogr.FieldDefn( 'Pixel', ogr.OFTReal )
+layer.CreateField( fd )
+
+fd = ogr.FieldDefn( 'Line', ogr.OFTReal )
+layer.CreateField( fd )
+
+fd = ogr.FieldDefn( 'Id', ogr.OFTString )
+layer.CreateField( fd )
+
+fd = ogr.FieldDefn( 'Info', ogr.OFTString )
+layer.CreateField( fd )
+
+# ----------------------------------------------------------------------------
+# Write GCPs. 
+# ----------------------------------------------------------------------------
+
+for gcp in gcps:
+
+    feat = ogr.Feature( layer.GetLayerDefn() )
+
+    geom = ogr.Geometry( ogr.wkbPoint25D )
+    geom.SetPoint( 0, gcp.GCPX, gcp.GCPY, gcp.GCPZ )
+    feat.SetGeometryDirectly( geom )
+
+    feat.SetField( 'Pixel', gcp.GCPPixel )
+    feat.SetField( 'Line',  gcp.GCPLine )
+    feat.SetField( 'Id',    gcp.Id )
+    feat.SetField( 'Info',  gcp.Info )
+
+    layer.CreateFeature( feat )
+
+feat = None
+ds.Destroy()
+
+