/ - Diff - Metacat - Ecoinformatics Redmine

« Previous | Next »

Revision 2639

Added by harris about 19 years ago

Added an attribute to the metacat shapefile called area_coef which is a function
of the areal extent of the datapackage -- so a datapackage which was related to
a global collection would have a large extent, whereas a 3-foot squared plot
would have a small extent. Also fixed a bug related to the ID field that was
causing errors in ESRI apps.

     /**
      *  '$RCSfile$'
      *  Copyright: 2000 Regents of the University of California and the
-...
     #include <stdlib.h>
     #include <string.h>
     #include <math.h>
     #include "iostream.h"
     #include <fstream>
     #include <string>
     #include <vector>
-...
     /************************************************************************
      * Method to write the points array to the shape file and the associated
      * dbf file using the shapelib-1.2.xx libraries
      * dbf file using the shapelib-1.2.ave_x libraries
+     *
      * @param nSHPType -- the type of shapefile (eg. poly vs point )
+     *
      * @param pszFilename -- the name of the shape file, sans the extension
+     *
      * @param gridx -- this is the float value of the x-location (assuming points)
      * @param min_pntx -- this is the float value of the min x-location
      * @param max_pntx -- this is the float value of the max x-location
+     *
      * @param gridy -- this is the float value of the y-location (assuming points)
      * @param min_pnty -- this is the float value of the min y-location
      * @param max_pnty -- this is the float value of the max y-location
+     *
      * @param gridz -- this is the float value of the optional p-value
      * @param pointz -- this is the float value of the optional p-value
      * (assuming points)
+     *
      * @param the number of lines in the coordinates arrays
      ************************************************************************/
     static void WritePointsArray ( int nSHPType, const char * pszFilename,
       float gridx[40000] [100], float gridy[40000] [100],float gridz[40000] [100],
       vector<string> urlvec, vector<string> docidvec, int number_data_lines ) {
     static void WritePointsArray ( int nSHPType,
       const char * pszFilename,
       float min_pntx[40000][100],
       float max_pntx[40000][100],
       float min_pnty[40000][100],
       float max_pnty[40000][100],
       float pointz[40000][100],
       vector<string> urlvec,
       vector<string> docidvec,
       int number_data_lines ) {
            SHPHandle hSHPHandle;
            SHPObject * psShape;
-...
+           }
            //add to the dbf file attributes for ID, pval(elevation), url
            if ( DBFAddField( hDBF, "ID", FTDouble, 10, 10 ) == -1 )
            if ( DBFAddField( hDBF, "ID", FTInteger, 10, 0 ) == -1 )
+           {
              printf( "DBFAddField(%s,FTDouble,%d,0) failed.\n" );
              printf( "DBFAddField(%s,FTInteger,%d,0) failed.\n" );
+           }
            if ( DBFAddField( hDBF, "pval", FTDouble, 10, 10 ) == -1 )
            if ( DBFAddField( hDBF, "area_coef", FTDouble, 10, 4 ) == -1 )
+           {
              printf( "DBFAddField(%s,FTDouble,%d,0) failed.\n" );
+           }
-...
           //cout<<" DBF Record Count --> "<<iRecord<<endl;
           int id_field = DBFGetFieldIndex( hDBF, "ID" );
           int pval_field = DBFGetFieldIndex( hDBF, "pval" );
           int url_field = DBFGetFieldIndex( hDBF, "url" );
           int docid_field = DBFGetFieldIndex( hDBF, "docid" );
           int area_field = DBFGetFieldIndex( hDBF, "area_coef" );
           //cout << "id field -> " << id_field << endl;
           //cout << "pval field -> " << pval_field << endl;
           //cout << "url field -> " << url_field << endl;
           hSHPHandle = SHPCreate( pszFilename, nSHPType );
            ii = 0;
            for ( ii = 0; ii < number_data_lines; ii++ )
+           {
              double xx = ( gridx[2] [ii] );
              double yy = ( gridy[2] [ii] );
              double zz = ( gridz[2] [ii] );
              double area_coef = sqrt( pow((max_pntx[2][ii]-min_pntx[2][ii]),2) + pow((max_pnty[2][ii]-min_pnty[2][ii]),2)  );
              double ave_x = ( (min_pntx[2][ii]+max_pntx[2][ii])/2 );
              double ave_y = ( (min_pnty[2][ii]+max_pnty[2][ii])/2 );
              double ave_z = ( pointz[2][ii] );
              //if (xx < 180 && xx > -180 && yy < 90 && yy >-90 ) {
              if (ave_x < 180 && ave_x > -180 && ave_y < 90 && ave_y >-90 ) {
              psShape = SHPCreateObject( SHPT_POINT, ii, 0, NULL, NULL, 1, &xx, &yy, &zz, NULL );
              psShape = SHPCreateObject( SHPT_POINT, ii, 0, NULL, NULL, 1, &ave_x, &ave_y, &ave_z, NULL );
              SHPWriteObject( hSHPHandle, -1, psShape );
              SHPDestroyObject( psShape );
              //do the update to the dbf file too
              iRecord = DBFGetRecordCount( hDBF );
             DBFWriteStringAttribute(hDBF, iRecord, url_field, urlvec[ii].c_str() );
             DBFWriteStringAttribute(hDBF, iRecord, docid_field, docidvec[ii].c_str() );
             DBFWriteDoubleAttribute(hDBF, iRecord, pval_field, 999 );
             DBFWriteDoubleAttribute(hDBF, iRecord, id_field, ii );
            //  }
             DBFWriteIntegerAttribute(hDBF, iRecord, id_field, int(ii) );
             DBFWriteDoubleAttribute(hDBF, iRecord, area_field, area_coef );
+            }
+           }
            //recompute the extents
            //SHPComputeExtents( psShape );
-...
       std::cout<<"  ##           harris@nceas.ucsb.edu                 ##\n";
       std::cout<<"  ##           Copyright 1999-2003                   ##\n";
       std::cout<<"  ##><><><><><><><><><><><><><><><><><><><><><><><><>##\n";
       std::cout << std::endl;
+    }
-...
+      }
       //these arrays are for the input x, y, z data
       float gridx[1000] [100];
       float gridy[1000] [100];
       float gridz[1000] [100];
       float min_pntx[1000] [100];
       float min_pnty[1000] [100];
       float max_pntx[1000] [100];
       float max_pnty[1000] [100];
       float pointz[1000] [100];
       vector<string> docidvec;
       vector<string> urlvec;
