shp 0.0.1

data/ext/shp/shapelib/contrib/shpcentrd.c

@@ -0,0 +1,159 @@
+/******************************************************************************
+ * Copyright (c) 1999, Carl Anderson
+ *
+ * this code is based in part on the earlier work of Frank Warmerdam
+ *
+ * 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.
+ ******************************************************************************
+ *
+ * shpcentrd.c  - compute XY centroid for complex shapes 
+ *			and create a new SHPT_PT file of then
+ * 			specifically undo compound objects but not complex ones
+ *
+ *
+ * $Log: shpcentrd.c,v $
+ * Revision 1.3  2011-07-24 03:17:46  fwarmerdam
+ * include string.h and stdlib.h where needed in contrib (#2146)
+ *
+ * Revision 1.2  1999-05-26 02:56:31  candrsn
+ * updates to shpdxf, dbfinfo, port from Shapelib 1.1.5 of dbfcat and shpinfo
+ *
+ *
+ * 
+ */
+
+
+/* the centroid is defined as
+ *
+ *      Cx = sum (x dArea ) / Total Area
+ *  and
+ *      Cy = sum (y dArea ) / Total Area
+ */           
+
+#include <stdlib.h>
+#include <string.h>
+#include "shapefil.h"
+#include "shpgeo.h"
+
+int main( int argc, char ** argv )
+
+{
+    SHPHandle	old_SHP, new_SHP;
+    DBFHandle   old_DBF, new_DBF;
+    int		nShapeType, nEntities, nVertices, nParts, *panParts, i, iPart;
+    double	*padVertices, adBounds[4];
+    const char 	*pszPlus;
+    DBFFieldType  idfld_type;
+    int		idfld, nflds;
+    char	kv[257] = "";
+    char	idfldName[120] = "";
+    char	fldName[120] = "";
+    char	shpFileName[120] = "";
+    char	dbfFileName[120] = "";
+    double	apeture[4];
+    char	*DBFRow = NULL;
+    int		Cpan[2] = { 0,0 };
+    int		byRing = 1;
+    PT		Centrd;
+    SHPObject	*psCShape, *cent_pt;
+
+
+    if( argc < 3 )
+    {
+	printf( "shpcentrd shp_file new_shp_file\n" );
+	exit( 1 );
+    }
+    
+    old_SHP = SHPOpen (argv[1], "rb" );
+    old_DBF = DBFOpen (argv[1], "rb");
+    if( old_SHP == NULL || old_DBF == NULL )
+    {
+	printf( "Unable to open old files:%s\n", argv[1] );
+	exit( 1 );
+    }
+
+    SHPGetInfo( old_SHP, &nEntities, &nShapeType, NULL, NULL );
+    new_SHP = SHPCreate ( argv[2], SHPT_POINT ); 
+    
+    new_DBF = DBFCloneEmpty (old_DBF, argv[2]);
+    if( new_SHP == NULL || new_DBF == NULL )
+    {
+	printf( "Unable to create new files:%s\n", argv[2] );
+	exit( 1 );
+    }
+
+    DBFRow = (char *) malloc ( (old_DBF->nRecordLength) + 15 );
+
+
+#ifdef 	DEBUG
+    printf ("ShpCentrd using shpgeo \n");
+#endif
+
+    for( i = 0; i < nEntities; i++ )
+    {
+	int		res ;
+ 
+	psCShape = SHPReadObject( old_SHP, i );
+
+        if ( byRing == 1 ) {
+          int 	ring;
+          for ( ring = 0; ring < psCShape->nParts; ring ++ ) {
+	    SHPObject 	*psO;
+	    psO = SHPClone ( psCShape, ring,  ring + 1 );
+
+            Centrd = SHPCentrd_2d ( psO ); 
+
+            cent_pt = SHPCreateSimpleObject ( SHPT_POINT, 1, 
+        	(double*) &(Centrd.x), (double*) &(Centrd.y), NULL ); 
+
+            SHPWriteObject ( new_SHP, -1, cent_pt ); 
+                        
+            memcpy ( DBFRow, DBFReadTuple ( old_DBF, i ),
+		 old_DBF->nRecordLength );
+            DBFWriteTuple ( new_DBF, new_DBF->nRecords, DBFRow );            
+
+            SHPDestroyObject ( cent_pt );
+
+	    SHPDestroyObject ( psO );
+           }
+
+          }
+        else {
+        
+          Centrd = SHPCentrd_2d ( psCShape ); 
+
+          cent_pt = SHPCreateSimpleObject ( SHPT_POINT, 1, 
+        	(double*) &(Centrd.x), (double*) &(Centrd.y), NULL ); 
+
+          SHPWriteObject ( new_SHP, -1, cent_pt ); 
+                        
+          memcpy ( DBFRow, DBFReadTuple ( old_DBF, i ),
+		 old_DBF->nRecordLength );
+          DBFWriteTuple ( new_DBF, new_DBF->nRecords, DBFRow );            
+
+          SHPDestroyObject ( cent_pt );
+         }
+    }
+
+    SHPClose( old_SHP );
+    SHPClose( new_SHP );
+    DBFClose( old_DBF );
+    DBFClose( new_DBF );
+    printf ("\n");
+}

