Changes between Initial Version and Version 1 of ProjectingWktFromStdin

Timestamp:

Feb 16, 2010, 8:59:25 AM (14 years ago)

Author:

hobu

Comment:

--

ProjectingWktFromStdin

@@ +1 @@
+{{{
+#!rst
+
+The following script takes in some WKT on stdin, projects it to the coordinate system of your choice, and emits it on stdout.  The first argument passed assigns the coordinate system to the WKT, and the second is your output coordinate system.
+
+::
+
+  project_wkt.py EPSG:26915 EPSG:4326 < input.wkt > output.wkt
+
+
+.. code-block:: python
+
+    #!/usr/bin/env python
+
+    import sys
+
+    from osgeo import ogr, osr
+
+    in_ref_text = sys.argv[1]
+    out_ref_text = sys.argv[2]
+
+    g = ogr.CreateGeometryFromWkt(sys.stdin.read())
+
+    in_ref = osr.SpatialReference()
+    in_ref.SetFromUserInput(in_ref_text)
+
+    out_ref = osr.SpatialReference()
+    out_ref.SetFromUserInput(out_ref_text)
+
+    g.AssignSpatialReference(in_ref)
+    g.TransformTo(out_ref)
+
+    sys.stdout.write(g.ExportToWkt())
+
+}}}