-...
       int cnt = 0;
       int number_lines = 0;
       char bufa[100];
       char buf[100];
       char buf1[100];
       char buf2[100];
       char buf_docid[100];
       char buf_xmin[100];
       char buf_xmax[100];
       char buf_ymin[100];
       char buf_ymax[100];
       char buf_url[100];
       SHPHandle hSHP;
       int nShapeType, nEntities, i, iPart;
-...
         exit( 2 );
+      }
       /* read the data into an array -- first the grid nodes */
       /* read the data into an array -- first the grid nodes
+       *
        * the data should have the following columns (from left to right)
+       *
        * nrs.239.1
        * -121.558333333333
        * -121.558333333333
        * 36.2083333333333
        * 36.2083333333333
        * http://metacat.nceas.ucsb.edu/knb/servlet/metacat?action=read&docid=nrs.239.1&qformat=knb
        */
       infile = fopen( argv[1], "r" );
       char const * out_shape = argv[2];
       while ( fgets( oneline, 199, infile ) ) {
         sscanf( oneline, "%s %s %s %s", bufa, buf, buf1, buf2 );
         sscanf( oneline, "%s %s %s %s %s %s %s", buf_docid, buf_xmin, buf_xmax, buf_ymin, buf_ymax, buf_url);
         /**** cout<<" docid -> " << bufa
           <<"\n cnt -> " << cnt
-...
           <<"\n p -> " << buf2
           <<endl; *******/
         gridx[2] [cnt] = atof( buf );
         gridy[2] [cnt] = atof( buf1 );
         gridz[2] [cnt] = atof( buf2 );
         //char const * sLine = buf2;
         urlvec.push_back(buf2);
         docidvec.push_back(bufa);
         min_pntx[2] [cnt] = atof(buf_xmin);
         max_pntx[2] [cnt] = atof(buf_xmax);
         min_pnty[2] [cnt] = atof(buf_ymin);
         max_pnty[2] [cnt] = atof(buf_ymax);
         urlvec.push_back(buf_url);
         docidvec.push_back(buf_docid);
         cnt++;
         number_lines = cnt;
+      }
       fclose( infile );
       cout<<"creating a new shape file, with "<< number_lines <<" records "<<endl;
       WritePointsArray( SHPT_POINT, out_shape, gridx, gridy, gridz,
                         urlvec, docidvec, number_lines );
       WritePointsArray( SHPT_POINT,
                         out_shape,
                         min_pntx,
                         max_pntx,
                         min_pnty,
                         max_pnty,
                         pointz,
                         urlvec,
                         docidvec,
                         number_lines );
     #ifdef USE_DBMALLOC
       malloc_dump( 2 );

          </java>
       </target>
       <target name="spatial_option" depends="init" description=" -- compiles code for the spatial option">
       <target name="spatial_option" depends="prepare" description=" -- compiles code for the spatial option">
         <echo message="Compiling the metatcat shapefile creator"/>
         <echo message="using the shapelib at: ${shapelib.dir}"/>
-...
         </exec>
         <exec dir="." executable="g++" os="Linux, Unix" failonerror="true">
          <arg line="build/metacat_shapefile.o ${shapelib.dir}/libshp.a -o bin/test.out" />
          <arg line="build/metacat_shapefile.o ${shapelib.dir}/.libs/libshp.a -o bin/metacat_shapefile" />
         </exec>
       </target>

Also available in: Unified diff

Project

General

Profile

Metacat

Revision 2639

Added by harris about 19 years ago