You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
284 lines
9.5 KiB
284 lines
9.5 KiB
/****************************************************************************** |
|
* $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 ); |
|
}
|
|
|