data/ext/shp/shapelib/contrib/shpdata.c

@@ -0,0 +1,129 @@
+/******************************************************************************
+ * Copyright (c) 1999, Carl Anderson
+ *
+ * this code is based in part on the earlier work of Frank Warmerdam
+ *
+ * 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.
+ ******************************************************************************
+ *
+ * shpdata.c  - utility program for testing elements of the libraries
+ *
+ *
+ * $Log: shpdata.c,v $
+ * Revision 1.2  1999-05-26 02:56:31  candrsn
+ * updates to shpdxf, dbfinfo, port from Shapelib 1.1.5 of dbfcat and shpinfo
+ *
+ *
+ * 
+ */
+      
+
+#include "shapefil.h"
+#include "shpgeo.h"
+
+int main( int argc, char ** argv )
+
+{
+    SHPHandle	old_SHP, new_SHP;
+    DBFHandle   old_DBF, new_DBF;
+    int		nShapeType, nEntities, nVertices, nParts, *panParts, i, iPart;
+    double	*padVertices, adBounds[4];
+    const char 	*pszPlus;
+    DBFFieldType  idfld_type;
+    int		idfld, nflds;
+    char	kv[257] = "";
+    char	idfldName[120] = "";
+    char	fldName[120] = "";
+    char	shpFileName[120] = "";
+    char	dbfFileName[120] = "";
+    char	*DBFRow = NULL;
+    int		Cpan[2] = { 0,0 };
+    int		byRing = 1;
+    PT		oCentrd, ringCentrd;
+    SHPObject	*psCShape, *cent_pt;
+    double	oArea = 0.0, oLen = 0.0;
+
+    if( argc < 2 )
+    {
+	printf( "shpdata shp_file \n" );
+	exit( 1 );
+    }
+    
+    old_SHP = SHPOpen (argv[1], "rb" );
+    old_DBF = DBFOpen (argv[1], "rb");
+    if( old_SHP == NULL || old_DBF == NULL )
+    {
+	printf( "Unable to open old files:%s\n", argv[1] );
+	exit( 1 );
+    }
+
+    SHPGetInfo( old_SHP, &nEntities, &nShapeType, NULL, NULL );
+    for( i = 0; i < nEntities; i++ )
+    {
+	int		res ;
+
+	psCShape = SHPReadObject( old_SHP, i );
+
+        if ( byRing == 1 ) {
+          int 	   ring, prevStart, ringDir;
+	  double   ringArea;
+
+          prevStart = psCShape->nVertices;
+          for ( ring = (psCShape->nParts - 1); ring >= 0; ring-- ) {
+	    SHPObject 	*psO;
+	    int		j, numVtx, rStart;
+            
+            rStart = psCShape->panPartStart[ring];
+            if ( ring == (psCShape->nParts -1) )
+              { numVtx = psCShape->nVertices - rStart; }
+             else
+              { numVtx = psCShape->panPartStart[ring+1] - rStart; }
+              
+            printf ("(shpdata) Ring(%d) (%d for %d) \n", ring, rStart, numVtx);              
+	    psO = SHPClone ( psCShape, ring,  ring + 1 );
+
+            ringDir = SHPRingDir_2d ( psO, 0 );
+            ringArea = RingArea_2d (psO->nVertices,(double*) psO->padfX,
+            	 (double*) psO->padfY);
+            RingCentroid_2d ( psO->nVertices, (double*) psO->padfX, 
+     		(double*) psO->padfY, &ringCentrd, &ringArea);  
+
+            	 
+            printf ("(shpdata)  Ring %d, %f Area %d dir \n", 
+           	ring, ringArea, ringDir );
+
+	    SHPDestroyObject ( psO );
+            printf ("(shpdata) End Ring \n");
+           }  /* (ring) [0,nParts  */
+
+          }  /* by ring   */
+
+	   oArea = SHPArea_2d ( psCShape );
+	   oLen = SHPLength_2d ( psCShape ); 
+	   oCentrd = SHPCentrd_2d ( psCShape );
+           printf ("(shpdata) Part (%d) %f Area  %f length, C (%f,%f)\n",
+           	 i, oArea, oLen, oCentrd.x, oCentrd.y );
+    }
+
+    SHPClose( old_SHP );
+
+    DBFClose( old_DBF );
+
+    printf ("\n");
+}

