#!/home/jason/perl/bin/perl -w
##############################
use strict;
use Geo::ShapeFile::Point comp_includes_z => 0, comp_includes_m => 0;
use Geo::ShapeFile;

# TODO: not sure if the real shpdump does z or m first
# TODO: documentation

my $obj = new Geo::ShapeFile(shift());

print "Shapefile Type: ".$obj->shape_type_text."   # of Shapes: ".$obj->shapes;
print "\n\n";
printf("File Bounds: ( %.3f,  %.3f,%d,%d)\n",
	$obj->x_min, $obj->y_min,
	($obj->m_min || 0), ($obj->z_min || 0),
);
printf("         to  ( %.3f,  %.3f,%d,%d)\n\n",
	$obj->x_max, $obj->y_max,
	($obj->m_max || 0), ($obj->z_max || 0),
);

for (1 .. $obj->shapes) {
	my $shape = $obj->get_shp_record($_);

	printf("Shape:%d (%s)  nVertices=%i, nParts=%i\n",
		$_-1, $shape->shape_type_text,$shape->num_points,$shape->num_parts,
	);
	printf("  Bounds:( %.3f,  %.3f,%d,%d)\n",
		$shape->x_min,$obj->y_min,($obj->m_min || 0),($obj->z_min || 0),
	);
	printf("      to ( %.3f,  %.3f, %d,%d)\n",
		$obj->x_max,$obj->y_max,($obj->m_max || 0),($obj->z_max || 0),
	);
	foreach my $p (1 .. $shape->num_parts) {
		my @part = $shape->get_part($p);
		my $labeled = 0;
		for(@part) {
			printf("   %1s ( %.3f,  %.3f, %d, %d) %s\n",
				((($p > 0) && (!$labeled))?"+":""),
				$_->X,$_->Y,($_->M || 0),($_->Z || 0),
				($labeled?"":"Ring"),
			);
			
			$labeled = 1;
		}
	}
	print "\n";
}