/****************************************************************************** * $Id: shpdump.c,v 1.19 2016-12-05 12:44:05 erouault Exp $ * * Project: Shapelib * Purpose: Sample application for dumping contents of a shapefile to * the terminal in human readable form. * Author: Frank Warmerdam, warmerdam@pobox.com * ****************************************************************************** * Copyright (c) 1999, Frank Warmerdam * * This software is available under the following "MIT Style" license, * or at the option of the licensee under the LGPL (see COPYING). This * option is discussed in more detail in shapelib.html. * * -- * * 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: shpdump.c,v $ * Revision 1.19 2016-12-05 12:44:05 erouault * * Major overhaul of Makefile build system to use autoconf/automake. * * * Warning fixes in contrib/ * * Revision 1.18 2011-07-24 03:05:14 fwarmerdam * use %.15g for formatting coordiantes in shpdump * * Revision 1.17 2010-07-01 07:33:04 fwarmerdam * do not crash in shpdump if null object returned * * Revision 1.16 2010-07-01 07:27:13 fwarmerdam * white space formatting adjustments * * Revision 1.15 2006-01-26 15:07:32 fwarmerdam * add bMeasureIsUsed flag from Craig Bruce: Bug 1249 * * Revision 1.14 2005/02/11 17:17:46 fwarmerdam * added panPartStart[0] validation * * Revision 1.13 2004/09/26 20:09:35 fwarmerdam * avoid rcsid warnings * * Revision 1.12 2004/01/27 18:05:35 fwarmerdam * Added the -ho (header only) switch. * * Revision 1.11 2004/01/09 16:39:49 fwarmerdam * include standard include files * * Revision 1.10 2002/04/10 16:59:29 warmerda * added -validate switch * * Revision 1.9 2002/01/15 14:36:07 warmerda * updated email address * * Revision 1.8 2000/07/07 13:39:45 warmerda * removed unused variables, and added system include files * * Revision 1.7 1999/11/05 14:12:04 warmerda * updated license terms * * Revision 1.6 1998/12/03 15:48:48 warmerda * Added report of shapefile type, and total number of shapes. * * Revision 1.5 1998/11/09 20:57:36 warmerda * use SHPObject. * * Revision 1.4 1995/10/21 03:14:49 warmerda * Changed to use binary file access. * * Revision 1.3 1995/08/23 02:25:25 warmerda * Added support for bounds. * * Revision 1.2 1995/08/04 03:18:11 warmerda * Added header. * */ #include #include #include "shapefil.h" SHP_CVSID("$Id: shpdump.c,v 1.19 2016-12-05 12:44:05 erouault Exp $") int main( int argc, char ** argv ) { SHPHandle hSHP; int nShapeType, nEntities, i, iPart, bValidate = 0,nInvalidCount=0; int bHeaderOnly = 0; const char *pszPlus; double adfMinBound[4], adfMaxBound[4]; int nPrecision = 15; if( argc > 1 && strcmp(argv[1],"-validate") == 0 ) { bValidate = 1; argv++; argc--; } if( argc > 1 && strcmp(argv[1],"-ho") == 0 ) { bHeaderOnly = 1; argv++; argc--; } if( argc > 2 && strcmp(argv[1],"-precision") == 0 ) { nPrecision = atoi(argv[2]); argv+=2; argc-=2; } /* -------------------------------------------------------------------- */ /* Display a usage message. */ /* -------------------------------------------------------------------- */ if( argc != 2 ) { printf( "shpdump [-validate] [-ho] [-precision number] shp_file\n" ); exit( 1 ); } /* -------------------------------------------------------------------- */ /* Open the passed shapefile. */ /* -------------------------------------------------------------------- */ hSHP = SHPOpen( argv[1], "rb" ); if( hSHP == NULL ) { printf( "Unable to open:%s\n", argv[1] ); exit( 1 ); } /* -------------------------------------------------------------------- */ /* Print out the file bounds. */ /* -------------------------------------------------------------------- */ SHPGetInfo( hSHP, &nEntities, &nShapeType, adfMinBound, adfMaxBound ); printf( "Shapefile Type: %s # of Shapes: %d\n\n", SHPTypeName( nShapeType ), nEntities ); printf( "File Bounds: (%.*g,%.*g,%.*g,%.*g)\n" " to (%.*g,%.*g,%.*g,%.*g)\n", nPrecision, adfMinBound[0], nPrecision, adfMinBound[1], nPrecision, adfMinBound[2], nPrecision, adfMinBound[3], nPrecision, adfMaxBound[0], nPrecision, adfMaxBound[1], nPrecision, adfMaxBound[2], nPrecision, adfMaxBound[3] ); /* -------------------------------------------------------------------- */ /* Skim over the list of shapes, printing all the vertices. */ /* -------------------------------------------------------------------- */ for( i = 0; i < nEntities && !bHeaderOnly; i++ ) { int j; SHPObject *psShape; psShape = SHPReadObject( hSHP, i ); if( psShape == NULL ) { fprintf( stderr, "Unable to read shape %d, terminating object reading.\n", i ); break; } if( psShape->bMeasureIsUsed ) printf( "\nShape:%d (%s) nVertices=%d, nParts=%d\n" " Bounds:(%.*g,%.*g, %.*g, %.*g)\n" " to (%.*g,%.*g, %.*g, %.*g)\n", i, SHPTypeName(psShape->nSHPType), psShape->nVertices, psShape->nParts, nPrecision, psShape->dfXMin, nPrecision, psShape->dfYMin, nPrecision, psShape->dfZMin, nPrecision, psShape->dfMMin, nPrecision, psShape->dfXMax, nPrecision, psShape->dfYMax, nPrecision, psShape->dfZMax, nPrecision, psShape->dfMMax ); else printf( "\nShape:%d (%s) nVertices=%d, nParts=%d\n" " Bounds:(%.*g,%.*g, %.*g)\n" " to (%.*g,%.*g, %.*g)\n", i, SHPTypeName(psShape->nSHPType), psShape->nVertices, psShape->nParts, nPrecision, psShape->dfXMin, nPrecision, psShape->dfYMin, nPrecision, psShape->dfZMin, nPrecision, psShape->dfXMax, nPrecision, psShape->dfYMax, nPrecision, psShape->dfZMax ); if( psShape->nParts > 0 && psShape->panPartStart[0] != 0 ) { fprintf( stderr, "panPartStart[0] = %d, not zero as expected.\n", psShape->panPartStart[0] ); } for( j = 0, iPart = 1; j < psShape->nVertices; j++ ) { const char *pszPartType = ""; if( j == 0 && psShape->nParts > 0 ) pszPartType = SHPPartTypeName( psShape->panPartType[0] ); if( iPart < psShape->nParts && psShape->panPartStart[iPart] == j ) { pszPartType = SHPPartTypeName( psShape->panPartType[iPart] ); iPart++; pszPlus = "+"; } else pszPlus = " "; if( psShape->bMeasureIsUsed ) printf(" %s (%.*g,%.*g, %.*g, %.*g) %s \n", pszPlus, nPrecision, psShape->padfX[j], nPrecision, psShape->padfY[j], nPrecision, psShape->padfZ[j], nPrecision, psShape->padfM[j], pszPartType ); else printf(" %s (%.*g,%.*g, %.*g) %s \n", pszPlus, nPrecision, psShape->padfX[j], nPrecision, psShape->padfY[j], nPrecision, psShape->padfZ[j], pszPartType ); } if( bValidate ) { int nAltered = SHPRewindObject( hSHP, psShape ); if( nAltered > 0 ) { printf( " %d rings wound in the wrong direction.\n", nAltered ); nInvalidCount++; } } SHPDestroyObject( psShape ); } SHPClose( hSHP ); if( bValidate ) { printf( "%d object has invalid ring orderings.\n", nInvalidCount ); } #ifdef USE_DBMALLOC malloc_dump(2); #endif exit( 0 ); }