data/ext/shp/shapelib/contrib/shpdxf.c

@@ -0,0 +1,340 @@
+/******************************************************************************
+ * Copyright (c) 1999, Carl Anderson
+ *
+ * This code is based in part on the earlier work of Frank Warmerdam
+ * 
+ * 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.
+ ******************************************************************************
+ *
+ * shp2dxf.c
+ *
+ * derived from a ESRI Avenue Script
+ * and DXF specification from AutoCad 3 (yes 1984)
+ *
+ * modifications Carl Andrson 11/96
+ * modifications Carl Andrson 3/97
+ *
+ * converted to C code 12/98
+ *
+ * requires shapelib 1.2
+ *   gcc shpdxf.c shpopen.o dbfopen.o -o shpdxf 
+ *
+ */
+
+#include <stdlib.h>
+#include <string.h>
+#include "shapefil.h"
+
+#define FLOAT_PREC "%16.5f\r\n"
+
+void dxf_hdr (x1,y1,x2,y2,df)
+double x1,y1,x2,y2;
+FILE  *df;
+{
+/* Create HEADER section */
+
+  fprintf( df, "  0\r\n");
+  fprintf( df, "SECTION\r\n");
+  fprintf( df, "  2\r\n" );
+  fprintf( df, "HEADER\r\n" );
+  fprintf( df, "  9\r\n" );
+  fprintf( df, "$EXTMAX\r\n" );
+  fprintf( df, " 10\r\n" );
+  fprintf( df, FLOAT_PREC, x2 );
+  fprintf( df, " 20\r\n" );
+  fprintf( df, FLOAT_PREC, y2 );
+  fprintf( df, "  9\r\n" );
+  fprintf( df, "$EXTMIN\r\n" );
+  fprintf( df, " 10\r\n" );
+  fprintf( df, FLOAT_PREC, x1 );
+  fprintf( df, " 20\r\n" );
+  fprintf( df, FLOAT_PREC, y1 );
+  fprintf( df, "  9\r\n" );
+  fprintf( df, "$LUPREC\r\n" );
+  fprintf( df, " 70\r\n" );
+  fprintf( df, "    14\r\n" );
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "ENDSEC\r\n" );
+
+/*  ' Create TABLES section */
+
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "SECTION\r\n" );
+  fprintf( df, "  2\r\n" );
+  fprintf( df, "TABLES\r\n" );
+/*  ' Table 1 - set up line type */
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "TABLE\r\n" );
+  fprintf( df, "  2\r\n" );
+  fprintf( df, "LTYPE\r\n" );
+  fprintf( df, " 70\r\n" );
+  fprintf( df, "2\r\n" );
+/*  ' Entry 1 of Table 1 */
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "LTYPE\r\n" );
+  fprintf( df, "  2\r\n" );
+  fprintf( df, "CONTINUOUS\r\n" );
+  fprintf( df, " 70\r\n" );
+  fprintf( df, "64\r\n" );
+  fprintf( df, "  3\r\n" );
+  fprintf( df, "Solid line\r\n" );
+  fprintf( df, " 72\r\n" );
+  fprintf( df, "65\r\n" );
+  fprintf( df, " 73\r\n" );
+  fprintf( df, "0\r\n" );
+  fprintf( df, " 40\r\n" );
+  fprintf( df, "0.0\r\n" );
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "ENDTAB\r\n" );
+ /* End of TABLES section */
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "ENDSEC\r\n" );
+
+  /* Create BLOCKS section  */
+  
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "SECTION\r\n" );
+  fprintf( df, "  2\r\n" );
+  fprintf( df, "BLOCKS\r\n" );
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "ENDSEC\r\n" );
+  fprintf( df, "  0\r\n" );
+  fprintf( df, "SECTION\r\n" );
+  fprintf( df, "  2\r\n" );
+  fprintf( df, "ENTITIES\r\n" );
+  
+}
+
+
+void
+dxf_ent_preamble (dxf_type, id, df)
+int dxf_type;
+char *id;
+FILE *df;
+{
+
+  fprintf( df, "  0\r\n" );
+
+  switch (dxf_type) {
+    case  SHPT_POLYGON:
+    case  SHPT_ARC:	  fprintf (df, "POLYLINE\r\n");
+          break;
+    default:  fprintf(df, "POINT\r\n");
+  }	   
+
+  fprintf( df, "  8\r\n", df);
+  fprintf( df, "%s\r\n", id );
+  switch ( dxf_type ) {
+    case SHPT_ARC:
+                    fprintf( df, "  6\r\n" );
+                    fprintf( df, "CONTINUOUS\r\n" );
+                    fprintf( df, " 66\r\n" );
+                    fprintf( df, "1\r\n" );
+                 break;
+    case SHPT_POLYGON:
+                    fprintf( df, "  6\r\n" );
+                    fprintf( df, "CONTINUOUS\r\n" );
+                    fprintf( df, " 66\r\n" );
+                    fprintf( df, "1\r\n" );
+                    fprintf( df, " 70\r\n");
+		    fprintf (df, "1\r\n");
+    default:     break;
+  } 
+
+}
+
+void
+dxf_ent (id, x, y, z, dxf_type, df)
+char *id;
+double x,y,z;
+int dxf_type;
+FILE *df;
+{
+  if ((dxf_type == SHPT_ARC) || ( dxf_type == SHPT_POLYGON))  {
+    fprintf( df, "  0\r\n");
+    fprintf( df, "VERTEX\r\n");
+    fprintf( df, "  8\r\n");
+    fprintf( df, "%s\r\n", id);
+  }  
+  fprintf( df, " 10\r\n" );
+  fprintf( df, FLOAT_PREC, x );
+  fprintf( df, " 20\r\n" );
+  fprintf( df, FLOAT_PREC, y );
+  fprintf( df, " 30\r\n" );
+  if ( z != 0 ) 
+    fprintf( df, FLOAT_PREC, z );
+  else
+    fprintf( df, "0.0\r\n" );	       
+}
+
+
+void
+dxf_ent_postamble (dxf_type, df)
+int dxf_type;
+FILE *df;
+{
+  if ((dxf_type == SHPT_ARC) || ( dxf_type == SHPT_POLYGON)) 
+    fprintf( df, "  0\r\nSEQEND\r\n  8\r\n0\r\n", df);
+}
+
+
+int
+main (int argc, char **argv)
+
+{
+    char shpFileName[80] = "", dbfFileName[80] = "";
+    char dxfFileName[80] = "";
+    char idfldName[15];
+    char zfldName[6] = "ELEV";
+    char fldName[15];
+    char id[255];
+    double elev;
+    int   parts, *panParts, nParts, nVertices;
+    FILE	*dxf;
+    SHPHandle shp;
+    DBFHandle dbf;
+    DBFFieldType  idfld_type;
+    double adfBoundsMin[4], adfBoundsMax[4];
+    int	vrtx, shp_type, shp_numrec, zfld, idfld, nflds, recNum, part;
+    unsigned int MaxElem = -1;
+ 
+    if ( argc < 2 )  {
+        printf ("usage: shpdxf shapefile {idfield}\r\n");
+        exit (-1);
+    }
+   
+    strcpy (shpFileName,argv[1]);
+    strncpy (dbfFileName, shpFileName, strlen(shpFileName)-3);
+    strcat (dbfFileName,"dbf"); 
+  
+    strncpy (dxfFileName, shpFileName,strlen(shpFileName)-3);
+    strcat( dxfFileName, "dxf");
+  
+    shp = SHPOpen (shpFileName, "rb");
+    dbf = DBFOpen (dbfFileName, "rb");
+    dxf = fopen( dxfFileName, "w");
+
+    printf("Starting conversion %s(%s) -> %s\r\n",
+           shpFileName,dbfFileName,dxfFileName);
+  
+    SHPGetInfo (shp, &shp_numrec, &shp_type, adfBoundsMin, adfBoundsMax );
+    printf ("file has %d objects\r\n", shp_numrec);
+    dxf_hdr(adfBoundsMin[0], adfBoundsMin[1], adfBoundsMax[0], adfBoundsMax[1],
+            dxf);
+      
+/* Before proceeding, allow the user to specify the ID field to use or default to the record number.... */
+
+    if ( argc > 3 )  MaxElem = atoi(argv[3]);
+  
+    nflds = DBFGetFieldCount(dbf);
+    if ( argc > 2 ) {
+        strcpy (idfldName, argv[2]);
+        for ( idfld=0; idfld < nflds; idfld++ ) {
+            idfld_type = DBFGetFieldInfo( dbf, idfld, fldName, NULL, NULL);
+            if (!strcmp (idfldName, fldName ))
+                break;
+        }
+        if ( idfld >= nflds ) {
+            printf ("Id field %s not found, using default\r\n",idfldName); 
+            idfld = -1;  
+        } else
+            printf ("proceeding with field %s for LayerNames\r\n",fldName);
+    }
+    else
+        idfld = -1;  
+    
+    for ( zfld=0; zfld < nflds; zfld++ ) {
+        DBFGetFieldInfo( dbf, zfld, fldName, NULL, NULL);
+        if (!strcmp (zfldName, fldName  ))
+            break;
+    }
+    if ( zfld >= nflds ) 
+        zfld = -1;
+//  printf ("proceeding with id = %d, elevation = %d\r\n",idfld, zfld);
+  
+/* Proceed to process data..... */
+  
+    for ( recNum = 0; (recNum < shp_numrec) && (recNum < MaxElem); recNum++)  {
+
+        SHPObject	*shape;
+
+        if ( idfld >= 0 )
+            switch (idfld_type) {
+              case FTString:  sprintf (id, "lvl_%s",DBFReadStringAttribute ( dbf, recNum, idfld ));
+                break;
+              default:    sprintf(id, "%-20.0lf", DBFReadDoubleAttribute (dbf, recNum, idfld));
+            }
+        else
+            sprintf (id,"lvl_%-20d",(recNum +1 )); 
+
+    
+        if ( zfld >= 0 ) 
+            elev = 0;
+        else  
+            elev = DBFReadDoubleAttribute ( dbf, recNum, zfld );
+  
+#ifdef DEBUG
+        printf("\r\nworking on obj %d", recNum);
+#endif
+
+        shape = SHPReadObject( shp, recNum );
+   
+        nVertices = shape->nVertices;
+        nParts = shape->nParts;
+        panParts = shape->panPartStart;
+        part = 0;
+        for (vrtx=0; vrtx < nVertices; vrtx ++ ) { 
+#ifdef DEBUG
+        printf("\rworking on part %d, vertex %d", part,vrtx);     
+#endif 
+            if ( panParts[part] == vrtx ) {
+#ifdef DEBUG
+       printf ("object preamble\r\n");
+#endif
+                dxf_ent_preamble (shp_type, id, dxf);
+            }
+     
+            dxf_ent (id, shape->padfX[vrtx], shape->padfY[vrtx],
+                     elev, shp_type, dxf);
+ 
+            if ((panParts[part] == (vrtx + 1))|| (vrtx == (nVertices -1)) ) {
+                dxf_ent_postamble (shp_type, dxf);
+                part ++;
+            }
+        }
+        SHPDestroyObject( shape );
+    }
+
+/* close out DXF file */
+    fprintf( dxf, "0\r\n" );
+    fprintf( dxf, "ENDSEC\r\n" );
+    fprintf( dxf, "0\r\n" );
+    fprintf( dxf, "EOF\r\n" );
+
+
+    SHPClose (shp);
+    DBFClose (dbf);
+    fclose (dxf);
+}
+
+
+
+
+
+