Changes between Initial Version and Version 1 of PerlMapScriptExamples35ex11

Timestamp:

Jan 29, 2009, 6:54:58 AM (15 years ago)

Author:

jmckenna

Comment:

--

PerlMapScriptExamples35ex11

@@ +1 @@
+= project.pl =
+{{{
+#!perl
+#!/usr/bin/perl
+
+# File: project.pl
+#
+# Purpose: This script takes a shapefile as input along with 2 Proj.4
+# projection descriptions.
+#
+# Syntax: project.pl [from projection] [to projection] [in shapefile] [out shapefile]
+#
+# Author: Stephen Lime (08/05/02)
+
+use mapscript;
+
+sub show_syntax() {
+  print "Syntax: project.pl [from projection] [to projection] [in shapefile] [out shapefile]\n";
+  exit 0;
+}
+
+&show_syntax() unless $#ARGV == 3;
+ 
+#
+# open the shapefiles
+#
+$in_shapefile = new shapefileObj($ARGV[2], -1) or die "Unable to open input shapefile.";
+$out_shapefile = new shapefileObj($ARGV[3], $in_shapefile->{type}) or die "Unable to open output shapefile.";
+ 
+#
+# make a copy of the dBase file
+#
+system("cp ". $ARGV[2] .".dbf ". $ARGV[3] .".dbf");
+ 
+#
+# create a couple of projection objects
+#
+$from_projection = new projectionObj($ARGV[0]) or die "Unable to initialize \"from\" projection.";
+$to_projection = new projectionObj($ARGV[1]) or die "Unable to initialize \"to\" projection.";
+ 
+#
+# and finally loop through the shapefile
+#
+$shape = new shapeObj($in_shapefile->{type});
+for($i=0; $i<$in_shapefile->{numshapes}; $i++) {
+  $status = $in_shapefile->get($i, $shape);
+  die "Error reading shape $i." unless $status == $mapscript::MS_SUCCESS;
+  $shape->project($from_projection, $to_projection);
+  $out_shapefile->add($shape);
+}
+ 
+undef $out_shapefile;
+exit;
+}}}
+----
+back to PerlMapScript