#!/usr/bin/perl
      use mapscript;
      use Getopt::Long;
      %shptypes = ( '1' => 'point',
                 '3' => 'arc',
                 '5' => 'polygon',
                 '8' => 'multipoint'
               );
      %typetypes = ( '0' => 'String',
                  '1' => 'Integer',
                  '2' => 'Double',
                  '3' => 'Invalid'
                );
      &GetOptions("file=s", \$file);
      if(!$file) {
        print "Syntax: shpinfo.pl -file=[filename]\n";
        exit 0;
      }
      $shapefile = new shapefileObj($file, -1) or die "Unable to open shapefile $file.";
      print "Shapefile $file:\n\n";
      print "\ttype: ". $shptypes{$shapefile->{type}} ."\n";
      print "\tnumber of features: ". $shapefile->{numshapes} ."\n";
      printf "\tbounds: (%f,%f) (%f,%f)\n", $shapefile->{bounds}->{minx}, $shapefile->{bounds}->{miny}, $shapefile->{bounds}->{maxx}, $shapefile->{bounds}->{maxy};
      #
      $table = mapscript::msDBFOpen("$file.dbf","r+b");
      print "\nDBFInfo table $file.dbf:\n\n";
      print "\tnumber of records: ". $table->{nRecords} ."\n";
      my $nofields = $table->{nFields};
      print "\tnumber of fields: ". $nofields ."\n\n";
      print "\tName             Type    Length Decimals\n";
      print "\t---------------- ------- ------ --------\n";
      for ($i=0; $i<$nofields; $i++) {
        printf "\t%-16s %-7s %6d %8d\n", $table->getFieldName($i), $typetypes{$table->getFieldType($i)}, $table->getFieldWidth($i); $table->getFieldDecimals($i)
      }

      for($j=0; $j<$numpoints; $j++) {
        $point = $line->get($j);
      }

      while ( $label = $map->nextLabel()) {
        if ( $label->{status} ) {
          print $label->{string} . "\n";
          $shape = $label->{poly};
          for ( $i=0; $i<$shape->{numlines}; i++ ) {
            $part = $shape->get($i);
            for ( $j=0; $j<$part->{numpoints}; $j++) {
              $point = $part->get($j);
              print $point->{x} . "," . $point->{y} . " ";
            }
            print "\n";
          }
        }
      }

      my $distance = $point->distanceToPoint($garagepoint);
    

      #
      # Open the map.
      my $map = new mapObj('your.map');
      #
      # Create the layer object to query.
      my $layer = new layerObj($map);
      #
      # Create a new shapefile for the selection set.
      my $newshpf = new shapefileObj("$newshpfname", 5);
      #
      # Set index value for selected shape.
      my $poly = 400;
      #
      # Create a shapefile object for getting queried shape.
      my $shpf = new shapefileObj("data_shapefile", -1);
      #
      # Create shape object for storing queried shape.
      my $shp = new shapeObj(-1);
      #
      # Retrieve shape into shape object.
      $shpf->get($poly, $shp);
      #
      # Put shape into new shapefile.
      $newshpf->add($shp);
      #
      # Create dbf to go with new.
      @toucherror = system("touch $newshpfname.dbf");
      #
      # Get the extent of the new shapefile.
      my $newrect = $newshpf->{bounds};
      my $newminx = $newrect->{minx};
      my $newminy = $newrect->{miny};
      my $newmaxx = $newrect->{maxx};
      my $newmaxy = $newrect->{maxy};
      #
      # Close new shapefile.
      undef $newshpf;
    

      #
      # Create shape object.
      $shp = new shapeObj(-1);
      #
      # Retrieve shape into shape object.
      $shpfile->get(1,$shp);
      #
      # How many lines in shape.
      my $num_lines = $shp->{numlines};
      #
      # Loop through each line.
      for ( $line_num=0; $line_num<$num_lines; $line_num++ ) {
        #
        # Get line.
        my $line = new lineObj();
        $line = $ishp->get($line_num);
        print "Got Line #$line_num\n";
        #
        # How many points in line.
        my $num_points = $line->{numpoints};
        #
        # Loop through each point.
        for ( $point_num=0; $point_num<$num_points; $point_num++ ) {
          #
          # Get the point.
          my $point = new pointObj();
          $point = $line->get($point_num);
          print "Got Point #$point_num";
          $px = $point->{x};
          $py = $point->{y};
          print " X=$px Y=$py\n";
        }
      }
      print "\n";
    

      for($i=0; $i<$numshapes; i++) {
        $shapefile->get($i, $shape);
        ...
      }

      #
      # Save the current map.
      $map->save("tobepdf.map");
      #
      # Convert the map to pdf.
      @topdferror = system("shp2pdf -m tobepdf.map -o map.pdf");