forked from openkylin/shapelib
285 lines
9.5 KiB
C
285 lines
9.5 KiB
C
/******************************************************************************
|
|
* $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 <string.h>
|
|
#include <stdlib.h>
|
|
#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 );
|
|
}
|