Skip to content
Snippets Groups Projects
shpdump.c 9.54 KiB
Newer Older
  • Learn to ignore specific revisions
  • /******************************************************************************
     * $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 );
    }