From ff34a41621c4fb659d0523e4e0a5c1c57e682a51 Mon Sep 17 00:00:00 2001
From: Emmanuel Christophe <emmanuel.christophe@orfeo-toolbox.org>
Date: Thu, 26 Feb 2009 14:32:24 +0800
Subject: [PATCH] OSSIM: update to 1.7.15
---
.../ossim/elevation/ossimElevationShapeIdx.h | 41 -
.../ossim/imaging/ossimEsriShapeFileCutter.h | 200 --
.../ossim/imaging/ossimEsriShapeFileFilter.h | 270 ---
.../ossim/imaging/ossimImageShapeFileIndex.h | 61 -
.../ossim/imaging/ossimImageSourceInterface.h | 109 -
.../include/ossim/vec/ossimShapeDatabase.h | 136 --
.../include/ossim/vec/ossimShapeFile.h | 312 ---
.../otbossim/include/ossim/vec/shapefil.h | 540 -----
.../elevation/ossimElevationShapeIdx.cpp | 275 ---
.../imaging/ossimEsriShapeFileCutter.cpp | 323 ---
.../imaging/ossimEsriShapeFileFilter.cpp | 813 -------
.../imaging/ossimImageShapeFileIndex.cpp | 197 --
.../imaging/ossimImageSourceInterface.cpp | 18 -
Utilities/otbossim/src/ossim/vec/dbfopen.c | 1552 -------------
.../src/ossim/vec/ossimShapeDatabase.cpp | 268 ---
.../otbossim/src/ossim/vec/ossimShapeFile.cpp | 220 --
Utilities/otbossim/src/ossim/vec/shpopen.c | 2021 -----------------
Utilities/otbossim/src/ossim/vec/shptree.c | 986 --------
18 files changed, 8342 deletions(-)
delete mode 100644 Utilities/otbossim/include/ossim/elevation/ossimElevationShapeIdx.h
delete mode 100644 Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileCutter.h
delete mode 100644 Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileFilter.h
delete mode 100644 Utilities/otbossim/include/ossim/imaging/ossimImageShapeFileIndex.h
delete mode 100644 Utilities/otbossim/include/ossim/imaging/ossimImageSourceInterface.h
delete mode 100644 Utilities/otbossim/include/ossim/vec/ossimShapeDatabase.h
delete mode 100644 Utilities/otbossim/include/ossim/vec/ossimShapeFile.h
delete mode 100644 Utilities/otbossim/include/ossim/vec/shapefil.h
delete mode 100644 Utilities/otbossim/src/ossim/elevation/ossimElevationShapeIdx.cpp
delete mode 100644 Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileCutter.cpp
delete mode 100644 Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileFilter.cpp
delete mode 100644 Utilities/otbossim/src/ossim/imaging/ossimImageShapeFileIndex.cpp
delete mode 100644 Utilities/otbossim/src/ossim/imaging/ossimImageSourceInterface.cpp
delete mode 100644 Utilities/otbossim/src/ossim/vec/dbfopen.c
delete mode 100644 Utilities/otbossim/src/ossim/vec/ossimShapeDatabase.cpp
delete mode 100644 Utilities/otbossim/src/ossim/vec/ossimShapeFile.cpp
delete mode 100644 Utilities/otbossim/src/ossim/vec/shpopen.c
delete mode 100644 Utilities/otbossim/src/ossim/vec/shptree.c
diff --git a/Utilities/otbossim/include/ossim/elevation/ossimElevationShapeIdx.h b/Utilities/otbossim/include/ossim/elevation/ossimElevationShapeIdx.h
deleted file mode 100644
index 6f55b6b5c3..0000000000
--- a/Utilities/otbossim/include/ossim/elevation/ossimElevationShapeIdx.h
+++ /dev/null
@@ -1,41 +0,0 @@
-#ifndef ossimElevationShapeIdx_HEADER
-#define ossimElevationShapeIdx_HEADER
-#include <exception>
-#include <ossim/base/ossimObject.h>
-#include <ossim/base/ossimRefPtr.h>
-#include <ossim/base/ossimFilename.h>
-
-class ossimElevationShapeIdxPrivate;
-
-class ossimElevationShapeIdx : public ossimObject
-{
-public:
- ossimElevationShapeIdx(const ossimFilename& shapeFilename = "",
- bool overWriteFlag = false);
- virtual ~ossimElevationShapeIdx();
- void setFilename(const ossimFilename& file, bool overWriteFlag=false);
-
- void add(const ossimFilename& filename,
- ossim_float64 minx, // decimal degrees
- ossim_float64 miny, // decimal degrees
- ossim_float64 minz, // height in meters meters
- ossim_float64 maxx, // decimal degrees
- ossim_float64 maxy, // decimal degrees
- ossim_float64 maxz); // height in meters meters
-
- void buildQuadTree(int maxDepth=0);
- static const ossimString& defaultName();
-
-protected:
- static ossimString theDefaultName;
- ossimElevationShapeIdxPrivate* thePrivateData;
-
- ossimFilename theShapeFilename;
- ossim_float64 theGlobalMinx;
- ossim_float64 theGlobalMiny;
- ossim_float64 theGlobalMaxx;
- ossim_float64 theGlobalMaxy;
-
-};
-
-#endif
diff --git a/Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileCutter.h b/Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileCutter.h
deleted file mode 100644
index 0d64ed4a23..0000000000
--- a/Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileCutter.h
+++ /dev/null
@@ -1,200 +0,0 @@
-//*******************************************************************
-// Copyright (C) 2000 ImageLinks Inc.
-//
-// License: LGPL
-//
-// See LICENSE.txt file in the top level directory for more details.
-//
-// Author: Garrett Potts
-//
-//*************************************************************************
-// $Id: ossimEsriShapeFileCutter.h 9094 2006-06-13 19:12:40Z dburken $
-#ifndef ossimEsriShapeFileCutter_HEADER
-#define ossimEsriShapeFileCutter_HEADER
-#include <ossim/imaging/ossimImageSourceFilter.h>
-#include <ossim/base/ossimViewInterface.h>
-#include <ossim/imaging/ossimMaskFilter.h>
-#include <ossim/imaging/ossimEsriShapeFileFilter.h>
-
-/**
- * This class wraps the ossimEsriShapeFileFilter and ossimMaskFilter to implement
- * a composite cutter.
- *
- * the ossimEsriShapeFileFilter is saved with a shape_file_filter. prefix and the
- * ossimMaskFilter is saved with a mask. prefix. Please see ossimEsriShapeFileFilter
- * and ossimMaskFilter classes for their keywords. theborder size has an optional units
- * associated with it.
- * <pre>
- *
- * The border size can be us-feet, feet, meters, or degrees and is only applied
- * to polygons.
- *
- * Example keyword list
- *
- * shape_file_filter.brush_color: 255 255 255
- * shape_file_filter.feature_name:
- * shape_file_filter.filename:
- * shape_file_filter.fill_flag: 0
- * shape_file_filter.max_quadtree_levels: 10
- * shape_file_filter.pen_color: 255 255 255
- * shape_file_filter.point_width_height: 1 1
- * shape_file_filter.thickness: 1
- * shape_file_filter.border_size: 50 meters
- * shape_file_filter.type: ossimEsriShapeFileFilter
- * mask.mask_type: select
- * mask.type: ossimMaskFilter
-
- * </pre>
- *
- */
-class ossimEsriShapeFileCutter : public ossimImageSourceFilter,
- public ossimViewInterface
-{
-public:
- /**
- * allocates the mask and esri loader
- */
- ossimEsriShapeFileCutter();
-
- /**
- * destroys the data members
- */
- virtual ~ossimEsriShapeFileCutter();
-
- /**
- * Sent to the input
- */
- virtual void getDecimationFactor(ossim_uint32 resLevel,
- ossimDpt& result)const;
-
- /**
- * The bounding rect call is passed to the mask filter
- */
- virtual ossimIrect getBoundingRect(ossim_uint32 resLevel=0)const;
-
- /**
- * The get tile call is passed to the mask filter
- */
- virtual ossimRefPtr<ossimImageData> getTile(const ossimIrect& rect,
- ossim_uint32 resLevel=0);
-
- /**
- * is passed to the filter and the shape file loader
- */
- virtual void initialize();
-
- /**
- * passed to the input source.
- */
- virtual void getDecimationFactors(vector<ossimDpt>& decimations)const;
-
- /**
- * passed to the input source.
- */
- virtual ossim_uint32 getNumberOfDecimationLevels()const;
-
- /**
- * passed to the input source.
- */
- virtual ossim_uint32 getNumberOfOutputBands() const;
-
- /**
- * passed to the input source.
- */
- virtual ossimScalarType getOutputScalarType() const;
-
- /**
- * passed to the input source.
- */
- virtual void getValidImageVertices(
- vector<ossimIpt>& validVertices,
- ossimVertexOrdering ordering=OSSIM_CLOCKWISE_ORDER,
- ossim_uint32 resLevel=0)const;
-
- /**
- * passed to the input source.
- */
- virtual ossim_uint32 getTileWidth() const;
-
- /**
- * passed to the input source.
- */
- virtual ossim_uint32 getTileHeight() const;
-
- /**
- * passed to the input source.
- */
- virtual double getNullPixelValue(ossim_uint32 band)const;
-
- /**
- * passed to the input source.
- */
- virtual double getMinPixelValue(ossim_uint32 band=0)const;
-
- /**
- * passed to the input source.
- */
- virtual double getMaxPixelValue(ossim_uint32 band=0)const;
-
- /**
- * passed to the shape file to load the shape file that will be used
- * by the mask filter as input.
- */
- void loadFile(const ossimFilename& filename);
-
- /**
- * this is here to support the view interface of the esri shape file.
- * it passes view calls to the esri filter
- */
- virtual bool setView(ossimObject* baseObject,
- bool ownsTheView = false);
-
- /**
- * this is here to support the view interface of the esri shape file.
- * Will return the view setting for the esri shape file.
- */
- virtual ossimObject* getView();
-
- /**
- * this is here to support the view interface of the esri shape file.
- * Will return the view setting for the esri shape file.
- */
- virtual const ossimObject* getView()const;
-
- /**
- * Return the current esri shape file being used
- */
- ossimFilename getFilename()const;
-
- /**
- * Sets the mask type. Passes the call to the Mask filter.
- */
- void setMaskType(ossimMaskFilter::ossimFileSelectionMaskType type);
-
- /**
- * returns the mask type from the maks filter.
- */
- ossimMaskFilter::ossimFileSelectionMaskType getMaskType()const;
-
- /**
- * Passes load state call to the esri shape file loader with prefix
- * shape_file_filter. Then passes the load state cll to the mask
- * filter with prefix = mask.
- */
- virtual bool loadState(const ossimKeywordlist& kwl,
- const char* prefix=0);
-
- /**
- * Passes save state call to the esri shape file loader with prefix
- * shape_file_filter. Then passes the load state cll to the mask
- * filter with prefix = mask.
- */
- virtual bool saveState(ossimKeywordlist& kwl,
- const char* prefix=0)const;
-protected:
- ossimEsriShapeFileFilter* theEsriShapeFile;
- ossimMaskFilter* theMaskFilter;
-TYPE_DATA
-};
-
-#endif
diff --git a/Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileFilter.h b/Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileFilter.h
deleted file mode 100644
index 18c3f85b05..0000000000
--- a/Utilities/otbossim/include/ossim/imaging/ossimEsriShapeFileFilter.h
+++ /dev/null
@@ -1,270 +0,0 @@
-//*******************************************************************
-//
-// License: LGPL
-//
-// See LICENSE.txt file in the top level directory for more details.
-//
-// Author: Garrett Potts
-//
-//*************************************************************************
-// $Id: ossimEsriShapeFileFilter.h 9094 2006-06-13 19:12:40Z dburken $
-
-#ifndef ossimEsriShapeFileFilter_HEADER
-#define ossimEsriShapeFileFilter_HEADER
-#include <map>
-#include <ossim/imaging/ossimGeoAnnotationSource.h>
-#include <ossim/base/ossimViewInterface.h>
-#include <ossim/vec/shapefil.h>
-#include <ossim/vec/ossimShapeFile.h>
-#include <ossim/base/ossimRgbVector.h>
-
-class ossimGeoAnnotationObject;
-class ossimAnnotationObject;
-
-/*!
- * class ossimEsriShapeFileFilter
- *
- * This class is used to render shape files. Iif this filter has
- * an input connection to an ossimImageSourceInterface then it will
- * draw any vectors over the input tile. If it is not connected
- * it will just render its vector data to a tile and return it.
- * <pre>
- * supported keywords:
- *
- * max_quadtree_levels: // number of levels for faster hit testing
- * // defaults to 10
- *
- * fill_flag: // 1 for true 0 for false. Default is false
- * // and specifies how all closed surfaces are
- * // to be drawn
- *
- * feature_name: // currently not used.
- *
- * pen_color: // specified in Red Green Blue and is the color used
- * // if the fill flag is not specified. Future changes
- * // might be to include this as an outline color for
- * // filled objects. Default is white. Example: 255 255 255
- *
- * brush_color: // Color used by objects that have the fill flag enabled
- * // format is Red Green Blue. Default is white.
- *
- * line_thickness: // line drawing thickness in pixels.
- *
- * border_size: // can be us feet, feet, meters, degrees,
- * // will automatically expand all polygons by the
- * // specified size. Uses the centroid to estimate
- * // non degree units when converting to degrees.
- *
- * point_width_height: // Specifies the width and height of a point object in
- * // pixels. Default is 1 and has the following format
- * // w h. Example 1 1
- *
- * filename: // The esri shape file to be used
- *
- * example Keyword list: See ossimAnnotationSource for any additional keywords
- *
- *
- * brush_color: 255 255 255
- * feature_name:
- * filename:
- * fill_flag: 0
- * max_quadtree_levels: 10
- * pen_color: 255 255 255
- * point_width_height: 1 1
- * thickness: 1
- * border_size: 25 meters
- * type: ossimEsriShapeFileFilter
- *
- * </pre>
- */
-class ossimEsriShapeFileFilter : public ossimAnnotationSource,
- public ossimViewInterface
-{
-public:
- ossimEsriShapeFileFilter(ossimImageSource* inputSource=NULL);
- ~ossimEsriShapeFileFilter();
-
- virtual bool setView(ossimObject* baseObject,
- bool ownsTheView = false);
-
- virtual ossimObject* getView();
- virtual const ossimObject* getView()const;
-
- virtual bool addObject(ossimAnnotationObject* anObject);
- virtual void transformObjects(ossimProjection* projection=0);
- virtual void setProjection(ossimProjection* projection,
- bool ownsProjectionFlag=false);
- virtual bool getImageGeometry(ossimKeywordlist& kwl,
- const char* prefix = 0);
- virtual ossimIrect getBoundingRect(ossim_uint32 resLevel=0)const;
- virtual void computeBoundingRect();
-
- virtual void drawAnnotations(ossimRefPtr<ossimImageData> tile);
- /*!
- * Will delete the current objects within the layer and add all
- * objects in the passed in file.
- */
- virtual bool loadShapeFile(const ossimFilename& shapeFile);
-
- virtual ossimFilename getFilename()const
- {
- return theShapeFile.getFilename();
- }
- virtual ossim_int32 getMaxQuadTreeLevels()const
- {
- return theMaxQuadTreeLevels;
- }
-
- virtual void setMaxQuadTreeLevels(ossim_int32 levels)
- {
- theMaxQuadTreeLevels = (levels>0?levels:10);
- }
-
- virtual void setBrushColor(const ossimRgbVector& brushColor)
- {
- theBrushColor = brushColor;
- }
-
- virtual void setPenColor(const ossimRgbVector& penColor)
- {
- thePenColor = penColor;
- }
-
- virtual ossimRgbVector getPenColor()const
- {
- return thePenColor;
- }
-
- virtual ossimRgbVector getBrushColor()const
- {
- return theBrushColor;
- }
-
- virtual bool getFillFlag()const
- {
- return theFillFlag;
- }
-
- virtual void setFillFlag(bool flag)
- {
- theFillFlag=flag;
- }
- virtual ossimString getFeatureName()const
- {
- return theFeatureName;
- }
- virtual void setFeatureName(const ossimString& name)
- {
- theFeatureName = name;
- }
-
- virtual void setThickness(ossim_int32 thickness)
- {
- theThickness = thickness >= 0? thickness:1;
- }
- virtual ossim_int32 getThickness()const
- {
- return theThickness;
- }
- virtual void setPointRadius(double r)
- {
- thePointWidthHeight = ossimDpt(fabs(r)*2, fabs(r)*2);
- }
- virtual double getPointRadius()const
- {
- return thePointWidthHeight.x/2.0;
- }
-
- virtual ossimAnnotationObject* nextObject(bool restart=false)
- {
- if(restart)
- {
- theCurrentObject = theShapeCache.begin();
- }
- else
- {
- if(theCurrentObject != theShapeCache.end())
- {
- ++theCurrentObject;
- }
-
- }
-
- if(theCurrentObject == theShapeCache.end())
- {
- return (ossimAnnotationObject*)NULL;
- }
- return (*theCurrentObject).second;
- }
-
-
- virtual bool saveState(ossimKeywordlist& kwl,
- const char* prefix=NULL)const;
-
- virtual bool loadState(const ossimKeywordlist& kwl,
- const char* prefix=NULL);
-protected:
-
- /*!
- * This holds the output projection. This is used to transform
- * the objects to the projection plane.
- */
- ossimProjection* theViewProjection;
-
- /*!
- * If the shape file is projected we must know the projector.
- * This will allow us to define the geographic objects correctly.
- */
- ossimProjection* theShapeFileProjection;
-
- /*!
- * Flag that indicates if we own the view projection.
- */
- bool theOwnsViewProjectionFlag;
-
- /*!
- * Specifies if its image/pixel space, Geographic, projected ...etc.
- */
- ossimCoordinateSystemType theCoordinateSystem;
-
- /*!
- * If its projected this spceifies the type of unit. We
- * need to know if the projection is in meters, feet, us survey feet
- * ... etc.
- */
- ossimUnitType theUnitType;
-
- ossim_SHPTree* theTree;
-
- ossimShapeFile theShapeFile;
-
- double theMinArray[4];
- double theMaxArray[4];
-
- int theMaxQuadTreeLevels;
- ossimRgbVector thePenColor;
- ossimRgbVector theBrushColor;
- bool theFillFlag;
- ossim_int32 theThickness;
- ossimString theFeatureName;
- ossimDpt thePointWidthHeight;
- double theBorderSize;
- ossimUnitType theBorderSizeUnits;
-
- mutable std::multimap<int, ossimAnnotationObject*>::iterator theCurrentObject;
-
- std::multimap<int, ossimAnnotationObject*> theShapeCache;
- ossimDrect theBoundingRect;
-
- void removeViewProjection();
- void deleteCache();
- void checkAndSetDefaultView();
-
- virtual void loadPolygon(ossimShapeObject& obj);
- virtual void loadPoint(ossimShapeObject& obj);
- virtual void loadArc(ossimShapeObject& obj);
-
-TYPE_DATA
-};
-
-#endif
diff --git a/Utilities/otbossim/include/ossim/imaging/ossimImageShapeFileIndex.h b/Utilities/otbossim/include/ossim/imaging/ossimImageShapeFileIndex.h
deleted file mode 100644
index 8973473fdd..0000000000
--- a/Utilities/otbossim/include/ossim/imaging/ossimImageShapeFileIndex.h
+++ /dev/null
@@ -1,61 +0,0 @@
-#ifndef ossimImageShapeFileIndex_HEADER
-#define ossimImageShapeFileIndex_HEADER
-#include <ossim/base/ossimObject.h>
-#include <ossim/base/ossimReferenced.h>
-#include <ossim/base/ossimRefPtr.h>
-#include <ossim/base/ossimFilename.h>
-
-class OSSIM_DLL ossimImageShapeFileIndex : public ossimObject
-{
-public:
- ossimImageShapeFileIndex(const ossimFilename& shapeFilename = "");
- void setFilename(const ossimFilename& shapeFilename);
-
- void add(const ossimFilename& filename,
- const ossim_int32 entry=-1, // default to all entries
- const ossimString& groupName="RASTER");
- /**
- *
- */
- void build(ossim_uint32 maxDepth=1)const;
-
-protected:
- class ossimImageShapeFileIndexRecord : public ossimReferenced
- {
- public:
- ossimImageShapeFileIndexRecord(const ossimFilename& filename="",
- ossim_uint32 entry=0,
- const ossimString& groupName="",
- ossim_float64 minx=0.0,
- ossim_float64 miny=0.0,
- ossim_float64 maxx=0.0,
- ossim_float64 maxy=0.0)
- :theFilename(filename),
- theEntry(entry),
- theGroupName(groupName),
- theMinx(minx),
- theMiny(miny),
- theMaxx(maxx),
- theMaxy(maxy)
- {
- }
-
- ossimFilename theFilename;
- ossim_uint32 theEntry;
- ossimString theGroupName;
- ossim_float64 theMinx;
- ossim_float64 theMiny;
- ossim_float64 theMaxx;
- ossim_float64 theMaxy;
- };
-
- std::vector<ossimRefPtr<ossimImageShapeFileIndexRecord> > theRecords;
-
- ossimFilename theFilename;
- ossim_float64 theGlobalMinx;
- ossim_float64 theGlobalMiny;
- ossim_float64 theGlobalMaxx;
- ossim_float64 theGlobalMaxy;
-};
-
-#endif
diff --git a/Utilities/otbossim/include/ossim/imaging/ossimImageSourceInterface.h b/Utilities/otbossim/include/ossim/imaging/ossimImageSourceInterface.h
deleted file mode 100644
index ea924e442a..0000000000
--- a/Utilities/otbossim/include/ossim/imaging/ossimImageSourceInterface.h
+++ /dev/null
@@ -1,109 +0,0 @@
-//*******************************************************************
-// Copyright (C) 2000 ImageLinks Inc.
-//
-// License: See top level LICENSE.txt file.
-//
-// Author: Garrett Potts
-//
-//*************************************************************************
-// $Id: ossimImageSourceInterface.h 9094 2006-06-13 19:12:40Z dburken $
-#ifndef ossimImageSourceInterface_HEADER
-#define ossimImageSourceInterface_HEADER
-
-#include <vector>
-
-#include <ossim/base/ossimErrorContext.h>
-#include <ossim/base/ossimConstants.h>
-#include <ossim/base/ossimRefPtr.h>
-#include <ossim/base/ossimIrect.h>
-#include <ossim/base/ossimKeywordlist.h>
-#include <ossim/base/ossimFilename.h>
-
-class ossimImageData;
-class ossimObject;
-
-class OSSIMDLLEXPORT ossimImageSourceInterface
-{
-public:
- virtual ~ossimImageSourceInterface();
-
- virtual ossimObject* getObject()=0;
- virtual const ossimObject* getObject()const=0;
-
- virtual ossimRefPtr<ossimImageData> getTile(const ossimIpt& origin,
- ossim_uint32 resLevel=0) = 0;
-
- virtual ossimRefPtr<ossimImageData> getTile(const ossimIrect& rect,
- ossim_uint32 resLevel=0) = 0;
-
- virtual void getDecimationFactor(ossim_uint32 resLevel,
- ossimDpt& result)const=0;
- virtual void getDecimationFactors(std::vector<ossimDpt>& decimations)const=0;
- virtual ossim_uint32 getNumberOfDecimationLevels()const=0;
-
- /*!
- * Returns the number of bands available from the input.
- */
- virtual ossim_uint32 getNumberOfInputBands() const = 0;
-
- /*!
- * Returns the number of bands in output.
- */
- virtual ossim_uint32 getNumberOfOutputBands() const = 0;
-
- /*!
- * Initializes bandList to the zero based order of output bands.
- */
- virtual void getOutputBandList(std::vector<ossim_uint32>& bandList) const=0;
-
- virtual ossimScalarType getOutputScalarType() const = 0;
-
- virtual ossim_uint32 getTileWidth() const = 0;
- virtual ossim_uint32 getTileHeight() const = 0;
-
- virtual double getNullPixelValue(ossim_uint32 band=0)const=0;
- virtual double getMinPixelValue(ossim_uint32 band=0)const=0;
- virtual double getMaxPixelValue(ossim_uint32 band=0)const=0;
-
- virtual ossimIrect getBoundingRect(ossim_uint32 resLevel=0)const=0;
-
- virtual void getValidImageVertices(std::vector<ossimIpt>& validVertices,
- ossimVertexOrdering ordering=OSSIM_CLOCKWISE_ORDER,
- ossim_uint32 resLevel=0)const=0;
-
- virtual bool getImageGeometry(ossimKeywordlist& kwl,
- const char* prefix=0)=0;
-
- /**
- * Sets the image geometry keyword list.
- *
- * @note Callers should note that this will override any existing
- * geometry.
- *
- * @note This does not write the geometry to disk, to do so call
- * one of the saveImageGeometry methods.
- */
- virtual void setImageGeometry(const ossimKeywordlist& kwl)=0;
-
- /**
- * Saves the image geometry to a file that will be made based on the image
- * file name.
- */
- virtual void saveImageGeometry()const=0;
-
- /**
- * Saves the image geomety to geom_file.
- */
- virtual void saveImageGeometry(const ossimFilename& geometry_file)const=0;
-
- virtual void initialize()=0;
-
- virtual bool loadState(const ossimKeywordlist& kwl,
- const char* prefix=0)=0;
- virtual bool saveState(ossimKeywordlist& kwl,
- const char* prefix=0)const=0;
-
-TYPE_DATA
-};
-
-#endif
diff --git a/Utilities/otbossim/include/ossim/vec/ossimShapeDatabase.h b/Utilities/otbossim/include/ossim/vec/ossimShapeDatabase.h
deleted file mode 100644
index 95465b7a9d..0000000000
--- a/Utilities/otbossim/include/ossim/vec/ossimShapeDatabase.h
+++ /dev/null
@@ -1,136 +0,0 @@
-#ifndef ossimShapeDatabase_HEADER
-#define ossimShapeDatabase_HEADER
-#include <iostream>
-#include <ossim/vec/shapefil.h>
-#include <ossim/base/ossimFilename.h>
-#include <ossim/base/ossimString.h>
-
-
-class OSSIM_DLL ossimShapeDatabaseField
-{
-public:
- ossimString fieldTypeAsString()const;
-
- ossimString theName;
- int theWidth;
- int theDecimals;
- ossim_DBFFieldType theFieldType;
- ossimString theValue;
-};
-
-class OSSIM_DLL ossimShapeDatabaseRecord
-{
-
-public:
- bool getField(ossimShapeDatabaseField& result,
- ossim_uint32 i)
- {
- if(i < theFieldArray.size())
- {
- result = theFieldArray[i];
- return true;
- }
- return false;
- }
- bool setField(const ossimShapeDatabaseField& field,
- ossim_uint32 i)
- {
- if(i < theFieldArray.size())
- {
- theFieldArray[i] = field;
- return true;
- }
- return false;
- }
- int getNumberOfFields()const
- {
- return theFieldArray.size();
- }
- void setNumberOfFields(int n)
- {
- if(n)
- {
- theFieldArray.resize(n);
- }
- else
- {
- theFieldArray.clear();
- }
- }
- ossim_int32 getFieldIdx(const ossimString& name,
- bool caseInsensitive=true)const;
-protected:
- std::vector<ossimShapeDatabaseField> theFieldArray;
-};
-
-class OSSIM_DLL ossimShapeDatabase
-{
-public:
- friend std::ostream& operator <<(std::ostream& out, const ossimShapeDatabase& rhs)
- {
- rhs.print(out);
-
- return out;
- }
-
- ossimShapeDatabase();
- virtual ~ossimShapeDatabase();
-
- virtual bool open(const ossimFilename& file,
- const ossimString& flags=ossimString("rb"));
-
- virtual void close();
-
- bool getRecord(ossimShapeDatabaseRecord& result);
-
- bool getRecord(ossimShapeDatabaseRecord& result,
- int recordNumber)
- {
- if(isOpen())
- {
- if(recordNumber < getNumberOfRecords())
- {
- theRecordNumber = recordNumber;
- return getRecord(result);
- }
- }
-
- return false;
- }
-
- bool getNextRecord(ossimShapeDatabaseRecord& result)
- {
- if(isOpen() && ((theRecordNumber+1) < getNumberOfRecords()))
- {
- ++theRecordNumber;
- return getRecord(result);
- }
-
- return false;
- }
-
- int getNumberOfRecords()const;
- bool isOpen()const
- {
- return (theHandle!=NULL);
- }
- virtual ossim_DBFHandle getHandle()
- {
- return theHandle;
- }
-
- virtual const ossim_DBFHandle getHandle()const
- {
- return theHandle;
- }
-
- virtual void print(std::ostream& out)const;
-
-protected:
- ossim_DBFHandle theHandle;
- ossimFilename theFilename;
-
- int theRecordNumber;
-};
-
-#endif
diff --git a/Utilities/otbossim/include/ossim/vec/ossimShapeFile.h b/Utilities/otbossim/include/ossim/vec/ossimShapeFile.h
deleted file mode 100644
index 02c45848c6..0000000000
--- a/Utilities/otbossim/include/ossim/vec/ossimShapeFile.h
+++ /dev/null
@@ -1,312 +0,0 @@
-//*******************************************************************
-//
-// License: See top level LICENSE.txt file.
-//
-// Author: Garrett Potts
-//
-//*************************************************************************
-// $Id: ossimShapeFile.h 11424 2007-07-27 17:01:23Z dburken $
-
-#ifndef ossimShapeFile_HEADER
-#define ossimShapeFile_HEADER
-#include <iostream>
-#include <ossim/base/ossimConstants.h>
-#include <ossim/base/ossimFilename.h>
-#include <ossim/base/ossimString.h>
-#include <ossim/base/ossimDrect.h>
-
-#include <ossim/vec/shapefil.h>
-
-class ossimShapeFile;
-class ossimAnnotationObject;
-
-class OSSIMDLLEXPORT ossimShapeObject
-{
-public:
- friend std::ostream& operator<<(std::ostream& out,
- const ossimShapeObject& rhs);
-
- ossimShapeObject()
- :theShape((ossim_SHPObject*)NULL),
- theIndex(-1)
- {
- }
-
- ~ossimShapeObject()
- {
- if(theShape)
- {
- ossim_SHPDestroyObject(theShape);
- theShape = NULL;
- }
- }
- void setShape(ossim_SHPObject* obj);
- bool isLoaded()const
- {
- return (theShape!=NULL);
- }
- long getIndex()const
- {
- return theIndex;
- }
- long getId()const
- {
- if(theShape)
- {
- return theShape->nShapeId;
- }
-
- return -1;
- }
-
- bool loadShape(const ossimShapeFile& shapeFile,
- long shapeRecord);
-
- ossim_uint32 getNumberOfParts()const
- {
- if(theShape)
- {
- return theShape->nParts;
- }
-
- return 0;
- }
- ossim_uint32 getNumberOfVertices()const
- {
- if(theShape)
- {
- return theShape->nVertices;
- }
-
- return 0;
- }
-
- void getBounds(double& minX, double& minY, double& minZ, double& minM,
- double& maxX, double& maxY, double& maxZ, double& maxM)const
- {
- if(theShape)
- {
- minX = theShape->dfXMin;
- minY = theShape->dfYMin;
- minZ = theShape->dfZMin;
- minM = theShape->dfMMin;
- maxX = theShape->dfXMax;
- maxY = theShape->dfYMax;
- maxZ = theShape->dfZMax;
- maxM = theShape->dfMMax;
- }
- else
- {
- minX = minY = minZ = minM =
- maxX = maxY = maxZ = maxM = ossim::nan();
- }
- }
-
- void getBounds(double& minX, double& minY,
- double& maxX, double& maxY)const
- {
- if(theShape)
- {
- minX = theShape->dfXMin;
- minY = theShape->dfYMin;
- maxX = theShape->dfXMax;
- maxY = theShape->dfYMax;
- }
- else
- {
- minX = minY =
- maxX = maxY = ossim::nan();
- }
- }
-
- void getBoundingRect(ossimDrect& result,
- ossimCoordSysOrientMode orient = OSSIM_RIGHT_HANDED)const;
-
- ossimDrect getBoundingRect(ossimCoordSysOrientMode orient = OSSIM_RIGHT_HANDED)const
- {
- ossimDrect result;
-
- getBoundingRect(result, orient);
-
- return result;
- }
-
- int getType()const
- {
- if(theShape)
- {
- return theShape->nSHPType;
- }
- return ossim_SHPT_NULL;
- }
- ossimString getTypeByName()const
- {
- if(theShape)
- {
- return ossimString(ossim_SHPTypeName(theShape->nSHPType));
- }
- return "unknown";
- }
- int getPartType(ossim_uint32 partIndex)const
- {
- if((partIndex > getNumberOfParts())||
- (!theShape))
- {
- return -1;
- }
-
- return theShape->panPartType[partIndex];
- }
- ossimString getPartByName(ossim_uint32 partIndex)const
- {
- if((partIndex > getNumberOfParts())||
- (!theShape))
- {
- }
- switch(theShape->panPartType[partIndex])
- {
- case ossim_SHPP_TRISTRIP:
- {
- ossimString("tristrip");
- break;
- }
- case ossim_SHPP_TRIFAN:
- {
- ossimString("trifan");
- break;
- }
- case ossim_SHPP_OUTERRING:
- {
- ossimString("outerring");
- break;
- }
- case ossim_SHPP_INNERRING:
- {
- ossimString("innerring");
- break;
- }
- case ossim_SHPP_RING:
- {
- ossimString("ring");
- break;
- }
- }
- return ossimString("unknown");
- }
- ossim_SHPObject* getShapeObject()
- {
- return theShape;
- }
- const ossim_SHPObject* getShapeObject()const
- {
- return theShape;
- }
-
-
-protected:
- ossim_SHPObject* theShape;
- long theIndex;
-};
-
-class OSSIMDLLEXPORT ossimShapeFile
-{
-public:
- friend std::ostream& operator <<(std::ostream& out, const ossimShapeFile& rhs)
- {
- rhs.print(out);
-
- return out;
- }
-
- ossimShapeFile();
- virtual ~ossimShapeFile();
-
- virtual bool open(const ossimFilename& file,
- const ossimString& flags=ossimString("rb"));
-
- virtual void close();
- bool isOpen()const
- {
- return (theHandle!=NULL);
- }
- virtual ossim_SHPHandle getHandle()
- {
- return theHandle;
- }
-
- virtual const ossim_SHPHandle getHandle()const
- {
- return theHandle;
- }
-
- virtual void print(std::ostream& out)const;
- virtual ossimString getShapeTypeString()const;
-
-
- virtual long getNumberOfShapes()const
- {
- if(theHandle)
- {
- return theHandle->nRecords;
- }
- return 0;
- }
- void getBounds(double& minX, double& minY, double& minZ, double& minM,
- double& maxX, double& maxY, double& maxZ, double& maxM)const
- {
- if(theHandle)
- {
- minX = theHandle->adBoundsMin[0];
- minY = theHandle->adBoundsMin[1];
- minZ = theHandle->adBoundsMin[2];
- minM = theHandle->adBoundsMin[3];
- maxX = theHandle->adBoundsMax[0];
- maxY = theHandle->adBoundsMax[1];
- maxZ = theHandle->adBoundsMax[2];
- maxM = theHandle->adBoundsMax[3];
- }
- else
- {
- minX = minY = minZ = minM =
- maxX = maxY = maxZ = maxM = ossim::nan();
- }
- }
-
- void getBounds(double& minX, double& minY,
- double& maxX, double& maxY)const
- {
- if(theHandle)
- {
- minX = theHandle->adBoundsMin[0];
- minY = theHandle->adBoundsMin[1];
- maxX = theHandle->adBoundsMax[0];
- maxY = theHandle->adBoundsMax[1];
- }
- else
- {
- minX = minY = maxX = maxY = ossim::nan();
- }
- }
-
- void getBoundingRect(ossimDrect& result,
- ossimCoordSysOrientMode orient = OSSIM_RIGHT_HANDED)const;
-
- ossimDrect getBoundingRect(ossimCoordSysOrientMode orient = OSSIM_RIGHT_HANDED)const
- {
- ossimDrect result;
-
- getBoundingRect(result, orient);
-
- return result;
- }
- const ossimFilename& getFilename()const
- {
- return theFilename;
- }
-protected:
- ossim_SHPHandle theHandle;
- ossimFilename theFilename;
-
-};
-
-#endif
diff --git a/Utilities/otbossim/include/ossim/vec/shapefil.h b/Utilities/otbossim/include/ossim/vec/shapefil.h
deleted file mode 100644
index 7dfa62d688..0000000000
--- a/Utilities/otbossim/include/ossim/vec/shapefil.h
+++ /dev/null
@@ -1,540 +0,0 @@
-#ifndef _SHAPEFILE_H_INCLUDED
-#define _SHAPEFILE_H_INCLUDED
-
-/******************************************************************************
- * $Id: shapefil.h 11039 2007-05-24 15:42:43Z gpotts $
- *
- * Project: Shapelib
- * Purpose: Primary include file for Shapelib.
- * 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 LICENSE.LGPL). 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$
- * Revision 1.5 2005/10/08 12:57:26 gpotts
- * Added adiitional patches
- *
- * Revision 1.14 2005/03/14 15:21:29 fwarmerdam
- * updated
- *
- * Revision 1.30 2005/01/03 22:30:13 fwarmerdam
- * added support for saved quadtrees
- *
- * Revision 1.29 2004/09/26 20:09:35 fwarmerdam
- * avoid rcsid warnings
- *
- * Revision 1.28 2003/12/29 06:02:18 fwarmerdam
- * added cpl_error.h option
- *
- * Revision 1.27 2003/04/21 18:30:37 warmerda
- * added header write/update public methods
- *
- * Revision 1.26 2002/09/29 00:00:08 warmerda
- * added FTLogical and logical attribute read/write calls
- *
- * Revision 1.25 2002/05/07 13:46:30 warmerda
- * added ossim_DBFWriteAttributeDirectly().
- *
- * Revision 1.24 2002/04/10 16:59:54 warmerda
- * added ossim_SHPRewindObject
- *
- * Revision 1.23 2002/01/15 14:36:07 warmerda
- * updated email address
- *
- * Revision 1.22 2002/01/15 14:32:00 warmerda
- * try to improve ossim_SHPAPI_CALL docs
- *
- * Revision 1.21 2001/11/01 16:29:55 warmerda
- * move pabyRec into ossim_SHPInfo for thread safety
- *
- * Revision 1.20 2001/07/20 13:06:02 warmerda
- * fixed ossim_SHPAPI attribute for ossim_SHPTreeFindLikelyShapes
- *
- * Revision 1.19 2001/05/31 19:20:13 warmerda
- * added ossim_DBFGetFieldIndex()
- *
- * Revision 1.18 2001/05/31 18:15:40 warmerda
- * Added support for NULL fields in ossim_DBF files
- *
- * Revision 1.17 2001/05/23 13:36:52 warmerda
- * added use of ossim_SHPAPI_CALL
- *
- * Revision 1.16 2000/09/25 14:15:59 warmerda
- * added ossim_DBFGetNativeFieldType()
- *
- * Revision 1.15 2000/02/16 16:03:51 warmerda
- * added null shape support
- *
- * Revision 1.14 1999/11/05 14:12:05 warmerda
- * updated license terms
- *
- * Revision 1.13 1999/06/02 18:24:21 warmerda
- * added trimming code
- *
- * Revision 1.12 1999/06/02 17:56:12 warmerda
- * added quad'' subnode support for trees
- *
- * Revision 1.11 1999/05/18 19:11:11 warmerda
- * Added example searching capability
- *
- * Revision 1.10 1999/05/18 17:49:38 warmerda
- * added initial quadtree support
- *
- * Revision 1.9 1999/05/11 03:19:28 warmerda
- * added new Tuple api, and improved extension handling - add from candrsn
- *
- * Revision 1.8 1999/03/23 17:22:27 warmerda
- * Added extern "C" protection for C++ users of shapefil.h.
- *
- * Revision 1.7 1998/12/31 15:31:07 warmerda
- * Added the TRIM_DBF_WHITESPACE and DISABLE_MULTIPATCH_MEASURE options.
- *
- * Revision 1.6 1998/12/03 15:48:15 warmerda
- * Added ossim_SHPCalculateExtents().
- *
- * Revision 1.5 1998/11/09 20:57:16 warmerda
- * Altered ossim_SHPGetInfo() call.
- *
- * Revision 1.4 1998/11/09 20:19:33 warmerda
- * Added 3D support, and use of ossim_SHPObject.
- *
- * Revision 1.3 1995/08/23 02:24:05 warmerda
- * Added support for reading bounds.
- *
- * Revision 1.2 1995/08/04 03:17:39 warmerda
- * Added header.
- *
- */
-
-#include <stdio.h>
-
-#ifdef USE_DBMALLOC
-#include <dbmalloc.h>
-#endif
-
-#ifdef USE_CPL
-#include "cpl_error.h"
-#endif
-
-#ifdef __cplusplus
-extern "C" {
-#endif
-
-/************************************************************************/
-/* Configuration options. */
-/************************************************************************/
-
-/* -------------------------------------------------------------------- */
-/* Should the ossim_DBFReadStringAttribute() strip leading and */
-/* trailing white space? */
-/* -------------------------------------------------------------------- */
-#define TRIM_DBF_WHITESPACE
-
-/* -------------------------------------------------------------------- */
-/* Should we write measure values to the Multipatch object? */
-/* Reportedly ArcView crashes if we do write it, so for now it */
-/* is disabled. */
-/* -------------------------------------------------------------------- */
-#define DISABLE_MULTIPATCH_MEASURE
-
-/* -------------------------------------------------------------------- */
-/* ossim_SHPAPI_CALL */
-/* */
-/* The following two macros are present to allow forcing */
-/* various calling conventions on the Shapelib API. */
-/* */
-/* To force __stdcall conventions (needed to call Shapelib */
-/* from Visual Basic and/or Dephi I believe) the makefile could */
-/* be modified to define: */
-/* */
-/* /DSHPAPI_CALL=__stdcall */
-/* */
-/* If it is desired to force export of the Shapelib API without */
-/* using the shapelib.def file, use the following definition. */
-/* */
-/* /DSHAPELIB_DLLEXPORT */
-/* */
-/* To get both at once it will be necessary to hack this */
-/* include file to define: */
-/* */
-/* #define ossim_SHPAPI_CALL __declspec(dllexport) __stdcall */
-/* #define ossim_SHPAPI_CALL1 __declspec(dllexport) * __stdcall */
-/* */
-/* The complexity of the situtation is partly caused by the */
-/* peculiar requirement of Visual C++ that __stdcall appear */
-/* after any "*"'s in the return value of a function while the */
-/* __declspec(dllexport) must appear before them. */
-/* -------------------------------------------------------------------- */
-#include <ossim/base/ossimConstants.h>
-
-#ifdef OSSIMMAKINGDLL
-# define ossim_SHPAPI_CALL OSSIM_DLL
-# define ossim_SHPAPI_CALL1(x) OSSIM_DLL_DATA(x)
-#endif
-#ifndef ossim_SHPAPI_CALL
-# define ossim_SHPAPI_CALL
-#endif
-#ifndef ossim_SHPAPI_CALL1
-# define ossim_SHPAPI_CALL1(x) x ossim_SHPAPI_CALL
-#endif
-
-#if 0
-#ifdef SHAPELIB_DLLEXPORT
-# define ossim_SHPAPI_CALL __declspec(dllexport)
-# define ossim_SHPAPI_CALL1(x) __declspec(dllexport) x
-#endif
-#ifndef ossim_SHPAPI_CALL
-# define ossim_SHPAPI_CALL
-#endif
-#ifndef ossim_SHPAPI_CALL1
-# define ossim_SHPAPI_CALL1(x) x ossim_SHPAPI_CALL
-#endif
-#endif
-
-/* -------------------------------------------------------------------- */
-/* Macros for controlling CVSID and ensuring they don't appear */
-/* as unreferenced variables resulting in lots of warnings. */
-/* -------------------------------------------------------------------- */
-#ifndef DISABLE_CVSID
-# define ossim_SHP_CVSID(string) static char cpl_cvsid[] = string; \
-static char *cvsid_aw() { return( cvsid_aw() ? ((char *) NULL) : cpl_cvsid ); }
-#else
-# define ossim_SHP_CVSID(string)
-#endif
-
-/************************************************************************/
-/* ossim_SHP Support. */
-/************************************************************************/
-typedef struct
-{
- FILE *fpSHP;
- FILE *fpSHX;
-
- int nShapeType; /* ossim_SHPT_* */
-
- int nFileSize; /* ossim_SHP file */
-
- int nRecords;
- int nMaxRecords;
- int *panRecOffset;
- int *panRecSize;
-
- double adBoundsMin[4];
- double adBoundsMax[4];
-
- int bUpdated;
-
- unsigned char *pabyRec;
- int nBufSize;
-} ossim_SHPInfo;
-
-typedef ossim_SHPInfo * ossim_SHPHandle;
-
-/* -------------------------------------------------------------------- */
-/* Shape types (nSHPType) */
-/* -------------------------------------------------------------------- */
-#define ossim_SHPT_NULL 0
-#define ossim_SHPT_POINT 1
-#define ossim_SHPT_ARC 3
-#define ossim_SHPT_POLYGON 5
-#define ossim_SHPT_MULTIPOINT 8
-#define ossim_SHPT_POINTZ 11
-#define ossim_SHPT_ARCZ 13
-#define ossim_SHPT_POLYGONZ 15
-#define ossim_SHPT_MULTIPOINTZ 18
-#define ossim_SHPT_POINTM 21
-#define ossim_SHPT_ARCM 23
-#define ossim_SHPT_POLYGONM 25
-#define ossim_SHPT_MULTIPOINTM 28
-#define ossim_SHPT_MULTIPATCH 31
-
-
-/* -------------------------------------------------------------------- */
-/* Part types - everything but ossim_SHPT_MULTIPATCH just uses */
-/* ossim_SHPP_RING. */
-/* -------------------------------------------------------------------- */
-
-#define ossim_SHPP_TRISTRIP 0
-#define ossim_SHPP_TRIFAN 1
-#define ossim_SHPP_OUTERRING 2
-#define ossim_SHPP_INNERRING 3
-#define ossim_SHPP_FIRSTRING 4
-#define ossim_SHPP_RING 5
-
-/* -------------------------------------------------------------------- */
-/* ossim_SHPObject - represents on shape (without attributes) read */
-/* from the .shp file. */
-/* -------------------------------------------------------------------- */
-typedef struct
-{
- int nSHPType;
-
- int nShapeId; /* -1 is unknown/unassigned */
-
- int nParts;
- int *panPartStart;
- int *panPartType;
-
- int nVertices;
- double *padfX;
- double *padfY;
- double *padfZ;
- double *padfM;
-
- double dfXMin;
- double dfYMin;
- double dfZMin;
- double dfMMin;
-
- double dfXMax;
- double dfYMax;
- double dfZMax;
- double dfMMax;
-} ossim_SHPObject;
-
-/* -------------------------------------------------------------------- */
-/* ossim_SHP API Prototypes */
-/* -------------------------------------------------------------------- */
-ossim_SHPHandle ossim_SHPAPI_CALL
- ossim_SHPOpen( const char * pszShapeFile, const char * pszAccess );
-ossim_SHPHandle ossim_SHPAPI_CALL
- ossim_SHPCreate( const char * pszShapeFile, int nShapeType );
-void ossim_SHPAPI_CALL
- ossim_SHPGetInfo( ossim_SHPHandle hSHP, int * pnEntities, int * pnShapeType,
- double * padfMinBound, double * padfMaxBound );
-
-ossim_SHPObject ossim_SHPAPI_CALL1(*)
- ossim_SHPReadObject( ossim_SHPHandle hSHP, int iShape );
-int ossim_SHPAPI_CALL
- ossim_SHPWriteObject( ossim_SHPHandle hSHP, int iShape, ossim_SHPObject * psObject );
-
-void ossim_SHPAPI_CALL
- ossim_SHPDestroyObject( ossim_SHPObject * psObject );
-void ossim_SHPAPI_CALL
- ossim_SHPComputeExtents( ossim_SHPObject * psObject );
-ossim_SHPObject ossim_SHPAPI_CALL1(*)
- ossim_SHPCreateObject( int nSHPType, int nShapeId, int nParts,
- const int * panPartStart, const int * panPartType,
- int nVertices,
- const double * padfX, const double * padfY,
- const double * padfZ, const double * padfM );
-ossim_SHPObject ossim_SHPAPI_CALL1(*)
- ossim_SHPCreateSimpleObject( int nSHPType, int nVertices,
- const double * padfX,
- const double * padfY,
- const double * padfZ );
-
-int ossim_SHPAPI_CALL
- ossim_SHPRewindObject( ossim_SHPHandle hSHP, ossim_SHPObject * psObject );
-
-void ossim_SHPAPI_CALL ossim_SHPClose( ossim_SHPHandle hSHP );
-void ossim_SHPAPI_CALL ossim_SHPWriteHeader( ossim_SHPHandle hSHP );
-
-const char ossim_SHPAPI_CALL1(*)
- ossim_SHPTypeName( int nSHPType );
-const char ossim_SHPAPI_CALL1(*)
- ossim_SHPPartTypeName( int nPartType );
-
-/* -------------------------------------------------------------------- */
-/* Shape quadtree indexing API. */
-/* -------------------------------------------------------------------- */
-
-/* this can be two or four for binary or quad tree */
-#define MAX_SUBNODE 4
-
-typedef struct shape_tree_node
-{
- /* region covered by this node */
- double adfBoundsMin[4];
- double adfBoundsMax[4];
-
- /* list of shapes stored at this node. The papsShapeObj pointers
- or the whole list can be NULL */
- int nShapeCount;
- int *panShapeIds;
- ossim_SHPObject **papsShapeObj;
-
- int nSubNodes;
- struct shape_tree_node *apsSubNode[MAX_SUBNODE];
-
-} ossim_SHPTreeNode;
-
-typedef struct
-{
- ossim_SHPHandle hSHP;
-
- int nMaxDepth;
- int nDimension;
- int nTotalCount;
-
- ossim_SHPTreeNode *psRoot;
-} ossim_SHPTree;
-
-ossim_SHPTree ossim_SHPAPI_CALL1(*)
- ossim_SHPCreateTree( ossim_SHPHandle hSHP, int nDimension, int nMaxDepth,
- double *padfBoundsMin, double *padfBoundsMax );
-void ossim_SHPAPI_CALL
- ossim_SHPDestroyTree( ossim_SHPTree * hTree );
-
-int ossim_SHPAPI_CALL
- ossim_SHPWriteTree( ossim_SHPTree *hTree, const char * pszFilename );
-ossim_SHPTree ossim_SHPAPI_CALL
- ossim_SHPReadTree( const char * pszFilename );
-
-int ossim_SHPAPI_CALL
- ossim_SHPTreeAddObject( ossim_SHPTree * hTree, ossim_SHPObject * psObject );
-int ossim_SHPAPI_CALL
- ossim_SHPTreeAddShapeId( ossim_SHPTree * hTree, ossim_SHPObject * psObject );
-int ossim_SHPAPI_CALL
- ossim_SHPTreeRemoveShapeId( ossim_SHPTree * hTree, int nShapeId );
-
-void ossim_SHPAPI_CALL
- ossim_SHPTreeTrimExtraNodes( ossim_SHPTree * hTree );
-
-int ossim_SHPAPI_CALL1(*)
- ossim_SHPTreeFindLikelyShapes( ossim_SHPTree * hTree,
- double * padfBoundsMin,
- double * padfBoundsMax,
- int * );
-int ossim_SHPAPI_CALL
- ossim_SHPCheckBoundsOverlap( double *, double *, double *, double *, int );
-
-int ossim_SHPAPI_CALL1(*)
- ossim_SHPSearchDiskTree( FILE *fp,
- double *padfBoundsMin, double *padfBoundsMax,
- int *pnShapeCount );
-
-/************************************************************************/
-/* ossim_DBF Support. */
-/************************************************************************/
-typedef struct
-{
- FILE *fp;
-
- int nRecords;
-
- int nRecordLength;
- int nHeaderLength;
- int nFields;
- int *panFieldOffset;
- int *panFieldSize;
- int *panFieldDecimals;
- char *pachFieldType;
-
- char *pszHeader;
-
- int nCurrentRecord;
- int bCurrentRecordModified;
- char *pszCurrentRecord;
-
- int bNoHeader;
- int bUpdated;
-} ossim_DBFInfo;
-
-typedef ossim_DBFInfo * ossim_DBFHandle;
-
-typedef enum {
- FTString,
- FTInteger,
- FTDouble,
- FTLogical,
- FTInvalid
-} ossim_DBFFieldType;
-
-#define XBASE_FLDHDR_SZ 32
-
-ossim_DBFHandle ossim_SHPAPI_CALL
- ossim_DBFOpen( const char * pszDBFFile, const char * pszAccess );
-ossim_DBFHandle ossim_SHPAPI_CALL
- ossim_DBFCreate( const char * pszDBFFile );
-
-int ossim_SHPAPI_CALL
- ossim_DBFGetFieldCount( ossim_DBFHandle psDBF );
-int ossim_SHPAPI_CALL
- ossim_DBFGetRecordCount( ossim_DBFHandle psDBF );
-int ossim_SHPAPI_CALL
- ossim_DBFAddField( ossim_DBFHandle hDBF, const char * pszFieldName,
- ossim_DBFFieldType eType, int nWidth, int nDecimals );
-
-ossim_DBFFieldType ossim_SHPAPI_CALL
- ossim_DBFGetFieldInfo( ossim_DBFHandle psDBF, int iField,
- char * pszFieldName, int * pnWidth, int * pnDecimals );
-
-int ossim_SHPAPI_CALL
- ossim_DBFGetFieldIndex(ossim_DBFHandle psDBF, const char *pszFieldName);
-
-int ossim_SHPAPI_CALL
- ossim_DBFReadIntegerAttribute( ossim_DBFHandle hDBF, int iShape, int iField );
-double ossim_SHPAPI_CALL
- ossim_DBFReadDoubleAttribute( ossim_DBFHandle hDBF, int iShape, int iField );
-const char ossim_SHPAPI_CALL1(*)
- ossim_DBFReadStringAttribute( ossim_DBFHandle hDBF, int iShape, int iField );
-const char ossim_SHPAPI_CALL1(*)
- ossim_DBFReadLogicalAttribute( ossim_DBFHandle hDBF, int iShape, int iField );
-int ossim_SHPAPI_CALL
- ossim_DBFIsAttributeNULL( ossim_DBFHandle hDBF, int iShape, int iField );
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteIntegerAttribute( ossim_DBFHandle hDBF, int iShape, int iField,
- int nFieldValue );
-int ossim_SHPAPI_CALL
- ossim_DBFWriteDoubleAttribute( ossim_DBFHandle hDBF, int iShape, int iField,
- double dFieldValue );
-int ossim_SHPAPI_CALL
- ossim_DBFWriteStringAttribute( ossim_DBFHandle hDBF, int iShape, int iField,
- const char * pszFieldValue );
-int ossim_SHPAPI_CALL
- ossim_DBFWriteNULLAttribute( ossim_DBFHandle hDBF, int iShape, int iField );
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteLogicalAttribute( ossim_DBFHandle hDBF, int iShape, int iField,
- const char lFieldValue);
-int ossim_SHPAPI_CALL
- ossim_DBFWriteAttributeDirectly(ossim_DBFHandle psDBF, int hEntity, int iField,
- void * pValue );
-const char ossim_SHPAPI_CALL1(*)
- ossim_DBFReadTuple(ossim_DBFHandle psDBF, int hEntity );
-int ossim_SHPAPI_CALL
- ossim_DBFWriteTuple(ossim_DBFHandle psDBF, int hEntity, void * pRawTuple );
-
-ossim_DBFHandle ossim_SHPAPI_CALL
- ossim_DBFCloneEmpty(ossim_DBFHandle psDBF, const char * pszFilename );
-
-void ossim_SHPAPI_CALL
- ossim_DBFClose( ossim_DBFHandle hDBF );
-void ossim_SHPAPI_CALL
- ossim_DBFUpdateHeader( ossim_DBFHandle hDBF );
-char ossim_SHPAPI_CALL
- ossim_DBFGetNativeFieldType( ossim_DBFHandle hDBF, int iField );
-
-#ifdef __cplusplus
-}
-#endif
-
-#endif /* ndef _SHAPEFILE_H_INCLUDED */
diff --git a/Utilities/otbossim/src/ossim/elevation/ossimElevationShapeIdx.cpp b/Utilities/otbossim/src/ossim/elevation/ossimElevationShapeIdx.cpp
deleted file mode 100644
index da29ea1366..0000000000
--- a/Utilities/otbossim/src/ossim/elevation/ossimElevationShapeIdx.cpp
+++ /dev/null
@@ -1,275 +0,0 @@
-#include <ossim/elevation/ossimElevationShapeIdx.h>
-#include <ossim/vec/shapefil.h>
-
-ossimString ossimElevationShapeIdx::theDefaultName = "elevation_idx.shp";
-
-class ossimElevationShapeIdxPrivate : public ossimReferenced
-{
-public:
- ossimElevationShapeIdxPrivate()
- :theShapefileHandle(0),
- theDbfHandle(0)
- {
- }
- virtual ~ossimElevationShapeIdxPrivate()
- {
- shpClose();
- }
- ossim_SHPHandle shpOpen(const ossimFilename& file, bool overWriteFlag)
- {
- shpClose();
- theShapeFilename = file;
- if(overWriteFlag)
- {
- theShapefileHandle = ossim_SHPCreate(file.c_str(),
- ossim_SHPT_POLYGON);
- }
- else
- {
- theShapefileHandle = ossim_SHPOpen(file.c_str(),
- "r+b");
- if(!theShapefileHandle)
- {
- theShapefileHandle = ossim_SHPCreate(file.c_str(),
- ossim_SHPT_POLYGON);
- }
- }
- dbfOpen(overWriteFlag);
-
- return theShapefileHandle;
- }
- void shpClose()
- {
- if(theShapefileHandle)
- {
- ossim_SHPClose(theShapefileHandle);
- theShapefileHandle = 0;
- }
- }
- ossim_SHPHandle shpHandle()
- {
- return theShapefileHandle;
- }
-
- ossim_DBFHandle dbfHandle()
- {
- return theDbfHandle;
- }
-
-protected:
- void dbfOpen(bool overWriteFlag)
- {
- dbfClose();
-
- if(overWriteFlag)
- {
- theDbfHandle = ossim_DBFCreate(theShapeFilename.c_str());
- if(theDbfHandle)
- {
- ossim_DBFAddField(theDbfHandle, "FILENAME", FTString, 256, 0);
- ossim_DBFAddField(theDbfHandle, "MINZ", FTDouble, 32, 17);
- ossim_DBFAddField(theDbfHandle, "MAXZ", FTDouble, 32, 17);
- }
- }
- else
- {
- theDbfHandle = ossim_DBFOpen(theShapeFilename.c_str(), "r+b");
- if(!theDbfHandle)
- {
- theDbfHandle = ossim_DBFCreate(theShapeFilename.c_str());
- if(theDbfHandle)
- {
- ossim_DBFAddField(theDbfHandle, "FILENAME", FTString, 256, 0);
- ossim_DBFAddField(theDbfHandle, "MINZ", FTDouble, 32, 17);
- ossim_DBFAddField(theDbfHandle, "MAXZ", FTDouble, 32, 17);
- }
- }
- }
- }
- void dbfClose()
- {
- if(theDbfHandle)
- {
- ossim_DBFClose(theDbfHandle);
- theDbfHandle = 0;
- }
- }
- ossimFilename theShapeFilename;
- ossim_SHPHandle theShapefileHandle;
- ossim_DBFHandle theDbfHandle;
-};
-
-ossimElevationShapeIdx::ossimElevationShapeIdx(const ossimFilename& shapeFilename,
- bool overWriteFlag)
- :thePrivateData(0)
-{
- setFilename(shapeFilename, overWriteFlag);
-}
-
-ossimElevationShapeIdx::~ossimElevationShapeIdx()
-{
- if(thePrivateData)
- {
- delete thePrivateData;
- thePrivateData = 0;
- }
-}
-
-const ossimString& ossimElevationShapeIdx::defaultName()
-{
- return theDefaultName;
-}
-
-void ossimElevationShapeIdx::setFilename(const ossimFilename& shapeFilename, bool overWriteFlag)
-{
- if(thePrivateData)
- {
- delete thePrivateData;
- thePrivateData = 0;
- }
- if(shapeFilename == "") return;
- if(shapeFilename.isDir())
- {
- theShapeFilename = shapeFilename.dirCat(theDefaultName);
- }
- else if(shapeFilename.isFile())
- {
- theShapeFilename = shapeFilename;
- theShapeFilename.setExtension("shp");
- }
- else
- {
- theShapeFilename = shapeFilename.path();
-
- theShapeFilename = theShapeFilename.dirCat(theDefaultName);
- theShapeFilename.setExtension("shp");
- }
- thePrivateData = new ossimElevationShapeIdxPrivate();
-
- if(!thePrivateData->shpOpen(theShapeFilename, overWriteFlag))
- {
- thePrivateData = 0;
- throw(("ossimElevationShapeIdx::setFilename: Unable to setup index file " + theShapeFilename).c_str());
- }
- double minBounds[4];
- double maxBounds[4];
- int nEntities = 0;
- int shapeType = ossim_SHPT_NULL;
- ossim_SHPGetInfo(thePrivateData->shpHandle(),
- &nEntities,
- &shapeType,
- minBounds,
- maxBounds);
- if(nEntities < 0)
- {
- // set it up for the adds
- //
- theGlobalMinx = 180.0;
- theGlobalMiny = 90.0;
- theGlobalMaxx = -180.0;
- theGlobalMaxy = -90.0;
- }
- else
- {
- theGlobalMinx = minBounds[0];
- theGlobalMiny = minBounds[1];
- theGlobalMaxx = maxBounds[0];
- theGlobalMaxy = maxBounds[1];
- }
-}
-
-void ossimElevationShapeIdx::add(const ossimFilename& filename,
- ossim_float64 minx, // decimal degrees
- ossim_float64 miny, // decimal degrees
- ossim_float64 minz, // meters
- ossim_float64 maxx, // decimal degrees
- ossim_float64 maxy, // decimal degrees
- ossim_float64 maxz )// meters
-{
- if(!thePrivateData) return;
- int iRecord = ossim_DBFGetRecordCount(thePrivateData->dbfHandle());
- double* x = new double[8];
- double* y = new double[8];
- int nVertices = 4;
- x[0] = minx;
- x[1] = maxx;
- x[2] = maxx;
- x[3] = minx;
-
- y[0] = maxy;
- y[1] = maxy;
- y[2] = miny;
- y[3] = miny;
-
-
- if(minx < theGlobalMinx)
- {
- theGlobalMinx = minx;
- }
- if(miny < theGlobalMiny)
- {
- theGlobalMiny = miny;
- }
- if(maxx > theGlobalMaxx)
- {
- theGlobalMaxx = maxx;
- }
- if(maxy > theGlobalMaxy)
- {
- theGlobalMaxy = maxy;
- }
- ossim_SHPObject* obj = ossim_SHPCreateSimpleObject(ossim_SHPT_POLYGON,
- nVertices,
- x,
- y,
- 0);
-
- ossim_SHPWriteObject( thePrivateData->shpHandle(),
- -1,
- obj );
- ossim_DBFWriteStringAttribute(thePrivateData->dbfHandle(),
- iRecord,
- 0,
- filename.c_str());
- ossim_DBFWriteDoubleAttribute(thePrivateData->dbfHandle(),
- iRecord,
- 0,
- (double)minz);
- ossim_DBFWriteDoubleAttribute(thePrivateData->dbfHandle(),
- iRecord,
- 0,
- (double)maxz);
-}
-
-void ossimElevationShapeIdx::buildQuadTree(int maxDepth)
-{
- if(!thePrivateData) return;
-
- ossimFilename qTree = theShapeFilename;
-
- qTree.setExtension("qtree");
-
- if(thePrivateData->shpHandle())
- {
- double minBounds[4] = {0.0,0.0,0.0,0.0};
- double maxBounds[4] = {0.0,0.0,0.0,0.0};
-
- minBounds[0] = theGlobalMinx;
- minBounds[1] = theGlobalMiny;
- maxBounds[0] = theGlobalMaxx;
- maxBounds[1] = theGlobalMaxy;
-
- ossim_SHPTree* tree = ossim_SHPCreateTree(thePrivateData->shpHandle(),
- 2,
- maxDepth,
- minBounds,
- maxBounds);
- if(tree)
- {
- ossim_SHPWriteTree(tree,
- qTree.c_str());
-
- ossim_SHPDestroyTree(tree);
- }
- }
-}
diff --git a/Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileCutter.cpp b/Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileCutter.cpp
deleted file mode 100644
index 08bb370054..0000000000
--- a/Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileCutter.cpp
+++ /dev/null
@@ -1,323 +0,0 @@
-//*******************************************************************
-// Copyright (C) 2000 ImageLinks Inc.
-//
-// License: LGPL
-//
-// See LICENSE.txt file in the top level directory for more details.
-//
-// Author: Garrett Potts
-//
-//*************************************************************************
-// $Id: ossimEsriShapeFileCutter.cpp 9094 2006-06-13 19:12:40Z dburken $
-#include <ossim/imaging/ossimEsriShapeFileCutter.h>
-
-RTTI_DEF2(ossimEsriShapeFileCutter,
- "ossimEsriShapeFileCutter",
- ossimImageSourceFilter,
- ossimViewInterface);
-
-ossimEsriShapeFileCutter::ossimEsriShapeFileCutter()
- :ossimImageSourceFilter(),
- ossimViewInterface(),
- theEsriShapeFile(NULL),
- theMaskFilter(NULL)
-{
- theObject = this;
- theEsriShapeFile = new ossimEsriShapeFileFilter;
- theEsriShapeFile->setFillFlag(true);
-
- theMaskFilter = new ossimMaskFilter;
-
- theMaskFilter->connectMyInputTo(1,
- theEsriShapeFile);
-}
-
-ossimEsriShapeFileCutter::~ossimEsriShapeFileCutter()
-{
- if(theEsriShapeFile)
- {
- delete theEsriShapeFile;
- theEsriShapeFile = NULL;
- }
-
- if(theMaskFilter)
- {
- delete theMaskFilter;
- theMaskFilter = NULL;
- }
-}
-
-ossimRefPtr<ossimImageData> ossimEsriShapeFileCutter::getTile(
- const ossimIrect& rect,
- ossim_uint32 resLevel)
-{
- if(isSourceEnabled())
- {
- if(theMaskFilter)
- {
- return theMaskFilter->getTile(rect, resLevel);
- }
- }
- else
- {
- if(getInput())
- {
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput());
-
- if(input)
- {
- return input->getTile(rect, resLevel);
- }
- }
- }
-
- return NULL;
-}
-
-void ossimEsriShapeFileCutter::initialize()
-{
- if(theEsriShapeFile)
- {
- theEsriShapeFile->initialize();
- }
-
-
- if(theMaskFilter)
- {
- theMaskFilter->disconnectMyInput((ossim_int32)0);
- theMaskFilter->disconnectMyInput((ossim_int32)1);
- theMaskFilter->connectMyInputTo((int)0,
- getInput(0));
- theMaskFilter->connectMyInputTo(1,
- theEsriShapeFile);
- }
-}
-
-bool ossimEsriShapeFileCutter::loadState(const ossimKeywordlist& kwl,
- const char* prefix)
-{
- ossimString shapeFilePrefix = ossimString(prefix) + "shape_file_filter.";
- ossimString maskPrefix = ossimString(prefix) + "mask.";
-
- theEsriShapeFile->loadState(kwl, shapeFilePrefix.c_str());
- theMaskFilter->loadState(kwl, maskPrefix.c_str());
-
- return ossimImageSourceFilter::loadState(kwl, prefix);
-}
-
-bool ossimEsriShapeFileCutter::saveState(ossimKeywordlist& kwl,
- const char* prefix)const
-{
- ossimString shapeFilePrefix = ossimString(prefix) + "shape_file_filter.";
- ossimString maskPrefix = ossimString(prefix) + "mask.";
-
- theEsriShapeFile->saveState(kwl, shapeFilePrefix.c_str());
- theMaskFilter->saveState(kwl, maskPrefix.c_str());
-
- return ossimImageSourceFilter::saveState(kwl, prefix);
-}
-
-void ossimEsriShapeFileCutter::getDecimationFactor(ossim_uint32 resLevel,
- ossimDpt& result)const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface, getInput(0));
- if(input)
- {
- input->getDecimationFactor(resLevel, result);
- }
-}
-
-ossimIrect ossimEsriShapeFileCutter::getBoundingRect(ossim_uint32 resLevel)const
-{
- ossimIrect result;
- result.makeNan();
-
- if(theMaskFilter)
- {
- result = theMaskFilter->getBoundingRect(resLevel);
- }
-
- return result;
-}
-
-void ossimEsriShapeFileCutter::getDecimationFactors(
- vector<ossimDpt>& decimations)const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- input->getDecimationFactors(decimations);
- }
-}
-
-ossim_uint32 ossimEsriShapeFileCutter::getNumberOfDecimationLevels()const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- return input->getNumberOfDecimationLevels();
- }
- return 0;
-}
-
-ossim_uint32 ossimEsriShapeFileCutter::getNumberOfOutputBands() const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- return input->getNumberOfOutputBands();
- }
- return 0;
-}
-
-ossimScalarType ossimEsriShapeFileCutter::getOutputScalarType() const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- return input->getOutputScalarType();
- }
-
- return OSSIM_SCALAR_UNKNOWN;
-}
-
-void ossimEsriShapeFileCutter::getValidImageVertices(
- vector<ossimIpt>& validVertices,
- ossimVertexOrdering ordering,
- ossim_uint32 resLevel) const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- input->getValidImageVertices(validVertices, ordering, resLevel);
- }
-}
-
-ossim_uint32 ossimEsriShapeFileCutter::getTileWidth() const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- return input->getTileWidth();
- }
-
- return 0;
-}
-
-ossim_uint32 ossimEsriShapeFileCutter::getTileHeight() const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- return input->getTileHeight();
- }
- return 0;
-}
-
-double ossimEsriShapeFileCutter::getNullPixelValue(ossim_uint32 band)const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- return input->getNullPixelValue(band);
- }
- return 0;
-}
-
-double ossimEsriShapeFileCutter::getMinPixelValue(ossim_uint32 band)const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput(0));
- if(input)
- {
- return input->getMinPixelValue(band);
- }
- return 0;
-}
-
-double ossimEsriShapeFileCutter::getMaxPixelValue(ossim_uint32 band)const
-{
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface, getInput(0));
- if(input)
- {
- return input->getMaxPixelValue(band);
- }
- return 0;
-}
-
-void ossimEsriShapeFileCutter::loadFile(const ossimFilename& filename)
-{
- if(theEsriShapeFile)
- {
- theEsriShapeFile->loadShapeFile(filename);
- }
-}
-
-bool ossimEsriShapeFileCutter::setView(ossimObject* baseObject,
- bool ownsTheView)
-{
- if(theEsriShapeFile)
- {
- return theEsriShapeFile->setView(baseObject, ownsTheView);
- }
-
- return false;
-}
-
-ossimObject* ossimEsriShapeFileCutter::getView()
-{
- if(theEsriShapeFile)
- {
- return theEsriShapeFile->getView();
- }
-
- return NULL;;
-}
-
-const ossimObject* ossimEsriShapeFileCutter::getView()const
-{
- if(theEsriShapeFile)
- {
- return theEsriShapeFile->getView();
- }
-
- return NULL;;
-}
-
-ossimFilename ossimEsriShapeFileCutter::getFilename()const
-{
- if(theEsriShapeFile)
- {
- return theEsriShapeFile->getFilename();
- }
-
- return ossimFilename();
-}
-
-void ossimEsriShapeFileCutter::setMaskType(
- ossimMaskFilter::ossimFileSelectionMaskType type)
-{
- if(theMaskFilter)
- {
- theMaskFilter->setMaskType(type);
- }
-}
-
-ossimMaskFilter::ossimFileSelectionMaskType ossimEsriShapeFileCutter::getMaskType()const
-{
- if(theMaskFilter)
- {
- return theMaskFilter->getMaskType();
- }
-
- return ossimMaskFilter::OSSIM_MASK_TYPE_SELECT;
-}
diff --git a/Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileFilter.cpp b/Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileFilter.cpp
deleted file mode 100644
index d47c07b7e0..0000000000
--- a/Utilities/otbossim/src/ossim/imaging/ossimEsriShapeFileFilter.cpp
+++ /dev/null
@@ -1,813 +0,0 @@
-//*******************************************************************
-// Copyright (C) 2000 ImageLinks Inc.
-//
-// License: See LICENSE.txt file in the top level directory.
-//
-// Author: Garrett Potts
-//
-//*************************************************************************
-// $Id: ossimEsriShapeFileFilter.cpp 11411 2007-07-27 13:53:51Z dburken $
-
-#include <cstdio>
-#include <cstdlib>
-#include <sstream>
-using namespace std;
-
-#include <ossim/imaging/ossimEsriShapeFileFilter.h>
-#include <ossim/vec/ossimShapeFile.h>
-#include <ossim/imaging/ossimAnnotationPolyObject.h>
-#include <ossim/imaging/ossimGeoAnnotationPolyLineObject.h>
-#include <ossim/imaging/ossimAnnotationObject.h>
-#include <ossim/imaging/ossimGeoAnnotationPolyObject.h>
-#include <ossim/imaging/ossimGeoAnnotationEllipseObject.h>
-#include <ossim/base/ossimTrace.h>
-#include <ossim/base/ossimCommon.h>
-#include <ossim/base/ossimConstants.h>
-#include <ossim/base/ossimKeywordNames.h>
-#include <ossim/base/ossimTrace.h>
-#include <ossim/base/ossimNotifyContext.h>
-#include <ossim/base/ossimFilename.h>
-#include <ossim/base/ossimGeoPolygon.h>
-#include <ossim/base/ossimUnitConversionTool.h>
-#include <ossim/projection/ossimProjection.h>
-#include <ossim/projection/ossimProjectionFactoryRegistry.h>
-#include <ossim/projection/ossimMapProjection.h>
-
-
-RTTI_DEF2(ossimEsriShapeFileFilter,
- "ossimEsriShapeFileFilter",
- ossimAnnotationSource,
- ossimViewInterface);
-
-static const ossimTrace traceDebug("ossimEsriShapeFileFilter:debug");
-
-ossimEsriShapeFileFilter::ossimEsriShapeFileFilter(ossimImageSource* inputSource)
- :ossimAnnotationSource(inputSource),
- ossimViewInterface(),
- theViewProjection(NULL),
- theShapeFileProjection(NULL),
- theOwnsViewProjectionFlag(true),
- theCoordinateSystem(OSSIM_GEOGRAPHIC_SPACE),
- theUnitType(OSSIM_METERS),
- theTree((ossim_SHPTree*)NULL),
- theMaxQuadTreeLevels(10),
- thePenColor(255,255,255),
- theBrushColor(255,255,255),
- theFillFlag(false),
- theThickness(1),
- thePointWidthHeight(1, 1),
- theBorderSize(0.0),
- theBorderSizeUnits(OSSIM_DEGREES)
-{
- ossimViewInterface::theObject = this;
- ossimAnnotationSource::setNumberOfBands(3);
- theBoundingRect.makeNan();
- theMinArray[0] = theMinArray[1] = theMinArray[2] = theMinArray[3] = ossim::nan();
- theMaxArray[0] = theMaxArray[1] = theMaxArray[2] = theMaxArray[3] = ossim::nan();
-}
-
-ossimEsriShapeFileFilter::~ossimEsriShapeFileFilter()
-{
- removeViewProjection();
-
- if(theTree)
- {
- ossim_SHPDestroyTree(theTree);
- }
-
- deleteCache();
-}
-
-bool ossimEsriShapeFileFilter::setView(ossimObject* baseObject,
- bool ownsTheView)
-{
- ossimProjection* proj = PTR_CAST(ossimProjection, baseObject);
- if(proj)
- {
- setProjection(proj, ownsTheView);
-
- return true;
- }
- else if(!baseObject)
- {
- setProjection((ossimProjection*)NULL, true);
- }
-
- return false;
-}
-
-ossimObject* ossimEsriShapeFileFilter::getView()
-{
- return theViewProjection;
-}
-
-const ossimObject* ossimEsriShapeFileFilter::getView()const
-{
- return theViewProjection;
-}
-
-bool ossimEsriShapeFileFilter::addObject(ossimAnnotationObject* anObject)
-{
- ossimNotify(ossimNotifyLevel_WARN)
- << "ossimEsriShapeFileFilter::addObject\n"
- << "Can't add objects to layer, must go through Esri loadShapeFile"
- <<endl;
-
- return false;
-}
-
-void ossimEsriShapeFileFilter::computeBoundingRect()
-{
-// ossimAnnotationSource::computeBoundingRect();
-
- std::multimap<int, ossimAnnotationObject*>::iterator iter = theShapeCache.begin();
-
- theBoundingRect.makeNan();
- while(iter != theShapeCache.end())
- {
- ossimDrect rect = (*iter).second->getBoundingRect();
- if(theBoundingRect.hasNans())
- {
- theBoundingRect = rect;
- }
- else
- {
- if(!rect.hasNans())
- {
- theBoundingRect = theBoundingRect.combine(rect);
- }
- }
-
- ++iter;
- }
-}
-
-ossimIrect ossimEsriShapeFileFilter::getBoundingRect(ossim_uint32 resLevel)const
-{
- if(!isSourceEnabled()||
- getInput())
- {
- if(getInput())
- {
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface,
- getInput());
- if(input)
- {
- return input->getBoundingRect(resLevel);
- }
- }
- }
- return theBoundingRect;
-}
-
-void ossimEsriShapeFileFilter::drawAnnotations(ossimRefPtr<ossimImageData> tile)
-{
- ossimAnnotationSource::drawAnnotations(tile);
-
- if (!theTree||!theShapeFile.isOpen()) return;
- if(theViewProjection)
- {
- ossimIrect rect = tile->getImageRectangle();
-
- rect = ossimIrect(rect.ul().x,
- rect.ul().y,
- rect.lr().x,
- rect.lr().y);
- double boundsMin[2];
- double boundsMax[2];
-
- ossimGpt gp1;
- ossimGpt gp2;
- ossimGpt gp3;
- ossimGpt gp4;
-
- theViewProjection->lineSampleToWorld(rect.ul(),
- gp1);
- theViewProjection->lineSampleToWorld(rect.ur(),
- gp2);
- theViewProjection->lineSampleToWorld(rect.lr(),
- gp3);
- theViewProjection->lineSampleToWorld(rect.ll(),
- gp4);
-
- ossimDrect boundsRect( ossimDpt(gp1.lond(),
- gp1.latd()),
- ossimDpt(gp2.lond(),
- gp2.latd()),
- ossimDpt(gp3.lond(),
- gp3.latd()),
- ossimDpt(gp4.lond(),
- gp4.latd()),
- OSSIM_RIGHT_HANDED);
-
- boundsMin[0] = boundsRect.ul().x;
- boundsMin[1] = boundsRect.lr().y;
-
- boundsMax[0] = boundsRect.lr().x;
- boundsMax[1] = boundsRect.ul().y;
-
- int n;
- int *array=(int*)NULL;
-
- array = ossim_SHPTreeFindLikelyShapes(theTree,
- boundsMin,
- boundsMax,
- &n);
-
- theImage->setCurrentImageData(tile);
- if(n&&array)
- {
- for(int i = 0; i < n; ++i)
- {
- std::multimap<int, ossimAnnotationObject*>::iterator iter = theShapeCache.find(array[i]);
- while( ((*iter).first == array[i]) && (iter != theShapeCache.end()) )
- {
- (*iter).second->draw(*theImage);
- ++iter;
- }
- }
-
- free(array);
- }
- }
-}
-
-void ossimEsriShapeFileFilter::transformObjects(ossimProjection* projection)
-{
- std::multimap<int, ossimAnnotationObject*>::iterator iter = theShapeCache.begin();
-
-
- ossimProjection* proj = theViewProjection;
- if(projection)
- {
- proj = projection;
- }
-
- if(!proj) return;
-
- while(iter != theShapeCache.end())
- {
- ossimGeoAnnotationObject* obj = PTR_CAST(ossimGeoAnnotationObject,
- (*iter).second);
- if(obj)
- {
- obj->transform(proj);
- }
- ++iter;
- }
-
- computeBoundingRect();
-}
-
-void ossimEsriShapeFileFilter::setProjection(ossimProjection* projection,
- bool ownsProjectionFlag)
-{
- removeViewProjection();
-
- theViewProjection = projection;
- theOwnsViewProjectionFlag = ownsProjectionFlag;
-
- transformObjects();
-}
-
-bool ossimEsriShapeFileFilter::getImageGeometry(ossimKeywordlist& kwl,
- const char* prefix)
-{
- if(theViewProjection)
- {
- return theViewProjection->saveState(kwl, prefix);
- }
- else
- {
- ossimImageSourceInterface* input = PTR_CAST(ossimImageSourceInterface, getInput(0));
- if(input)
- {
- return input->getImageGeometry(kwl, prefix);
- }
- else
- {
- kwl.add(prefix,
- ossimKeywordNames::TYPE_KW,
- "ossimEquDistCylProjection",
- true);
- kwl.add(prefix,
- ossimKeywordNames::DECIMAL_DEGREES_PER_PIXEL_LAT,
- .1,
- true);
- kwl.add(prefix,
- ossimKeywordNames::DECIMAL_DEGREES_PER_PIXEL_LON,
- .1,
- true);
- kwl.add(prefix,
- ossimKeywordNames::CENTRAL_MERIDIAN_KW,
- 0.0,
- true);
- kwl.add(prefix,
- ossimKeywordNames::ORIGIN_LATITUDE_KW,
- 0.0,
- true);
- kwl.add(prefix,
- ossimKeywordNames::DATUM_KW,
- "WGE",
- true);
- return true;
- }
- }
-
- return false;
-}
-
-void ossimEsriShapeFileFilter::removeViewProjection()
-{
- if(theViewProjection&&theOwnsViewProjectionFlag)
- {
- delete theViewProjection;
- }
-
- theViewProjection = NULL;
-}
-
-void ossimEsriShapeFileFilter::deleteCache()
-{
- std::multimap<int, ossimAnnotationObject*>::iterator iter = theShapeCache.begin();
-
-
- while(iter != theShapeCache.end())
- {
- if((*iter).second)
- {
- delete (*iter).second;
- }
- ++iter;
- }
-
- theShapeCache.clear();
-}
-
-void ossimEsriShapeFileFilter::checkAndSetDefaultView()
-{
- if(!theViewProjection)
- {
- ossimKeywordlist tempKwl;
- getImageGeometry(tempKwl);
-
- theViewProjection = ossimProjectionFactoryRegistry::instance()->createProjection(tempKwl);
-
- if(theViewProjection)
- {
- transformObjects();
- }
- }
-}
-
-bool ossimEsriShapeFileFilter::loadShapeFile(const ossimFilename& shapeFile)
-{
- if(theTree)
- {
- ossim_SHPDestroyTree(theTree);
- theTree = (ossim_SHPTree*)NULL;
- }
- theShapeFile.open(shapeFile);
- deleteCache();
- deleteAll();
-
- if(theShapeFile.isOpen())
- {
- theShapeFile.getBounds(theMinArray[0],theMinArray[1],theMinArray[2],theMinArray[3],
- theMaxArray[0],theMaxArray[1],theMaxArray[2],theMaxArray[3]);
-
- theTree = ossim_SHPCreateTree(theShapeFile.getHandle(),
- 2,
- theMaxQuadTreeLevels,
- theMinArray,
- theMaxArray);
-
- ossimShapeObject obj;
- for(int index = 0 ; index < theShapeFile.getNumberOfShapes(); ++index)
- {
- obj.loadShape(theShapeFile,
- index);
-
- if(obj.isLoaded())
- {
- switch(obj.getType())
- {
- case ossim_SHPT_POLYGON:
- case ossim_SHPT_POLYGONZ:
- {
- loadPolygon(obj);
- break;
- }
- case ossim_SHPT_POINT:
- case ossim_SHPT_POINTZ:
- {
- loadPoint(obj);
- break;
- }
- case ossim_SHPT_ARC:
- case ossim_SHPT_ARCZ:
- {
- loadArc(obj);
- break;
- }
- case ossim_SHPT_NULL:
- {
- break;
- }
- default:
- {
- ossimNotify(ossimNotifyLevel_WARN)
- << "ossimEsriShapeFileFilter::loadShapeFile\n"
- << "SHAPE " << obj.getTypeByName()
- << " Not supported" << endl;
- break;
- }
- }
- }
- }
-
- theCurrentObject = theShapeCache .begin();
- if(theViewProjection)
- {
- transformObjects();
- }
- else
- {
- checkAndSetDefaultView();
- }
- }
-
- return true;
-}
-
-void ossimEsriShapeFileFilter::loadPolygon(ossimShapeObject& obj)
-{
- int starti = 0;
- int endi = 0;
- if(obj.getNumberOfParts() > 1)
- {
- starti = obj.getShapeObject()->panPartStart[0];
- endi = obj.getShapeObject()->panPartStart[1];
- }
- else
- {
- starti = 0;
- endi = obj.getShapeObject()->nVertices;
- }
-
- vector<ossimGpt> groundPolygon;
- for(ossim_uint32 part = 0; part < obj.getNumberOfParts(); ++part)
- {
- if(obj.getPartType(part) != ossim_SHPP_RING)
- {
- ossimNotify(ossimNotifyLevel_WARN)
- << "ossimEsriShapeFileFilter::loadPolygon\n"
- << "Part = " << obj.getPartByName(part)
- << " not supported for shape = "
- << obj.getTypeByName() << endl;
- break;
- }
- groundPolygon.clear();
- for(ossim_int32 vertexNumber = starti; vertexNumber < endi; ++vertexNumber)
- {
- groundPolygon.push_back(ossimGpt(obj.getShapeObject()->padfY[vertexNumber],
- obj.getShapeObject()->padfX[vertexNumber]));
-
- }
- starti = endi;
- if((part + 2) < obj.getNumberOfParts())
- {
- endi = obj.getShapeObject()->panPartStart[part+2];
- }
- else
- {
- endi = obj.getShapeObject()->nVertices;
- }
-
- ossimRgbVector color;
-
- if(theFillFlag)
- {
- color = theBrushColor;
- }
- else
- {
- color = thePenColor;
- }
-
- if(theBorderSize != 0.0)
- {
- ossimGeoPolygon tempPoly(groundPolygon);
- ossimGeoPolygon tempPoly2;
-
- tempPoly.stretchOut(tempPoly2,
- theBorderSize);
- groundPolygon = tempPoly2.getVertexList();
-
- }
-
- ossimGeoAnnotationObject *newGeoObj = new ossimGeoAnnotationPolyObject(groundPolygon,
- theFillFlag,
- color.getR(),
- color.getG(),
- color.getB(),
- theThickness);
- newGeoObj->setName(theFeatureName);
- theShapeCache.insert(make_pair(obj.getId(),
- newGeoObj));
- }
-}
-
-void ossimEsriShapeFileFilter::loadPoint(ossimShapeObject& obj)
-{
- int n = obj.getNumberOfVertices();
-
- if(n)
- {
- ossimGpt gpt(obj.getShapeObject()->padfY[0],
- obj.getShapeObject()->padfX[0]);
-
- ossimRgbVector color;
-
- if(theFillFlag)
- {
- color = theBrushColor;
- }
- else
- {
- color = thePenColor;
- }
- ossimGeoAnnotationEllipseObject *newGeoObj =
- new ossimGeoAnnotationEllipseObject(gpt,
- thePointWidthHeight,
- theFillFlag,
- color.getR(),
- color.getG(),
- color.getB(),
- theThickness);
- newGeoObj->setEllipseWidthHeightUnitType(OSSIM_PIXEL);
- newGeoObj->setName(theFeatureName);
- theShapeCache.insert(make_pair(obj.getId(),
- newGeoObj));
-
- }
-}
-
-void ossimEsriShapeFileFilter::loadArc(ossimShapeObject& obj)
-{
- int starti = 0;
- int endi = 0;
- if(obj.getNumberOfParts() > 1)
- {
- starti = obj.getShapeObject()->panPartStart[0];
- endi = obj.getShapeObject()->panPartStart[1];
- }
- else
- {
- starti = 0;
- endi = obj.getShapeObject()->nVertices;
- }
-
- vector<ossimGpt> groundPolygon;
- for(ossim_uint32 part = 0; part < obj.getNumberOfParts(); ++part)
- {
- if(obj.getPartType(part) != ossim_SHPP_RING)
- {
- ossimNotify(ossimNotifyLevel_WARN)
- << "ossimEsriShapeFileFilter::loadArc\n"
- << "Part = " << obj.getPartByName(part)
- << " not supported for shape = "
- << obj.getTypeByName() << endl;
- break;
- }
- groundPolygon.clear();
- for(ossim_int32 vertexNumber = starti; vertexNumber < endi; ++vertexNumber)
- {
- groundPolygon.push_back(ossimGpt(obj.getShapeObject()->padfY[vertexNumber],
- obj.getShapeObject()->padfX[vertexNumber]));
-
- }
- starti = endi;
- if((part + 2) < obj.getNumberOfParts())
- {
- endi = obj.getShapeObject()->panPartStart[part+2];
- }
- else
- {
- endi = obj.getShapeObject()->nVertices;
- }
-
- ossimRgbVector color;
-
- if(theFillFlag)
- {
- color = theBrushColor;
- }
- else
- {
- color = thePenColor;
- }
-
- ossimGeoAnnotationObject *newGeoObj = new ossimGeoAnnotationPolyLineObject(groundPolygon,
- color.getR(),
- color.getG(),
- color.getB(),
- theThickness);
- newGeoObj->setName(theFeatureName);
- theShapeCache.insert(make_pair(obj.getId(),
- newGeoObj));
- }
-}
-
-bool ossimEsriShapeFileFilter::saveState(ossimKeywordlist& kwl,
- const char* prefix)const
-{
- ossimString s;
-
- kwl.add(prefix,
- ossimKeywordNames::FILENAME_KW,
- theShapeFile.getFilename(),
- true);
-
- kwl.add(prefix,
- ossimKeywordNames::MAX_QUADTREE_LEVELS_KW,
- theMaxQuadTreeLevels,
- true);
-
- s = ossimString::toString((int)thePenColor.getR()) + " " +
- ossimString::toString((int)thePenColor.getG()) + " " +
- ossimString::toString((int)thePenColor.getB());
-
- kwl.add(prefix,
- ossimKeywordNames::PEN_COLOR_KW,
- s.c_str(),
- true);
-
- s = ossimString::toString((int)theBrushColor.getR()) + " " +
- ossimString::toString((int)theBrushColor.getG()) + " " +
- ossimString::toString((int)theBrushColor.getB());
-
- kwl.add(prefix,
- ossimKeywordNames::BRUSH_COLOR_KW,
- s.c_str(),
- true);
-
- kwl.add(prefix,
- ossimKeywordNames::FILL_FLAG_KW,
- (int)theFillFlag,
- true);
-
- kwl.add(prefix,
- ossimKeywordNames::FEATURE_NAME_KW,
- theFeatureName.c_str(),
- true);
-
- kwl.add(prefix,
- ossimKeywordNames::THICKNESS_KW,
- theThickness,
- true);
-
- ossimString border;
- border = ossimString::toString(theBorderSize);
- border += " degrees";
- kwl.add(prefix,
- ossimKeywordNames::BORDER_SIZE_KW,
- border,
- true);
-
- s = ossimString::toString((int)thePointWidthHeight.x) + " " +
- ossimString::toString((int)thePointWidthHeight.y) + " ";
-
- kwl.add(prefix,
- ossimKeywordNames::POINT_WIDTH_HEIGHT_KW,
- s.c_str(),
- true);
-
- if(theViewProjection)
- {
- ossimString newPrefix = prefix;
- newPrefix += "view_proj.";
- theViewProjection->saveState(kwl,
- newPrefix.c_str());
- }
-
- return ossimAnnotationSource::saveState(kwl, prefix);
-}
-
-bool ossimEsriShapeFileFilter::loadState(const ossimKeywordlist& kwl,
- const char* prefix)
-{
-
- const char* quadLevels = kwl.find(prefix, ossimKeywordNames::MAX_QUADTREE_LEVELS_KW);
- const char* filename = kwl.find(prefix, ossimKeywordNames::FILENAME_KW);
- const char* penColor = kwl.find(prefix, ossimKeywordNames::PEN_COLOR_KW);
- const char* brushColor = kwl.find(prefix, ossimKeywordNames::BRUSH_COLOR_KW);
- const char* featureName = kwl.find(prefix, ossimKeywordNames::FEATURE_NAME_KW);
- const char* fillFlag = kwl.find(prefix, ossimKeywordNames::FILL_FLAG_KW);
- const char* thickness = kwl.find(prefix, ossimKeywordNames::THICKNESS_KW);
- const char* pointWh = kwl.find(prefix, ossimKeywordNames::POINT_WIDTH_HEIGHT_KW);
- const char* border_size = kwl.find(prefix, ossimKeywordNames::BORDER_SIZE_KW);
-
- deleteCache();
-
- if(thickness)
- {
- theThickness = ossimString(thickness).toLong();
- }
- if(quadLevels)
- {
- theMaxQuadTreeLevels = ossimString(quadLevels).toLong();
- }
-
- if(penColor)
- {
- int r, g, b;
- istringstream s(penColor);
- s>>r>>g>>b;
- thePenColor = ossimRgbVector((ossim_uint8)r, (ossim_uint8)g, (ossim_uint8)b);
- }
-
- if(brushColor)
- {
- int r, g, b;
- istringstream s(brushColor);
- s>>r>>g>>b;
- theBrushColor = ossimRgbVector((ossim_uint8)r, (ossim_uint8)g, (ossim_uint8)b);
- }
- if(pointWh)
- {
- double w, h;
- istringstream s(pointWh);
- s>>w>>h;
- thePointWidthHeight = ossimDpt(w, h);
- }
-
- if(fillFlag)
- {
- theFillFlag = ossimString(fillFlag).toBool();
- }
-
- if(border_size)
- {
- istringstream input(border_size);
-
- ossimString s;
- input >> s;
-
- theBorderSize = s.toDouble();
-
- ossimString s2;
-
- input >> s2;
-
- s2 = s2.upcase();
-
- if(s2 == "US")
- {
- theBorderSizeUnits = OSSIM_US_SURVEY_FEET;
- }
- else if(s2 == "METERS")
- {
- theBorderSizeUnits = OSSIM_METERS;
- }
- else if(s2 == "FEET")
- {
- theBorderSizeUnits = OSSIM_FEET;
- }
- else
- {
- theBorderSizeUnits = OSSIM_DEGREES;
- }
- ossimUnitConversionTool unitConvert(theBorderSize,
- theBorderSizeUnits);
-
- theBorderSize = unitConvert.getValue(OSSIM_DEGREES);
- theBorderSizeUnits = OSSIM_DEGREES;
- }
- else
- {
- theBorderSize = 0.0;
- theBorderSizeUnits = OSSIM_DEGREES;
- }
-
- theFeatureName = featureName;
-
- removeViewProjection();
-
- ossimString newPrefix = prefix;
- newPrefix += "view_proj.";
-
- theViewProjection = ossimProjectionFactoryRegistry::instance()->createProjection(kwl,
- newPrefix.c_str());
- theOwnsViewProjectionFlag = true;
-
- if(filename)
- {
- loadShapeFile(ossimFilename(filename));
- }
-
- checkAndSetDefaultView();
-
-
- return ossimAnnotationSource::loadState(kwl, prefix);
-}
diff --git a/Utilities/otbossim/src/ossim/imaging/ossimImageShapeFileIndex.cpp b/Utilities/otbossim/src/ossim/imaging/ossimImageShapeFileIndex.cpp
deleted file mode 100644
index 90dd48d374..0000000000
--- a/Utilities/otbossim/src/ossim/imaging/ossimImageShapeFileIndex.cpp
+++ /dev/null
@@ -1,197 +0,0 @@
-#include <ossim/imaging/ossimImageShapeFileIndex.h>
-#include <ossim/imaging/ossimImageHandlerRegistry.h>
-#include <ossim/imaging/ossimImageHandler.h>
-#include <ossim/projection/ossimProjectionFactoryRegistry.h>
-#include <ossim/projection/ossimProjection.h>
-#include <ossim/base/ossimTrace.h>
-
-#include <ossim/vec/shapefil.h>
-
-static ossimTrace traceDebug("ossimImageShapeFileIndex:debug");
-
-ossimImageShapeFileIndex::ossimImageShapeFileIndex(const ossimFilename& shapeFilename)
- :theGlobalMinx(1.0/DBL_EPSILON),
- theGlobalMiny(1.0/DBL_EPSILON),
- theGlobalMaxx(-1.0/DBL_EPSILON),
- theGlobalMaxy(-1.0/DBL_EPSILON)
-{
- setFilename(shapeFilename);
-}
-
-void ossimImageShapeFileIndex::setFilename(const ossimFilename& shapeFilename)
-{
- theFilename = shapeFilename;
- theFilename = theFilename.setExtension("");
-}
-
-void ossimImageShapeFileIndex::add(const ossimFilename& filename,
- const ossim_int32 entry, // default to all entries
- const ossimString& groupName)
-{
- if(traceDebug())
- {
- ossimNotify(ossimNotifyLevel_DEBUG) << "ossimImageShapeFileIndex::add: Entered ...." << std::endl;
- }
- ossimRefPtr<ossimImageHandler> h = ossimImageHandlerRegistry::instance()->open(filename);
- ossimRefPtr<ossimProjection> proj;
-
- if(h.valid())
- {
- std::vector<ossim_uint32> entryList;
- if(entry < 0)
- {
- h->getEntryList(entryList);
- }
- else
- {
- entryList.push_back(static_cast<ossim_uint32>(entry));
- }
-
- ossim_uint32 idx = 0;
- for(idx = 0; idx < entryList.size(); ++idx)
- {
- ossimKeywordlist kwl;
- if(h->setCurrentEntry(entryList[idx]))
- {
- ossimIrect bounds = h->getBoundingRect();
-
- h->getImageGeometry(kwl);
- proj = ossimProjectionFactoryRegistry::instance()->createProjection(kwl);
-
- if(proj.valid())
- {
- if(traceDebug())
- {
- ossimNotify(ossimNotifyLevel_DEBUG) << "Adding file: " << filename << std::endl
- << " entry: " << entryList[idx] << std::endl;
- }
- std::vector<ossimDpt> points;
- ossim_float64 minx;
- ossim_float64 miny;
- ossim_float64 maxx;
- ossim_float64 maxy;
- ossimGpt ulgpt;
- ossimGpt urgpt;
- ossimGpt lrgpt;
- ossimGpt llgpt;
- const ossimDatum* datum = ossimDatumFactory::instance()->wgs84();
- proj->lineSampleToWorld(bounds.ul(), ulgpt);
- proj->lineSampleToWorld(bounds.ur(), urgpt);
- proj->lineSampleToWorld(bounds.lr(), lrgpt);
- proj->lineSampleToWorld(bounds.ll(), llgpt);
-
- ulgpt.changeDatum(datum);
- urgpt.changeDatum(datum);
- lrgpt.changeDatum(datum);
- llgpt.changeDatum(datum);
-
- points.push_back(ossimDpt(ulgpt));
- points.push_back(ossimDpt(urgpt));
- points.push_back(ossimDpt(lrgpt));
- points.push_back(ossimDpt(llgpt));
-
- ossimDrect rect(points, OSSIM_RIGHT_HANDED);
-
- minx = rect.ll().x;
- miny = rect.ll().y;
- maxx = rect.ur().x;
- maxy = rect.ur().y;
-
- if(minx < theGlobalMinx) theGlobalMinx = minx;
- if(maxx > theGlobalMaxx) theGlobalMaxx = maxx;
- if(miny < theGlobalMiny) theGlobalMiny = miny;
- if(maxy > theGlobalMaxy) theGlobalMaxy = maxy;
-
- theRecords.push_back(new ossimImageShapeFileIndexRecord(filename,
- entryList[idx],
- groupName,
- minx,
- miny,
- maxx,
- maxy));
- }
- }
- }
- }
- if(traceDebug())
- {
- ossimNotify(ossimNotifyLevel_DEBUG) << "ossimImageShapeFileIndex::add: Leaving ...." << std::endl;
- }
-}
-
-void ossimImageShapeFileIndex::build(ossim_uint32 maxDepth)const
-{
- if(theRecords.size())
- {
- ossim_SHPHandle hSHP = ossim_SHPCreate( theFilename.c_str(), ossim_SHPT_POLYGON);
- ossim_DBFHandle hDBF = ossim_DBFCreate( theFilename.c_str());
-
- ossim_DBFAddField(hDBF, "GROUP", FTString, 256, 0);
- ossim_DBFAddField(hDBF, "FILENAME", FTString, 256, 0);
- ossim_DBFAddField(hDBF, "ENTRY", FTInteger, 10, 0);
- ossim_DBFAddField(hDBF, "DATUM", FTString, 16, 0);
-
- double* x = new double[4];
- double* y = new double[4];
- int nVertices = 4;
- ossim_uint32 idx = 0;
- int iRecord = 0;
- for(idx = 0; idx < theRecords.size(); ++idx)
- {
- x[0] = theRecords[idx]->theMinx;
- x[1] = theRecords[idx]->theMaxx;
- x[2] = theRecords[idx]->theMaxx;
- x[3] = theRecords[idx]->theMinx;
-
- y[0] = theRecords[idx]->theMaxy;
- y[1] = theRecords[idx]->theMaxy;
- y[2] = theRecords[idx]->theMiny;
- y[3] = theRecords[idx]->theMiny;
-
- ossim_SHPObject* obj = ossim_SHPCreateSimpleObject(ossim_SHPT_POLYGON,
- nVertices,
- x,
- y,
- 0);
- ossim_SHPWriteObject( hSHP, -1, obj );
- iRecord = ossim_DBFGetRecordCount( hDBF );
- ossim_DBFWriteStringAttribute(hDBF, iRecord, 0, theRecords[idx]->theGroupName.c_str());
- ossim_DBFWriteStringAttribute(hDBF, iRecord, 1, theRecords[idx]->theFilename.c_str());
- ossim_DBFWriteIntegerAttribute(hDBF, iRecord, 2, (int)(theRecords[idx]->theEntry));
- ossim_DBFWriteStringAttribute(hDBF, iRecord, 3, "WGE");
- ossim_SHPDestroyObject(obj);
- }
-
- ossim_DBFClose(hDBF);
- ossim_SHPTree * shpTree = 0;
- double minV[4] = {0.0,0.0,0.0,0.0};
- double maxV[4] = {0.0,0.0,0.0,0.0};
-
- minV[0] = theGlobalMinx;
- minV[1] = theGlobalMiny;
-
- maxV[0] = theGlobalMaxx;
- maxV[1] = theGlobalMaxy;
- if(maxDepth>1)
- {
- if(traceDebug())
- {
- ossimNotify(ossimNotifyLevel_DEBUG) << "building quadtree with " << maxDepth << " levels." << std::endl
- << "minx = " << theGlobalMinx << std::endl
- << "miny = " << theGlobalMiny << std::endl
- << "maxx = " << theGlobalMaxx << std::endl
- << "maxy = " << theGlobalMaxy << std::endl;
- }
- shpTree = ossim_SHPCreateTree(hSHP, 2, maxDepth, minV, maxV);
- }
-
- if(shpTree)
- {
- ossim_SHPWriteTree(shpTree, (theFilename+".qtree").c_str());
- ossim_SHPDestroyTree(shpTree);
- }
- ossim_SHPClose(hSHP);
- }
-
-}
-
diff --git a/Utilities/otbossim/src/ossim/imaging/ossimImageSourceInterface.cpp b/Utilities/otbossim/src/ossim/imaging/ossimImageSourceInterface.cpp
deleted file mode 100644
index 51a2e25dfe..0000000000
--- a/Utilities/otbossim/src/ossim/imaging/ossimImageSourceInterface.cpp
+++ /dev/null
@@ -1,18 +0,0 @@
-//*******************************************************************
-// Copyright (C) 2000 ImageLinks Inc.
-//
-// License: LGPL
-//
-// See LICENSE.txt file in the top level directory for more details.
-//
-// Author: Garrett Potts
-//
-//*************************************************************************
-// $Id: ossimImageSourceInterface.cpp 9094 2006-06-13 19:12:40Z dburken $
-#include <ossim/imaging/ossimImageSourceInterface.h>
-
-RTTI_DEF(ossimImageSourceInterface, "ossimImageSourceInterface");
-
-ossimImageSourceInterface::~ossimImageSourceInterface()
-{
-}
diff --git a/Utilities/otbossim/src/ossim/vec/dbfopen.c b/Utilities/otbossim/src/ossim/vec/dbfopen.c
deleted file mode 100644
index 1ec330cd4f..0000000000
--- a/Utilities/otbossim/src/ossim/vec/dbfopen.c
+++ /dev/null
@@ -1,1552 +0,0 @@
-/******************************************************************************
- * $Id: dbfopen.c 9094 2006-06-13 19:12:40Z dburken $
- *
- * Project: Shapelib
- * Purpose: Implementation of .dbf access API documented in dbf_api.html.
- * 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 LICENSE.LGPL). 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$
- * Revision 1.6 2006/06/13 19:11:32 dburken
- * Global include path change.
- *
- * Revision 1.5 2005/10/08 12:57:26 gpotts
- * Added adiitional patches
- *
- * Revision 1.59 2005/03/14 15:20:28 fwarmerdam
- * Fixed last change.
- *
- * Revision 1.58 2005/03/14 15:18:54 fwarmerdam
- * Treat very wide fields with no decimals as double. This is
- * more than 32bit integer fields.
- *
- * Revision 1.57 2005/02/10 20:16:54 fwarmerdam
- * Make the pszStringField buffer for ossim_DBFReadAttribute() static char [256]
- * as per bug 306.
- *
- * Revision 1.56 2005/02/10 20:07:56 fwarmerdam
- * Fixed bug 305 in ossim_DBFCloneEmpty() - header length problem.
- *
- * Revision 1.55 2004/09/26 20:23:46 fwarmerdam
- * avoid warnings with rcsid and signed/unsigned stuff
- *
- * Revision 1.54 2004/09/15 16:26:10 fwarmerdam
- * Treat all blank numeric fields as null too.
- *
- * Revision 1.53 2003/12/29 00:00:30 fwarmerdam
- * mark ossim_DBFWriteAttributeDirectly as ossim_SHPAPI_CALL
- *
- * Revision 1.52 2003/07/08 15:20:03 warmerda
- * avoid warnings about downcasting to unsigned char
- *
- * Revision 1.51 2003/07/08 13:50:15 warmerda
- * ossim_DBFIsAttributeNULL check for pszValue==NULL - bug 360
- *
- * Revision 1.50 2003/04/21 18:58:25 warmerda
- * ensure current record is flushed at same time as header is updated
- *
- * Revision 1.49 2003/04/21 18:30:37 warmerda
- * added header write/update public methods
- *
- * Revision 1.48 2003/03/10 14:51:27 warmerda
- * ossim_DBFWrite* calls now return FALSE if they have to truncate
- *
- * Revision 1.47 2002/11/20 03:32:22 warmerda
- * Ensure field name in ossim_DBFGetFieldIndex() is properly terminated.
- *
- * Revision 1.46 2002/10/09 13:10:21 warmerda
- * Added check that width is positive.
- *
- * Revision 1.45 2002/09/29 00:00:08 warmerda
- * added FTLogical and logical attribute read/write calls
- *
- * Revision 1.44 2002/05/07 13:46:11 warmerda
- * Added ossim_DBFWriteAttributeDirectly().
- *
- * Revision 1.43 2002/02/13 19:39:21 warmerda
- * Fix casting issues in ossim_DBFCloneEmpty().
- *
- * Revision 1.42 2002/01/15 14:36:07 warmerda
- * updated email address
- *
- * Revision 1.41 2002/01/15 14:31:49 warmerda
- * compute rather than copying nHeaderLength in ossim_DBFCloneEmpty()
- *
- * Revision 1.40 2002/01/09 04:32:35 warmerda
- * fixed to read correct amount of header
- *
- * Revision 1.39 2001/12/11 22:41:03 warmerda
- * improve io related error checking when reading header
- *
- * Revision 1.38 2001/11/28 16:07:31 warmerda
- * Cleanup to avoid compiler warnings as suggested by Richard Hash.
- *
- * Revision 1.37 2001/07/04 05:18:09 warmerda
- * do last fix properly
- *
- * Revision 1.36 2001/07/04 05:16:09 warmerda
- * fixed fieldname comparison in ossim_DBFGetFieldIndex
- *
- * Revision 1.35 2001/06/22 02:10:06 warmerda
- * fixed NULL shape support with help from Jim Matthews
- *
- * Revision 1.33 2001/05/31 19:20:13 warmerda
- * added ossim_DBFGetFieldIndex()
- *
- * Revision 1.32 2001/05/31 18:15:40 warmerda
- * Added support for NULL fields in ossim_DBF files
- *
- * Revision 1.31 2001/05/23 13:36:52 warmerda
- * added use of ossim_SHPAPI_CALL
- *
- * Revision 1.30 2000/12/05 14:43:38 warmerda
- * DBReadAttribute() white space trimming bug fix
- *
- * Revision 1.29 2000/10/05 14:36:44 warmerda
- * fix bug with writing very wide numeric fields
- *
- * Revision 1.28 2000/09/25 14:18:07 warmerda
- * Added some casts of strlen() return result to fix warnings on some
- * systems, as submitted by Daniel.
- *
- * Revision 1.27 2000/09/25 14:15:51 warmerda
- * added ossim_DBFGetNativeFieldType()
- *
- * Revision 1.26 2000/07/07 13:39:45 warmerda
- * removed unused variables, and added system include files
- *
- * Revision 1.25 2000/05/29 18:19:13 warmerda
- * avoid use of uchar, and adding casting fix
- *
- * Revision 1.24 2000/05/23 13:38:27 warmerda
- * Added error checks on return results of fread() and fseek().
- *
- * Revision 1.23 2000/05/23 13:25:49 warmerda
- * Avoid crashing if field or record are out of range in dbfread*attribute().
- *
- * Revision 1.22 1999/12/15 13:47:24 warmerda
- * Added stdlib.h to ensure that atof() is prototyped.
- *
- * Revision 1.21 1999/12/13 17:25:46 warmerda
- * Added support for upper case .DBF extention.
- *
- * Revision 1.20 1999/11/30 16:32:11 warmerda
- * Use atof() instead of sscanf().
- *
- * Revision 1.19 1999/11/05 14:12:04 warmerda
- * updated license terms
- *
- * Revision 1.18 1999/07/27 00:53:28 warmerda
- * ensure that whole old field value clear on write of string
- *
- * Revision 1.1 1999/07/05 18:58:07 warmerda
- * New
- *
- * Revision 1.17 1999/06/11 19:14:12 warmerda
- * Fixed some memory leaks.
- *
- * Revision 1.16 1999/06/11 19:04:11 warmerda
- * Remoted some unused variables.
- *
- * Revision 1.15 1999/05/11 03:19:28 warmerda
- * added new Tuple api, and improved extension handling - add from candrsn
- *
- * Revision 1.14 1999/05/04 15:01:48 warmerda
- * Added 'F' support.
- *
- * Revision 1.13 1999/03/23 17:38:59 warmerda
- * ossim_DBFAddField() now actually does return the new field number, or -1 if
- * it fails.
- *
- * Revision 1.12 1999/03/06 02:54:46 warmerda
- * Added logic to convert shapefile name to dbf filename in ossim_DBFOpen()
- * for convenience.
- *
- * Revision 1.11 1998/12/31 15:30:34 warmerda
- * Improved the interchangability of numeric and string attributes. Add
- * white space trimming option for attributes.
- *
- * Revision 1.10 1998/12/03 16:36:44 warmerda
- * Use r+b instead of rb+ for binary access.
- *
- * Revision 1.9 1998/12/03 15:34:23 warmerda
- * Updated copyright message.
- *
- * Revision 1.8 1997/12/04 15:40:15 warmerda
- * Added newline character after field definitions.
- *
- * Revision 1.7 1997/03/06 14:02:10 warmerda
- * Ensure bUpdated is initialized.
- *
- * Revision 1.6 1996/02/12 04:54:41 warmerda
- * Ensure that ossim_DBFWriteAttribute() returns TRUE if it succeeds.
- *
- * Revision 1.5 1995/10/21 03:15:12 warmerda
- * Changed to use binary file access, and ensure that the
- * field name field is zero filled, and limited to 10 chars.
- *
- * Revision 1.4 1995/08/24 18:10:42 warmerda
- * Added use of SfRealloc() to avoid pre-ANSI realloc() functions such
- * as on the Sun.
- *
- * Revision 1.3 1995/08/04 03:15:16 warmerda
- * Fixed up header.
- *
- * Revision 1.2 1995/08/04 03:14:43 warmerda
- * Added header.
- */
-
-#include <ossim/vec/shapefil.h>
-
-#include <math.h>
-#include <stdlib.h>
-#include <ctype.h>
-#include <string.h>
-
-ossim_SHP_CVSID("$Id: dbfopen.c 9094 2006-06-13 19:12:40Z dburken $")
-
-#ifndef FALSE
-# define FALSE 0
-# define TRUE 1
-#endif
-
-/************************************************************************/
-/* SfRealloc() */
-/* */
-/* A realloc cover function that will access a NULL pointer as */
-/* a valid input. */
-/************************************************************************/
-
-static void * SfRealloc( void * pMem, int nNewSize )
-
-{
- if( pMem == NULL )
- return( (void *) malloc(nNewSize) );
- else
- return( (void *) realloc(pMem,nNewSize) );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteHeader() */
-/* */
-/* This is called to write out the file header, and field */
-/* descriptions before writing any actual data records. This */
-/* also computes all the ossim_DBFDataSet field offset/size/decimals */
-/* and so forth values. */
-/************************************************************************/
-
-static void ossim_DBFWriteHeader(ossim_DBFHandle psDBF)
-
-{
- unsigned char abyHeader[XBASE_FLDHDR_SZ];
- int i;
-
- if( !psDBF->bNoHeader )
- return;
-
- psDBF->bNoHeader = FALSE;
-
-/* -------------------------------------------------------------------- */
-/* Initialize the file header information. */
-/* -------------------------------------------------------------------- */
- for( i = 0; i < XBASE_FLDHDR_SZ; i++ )
- abyHeader[i] = 0;
-
- abyHeader[0] = 0x03; /* memo field? - just copying */
-
- /* write out a dummy date */
- abyHeader[1] = 95; /* YY */
- abyHeader[2] = 7; /* MM */
- abyHeader[3] = 26; /* DD */
-
- /* record count preset at zero */
-
- abyHeader[8] = (unsigned char) (psDBF->nHeaderLength % 256);
- abyHeader[9] = (unsigned char) (psDBF->nHeaderLength / 256);
-
- abyHeader[10] = (unsigned char) (psDBF->nRecordLength % 256);
- abyHeader[11] = (unsigned char) (psDBF->nRecordLength / 256);
-
-/* -------------------------------------------------------------------- */
-/* Write the initial 32 byte file header, and all the field */
-/* descriptions. */
-/* -------------------------------------------------------------------- */
- fseek( psDBF->fp, 0, 0 );
- fwrite( abyHeader, XBASE_FLDHDR_SZ, 1, psDBF->fp );
- fwrite( psDBF->pszHeader, XBASE_FLDHDR_SZ, psDBF->nFields, psDBF->fp );
-
-/* -------------------------------------------------------------------- */
-/* Write out the newline character if there is room for it. */
-/* -------------------------------------------------------------------- */
- if( psDBF->nHeaderLength > 32*psDBF->nFields + 32 )
- {
- char cNewline;
-
- cNewline = 0x0d;
- fwrite( &cNewline, 1, 1, psDBF->fp );
- }
-}
-
-/************************************************************************/
-/* ossim_DBFFlushRecord() */
-/* */
-/* Write out the current record if there is one. */
-/************************************************************************/
-
-static void ossim_DBFFlushRecord( ossim_DBFHandle psDBF )
-
-{
- int nRecordOffset;
-
- if( psDBF->bCurrentRecordModified && psDBF->nCurrentRecord > -1 )
- {
- psDBF->bCurrentRecordModified = FALSE;
-
- nRecordOffset = psDBF->nRecordLength * psDBF->nCurrentRecord
- + psDBF->nHeaderLength;
-
- fseek( psDBF->fp, nRecordOffset, 0 );
- fwrite( psDBF->pszCurrentRecord, psDBF->nRecordLength, 1, psDBF->fp );
- }
-}
-
-/************************************************************************/
-/* ossim_DBFUpdateHeader() */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_DBFUpdateHeader( ossim_DBFHandle psDBF )
-
-{
- unsigned char abyFileHeader[32];
-
- if( psDBF->bNoHeader )
- ossim_DBFWriteHeader( psDBF );
-
- ossim_DBFFlushRecord( psDBF );
-
- fseek( psDBF->fp, 0, 0 );
- fread( abyFileHeader, 32, 1, psDBF->fp );
-
- abyFileHeader[4] = (unsigned char) (psDBF->nRecords % 256);
- abyFileHeader[5] = (unsigned char) ((psDBF->nRecords/256) % 256);
- abyFileHeader[6] = (unsigned char) ((psDBF->nRecords/(256*256)) % 256);
- abyFileHeader[7] = (unsigned char) ((psDBF->nRecords/(256*256*256)) % 256);
-
- fseek( psDBF->fp, 0, 0 );
- fwrite( abyFileHeader, 32, 1, psDBF->fp );
-
- fflush( psDBF->fp );
-}
-
-/************************************************************************/
-/* ossim_DBFOpen() */
-/* */
-/* Open a .dbf file. */
-/************************************************************************/
-
-ossim_DBFHandle ossim_SHPAPI_CALL
- ossim_DBFOpen( const char * pszFilename, const char * pszAccess )
-
-{
- ossim_DBFHandle psDBF;
- unsigned char *pabyBuf;
- int nFields, nHeadLen, nRecLen, iField, i;
- char *pszBasename, *pszFullname;
-
-/* -------------------------------------------------------------------- */
-/* We only allow the access strings "rb" and "r+". */
-/* -------------------------------------------------------------------- */
- if( strcmp(pszAccess,"r") != 0 && strcmp(pszAccess,"r+") != 0
- && strcmp(pszAccess,"rb") != 0 && strcmp(pszAccess,"rb+") != 0
- && strcmp(pszAccess,"r+b") != 0 )
- return( NULL );
-
- if( strcmp(pszAccess,"r") == 0 )
- pszAccess = "rb";
-
- if( strcmp(pszAccess,"r+") == 0 )
- pszAccess = "rb+";
-
-/* -------------------------------------------------------------------- */
-/* Compute the base (layer) name. If there is any extension */
-/* on the passed in filename we will strip it off. */
-/* -------------------------------------------------------------------- */
- pszBasename = (char *) malloc(strlen(pszFilename)+5);
- strcpy( pszBasename, pszFilename );
- for( i = strlen(pszBasename)-1;
- i > 0 && pszBasename[i] != '.' && pszBasename[i] != '/'
- && pszBasename[i] != '\\';
- i-- ) {}
-
- if( pszBasename[i] == '.' )
- pszBasename[i] = '\0';
-
- pszFullname = (char *) malloc(strlen(pszBasename) + 5);
- sprintf( pszFullname, "%s.dbf", pszBasename );
-
- psDBF = (ossim_DBFHandle) calloc( 1, sizeof(ossim_DBFInfo) );
- psDBF->fp = fopen( pszFullname, pszAccess );
-
- if( psDBF->fp == NULL )
- {
- sprintf( pszFullname, "%s.DBF", pszBasename );
- psDBF->fp = fopen(pszFullname, pszAccess );
- }
-
- free( pszBasename );
- free( pszFullname );
-
- if( psDBF->fp == NULL )
- {
- free( psDBF );
- return( NULL );
- }
-
- psDBF->bNoHeader = FALSE;
- psDBF->nCurrentRecord = -1;
- psDBF->bCurrentRecordModified = FALSE;
-
-/* -------------------------------------------------------------------- */
-/* Read Table Header info */
-/* -------------------------------------------------------------------- */
- pabyBuf = (unsigned char *) malloc(500);
- if( fread( pabyBuf, 32, 1, psDBF->fp ) != 1 )
- {
- fclose( psDBF->fp );
- free( pabyBuf );
- free( psDBF );
- return NULL;
- }
-
- psDBF->nRecords =
- pabyBuf[4] + pabyBuf[5]*256 + pabyBuf[6]*256*256 + pabyBuf[7]*256*256*256;
-
- psDBF->nHeaderLength = nHeadLen = pabyBuf[8] + pabyBuf[9]*256;
- psDBF->nRecordLength = nRecLen = pabyBuf[10] + pabyBuf[11]*256;
-
- psDBF->nFields = nFields = (nHeadLen - 32) / 32;
-
- psDBF->pszCurrentRecord = (char *) malloc(nRecLen);
-
-/* -------------------------------------------------------------------- */
-/* Read in Field Definitions */
-/* -------------------------------------------------------------------- */
-
- pabyBuf = (unsigned char *) SfRealloc(pabyBuf,nHeadLen);
- psDBF->pszHeader = (char *) pabyBuf;
-
- fseek( psDBF->fp, 32, 0 );
- if( fread( pabyBuf, nHeadLen-32, 1, psDBF->fp ) != 1 )
- {
- fclose( psDBF->fp );
- free( pabyBuf );
- free( psDBF );
- return NULL;
- }
-
- psDBF->panFieldOffset = (int *) malloc(sizeof(int) * nFields);
- psDBF->panFieldSize = (int *) malloc(sizeof(int) * nFields);
- psDBF->panFieldDecimals = (int *) malloc(sizeof(int) * nFields);
- psDBF->pachFieldType = (char *) malloc(sizeof(char) * nFields);
-
- for( iField = 0; iField < nFields; iField++ )
- {
- unsigned char *pabyFInfo;
-
- pabyFInfo = pabyBuf+iField*32;
-
- if( pabyFInfo[11] == 'N' || pabyFInfo[11] == 'F' )
- {
- psDBF->panFieldSize[iField] = pabyFInfo[16];
- psDBF->panFieldDecimals[iField] = pabyFInfo[17];
- }
- else
- {
- psDBF->panFieldSize[iField] = pabyFInfo[16] + pabyFInfo[17]*256;
- psDBF->panFieldDecimals[iField] = 0;
- }
-
- psDBF->pachFieldType[iField] = (char) pabyFInfo[11];
- if( iField == 0 )
- psDBF->panFieldOffset[iField] = 1;
- else
- psDBF->panFieldOffset[iField] =
- psDBF->panFieldOffset[iField-1] + psDBF->panFieldSize[iField-1];
- }
-
- return( psDBF );
-}
-
-/************************************************************************/
-/* ossim_DBFClose() */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_DBFClose(ossim_DBFHandle psDBF)
-{
-/* -------------------------------------------------------------------- */
-/* Write out header if not already written. */
-/* -------------------------------------------------------------------- */
- if( psDBF->bNoHeader )
- ossim_DBFWriteHeader( psDBF );
-
- ossim_DBFFlushRecord( psDBF );
-
-/* -------------------------------------------------------------------- */
-/* Update last access date, and number of records if we have */
-/* write access. */
-/* -------------------------------------------------------------------- */
- if( psDBF->bUpdated )
- ossim_DBFUpdateHeader( psDBF );
-
-/* -------------------------------------------------------------------- */
-/* Close, and free resources. */
-/* -------------------------------------------------------------------- */
- fclose( psDBF->fp );
-
- if( psDBF->panFieldOffset != NULL )
- {
- free( psDBF->panFieldOffset );
- free( psDBF->panFieldSize );
- free( psDBF->panFieldDecimals );
- free( psDBF->pachFieldType );
- }
-
- free( psDBF->pszHeader );
- free( psDBF->pszCurrentRecord );
-
- free( psDBF );
-}
-
-/************************************************************************/
-/* ossim_DBFCreate() */
-/* */
-/* Create a new .dbf file. */
-/************************************************************************/
-
-ossim_DBFHandle ossim_SHPAPI_CALL
- ossim_DBFCreate( const char * pszFilename )
-
-{
- ossim_DBFHandle psDBF;
- FILE *fp;
- char *pszFullname, *pszBasename;
- int i;
-
-/* -------------------------------------------------------------------- */
-/* Compute the base (layer) name. If there is any extension */
-/* on the passed in filename we will strip it off. */
-/* -------------------------------------------------------------------- */
- pszBasename = (char *) malloc(strlen(pszFilename)+5);
- strcpy( pszBasename, pszFilename );
- for( i = strlen(pszBasename)-1;
- i > 0 && pszBasename[i] != '.' && pszBasename[i] != '/'
- && pszBasename[i] != '\\';
- i-- ) {}
-
- if( pszBasename[i] == '.' )
- pszBasename[i] = '\0';
-
- pszFullname = (char *) malloc(strlen(pszBasename) + 5);
- sprintf( pszFullname, "%s.dbf", pszBasename );
- free( pszBasename );
-
-/* -------------------------------------------------------------------- */
-/* Create the file. */
-/* -------------------------------------------------------------------- */
- fp = fopen( pszFullname, "wb" );
- if( fp == NULL )
- return( NULL );
-
- fputc( 0, fp );
- fclose( fp );
-
- fp = fopen( pszFullname, "rb+" );
- if( fp == NULL )
- return( NULL );
-
- free( pszFullname );
-
-/* -------------------------------------------------------------------- */
-/* Create the info structure. */
-/* -------------------------------------------------------------------- */
- psDBF = (ossim_DBFHandle) malloc(sizeof(ossim_DBFInfo));
-
- psDBF->fp = fp;
- psDBF->nRecords = 0;
- psDBF->nFields = 0;
- psDBF->nRecordLength = 1;
- psDBF->nHeaderLength = 33;
-
- psDBF->panFieldOffset = NULL;
- psDBF->panFieldSize = NULL;
- psDBF->panFieldDecimals = NULL;
- psDBF->pachFieldType = NULL;
- psDBF->pszHeader = NULL;
-
- psDBF->nCurrentRecord = -1;
- psDBF->bCurrentRecordModified = FALSE;
- psDBF->pszCurrentRecord = NULL;
-
- psDBF->bNoHeader = TRUE;
-
- return( psDBF );
-}
-
-/************************************************************************/
-/* ossim_DBFAddField() */
-/* */
-/* Add a field to a newly created .dbf file before any records */
-/* are written. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFAddField(ossim_DBFHandle psDBF, const char * pszFieldName,
- ossim_DBFFieldType eType, int nWidth, int nDecimals )
-
-{
- char *pszFInfo;
- int i;
-
-/* -------------------------------------------------------------------- */
-/* Do some checking to ensure we can add records to this file. */
-/* -------------------------------------------------------------------- */
- if( psDBF->nRecords > 0 )
- return( -1 );
-
- if( !psDBF->bNoHeader )
- return( -1 );
-
- if( eType != FTDouble && nDecimals != 0 )
- return( -1 );
-
- if( nWidth < 1 )
- return -1;
-
-/* -------------------------------------------------------------------- */
-/* SfRealloc all the arrays larger to hold the additional field */
-/* information. */
-/* -------------------------------------------------------------------- */
- psDBF->nFields++;
-
- psDBF->panFieldOffset = (int *)
- SfRealloc( psDBF->panFieldOffset, sizeof(int) * psDBF->nFields );
-
- psDBF->panFieldSize = (int *)
- SfRealloc( psDBF->panFieldSize, sizeof(int) * psDBF->nFields );
-
- psDBF->panFieldDecimals = (int *)
- SfRealloc( psDBF->panFieldDecimals, sizeof(int) * psDBF->nFields );
-
- psDBF->pachFieldType = (char *)
- SfRealloc( psDBF->pachFieldType, sizeof(char) * psDBF->nFields );
-
-/* -------------------------------------------------------------------- */
-/* Assign the new field information fields. */
-/* -------------------------------------------------------------------- */
- psDBF->panFieldOffset[psDBF->nFields-1] = psDBF->nRecordLength;
- psDBF->nRecordLength += nWidth;
- psDBF->panFieldSize[psDBF->nFields-1] = nWidth;
- psDBF->panFieldDecimals[psDBF->nFields-1] = nDecimals;
-
- if( eType == FTLogical )
- psDBF->pachFieldType[psDBF->nFields-1] = 'L';
- else if( eType == FTString )
- psDBF->pachFieldType[psDBF->nFields-1] = 'C';
- else
- psDBF->pachFieldType[psDBF->nFields-1] = 'N';
-
-/* -------------------------------------------------------------------- */
-/* Extend the required header information. */
-/* -------------------------------------------------------------------- */
- psDBF->nHeaderLength += 32;
- psDBF->bUpdated = FALSE;
-
- psDBF->pszHeader = (char *) SfRealloc(psDBF->pszHeader,psDBF->nFields*32);
-
- pszFInfo = psDBF->pszHeader + 32 * (psDBF->nFields-1);
-
- for( i = 0; i < 32; i++ )
- pszFInfo[i] = '\0';
-
- if( (int) strlen(pszFieldName) < 10 )
- strncpy( pszFInfo, pszFieldName, strlen(pszFieldName));
- else
- strncpy( pszFInfo, pszFieldName, 10);
-
- pszFInfo[11] = psDBF->pachFieldType[psDBF->nFields-1];
-
- if( eType == FTString )
- {
- pszFInfo[16] = (unsigned char) (nWidth % 256);
- pszFInfo[17] = (unsigned char) (nWidth / 256);
- }
- else
- {
- pszFInfo[16] = (unsigned char) nWidth;
- pszFInfo[17] = (unsigned char) nDecimals;
- }
-
-/* -------------------------------------------------------------------- */
-/* Make the current record buffer appropriately larger. */
-/* -------------------------------------------------------------------- */
- psDBF->pszCurrentRecord = (char *) SfRealloc(psDBF->pszCurrentRecord,
- psDBF->nRecordLength);
-
- return( psDBF->nFields-1 );
-}
-
-/************************************************************************/
-/* ossim_DBFReadAttribute() */
-/* */
-/* Read one of the attribute fields of a record. */
-/************************************************************************/
-
-static void *ossim_DBFReadAttribute(ossim_DBFHandle psDBF, int hEntity, int iField,
- char chReqType )
-
-{
- int nRecordOffset;
- unsigned char *pabyRec;
- void *pReturnField = NULL;
-
- static double dDoubleField;
- static char szStringField[257];
-
-/* -------------------------------------------------------------------- */
-/* Verify selection. */
-/* -------------------------------------------------------------------- */
- if( hEntity < 0 || hEntity >= psDBF->nRecords )
- return( NULL );
-
- if( iField < 0 || iField >= psDBF->nFields )
- return( NULL );
-
-/* -------------------------------------------------------------------- */
-/* Have we read the record? */
-/* -------------------------------------------------------------------- */
- if( psDBF->nCurrentRecord != hEntity )
- {
- ossim_DBFFlushRecord( psDBF );
-
- nRecordOffset = psDBF->nRecordLength * hEntity + psDBF->nHeaderLength;
-
- if( fseek( psDBF->fp, nRecordOffset, 0 ) != 0 )
- {
- fprintf( stderr, "fseek(%d) failed on ossim_DBF file.\n",
- nRecordOffset );
- return NULL;
- }
-
- if( fread( psDBF->pszCurrentRecord, psDBF->nRecordLength,
- 1, psDBF->fp ) != 1 )
- {
- fprintf( stderr, "fread(%d) failed on ossim_DBF file.\n",
- psDBF->nRecordLength );
- return NULL;
- }
-
- psDBF->nCurrentRecord = hEntity;
- }
-
- pabyRec = (unsigned char *) psDBF->pszCurrentRecord;
-
-/* -------------------------------------------------------------------- */
-/* Extract the requested field. */
-/* -------------------------------------------------------------------- */
- strncpy( szStringField,
- ((const char *) pabyRec) + psDBF->panFieldOffset[iField],
- psDBF->panFieldSize[iField] );
- szStringField[psDBF->panFieldSize[iField]] = '\0';
-
- pReturnField = szStringField;
-
-/* -------------------------------------------------------------------- */
-/* Decode the field. */
-/* -------------------------------------------------------------------- */
- if( chReqType == 'N' )
- {
- dDoubleField = atof(szStringField);
-
- pReturnField = &dDoubleField;
- }
-
-/* -------------------------------------------------------------------- */
-/* Should we trim white space off the string attribute value? */
-/* -------------------------------------------------------------------- */
-#ifdef TRIM_DBF_WHITESPACE
- else
- {
- char *pchSrc, *pchDst;
-
- pchDst = pchSrc = szStringField;
- while( *pchSrc == ' ' )
- pchSrc++;
-
- while( *pchSrc != '\0' )
- *(pchDst++) = *(pchSrc++);
- *pchDst = '\0';
-
- while( pchDst != szStringField && *(--pchDst) == ' ' )
- *pchDst = '\0';
- }
-#endif
-
- return( pReturnField );
-}
-
-/************************************************************************/
-/* ossim_DBFReadIntAttribute() */
-/* */
-/* Read an integer attribute. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFReadIntegerAttribute( ossim_DBFHandle psDBF, int iRecord, int iField )
-
-{
- double *pdValue;
-
- pdValue = (double *) ossim_DBFReadAttribute( psDBF, iRecord, iField, 'N' );
-
- if( pdValue == NULL )
- return 0;
- else
- return( (int) *pdValue );
-}
-
-/************************************************************************/
-/* ossim_DBFReadDoubleAttribute() */
-/* */
-/* Read a double attribute. */
-/************************************************************************/
-
-double ossim_SHPAPI_CALL
- ossim_DBFReadDoubleAttribute( ossim_DBFHandle psDBF, int iRecord, int iField )
-
-{
- double *pdValue;
-
- pdValue = (double *) ossim_DBFReadAttribute( psDBF, iRecord, iField, 'N' );
-
- if( pdValue == NULL )
- return 0.0;
- else
- return( *pdValue );
-}
-
-/************************************************************************/
-/* ossim_DBFReadStringAttribute() */
-/* */
-/* Read a string attribute. */
-/************************************************************************/
-
-const char ossim_SHPAPI_CALL1(*)
- ossim_DBFReadStringAttribute( ossim_DBFHandle psDBF, int iRecord, int iField )
-
-{
- return( (const char *) ossim_DBFReadAttribute( psDBF, iRecord, iField, 'C' ) );
-}
-
-/************************************************************************/
-/* ossim_DBFReadLogicalAttribute() */
-/* */
-/* Read a logical attribute. */
-/************************************************************************/
-
-const char ossim_SHPAPI_CALL1(*)
- ossim_DBFReadLogicalAttribute( ossim_DBFHandle psDBF, int iRecord, int iField )
-
-{
- return( (const char *) ossim_DBFReadAttribute( psDBF, iRecord, iField, 'L' ) );
-}
-
-/************************************************************************/
-/* ossim_DBFIsAttributeNULL() */
-/* */
-/* Return TRUE if value for field is NULL. */
-/* */
-/* Contributed by Jim Matthews. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFIsAttributeNULL( ossim_DBFHandle psDBF, int iRecord, int iField )
-
-{
- const char *pszValue;
- int i;
-
- pszValue = ossim_DBFReadStringAttribute( psDBF, iRecord, iField );
-
- if( pszValue == NULL )
- return TRUE;
-
- switch(psDBF->pachFieldType[iField])
- {
- case 'N':
- case 'F':
- /*
- ** We accept all asterisks or all blanks as NULL
- ** though according to the spec I think it should be all
- ** asterisks.
- */
- if( pszValue[0] == '*' )
- return TRUE;
-
- for( i = 0; pszValue[i] != '\0'; i++ )
- {
- if( pszValue[i] != ' ' )
- return FALSE;
- }
- return TRUE;
-
- case 'D':
- /* NULL date fields have value "00000000" */
- return strncmp(pszValue,"00000000",8) == 0;
-
- case 'L':
- /* NULL boolean fields have value "?" */
- return pszValue[0] == '?';
-
- default:
- /* empty string fields are considered NULL */
- return strlen(pszValue) == 0;
- }
-}
-
-/************************************************************************/
-/* ossim_DBFGetFieldCount() */
-/* */
-/* Return the number of fields in this table. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFGetFieldCount( ossim_DBFHandle psDBF )
-
-{
- return( psDBF->nFields );
-}
-
-/************************************************************************/
-/* ossim_DBFGetRecordCount() */
-/* */
-/* Return the number of records in this table. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFGetRecordCount( ossim_DBFHandle psDBF )
-
-{
- return( psDBF->nRecords );
-}
-
-/************************************************************************/
-/* ossim_DBFGetFieldInfo() */
-/* */
-/* Return any requested information about the field. */
-/************************************************************************/
-
-ossim_DBFFieldType ossim_SHPAPI_CALL
- ossim_DBFGetFieldInfo( ossim_DBFHandle psDBF, int iField, char * pszFieldName,
- int * pnWidth, int * pnDecimals )
-
-{
- if( iField < 0 || iField >= psDBF->nFields )
- return( FTInvalid );
-
- if( pnWidth != NULL )
- *pnWidth = psDBF->panFieldSize[iField];
-
- if( pnDecimals != NULL )
- *pnDecimals = psDBF->panFieldDecimals[iField];
-
- if( pszFieldName != NULL )
- {
- int i;
-
- strncpy( pszFieldName, (char *) psDBF->pszHeader+iField*32, 11 );
- pszFieldName[11] = '\0';
- for( i = 10; i > 0 && pszFieldName[i] == ' '; i-- )
- pszFieldName[i] = '\0';
- }
-
- if ( psDBF->pachFieldType[iField] == 'L' )
- return( FTLogical);
-
- else if( psDBF->pachFieldType[iField] == 'N'
- || psDBF->pachFieldType[iField] == 'F'
- || psDBF->pachFieldType[iField] == 'D' )
- {
- if( psDBF->panFieldDecimals[iField] > 0
- || psDBF->panFieldSize[iField] > 10 )
- return( FTDouble );
- else
- return( FTInteger );
- }
- else
- {
- return( FTString );
- }
-}
-
-/************************************************************************/
-/* ossim_DBFWriteAttribute() */
-/* */
-/* Write an attribute record to the file. */
-/************************************************************************/
-
-static int ossim_DBFWriteAttribute(ossim_DBFHandle psDBF, int hEntity, int iField,
- void * pValue )
-
-{
- int nRecordOffset, i, j, nRetResult = TRUE;
- unsigned char *pabyRec;
- char szSField[400], szFormat[20];
-
-/* -------------------------------------------------------------------- */
-/* Is this a valid record? */
-/* -------------------------------------------------------------------- */
- if( hEntity < 0 || hEntity > psDBF->nRecords )
- return( FALSE );
-
- if( psDBF->bNoHeader )
- ossim_DBFWriteHeader(psDBF);
-
-/* -------------------------------------------------------------------- */
-/* Is this a brand new record? */
-/* -------------------------------------------------------------------- */
- if( hEntity == psDBF->nRecords )
- {
- ossim_DBFFlushRecord( psDBF );
-
- psDBF->nRecords++;
- for( i = 0; i < psDBF->nRecordLength; i++ )
- psDBF->pszCurrentRecord[i] = ' ';
-
- psDBF->nCurrentRecord = hEntity;
- }
-
-/* -------------------------------------------------------------------- */
-/* Is this an existing record, but different than the last one */
-/* we accessed? */
-/* -------------------------------------------------------------------- */
- if( psDBF->nCurrentRecord != hEntity )
- {
- ossim_DBFFlushRecord( psDBF );
-
- nRecordOffset = psDBF->nRecordLength * hEntity + psDBF->nHeaderLength;
-
- fseek( psDBF->fp, nRecordOffset, 0 );
- fread( psDBF->pszCurrentRecord, psDBF->nRecordLength, 1, psDBF->fp );
-
- psDBF->nCurrentRecord = hEntity;
- }
-
- pabyRec = (unsigned char *) psDBF->pszCurrentRecord;
-
- psDBF->bCurrentRecordModified = TRUE;
- psDBF->bUpdated = TRUE;
-
-/* -------------------------------------------------------------------- */
-/* Translate NULL value to valid ossim_DBF file representation. */
-/* */
-/* Contributed by Jim Matthews. */
-/* -------------------------------------------------------------------- */
- if( pValue == NULL )
- {
- switch(psDBF->pachFieldType[iField])
- {
- case 'N':
- case 'F':
- /* NULL numeric fields have value "****************" */
- memset( (char *) (pabyRec+psDBF->panFieldOffset[iField]), '*',
- psDBF->panFieldSize[iField] );
- break;
-
- case 'D':
- /* NULL date fields have value "00000000" */
- memset( (char *) (pabyRec+psDBF->panFieldOffset[iField]), '0',
- psDBF->panFieldSize[iField] );
- break;
-
- case 'L':
- /* NULL boolean fields have value "?" */
- memset( (char *) (pabyRec+psDBF->panFieldOffset[iField]), '?',
- psDBF->panFieldSize[iField] );
- break;
-
- default:
- /* empty string fields are considered NULL */
- memset( (char *) (pabyRec+psDBF->panFieldOffset[iField]), '\0',
- psDBF->panFieldSize[iField] );
- break;
- }
- return TRUE;
- }
-
-/* -------------------------------------------------------------------- */
-/* Assign all the record fields. */
-/* -------------------------------------------------------------------- */
- switch( psDBF->pachFieldType[iField] )
- {
- case 'D':
- case 'N':
- case 'F':
- if( psDBF->panFieldDecimals[iField] == 0 )
- {
- int nWidth = psDBF->panFieldSize[iField];
-
- if( (int) sizeof(szSField)-2 < nWidth )
- nWidth = sizeof(szSField)-2;
-
- sprintf( szFormat, "%%%dd", nWidth );
- sprintf(szSField, szFormat, (int) *((double *) pValue) );
- if( (int)strlen(szSField) > psDBF->panFieldSize[iField] )
- {
- szSField[psDBF->panFieldSize[iField]] = '\0';
- nRetResult = FALSE;
- }
-
- strncpy((char *) (pabyRec+psDBF->panFieldOffset[iField]),
- szSField, strlen(szSField) );
- }
- else
- {
- int nWidth = psDBF->panFieldSize[iField];
-
- if( (int) sizeof(szSField)-2 < nWidth )
- nWidth = sizeof(szSField)-2;
-
- sprintf( szFormat, "%%%d.%df",
- nWidth, psDBF->panFieldDecimals[iField] );
- sprintf(szSField, szFormat, *((double *) pValue) );
- if( (int) strlen(szSField) > psDBF->panFieldSize[iField] )
- {
- szSField[psDBF->panFieldSize[iField]] = '\0';
- nRetResult = FALSE;
- }
- strncpy((char *) (pabyRec+psDBF->panFieldOffset[iField]),
- szSField, strlen(szSField) );
- }
- break;
-
- case 'L':
- if (psDBF->panFieldSize[iField] >= 1 &&
- (*(char*)pValue == 'F' || *(char*)pValue == 'T'))
- *(pabyRec+psDBF->panFieldOffset[iField]) = *(char*)pValue;
- break;
-
- default:
- if( (int) strlen((char *) pValue) > psDBF->panFieldSize[iField] )
- {
- j = psDBF->panFieldSize[iField];
- nRetResult = FALSE;
- }
- else
- {
- memset( pabyRec+psDBF->panFieldOffset[iField], ' ',
- psDBF->panFieldSize[iField] );
- j = strlen((char *) pValue);
- }
-
- strncpy((char *) (pabyRec+psDBF->panFieldOffset[iField]),
- (char *) pValue, j );
- break;
- }
-
- return( nRetResult );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteAttributeDirectly() */
-/* */
-/* Write an attribute record to the file, but without any */
-/* reformatting based on type. The provided buffer is written */
-/* as is to the field position in the record. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteAttributeDirectly(ossim_DBFHandle psDBF, int hEntity, int iField,
- void * pValue )
-
-{
- int nRecordOffset, i, j;
- unsigned char *pabyRec;
-
-/* -------------------------------------------------------------------- */
-/* Is this a valid record? */
-/* -------------------------------------------------------------------- */
- if( hEntity < 0 || hEntity > psDBF->nRecords )
- return( FALSE );
-
- if( psDBF->bNoHeader )
- ossim_DBFWriteHeader(psDBF);
-
-/* -------------------------------------------------------------------- */
-/* Is this a brand new record? */
-/* -------------------------------------------------------------------- */
- if( hEntity == psDBF->nRecords )
- {
- ossim_DBFFlushRecord( psDBF );
-
- psDBF->nRecords++;
- for( i = 0; i < psDBF->nRecordLength; i++ )
- psDBF->pszCurrentRecord[i] = ' ';
-
- psDBF->nCurrentRecord = hEntity;
- }
-
-/* -------------------------------------------------------------------- */
-/* Is this an existing record, but different than the last one */
-/* we accessed? */
-/* -------------------------------------------------------------------- */
- if( psDBF->nCurrentRecord != hEntity )
- {
- ossim_DBFFlushRecord( psDBF );
-
- nRecordOffset = psDBF->nRecordLength * hEntity + psDBF->nHeaderLength;
-
- fseek( psDBF->fp, nRecordOffset, 0 );
- fread( psDBF->pszCurrentRecord, psDBF->nRecordLength, 1, psDBF->fp );
-
- psDBF->nCurrentRecord = hEntity;
- }
-
- pabyRec = (unsigned char *) psDBF->pszCurrentRecord;
-
-/* -------------------------------------------------------------------- */
-/* Assign all the record fields. */
-/* -------------------------------------------------------------------- */
- if( (int)strlen((char *) pValue) > psDBF->panFieldSize[iField] )
- j = psDBF->panFieldSize[iField];
- else
- {
- memset( pabyRec+psDBF->panFieldOffset[iField], ' ',
- psDBF->panFieldSize[iField] );
- j = strlen((char *) pValue);
- }
-
- strncpy((char *) (pabyRec+psDBF->panFieldOffset[iField]),
- (char *) pValue, j );
-
- psDBF->bCurrentRecordModified = TRUE;
- psDBF->bUpdated = TRUE;
-
- return( TRUE );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteDoubleAttribute() */
-/* */
-/* Write a double attribute. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteDoubleAttribute( ossim_DBFHandle psDBF, int iRecord, int iField,
- double dValue )
-
-{
- return( ossim_DBFWriteAttribute( psDBF, iRecord, iField, (void *) &dValue ) );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteIntegerAttribute() */
-/* */
-/* Write a integer attribute. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteIntegerAttribute( ossim_DBFHandle psDBF, int iRecord, int iField,
- int nValue )
-
-{
- double dValue = nValue;
-
- return( ossim_DBFWriteAttribute( psDBF, iRecord, iField, (void *) &dValue ) );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteStringAttribute() */
-/* */
-/* Write a string attribute. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteStringAttribute( ossim_DBFHandle psDBF, int iRecord, int iField,
- const char * pszValue )
-
-{
- return( ossim_DBFWriteAttribute( psDBF, iRecord, iField, (void *) pszValue ) );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteNULLAttribute() */
-/* */
-/* Write a string attribute. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteNULLAttribute( ossim_DBFHandle psDBF, int iRecord, int iField )
-
-{
- return( ossim_DBFWriteAttribute( psDBF, iRecord, iField, NULL ) );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteLogicalAttribute() */
-/* */
-/* Write a logical attribute. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteLogicalAttribute( ossim_DBFHandle psDBF, int iRecord, int iField,
- const char lValue)
-
-{
- return( ossim_DBFWriteAttribute( psDBF, iRecord, iField, (void *) (&lValue) ) );
-}
-
-/************************************************************************/
-/* ossim_DBFWriteTuple() */
-/* */
-/* Write an attribute record to the file. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFWriteTuple(ossim_DBFHandle psDBF, int hEntity, void * pRawTuple )
-
-{
- int nRecordOffset, i;
- unsigned char *pabyRec;
-
-/* -------------------------------------------------------------------- */
-/* Is this a valid record? */
-/* -------------------------------------------------------------------- */
- if( hEntity < 0 || hEntity > psDBF->nRecords )
- return( FALSE );
-
- if( psDBF->bNoHeader )
- ossim_DBFWriteHeader(psDBF);
-
-/* -------------------------------------------------------------------- */
-/* Is this a brand new record? */
-/* -------------------------------------------------------------------- */
- if( hEntity == psDBF->nRecords )
- {
- ossim_DBFFlushRecord( psDBF );
-
- psDBF->nRecords++;
- for( i = 0; i < psDBF->nRecordLength; i++ )
- psDBF->pszCurrentRecord[i] = ' ';
-
- psDBF->nCurrentRecord = hEntity;
- }
-
-/* -------------------------------------------------------------------- */
-/* Is this an existing record, but different than the last one */
-/* we accessed? */
-/* -------------------------------------------------------------------- */
- if( psDBF->nCurrentRecord != hEntity )
- {
- ossim_DBFFlushRecord( psDBF );
-
- nRecordOffset = psDBF->nRecordLength * hEntity + psDBF->nHeaderLength;
-
- fseek( psDBF->fp, nRecordOffset, 0 );
- fread( psDBF->pszCurrentRecord, psDBF->nRecordLength, 1, psDBF->fp );
-
- psDBF->nCurrentRecord = hEntity;
- }
-
- pabyRec = (unsigned char *) psDBF->pszCurrentRecord;
-
- memcpy ( pabyRec, pRawTuple, psDBF->nRecordLength );
-
- psDBF->bCurrentRecordModified = TRUE;
- psDBF->bUpdated = TRUE;
-
- return( TRUE );
-}
-
-/************************************************************************/
-/* ossim_DBFReadTuple() */
-/* */
-/* Read one of the attribute fields of a record. */
-/************************************************************************/
-
-const char ossim_SHPAPI_CALL1(*)
- ossim_DBFReadTuple(ossim_DBFHandle psDBF, int hEntity )
-
-{
- int nRecordOffset;
- unsigned char *pabyRec;
- static char *pReturnTuple = NULL;
-
- static int nTupleLen = 0;
-
-/* -------------------------------------------------------------------- */
-/* Have we read the record? */
-/* -------------------------------------------------------------------- */
- if( hEntity < 0 || hEntity >= psDBF->nRecords )
- return( NULL );
-
- if( psDBF->nCurrentRecord != hEntity )
- {
- ossim_DBFFlushRecord( psDBF );
-
- nRecordOffset = psDBF->nRecordLength * hEntity + psDBF->nHeaderLength;
-
- fseek( psDBF->fp, nRecordOffset, 0 );
- fread( psDBF->pszCurrentRecord, psDBF->nRecordLength, 1, psDBF->fp );
-
- psDBF->nCurrentRecord = hEntity;
- }
-
- pabyRec = (unsigned char *) psDBF->pszCurrentRecord;
-
- if ( nTupleLen < psDBF->nRecordLength) {
- nTupleLen = psDBF->nRecordLength;
- pReturnTuple = (char *) SfRealloc(pReturnTuple, psDBF->nRecordLength);
- }
-
- memcpy ( pReturnTuple, pabyRec, psDBF->nRecordLength );
-
- return( pReturnTuple );
-}
-
-/************************************************************************/
-/* ossim_DBFCloneEmpty() */
-/* */
-/* Read one of the attribute fields of a record. */
-/************************************************************************/
-
-ossim_DBFHandle ossim_SHPAPI_CALL
- ossim_DBFCloneEmpty(ossim_DBFHandle psDBF, const char * pszFilename )
-{
- ossim_DBFHandle newDBF;
-
- newDBF = ossim_DBFCreate ( pszFilename );
- if ( newDBF == NULL ) return ( NULL );
-
- newDBF->nFields = psDBF->nFields;
- newDBF->nRecordLength = psDBF->nRecordLength;
- newDBF->nHeaderLength = psDBF->nHeaderLength;
-
- newDBF->pszHeader = (char *) malloc ( newDBF->nHeaderLength );
- memcpy ( newDBF->pszHeader, psDBF->pszHeader, newDBF->nHeaderLength );
-
- newDBF->panFieldOffset = (int *) malloc ( sizeof(int) * psDBF->nFields );
- memcpy ( newDBF->panFieldOffset, psDBF->panFieldOffset, sizeof(int) * psDBF->nFields );
- newDBF->panFieldSize = (int *) malloc ( sizeof(int) * psDBF->nFields );
- memcpy ( newDBF->panFieldSize, psDBF->panFieldSize, sizeof(int) * psDBF->nFields );
- newDBF->panFieldDecimals = (int *) malloc ( sizeof(int) * psDBF->nFields );
- memcpy ( newDBF->panFieldDecimals, psDBF->panFieldDecimals, sizeof(int) * psDBF->nFields );
- newDBF->pachFieldType = (char *) malloc ( sizeof(int) * psDBF->nFields );
- memcpy ( newDBF->pachFieldType, psDBF->pachFieldType, sizeof(int) * psDBF->nFields );
-
- newDBF->bNoHeader = TRUE;
- newDBF->bUpdated = TRUE;
-
- ossim_DBFWriteHeader ( newDBF );
- ossim_DBFClose ( newDBF );
-
- newDBF = ossim_DBFOpen ( pszFilename, "rb+" );
-
- return ( newDBF );
-}
-
-/************************************************************************/
-/* ossim_DBFGetNativeFieldType() */
-/* */
-/* Return the DBase field type for the specified field. */
-/* */
-/* Value can be one of: 'C' (String), 'D' (Date), 'F' (Float), */
-/* 'N' (Numeric, with or without decimal), */
-/* 'L' (Logical), */
-/* 'M' (Memo: 10 digits .DBT block ptr) */
-/************************************************************************/
-
-char ossim_SHPAPI_CALL
- ossim_DBFGetNativeFieldType( ossim_DBFHandle psDBF, int iField )
-
-{
- if( iField >=0 && iField < psDBF->nFields )
- return psDBF->pachFieldType[iField];
-
- return ' ';
-}
-
-/************************************************************************/
-/* str_to_upper() */
-/************************************************************************/
-
-static void str_to_upper (char *string)
-{
- int len;
- short i = -1;
-
- len = strlen (string);
-
- while (++i < len)
- if (isalpha(string[i]) && islower(string[i]))
- string[i] = (char) toupper ((int)string[i]);
-}
-
-/************************************************************************/
-/* ossim_DBFGetFieldIndex() */
-/* */
-/* Get the index number for a field in a .dbf file. */
-/* */
-/* Contributed by Jim Matthews. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_DBFGetFieldIndex(ossim_DBFHandle psDBF, const char *pszFieldName)
-
-{
- char name[12], name1[12], name2[12];
- int i;
-
- strncpy(name1, pszFieldName,11);
- name1[11] = '\0';
- str_to_upper(name1);
-
- for( i = 0; i < ossim_DBFGetFieldCount(psDBF); i++ )
- {
- ossim_DBFGetFieldInfo( psDBF, i, name, NULL, NULL );
- strncpy(name2,name,11);
- str_to_upper(name2);
-
- if(!strncmp(name1,name2,10))
- return(i);
- }
- return(-1);
-}
diff --git a/Utilities/otbossim/src/ossim/vec/ossimShapeDatabase.cpp b/Utilities/otbossim/src/ossim/vec/ossimShapeDatabase.cpp
deleted file mode 100644
index 5ed5a2831d..0000000000
--- a/Utilities/otbossim/src/ossim/vec/ossimShapeDatabase.cpp
+++ /dev/null
@@ -1,268 +0,0 @@
-
-#include <iomanip>
-
-#include <ossim/vec/ossimShapeDatabase.h>
-#include <ossim/vec/ossimShapeFile.h>
-
-
-
-ossimString ossimShapeDatabaseField::fieldTypeAsString()const
-{
- switch(theFieldType)
- {
- case FTString:
- {
- return "String";
- }
- case FTInteger:
- {
- return "Integer";
- }
- case FTDouble:
- {
- return "Double";
- }
- default:
- {
- return "Unknown";
- }
- };
-
- return "";
-
-}
-
-
-ossim_int32 ossimShapeDatabaseRecord::getFieldIdx(const ossimString& name,
- bool caseInsensitive)const
-{
- ossimString searchString = name;
- if(caseInsensitive) searchString = searchString.downcase();
- ossim_int32 idx = 0;
- for(idx = 0; idx < (int)theFieldArray.size(); ++idx)
- {
- if(caseInsensitive)
- {
- if(ossimString(theFieldArray[idx].theName).downcase() == searchString)
- {
- return idx;
- }
- }
- else
- {
- if(theFieldArray[idx].theName == searchString)
- {
- return idx;
- }
- }
- }
- return -1;
-}
-
-ossimShapeDatabase::ossimShapeDatabase()
- :theHandle(NULL),
- theFilename("")
-{
- theRecordNumber = -1;
-}
-
-ossimShapeDatabase::~ossimShapeDatabase()
-{
- close();
-}
-
-bool ossimShapeDatabase::open(const ossimFilename& file,
- const ossimString& flags)
-{
- if(isOpen()) close();
-
- theHandle = ossim_DBFOpen(file.c_str(), flags.c_str());
- if(theHandle)
- {
- theFilename = file;
- theRecordNumber = -1;
- }
-
- return (theHandle != NULL);
-}
-
-void ossimShapeDatabase::close()
-{
- if(isOpen())
- {
- ossim_DBFClose(theHandle);
- theHandle = NULL;
- theRecordNumber = -1;
- }
-}
-int ossimShapeDatabase::getNumberOfRecords()const
-{
- if(isOpen())
- {
- return theHandle->nRecords;
- }
-
- return 0;
-}
-bool ossimShapeDatabase::getRecord(ossimShapeDatabaseRecord& result)
-{
- if(isOpen()&&( (theRecordNumber < theHandle->nRecords) ))
- {
- if(result.getNumberOfFields() != theHandle->nFields)
- {
- result.setNumberOfFields(theHandle->nFields);
- }
-
- char name[1024] = {'\0'};
- int width = 0;
- int decimals = 0;
- int iField = 0;
- std::vector<int> fieldWidths;
-
- for(iField = 0; iField < theHandle->nFields; ++iField)
- {
- ossim_DBFFieldType fieldType = ossim_DBFGetFieldInfo(theHandle,
- iField,
- name,
- &width,
- &decimals);
- ossimShapeDatabaseField field;
- field.theName = name;
- field.theWidth = width;
- field.theDecimals = decimals;
- field.theFieldType = fieldType;
-
- ossimString key = "field";
- key+=ossimString::toString(iField+1);
- key+=(ossimString(".")+name+":");
-
- switch(fieldType)
- {
- case FTString:
- {
- field.theValue = ossim_DBFReadStringAttribute(theHandle, theRecordNumber, iField);
- break;
- }
- case FTInteger:
- {
- field.theValue = ossimString::toString(ossim_DBFReadIntegerAttribute(theHandle, theRecordNumber, iField));
- break;
- }
- case FTDouble:
- {
- field.theValue = ossimString::toString(ossim_DBFReadDoubleAttribute(theHandle, theRecordNumber, iField));
- break;
- }
- case FTLogical:
- {
- break;
- }
- case FTInvalid:
- {
- break;
- }
- }
-
- result.setField(field,
- iField);
- }
- return true;
- }
-
- return false;
-}
-
-void ossimShapeDatabase::print(std::ostream& out)const
-{
- if(isOpen())
- {
- out << std::setw(15)<<setiosflags(std::ios::left)
- <<"DB filename:" << theFilename << std::endl
- << std::setw(15)<<setiosflags(std::ios::left)
- <<"records:" << theHandle->nRecords << std::endl
- << std::setw(15)<<setiosflags(std::ios::left)
- <<"fields:" << theHandle->nFields << std::endl;
-
- char name[1024] = {'\0'};
- int width = 0;
- int decimals = 0;
- int iField = 0;
- std::vector<int> fieldWidths;
-
-
- for(iField = 0; iField < theHandle->nFields; ++iField)
- {
- ossim_DBFFieldType fieldType = ossim_DBFGetFieldInfo(theHandle,
- iField,
- name,
- &width,
- &decimals);
- ossimString s = "field " + ossimString::toString(iField+1) + " name:";
- switch(fieldType)
- {
- case FTString:
- case FTInteger:
- case FTDouble:
- {
- out << std::setw(15) << setiosflags(std::ios::left) << s.c_str() << name << std::endl;
- break;
- }
- default:
- {
- out << std::setw(15) << setiosflags(std::ios::left) << s.c_str() << "INVALID"<<std::endl;
- break;
- }
- }
- }
- for(int iShape = 0; iShape < theHandle->nRecords; ++iShape)
- {
- for(iField = 0; iField < theHandle->nFields; ++iField)
- {
- ossim_DBFFieldType fieldType = ossim_DBFGetFieldInfo(theHandle,
- iField,
- name,
- &width,
- &decimals);
-
- ossimString key = "field";
- key+=ossimString::toString(iField+1);
- key+=(ossimString(".")+name+":");
-
- switch(fieldType)
- {
- case FTString:
- {
-
- out << std::setw(25) << setiosflags(std::ios::left) << key.c_str()
- << ossim_DBFReadStringAttribute(theHandle, iShape, iField) <<std::endl;
-
- break;
- }
- case FTInteger:
- {
- out << std::setw(25) << setiosflags(std::ios::left) << key.c_str()
- << ossim_DBFReadIntegerAttribute(theHandle, iShape, iField) << std::endl;
-
- break;
- }
- case FTDouble:
- {
- out << std::setw(25) << setiosflags(std::ios::left) << key.c_str()
- << ossim_DBFReadDoubleAttribute(theHandle, iShape, iField) << std::endl;
-
- break;
- }
- case FTLogical:
- {
- break;
- }
- case FTInvalid:
- {
- break;
- }
- }
- }
- out << "_________________________________"<<std::endl;
- }
- }
-}
diff --git a/Utilities/otbossim/src/ossim/vec/ossimShapeFile.cpp b/Utilities/otbossim/src/ossim/vec/ossimShapeFile.cpp
deleted file mode 100644
index 2b15932ea2..0000000000
--- a/Utilities/otbossim/src/ossim/vec/ossimShapeFile.cpp
+++ /dev/null
@@ -1,220 +0,0 @@
-//*******************************************************************
-// License: See top level LICENSE.txt file.
-//
-// Author: Garrett Potts
-//
-// Description: This class extends the stl's string class.
-//
-//********************************************************************
-// $Id: ossimShapeFile.cpp 9963 2006-11-28 21:11:01Z gpotts $
-//
-
-#include <iomanip>
-#include <ossim/vec/ossimShapeFile.h>
-
-
-
-std::ostream& operator<<(std::ostream& out,
- const ossimShapeObject& rhs)
-{
- if(rhs.isLoaded())
- {
- double minx, miny, minz, minm;
- double maxx, maxy, maxz, maxm;
- ossim_uint32 i = 0;
-
- rhs.getBounds(minx, miny, minz, minm,
- maxx, maxy, maxz, maxm);
-
- out << std::setw(15) << setiosflags(std::ios::left)<<"type:"<<rhs.getTypeByName() << std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"id:"<<rhs.getId()<<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"minx:"<<minx <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"miny:"<<miny <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"minz:"<<minz <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"minm:"<<minm <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"maxx:"<<maxx <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"maxy:"<<maxy <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"maxz:"<<maxz <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"maxm:"<<maxm <<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"parts:"<<rhs.getNumberOfParts()<<std::endl
- << std::setw(15) << setiosflags(std::ios::left)<<"vertices:"<<rhs.getNumberOfVertices();
-
- if(rhs.getNumberOfParts())
- {
- out << std::endl;
- for(i = 0; i < rhs.getNumberOfParts()-1; ++i)
- {
- ossimString s1 = "part start ";
- ossimString s2 = "part type ";
- s1 += (ossimString::toString(i+1)+":");
- s2 += (ossimString::toString(i+1)+":");
-
- out << std::setw(15) << setiosflags(std::ios::left) << s1.c_str() << rhs.theShape->panPartStart[i]<<std::endl;
- out << std::setw(15) << setiosflags(std::ios::left) << s2.c_str() << ossim_SHPPartTypeName(rhs.theShape->panPartType[i])<<std::endl;
- }
- ossimString s1 = "part start ";
- ossimString s2 = "part type ";
-
- s1 += (ossimString::toString(rhs.getNumberOfParts())+":");
- out << std::setw(15) << setiosflags(std::ios::left) << s1.c_str() << rhs.theShape->panPartStart[i]<<std::endl;
- out << std::setw(15) << setiosflags(std::ios::left) << s2.c_str() << ossim_SHPPartTypeName(rhs.theShape->panPartType[i]);
- }
-
- out << std::setw(0);
-
- }
- return out;
-}
-
-void ossimShapeObject::getBoundingRect(ossimDrect& result,
- ossimCoordSysOrientMode orient)const
-{
- double minx, miny, maxx, maxy;
-
- if(theShape)
- {
- getBounds(minx, miny, maxx, maxy);
-
- if(orient == OSSIM_RIGHT_HANDED)
- {
- result = ossimDrect(minx, maxy, maxx, miny, orient);
- }
- else
- {
- result = ossimDrect(minx, miny, maxx, maxy, orient);
- }
- }
- else
- {
- result = ossimDrect(0,0,0,0,orient);
- result.makeNan();
- }
-}
-
-void ossimShapeObject::setShape(ossim_SHPObject* obj)
-{
- if(theShape)
- {
- ossim_SHPDestroyObject(theShape);
- theShape = NULL;
- }
-
- theShape = obj;
-}
-
-bool ossimShapeObject::loadShape(const ossimShapeFile& shapeFile,
- long shapeRecord)
-{
- if(theShape)
- {
- ossim_SHPDestroyObject(theShape);
- theShape = NULL;
- }
- if(shapeFile.isOpen())
- {
- theShape = ossim_SHPReadObject(shapeFile.getHandle(),
- shapeRecord);
- theIndex = shapeRecord;
- }
- else
- {
- theIndex = -1;
- }
- return (theShape != (ossim_SHPObject*)NULL);
-}
-
-
-ossimShapeFile::ossimShapeFile()
- :theHandle(NULL)
-{
-}
-
-ossimShapeFile::~ossimShapeFile()
-{
- close();
-}
-
-void ossimShapeFile::close()
-{
- if(theHandle)
- {
- ossim_SHPClose(theHandle);
- theHandle = NULL;
- }
-}
-
-bool ossimShapeFile::open(const ossimFilename& file,
- const ossimString& flags)
-{
- if(isOpen()) close();
-
- theHandle = ossim_SHPOpen( file.c_str(),
- flags.c_str());
- if(isOpen())
- {
- theFilename = file;
- }
-
- return (theHandle != NULL);
-}
-
-
-
-void ossimShapeFile::print(std::ostream& out)const
-{
- if(isOpen())
- {
- out << std::setw(15) << setiosflags(std::ios::left)<<"Shp filename:" << theFilename << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"Record count:" << theHandle->nRecords << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"File type:" << getShapeTypeString().c_str() << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"minx:" << theHandle->adBoundsMin[0] << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"miny:" << theHandle->adBoundsMin[1] << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"minz:" << theHandle->adBoundsMin[2] << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"minm:" << theHandle->adBoundsMin[3] << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"maxx:" << theHandle->adBoundsMax[0] << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"maxy:" << theHandle->adBoundsMax[1] << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"maxz:" << theHandle->adBoundsMax[2] << std::endl;
- out << std::setw(15) << setiosflags(std::ios::left)<<"maxm:" << theHandle->adBoundsMax[3] << std::endl;
-
- ossimShapeObject shape;
-
- if(theHandle->nRecords)
- {
- out << std::setw(30) << std::setfill('_') << "" << std::setfill(' ')<<std::endl;
- }
- for(int i=0; i < theHandle->nRecords; ++i)
- {
- if(shape.loadShape(*this, i))
- {
- out << shape << std::endl;
- out << std::setw(30) << std::setfill('_') << ""<<std::setfill(' ')<<std::endl;
- }
- }
- }
-}
-
-
-ossimString ossimShapeFile::getShapeTypeString()const
-{
- if(theHandle)
- {
- return ossim_SHPTypeName(theHandle->nShapeType);
- }
- return "";
-}
-
-void ossimShapeFile::getBoundingRect(ossimDrect& result,
- ossimCoordSysOrientMode orient)const
-{
- double minx, miny, maxx, maxy;
-
- getBounds(minx, miny, maxx, maxy);
- if(orient == OSSIM_RIGHT_HANDED)
- {
- result = ossimDrect(minx, maxy, maxx, miny, orient);
- }
- else
- {
- result = ossimDrect(minx, miny, maxx, maxy, orient);
- }
-}
diff --git a/Utilities/otbossim/src/ossim/vec/shpopen.c b/Utilities/otbossim/src/ossim/vec/shpopen.c
deleted file mode 100644
index a8aa5a1e10..0000000000
--- a/Utilities/otbossim/src/ossim/vec/shpopen.c
+++ /dev/null
@@ -1,2021 +0,0 @@
-/******************************************************************************
- * $Id: shpopen.c 9094 2006-06-13 19:12:40Z dburken $
- *
- * Project: Shapelib
- * Purpose: Implementation of core Shapefile read/write functions.
- * Author: Frank Warmerdam, warmerdam@pobox.com
- *
- ******************************************************************************
- * Copyright (c) 1999, 2001, Frank Warmerdam
- *
- * This software is available under the following "MIT Style" license,
- * or at the option of the licensee under the LGPL (see LICENSE.LGPL). 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$
- * Revision 1.6 2006/06/13 19:11:33 dburken
- * Global include path change.
- *
- * Revision 1.5 2005/10/08 12:57:26 gpotts
- * Added adiitional patches
- *
- * Revision 1.46 2005/02/11 17:17:46 fwarmerdam
- * added panPartStart[0] validation
- *
- * Revision 1.45 2004/09/26 20:09:48 fwarmerdam
- * const correctness changes
- *
- * Revision 1.44 2003/12/29 00:18:39 fwarmerdam
- * added error checking for failed IO and optional CPL error reporting
- *
- * Revision 1.43 2003/12/01 16:20:08 warmerda
- * be careful of zero vertex shapes
- *
- * Revision 1.42 2003/12/01 14:58:27 warmerda
- * added degenerate object check in ossim_SHPRewindObject()
- *
- * Revision 1.41 2003/07/08 15:22:43 warmerda
- * avoid warning
- *
- * Revision 1.40 2003/04/21 18:30:37 warmerda
- * added header write/update public methods
- *
- * Revision 1.39 2002/08/26 06:46:56 warmerda
- * avoid c++ comments
- *
- * Revision 1.38 2002/05/07 16:43:39 warmerda
- * Removed debugging printf.
- *
- * Revision 1.37 2002/04/10 17:35:22 warmerda
- * fixed bug in ring reversal code
- *
- * Revision 1.36 2002/04/10 16:59:54 warmerda
- * added ossim_SHPRewindObject
- *
- * Revision 1.35 2001/12/07 15:10:44 warmerda
- * fix if .shx fails to open
- *
- * Revision 1.34 2001/11/01 16:29:55 warmerda
- * move pabyRec into ossim_SHPInfo for thread safety
- *
- * Revision 1.33 2001/07/03 12:18:15 warmerda
- * Improved cleanup if SHX not found, provied by Riccardo Cohen.
- *
- * Revision 1.32 2001/06/22 01:58:07 warmerda
- * be more careful about establishing initial bounds in face of NULL shapes
- *
- * Revision 1.31 2001/05/31 19:35:29 warmerda
- * added support for writing null shapes
- *
- * Revision 1.30 2001/05/28 12:46:29 warmerda
- * Add some checking on reasonableness of record count when opening.
- *
- * Revision 1.29 2001/05/23 13:36:52 warmerda
- * added use of ossim_SHPAPI_CALL
- *
- * Revision 1.28 2001/02/06 22:25:06 warmerda
- * fixed memory leaks when ossim_SHPOpen() fails
- *
- * Revision 1.27 2000/07/18 15:21:33 warmerda
- * added better enforcement of -1 for append in ossim_SHPWriteObject
- *
- * Revision 1.26 2000/02/16 16:03:51 warmerda
- * added null shape support
- *
- * Revision 1.25 1999/12/15 13:47:07 warmerda
- * Fixed record size settings in .shp file (was 4 words too long)
- * Added stdlib.h.
- *
- * Revision 1.24 1999/11/05 14:12:04 warmerda
- * updated license terms
- *
- * Revision 1.23 1999/07/27 00:53:46 warmerda
- * added support for rewriting shapes
- *
- * Revision 1.22 1999/06/11 19:19:11 warmerda
- * Cleanup pabyRec static buffer on ossim_SHPClose().
- *
- * Revision 1.21 1999/06/02 14:57:56 kshih
- * Remove unused variables
- *
- * Revision 1.20 1999/04/19 21:04:17 warmerda
- * Fixed syntax error.
- *
- * Revision 1.19 1999/04/19 21:01:57 warmerda
- * Force access string to binary in ossim_SHPOpen().
- *
- * Revision 1.18 1999/04/01 18:48:07 warmerda
- * Try upper case extensions if lower case doesn't work.
- *
- * Revision 1.17 1998/12/31 15:29:39 warmerda
- * Disable writing measure values to multipatch objects if
- * DISABLE_MULTIPATCH_MEASURE is defined.
- *
- * Revision 1.16 1998/12/16 05:14:33 warmerda
- * Added support to write MULTIPATCH. Fixed reading Z coordinate of
- * MULTIPATCH. Fixed record size written for all feature types.
- *
- * Revision 1.15 1998/12/03 16:35:29 warmerda
- * r+b is proper binary access string, not rb+.
- *
- * Revision 1.14 1998/12/03 15:47:56 warmerda
- * Fixed setting of nVertices in ossim_SHPCreateObject().
- *
- * Revision 1.13 1998/12/03 15:33:54 warmerda
- * Made ossim_SHPCalculateExtents() separately callable.
- *
- * Revision 1.12 1998/11/11 20:01:50 warmerda
- * Fixed bug writing ArcM/Z, and PolygonM/Z for big endian machines.
- *
- * Revision 1.11 1998/11/09 20:56:44 warmerda
- * Fixed up handling of file wide bounds.
- *
- * Revision 1.10 1998/11/09 20:18:51 warmerda
- * Converted to support 3D shapefiles, and use of ossim_SHPObject.
- *
- * Revision 1.9 1998/02/24 15:09:05 warmerda
- * Fixed memory leak.
- *
- * Revision 1.8 1997/12/04 15:40:29 warmerda
- * Fixed byte swapping of record number, and record length fields in the
- * .shp file.
- *
- * Revision 1.7 1995/10/21 03:15:58 warmerda
- * Added support for binary file access, the magic cookie 9997
- * and tried to improve the int32 selection logic for 16bit systems.
- *
- * Revision 1.6 1995/09/04 04:19:41 warmerda
- * Added fix for file bounds.
- *
- * Revision 1.5 1995/08/25 15:16:44 warmerda
- * Fixed a couple of problems with big endian systems ... one with bounds
- * and the other with multipart polygons.
- *
- * Revision 1.4 1995/08/24 18:10:17 warmerda
- * Switch to use SfRealloc() to avoid problems with pre-ANSI realloc()
- * functions (such as on the Sun).
- *
- * Revision 1.3 1995/08/23 02:23:15 warmerda
- * Added support for reading bounds, and fixed up problems in setting the
- * file wide bounds.
- *
- * Revision 1.2 1995/08/04 03:16:57 warmerda
- * Added header.
- *
- */
-
-#include <ossim/vec/shapefil.h>
-
-#include <math.h>
-#include <limits.h>
-#include <assert.h>
-#include <stdlib.h>
-#include <string.h>
-#include <stdio.h>
-
-ossim_SHP_CVSID("$Id: shpopen.c 9094 2006-06-13 19:12:40Z dburken $")
-
-typedef unsigned char uchar;
-
-#if UINT_MAX == 65535
-typedef long int32;
-#else
-typedef int int32;
-#endif
-
-#ifndef FALSE
-# define FALSE 0
-# define TRUE 1
-#endif
-
-#define ByteCopy( a, b, c ) memcpy( b, a, c )
-#ifndef MAX
-# define MIN(a,b) ((a<b) ? a : b)
-# define MAX(a,b) ((a>b) ? a : b)
-#endif
-
-static int bBigEndian;
-
-
-/************************************************************************/
-/* SwapWord() */
-/* */
-/* Swap a 2, 4 or 8 byte word. */
-/************************************************************************/
-
-static void SwapWord( int length, void * wordP )
-
-{
- int i;
- uchar temp;
-
- for( i=0; i < length/2; i++ )
- {
- temp = ((uchar *) wordP)[i];
- ((uchar *)wordP)[i] = ((uchar *) wordP)[length-i-1];
- ((uchar *) wordP)[length-i-1] = temp;
- }
-}
-
-/************************************************************************/
-/* SfRealloc() */
-/* */
-/* A realloc cover function that will access a NULL pointer as */
-/* a valid input. */
-/************************************************************************/
-
-static void * SfRealloc( void * pMem, int nNewSize )
-
-{
- if( pMem == NULL )
- return( (void *) malloc(nNewSize) );
- else
- return( (void *) realloc(pMem,nNewSize) );
-}
-
-/************************************************************************/
-/* ossim_SHPWriteHeader() */
-/* */
-/* Write out a header for the .shp and .shx files as well as the */
-/* contents of the index (.shx) file. */
-/************************************************************************/
-
-void ossim_SHPWriteHeader( ossim_SHPHandle psSHP )
-
-{
- uchar abyHeader[100];
- int i;
- int32 i32;
- double dValue;
- int32 *panSHX;
-
-/* -------------------------------------------------------------------- */
-/* Prepare header block for .shp file. */
-/* -------------------------------------------------------------------- */
- for( i = 0; i < 100; i++ )
- abyHeader[i] = 0;
-
- abyHeader[2] = 0x27; /* magic cookie */
- abyHeader[3] = 0x0a;
-
- i32 = psSHP->nFileSize/2; /* file size */
- ByteCopy( &i32, abyHeader+24, 4 );
- if( !bBigEndian ) SwapWord( 4, abyHeader+24 );
-
- i32 = 1000; /* version */
- ByteCopy( &i32, abyHeader+28, 4 );
- if( bBigEndian ) SwapWord( 4, abyHeader+28 );
-
- i32 = psSHP->nShapeType; /* shape type */
- ByteCopy( &i32, abyHeader+32, 4 );
- if( bBigEndian ) SwapWord( 4, abyHeader+32 );
-
- dValue = psSHP->adBoundsMin[0]; /* set bounds */
- ByteCopy( &dValue, abyHeader+36, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+36 );
-
- dValue = psSHP->adBoundsMin[1];
- ByteCopy( &dValue, abyHeader+44, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+44 );
-
- dValue = psSHP->adBoundsMax[0];
- ByteCopy( &dValue, abyHeader+52, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+52 );
-
- dValue = psSHP->adBoundsMax[1];
- ByteCopy( &dValue, abyHeader+60, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+60 );
-
- dValue = psSHP->adBoundsMin[2]; /* z */
- ByteCopy( &dValue, abyHeader+68, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+68 );
-
- dValue = psSHP->adBoundsMax[2];
- ByteCopy( &dValue, abyHeader+76, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+76 );
-
- dValue = psSHP->adBoundsMin[3]; /* m */
- ByteCopy( &dValue, abyHeader+84, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+84 );
-
- dValue = psSHP->adBoundsMax[3];
- ByteCopy( &dValue, abyHeader+92, 8 );
- if( bBigEndian ) SwapWord( 8, abyHeader+92 );
-
-/* -------------------------------------------------------------------- */
-/* Write .shp file header. */
-/* -------------------------------------------------------------------- */
- if( fseek( psSHP->fpSHP, 0, 0 ) != 0
- || fwrite( abyHeader, 100, 1, psSHP->fpSHP ) != 1 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_OpenFailed,
- "Failure writing .shp header." );
-#endif
- return;
- }
-
-/* -------------------------------------------------------------------- */
-/* Prepare, and write .shx file header. */
-/* -------------------------------------------------------------------- */
- i32 = (psSHP->nRecords * 2 * sizeof(int32) + 100)/2; /* file size */
- ByteCopy( &i32, abyHeader+24, 4 );
- if( !bBigEndian ) SwapWord( 4, abyHeader+24 );
-
- if( fseek( psSHP->fpSHX, 0, 0 ) != 0
- || fwrite( abyHeader, 100, 1, psSHP->fpSHX ) != 1 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_OpenFailed,
- "Failure writing .shx header." );
-#endif
- return;
- }
-
-/* -------------------------------------------------------------------- */
-/* Write out the .shx contents. */
-/* -------------------------------------------------------------------- */
- panSHX = (int32 *) malloc(sizeof(int32) * 2 * psSHP->nRecords);
-
- for( i = 0; i < psSHP->nRecords; i++ )
- {
- panSHX[i*2 ] = psSHP->panRecOffset[i]/2;
- panSHX[i*2+1] = psSHP->panRecSize[i]/2;
- if( !bBigEndian ) SwapWord( 4, panSHX+i*2 );
- if( !bBigEndian ) SwapWord( 4, panSHX+i*2+1 );
- }
-
- if( (int)fwrite( panSHX, sizeof(int32)*2, psSHP->nRecords, psSHP->fpSHX )
- != psSHP->nRecords )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_OpenFailed,
- "Failure writing .shx contents." );
-#endif
- }
-
- free( panSHX );
-
-/* -------------------------------------------------------------------- */
-/* Flush to disk. */
-/* -------------------------------------------------------------------- */
- fflush( psSHP->fpSHP );
- fflush( psSHP->fpSHX );
-}
-
-/************************************************************************/
-/* shpopen() */
-/* */
-/* Open the .shp and .shx files based on the basename of the */
-/* files or either file name. */
-/************************************************************************/
-
-ossim_SHPHandle ossim_SHPAPI_CALL
- ossim_SHPOpen( const char * pszLayer, const char * pszAccess )
-
-{
- char *pszFullname, *pszBasename;
- ossim_SHPHandle psSHP;
-
- uchar *pabyBuf;
- int i;
- double dValue;
-
-/* -------------------------------------------------------------------- */
-/* Ensure the access string is one of the legal ones. We */
-/* ensure the result string indicates binary to avoid common */
-/* problems on Windows. */
-/* -------------------------------------------------------------------- */
- if( strcmp(pszAccess,"rb+") == 0 || strcmp(pszAccess,"r+b") == 0
- || strcmp(pszAccess,"r+") == 0 )
- pszAccess = "r+b";
- else
- pszAccess = "rb";
-
-/* -------------------------------------------------------------------- */
-/* Establish the byte order on this machine. */
-/* -------------------------------------------------------------------- */
- i = 1;
- if( *((uchar *) &i) == 1 )
- bBigEndian = FALSE;
- else
- bBigEndian = TRUE;
-
-/* -------------------------------------------------------------------- */
-/* Initialize the info structure. */
-/* -------------------------------------------------------------------- */
- psSHP = (ossim_SHPHandle) calloc(sizeof(ossim_SHPInfo),1);
-
- psSHP->bUpdated = FALSE;
-
-/* -------------------------------------------------------------------- */
-/* Compute the base (layer) name. If there is any extension */
-/* on the passed in filename we will strip it off. */
-/* -------------------------------------------------------------------- */
- pszBasename = (char *) malloc(strlen(pszLayer)+5);
- strcpy( pszBasename, pszLayer );
- for( i = strlen(pszBasename)-1;
- i > 0 && pszBasename[i] != '.' && pszBasename[i] != '/'
- && pszBasename[i] != '\\';
- i-- ) {}
-
- if( pszBasename[i] == '.' )
- pszBasename[i] = '\0';
-
-/* -------------------------------------------------------------------- */
-/* Open the .shp and .shx files. Note that files pulled from */
-/* a PC to Unix with upper case filenames won't work! */
-/* -------------------------------------------------------------------- */
- pszFullname = (char *) malloc(strlen(pszBasename) + 5);
- sprintf( pszFullname, "%s.shp", pszBasename );
- psSHP->fpSHP = fopen(pszFullname, pszAccess );
- if( psSHP->fpSHP == NULL )
- {
- sprintf( pszFullname, "%s.SHP", pszBasename );
- psSHP->fpSHP = fopen(pszFullname, pszAccess );
- }
-
- if( psSHP->fpSHP == NULL )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_OpenFailed,
- "Unable to open %s.shp or %s.SHP.",
- pszBasename, pszBasename );
-#endif
- free( psSHP );
- free( pszBasename );
- free( pszFullname );
- return( NULL );
- }
-
- sprintf( pszFullname, "%s.shx", pszBasename );
- psSHP->fpSHX = fopen(pszFullname, pszAccess );
- if( psSHP->fpSHX == NULL )
- {
- sprintf( pszFullname, "%s.SHX", pszBasename );
- psSHP->fpSHX = fopen(pszFullname, pszAccess );
- }
-
- if( psSHP->fpSHX == NULL )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_OpenFailed,
- "Unable to open %s.shx or %s.SHX.",
- pszBasename, pszBasename );
-#endif
- fclose( psSHP->fpSHP );
- free( psSHP );
- free( pszBasename );
- free( pszFullname );
- return( NULL );
- }
-
- free( pszFullname );
- free( pszBasename );
-
-/* -------------------------------------------------------------------- */
-/* Read the file size from the ossim_SHP file. */
-/* -------------------------------------------------------------------- */
- pabyBuf = (uchar *) malloc(100);
- fread( pabyBuf, 100, 1, psSHP->fpSHP );
-
- psSHP->nFileSize = (pabyBuf[24] * 256 * 256 * 256
- + pabyBuf[25] * 256 * 256
- + pabyBuf[26] * 256
- + pabyBuf[27]) * 2;
-
-/* -------------------------------------------------------------------- */
-/* Read SHX file Header info */
-/* -------------------------------------------------------------------- */
- if( fread( pabyBuf, 100, 1, psSHP->fpSHX ) != 1
- || pabyBuf[0] != 0
- || pabyBuf[1] != 0
- || pabyBuf[2] != 0x27
- || (pabyBuf[3] != 0x0a && pabyBuf[3] != 0x0d) )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- ".shx file is unreadable, or corrupt." );
-#endif
- fclose( psSHP->fpSHP );
- fclose( psSHP->fpSHX );
- free( psSHP );
-
- return( NULL );
- }
-
- psSHP->nRecords = pabyBuf[27] + pabyBuf[26] * 256
- + pabyBuf[25] * 256 * 256 + pabyBuf[24] * 256 * 256 * 256;
- psSHP->nRecords = (psSHP->nRecords*2 - 100) / 8;
-
- psSHP->nShapeType = pabyBuf[32];
-
- if( psSHP->nRecords < 0 || psSHP->nRecords > 256000000 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- "Record count in .shp header is %d, which seems\n"
- "unreasonable. Assuming header is corrupt.",
- psSHP->nRecords );
-#endif
- fclose( psSHP->fpSHP );
- fclose( psSHP->fpSHX );
- free( psSHP );
-
- return( NULL );
- }
-
-/* -------------------------------------------------------------------- */
-/* Read the bounds. */
-/* -------------------------------------------------------------------- */
- if( bBigEndian ) SwapWord( 8, pabyBuf+36 );
- memcpy( &dValue, pabyBuf+36, 8 );
- psSHP->adBoundsMin[0] = dValue;
-
- if( bBigEndian ) SwapWord( 8, pabyBuf+44 );
- memcpy( &dValue, pabyBuf+44, 8 );
- psSHP->adBoundsMin[1] = dValue;
-
- if( bBigEndian ) SwapWord( 8, pabyBuf+52 );
- memcpy( &dValue, pabyBuf+52, 8 );
- psSHP->adBoundsMax[0] = dValue;
-
- if( bBigEndian ) SwapWord( 8, pabyBuf+60 );
- memcpy( &dValue, pabyBuf+60, 8 );
- psSHP->adBoundsMax[1] = dValue;
-
- if( bBigEndian ) SwapWord( 8, pabyBuf+68 ); /* z */
- memcpy( &dValue, pabyBuf+68, 8 );
- psSHP->adBoundsMin[2] = dValue;
-
- if( bBigEndian ) SwapWord( 8, pabyBuf+76 );
- memcpy( &dValue, pabyBuf+76, 8 );
- psSHP->adBoundsMax[2] = dValue;
-
- if( bBigEndian ) SwapWord( 8, pabyBuf+84 ); /* z */
- memcpy( &dValue, pabyBuf+84, 8 );
- psSHP->adBoundsMin[3] = dValue;
-
- if( bBigEndian ) SwapWord( 8, pabyBuf+92 );
- memcpy( &dValue, pabyBuf+92, 8 );
- psSHP->adBoundsMax[3] = dValue;
-
- free( pabyBuf );
-
-/* -------------------------------------------------------------------- */
-/* Read the .shx file to get the offsets to each record in */
-/* the .shp file. */
-/* -------------------------------------------------------------------- */
- psSHP->nMaxRecords = psSHP->nRecords;
-
- psSHP->panRecOffset =
- (int *) malloc(sizeof(int) * MAX(1,psSHP->nMaxRecords) );
- psSHP->panRecSize =
- (int *) malloc(sizeof(int) * MAX(1,psSHP->nMaxRecords) );
-
- pabyBuf = (uchar *) malloc(8 * MAX(1,psSHP->nRecords) );
- if( (int) fread( pabyBuf, 8, psSHP->nRecords, psSHP->fpSHX )
- != psSHP->nRecords )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- "Failed to read all values for %d records in .shx file.",
- psSHP->nRecords );
-#endif
- /* SHX is short or unreadable for some reason. */
- fclose( psSHP->fpSHP );
- fclose( psSHP->fpSHX );
- free( psSHP->panRecOffset );
- free( psSHP->panRecSize );
- free( psSHP );
-
- return( NULL );
- }
-
- for( i = 0; i < psSHP->nRecords; i++ )
- {
- int32 nOffset, nLength;
-
- memcpy( &nOffset, pabyBuf + i * 8, 4 );
- if( !bBigEndian ) SwapWord( 4, &nOffset );
-
- memcpy( &nLength, pabyBuf + i * 8 + 4, 4 );
- if( !bBigEndian ) SwapWord( 4, &nLength );
-
- psSHP->panRecOffset[i] = nOffset*2;
- psSHP->panRecSize[i] = nLength*2;
- }
- free( pabyBuf );
-
- return( psSHP );
-}
-
-/************************************************************************/
-/* ossim_SHPClose() */
-/* */
-/* Close the .shp and .shx files. */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPClose(ossim_SHPHandle psSHP )
-
-{
- if( psSHP == NULL )
- return;
-
-/* -------------------------------------------------------------------- */
-/* Update the header if we have modified anything. */
-/* -------------------------------------------------------------------- */
- if( psSHP->bUpdated )
- ossim_SHPWriteHeader( psSHP );
-
-/* -------------------------------------------------------------------- */
-/* Free all resources, and close files. */
-/* -------------------------------------------------------------------- */
- free( psSHP->panRecOffset );
- free( psSHP->panRecSize );
-
- fclose( psSHP->fpSHX );
- fclose( psSHP->fpSHP );
-
- if( psSHP->pabyRec != NULL )
- {
- free( psSHP->pabyRec );
- }
-
- free( psSHP );
-}
-
-/************************************************************************/
-/* ossim_SHPGetInfo() */
-/* */
-/* Fetch general information about the shape file. */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPGetInfo(ossim_SHPHandle psSHP, int * pnEntities, int * pnShapeType,
- double * padfMinBound, double * padfMaxBound )
-
-{
- int i;
-
- if( psSHP == NULL )
- return;
-
- if( pnEntities != NULL )
- *pnEntities = psSHP->nRecords;
-
- if( pnShapeType != NULL )
- *pnShapeType = psSHP->nShapeType;
-
- for( i = 0; i < 4; i++ )
- {
- if( padfMinBound != NULL )
- padfMinBound[i] = psSHP->adBoundsMin[i];
- if( padfMaxBound != NULL )
- padfMaxBound[i] = psSHP->adBoundsMax[i];
- }
-}
-
-/************************************************************************/
-/* ossim_SHPCreate() */
-/* */
-/* Create a new shape file and return a handle to the open */
-/* shape file with read/write access. */
-/************************************************************************/
-
-ossim_SHPHandle ossim_SHPAPI_CALL
- ossim_SHPCreate( const char * pszLayer, int nShapeType )
-
-{
- char *pszBasename, *pszFullname;
- int i;
- FILE *fpSHP, *fpSHX;
- uchar abyHeader[100];
- int32 i32;
- double dValue;
-
-/* -------------------------------------------------------------------- */
-/* Establish the byte order on this system. */
-/* -------------------------------------------------------------------- */
- i = 1;
- if( *((uchar *) &i) == 1 )
- bBigEndian = FALSE;
- else
- bBigEndian = TRUE;
-
-/* -------------------------------------------------------------------- */
-/* Compute the base (layer) name. If there is any extension */
-/* on the passed in filename we will strip it off. */
-/* -------------------------------------------------------------------- */
- pszBasename = (char *) malloc(strlen(pszLayer)+5);
- strcpy( pszBasename, pszLayer );
- for( i = strlen(pszBasename)-1;
- i > 0 && pszBasename[i] != '.' && pszBasename[i] != '/'
- && pszBasename[i] != '\\';
- i-- ) {}
-
- if( pszBasename[i] == '.' )
- pszBasename[i] = '\0';
-
-/* -------------------------------------------------------------------- */
-/* Open the two files so we can write their headers. */
-/* -------------------------------------------------------------------- */
- pszFullname = (char *) malloc(strlen(pszBasename) + 5);
- sprintf( pszFullname, "%s.shp", pszBasename );
- fpSHP = fopen(pszFullname, "wb" );
- if( fpSHP == NULL )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- "Failed to create file %s.",
- pszFullname );
-#endif
- return( NULL );
- }
-
- sprintf( pszFullname, "%s.shx", pszBasename );
- fpSHX = fopen(pszFullname, "wb" );
- if( fpSHX == NULL )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- "Failed to create file %s.",
- pszFullname );
-#endif
- return( NULL );
- }
-
- free( pszFullname );
- free( pszBasename );
-
-/* -------------------------------------------------------------------- */
-/* Prepare header block for .shp file. */
-/* -------------------------------------------------------------------- */
- for( i = 0; i < 100; i++ )
- abyHeader[i] = 0;
-
- abyHeader[2] = 0x27; /* magic cookie */
- abyHeader[3] = 0x0a;
-
- i32 = 50; /* file size */
- ByteCopy( &i32, abyHeader+24, 4 );
- if( !bBigEndian ) SwapWord( 4, abyHeader+24 );
-
- i32 = 1000; /* version */
- ByteCopy( &i32, abyHeader+28, 4 );
- if( bBigEndian ) SwapWord( 4, abyHeader+28 );
-
- i32 = nShapeType; /* shape type */
- ByteCopy( &i32, abyHeader+32, 4 );
- if( bBigEndian ) SwapWord( 4, abyHeader+32 );
-
- dValue = 0.0; /* set bounds */
- ByteCopy( &dValue, abyHeader+36, 8 );
- ByteCopy( &dValue, abyHeader+44, 8 );
- ByteCopy( &dValue, abyHeader+52, 8 );
- ByteCopy( &dValue, abyHeader+60, 8 );
-
-/* -------------------------------------------------------------------- */
-/* Write .shp file header. */
-/* -------------------------------------------------------------------- */
- if( fwrite( abyHeader, 100, 1, fpSHP ) != 1 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- "Failed to write .shp header." );
-#endif
- return NULL;
- }
-
-/* -------------------------------------------------------------------- */
-/* Prepare, and write .shx file header. */
-/* -------------------------------------------------------------------- */
- i32 = 50; /* file size */
- ByteCopy( &i32, abyHeader+24, 4 );
- if( !bBigEndian ) SwapWord( 4, abyHeader+24 );
-
- if( fwrite( abyHeader, 100, 1, fpSHX ) != 1 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- "Failed to write .shx header." );
-#endif
- return NULL;
- }
-
-/* -------------------------------------------------------------------- */
-/* Close the files, and then open them as regular existing files. */
-/* -------------------------------------------------------------------- */
- fclose( fpSHP );
- fclose( fpSHX );
-
- return( ossim_SHPOpen( pszLayer, "r+b" ) );
-}
-
-/************************************************************************/
-/* _SHPSetBounds() */
-/* */
-/* Compute a bounds rectangle for a shape, and set it into the */
-/* indicated location in the record. */
-/************************************************************************/
-
-static void _SHPSetBounds( uchar * pabyRec, ossim_SHPObject * psShape )
-
-{
- ByteCopy( &(psShape->dfXMin), pabyRec + 0, 8 );
- ByteCopy( &(psShape->dfYMin), pabyRec + 8, 8 );
- ByteCopy( &(psShape->dfXMax), pabyRec + 16, 8 );
- ByteCopy( &(psShape->dfYMax), pabyRec + 24, 8 );
-
- if( bBigEndian )
- {
- SwapWord( 8, pabyRec + 0 );
- SwapWord( 8, pabyRec + 8 );
- SwapWord( 8, pabyRec + 16 );
- SwapWord( 8, pabyRec + 24 );
- }
-}
-
-/************************************************************************/
-/* ossim_SHPComputeExtents() */
-/* */
-/* Recompute the extents of a shape. Automatically done by */
-/* ossim_SHPCreateObject(). */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPComputeExtents( ossim_SHPObject * psObject )
-
-{
- int i;
-
-/* -------------------------------------------------------------------- */
-/* Build extents for this object. */
-/* -------------------------------------------------------------------- */
- if( psObject->nVertices > 0 )
- {
- psObject->dfXMin = psObject->dfXMax = psObject->padfX[0];
- psObject->dfYMin = psObject->dfYMax = psObject->padfY[0];
- psObject->dfZMin = psObject->dfZMax = psObject->padfZ[0];
- psObject->dfMMin = psObject->dfMMax = psObject->padfM[0];
- }
-
- for( i = 0; i < psObject->nVertices; i++ )
- {
- psObject->dfXMin = MIN(psObject->dfXMin, psObject->padfX[i]);
- psObject->dfYMin = MIN(psObject->dfYMin, psObject->padfY[i]);
- psObject->dfZMin = MIN(psObject->dfZMin, psObject->padfZ[i]);
- psObject->dfMMin = MIN(psObject->dfMMin, psObject->padfM[i]);
-
- psObject->dfXMax = MAX(psObject->dfXMax, psObject->padfX[i]);
- psObject->dfYMax = MAX(psObject->dfYMax, psObject->padfY[i]);
- psObject->dfZMax = MAX(psObject->dfZMax, psObject->padfZ[i]);
- psObject->dfMMax = MAX(psObject->dfMMax, psObject->padfM[i]);
- }
-}
-
-/************************************************************************/
-/* ossim_SHPCreateObject() */
-/* */
-/* Create a shape object. It should be freed with */
-/* ossim_SHPDestroyObject(). */
-/************************************************************************/
-
-ossim_SHPObject ossim_SHPAPI_CALL1(*)
- ossim_SHPCreateObject( int nSHPType, int nShapeId, int nParts,
- const int * panPartStart, const int * panPartType,
- int nVertices, const double *padfX, const double *padfY,
- const double * padfZ, const double * padfM )
-
-{
- ossim_SHPObject *psObject;
- int i, bHasM, bHasZ;
-
- psObject = (ossim_SHPObject *) calloc(1,sizeof(ossim_SHPObject));
- psObject->nSHPType = nSHPType;
- psObject->nShapeId = nShapeId;
-
-/* -------------------------------------------------------------------- */
-/* Establish whether this shape type has M, and Z values. */
-/* -------------------------------------------------------------------- */
- if( nSHPType == ossim_SHPT_ARCM
- || nSHPType == ossim_SHPT_POINTM
- || nSHPType == ossim_SHPT_POLYGONM
- || nSHPType == ossim_SHPT_MULTIPOINTM )
- {
- bHasM = TRUE;
- bHasZ = FALSE;
- }
- else if( nSHPType == ossim_SHPT_ARCZ
- || nSHPType == ossim_SHPT_POINTZ
- || nSHPType == ossim_SHPT_POLYGONZ
- || nSHPType == ossim_SHPT_MULTIPOINTZ
- || nSHPType == ossim_SHPT_MULTIPATCH )
- {
- bHasM = TRUE;
- bHasZ = TRUE;
- }
- else
- {
- bHasM = FALSE;
- bHasZ = FALSE;
- }
-
-/* -------------------------------------------------------------------- */
-/* Capture parts. Note that part type is optional, and */
-/* defaults to ring. */
-/* -------------------------------------------------------------------- */
- if( nSHPType == ossim_SHPT_ARC || nSHPType == ossim_SHPT_POLYGON
- || nSHPType == ossim_SHPT_ARCM || nSHPType == ossim_SHPT_POLYGONM
- || nSHPType == ossim_SHPT_ARCZ || nSHPType == ossim_SHPT_POLYGONZ
- || nSHPType == ossim_SHPT_MULTIPATCH )
- {
- psObject->nParts = MAX(1,nParts);
-
- psObject->panPartStart = (int *)
- malloc(sizeof(int) * psObject->nParts);
- psObject->panPartType = (int *)
- malloc(sizeof(int) * psObject->nParts);
-
- psObject->panPartStart[0] = 0;
- psObject->panPartType[0] = ossim_SHPP_RING;
-
- for( i = 0; i < nParts; i++ )
- {
- psObject->panPartStart[i] = panPartStart[i];
-
- if( panPartType != NULL )
- psObject->panPartType[i] = panPartType[i];
- else
- psObject->panPartType[i] = ossim_SHPP_RING;
- }
-
- if( psObject->panPartStart[0] != 0 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_AppDefined,
- "panPartStart[0] != 0, patching internally. Please fix your code!\n" );
-#else
- fprintf( stderr, "panPartStart[0] != 0, patching internally. Please fix your code!\n" );
-#endif
- psObject->panPartStart[0] = 0;
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Capture vertices. Note that Z and M are optional, but X and */
-/* Y are not. */
-/* -------------------------------------------------------------------- */
- if( nVertices > 0 )
- {
- psObject->padfX = (double *) calloc(sizeof(double),nVertices);
- psObject->padfY = (double *) calloc(sizeof(double),nVertices);
- psObject->padfZ = (double *) calloc(sizeof(double),nVertices);
- psObject->padfM = (double *) calloc(sizeof(double),nVertices);
-
- assert( padfX != NULL );
- assert( padfY != NULL );
-
- for( i = 0; i < nVertices; i++ )
- {
- psObject->padfX[i] = padfX[i];
- psObject->padfY[i] = padfY[i];
- if( padfZ != NULL && bHasZ )
- psObject->padfZ[i] = padfZ[i];
- if( padfM != NULL && bHasM )
- psObject->padfM[i] = padfM[i];
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Compute the extents. */
-/* -------------------------------------------------------------------- */
- psObject->nVertices = nVertices;
- ossim_SHPComputeExtents( psObject );
-
- return( psObject );
-}
-
-/************************************************************************/
-/* ossim_SHPCreateSimpleObject() */
-/* */
-/* Create a simple (common) shape object. Destroy with */
-/* ossim_SHPDestroyObject(). */
-/************************************************************************/
-
-ossim_SHPObject ossim_SHPAPI_CALL1(*)
- ossim_SHPCreateSimpleObject( int nSHPType, int nVertices,
- const double * padfX, const double * padfY,
- const double * padfZ )
-
-{
- return( ossim_SHPCreateObject( nSHPType, -1, 0, NULL, NULL,
- nVertices, padfX, padfY, padfZ, NULL ) );
-}
-
-/************************************************************************/
-/* ossim_SHPWriteObject() */
-/* */
-/* Write out the vertices of a new structure. Note that it is */
-/* only possible to write vertices at the end of the file. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_SHPWriteObject(ossim_SHPHandle psSHP, int nShapeId, ossim_SHPObject * psObject )
-
-{
- int nRecordOffset, i, nRecordSize=0;
- uchar *pabyRec;
- int32 i32;
-
- psSHP->bUpdated = TRUE;
-
-/* -------------------------------------------------------------------- */
-/* Ensure that shape object matches the type of the file it is */
-/* being written to. */
-/* -------------------------------------------------------------------- */
- assert( psObject->nSHPType == psSHP->nShapeType
- || psObject->nSHPType == ossim_SHPT_NULL );
-
-/* -------------------------------------------------------------------- */
-/* Ensure that -1 is used for appends. Either blow an */
-/* assertion, or if they are disabled, set the shapeid to -1 */
-/* for appends. */
-/* -------------------------------------------------------------------- */
- assert( nShapeId == -1
- || (nShapeId >= 0 && nShapeId < psSHP->nRecords) );
-
- if( nShapeId != -1 && nShapeId >= psSHP->nRecords )
- nShapeId = -1;
-
-/* -------------------------------------------------------------------- */
-/* Add the new entity to the in memory index. */
-/* -------------------------------------------------------------------- */
- if( nShapeId == -1 && psSHP->nRecords+1 > psSHP->nMaxRecords )
- {
- psSHP->nMaxRecords =(int) ( psSHP->nMaxRecords * 1.3 + 100);
-
- psSHP->panRecOffset = (int *)
- SfRealloc(psSHP->panRecOffset,sizeof(int) * psSHP->nMaxRecords );
- psSHP->panRecSize = (int *)
- SfRealloc(psSHP->panRecSize,sizeof(int) * psSHP->nMaxRecords );
- }
-
-/* -------------------------------------------------------------------- */
-/* Initialize record. */
-/* -------------------------------------------------------------------- */
- pabyRec = (uchar *) malloc(psObject->nVertices * 4 * sizeof(double)
- + psObject->nParts * 8 + 128);
-
-/* -------------------------------------------------------------------- */
-/* Extract vertices for a Polygon or Arc. */
-/* -------------------------------------------------------------------- */
- if( psObject->nSHPType == ossim_SHPT_POLYGON
- || psObject->nSHPType == ossim_SHPT_POLYGONZ
- || psObject->nSHPType == ossim_SHPT_POLYGONM
- || psObject->nSHPType == ossim_SHPT_ARC
- || psObject->nSHPType == ossim_SHPT_ARCZ
- || psObject->nSHPType == ossim_SHPT_ARCM
- || psObject->nSHPType == ossim_SHPT_MULTIPATCH )
- {
- int32 nPoints, nParts;
- int i;
-
- nPoints = psObject->nVertices;
- nParts = psObject->nParts;
-
- _SHPSetBounds( pabyRec + 12, psObject );
-
- if( bBigEndian ) SwapWord( 4, &nPoints );
- if( bBigEndian ) SwapWord( 4, &nParts );
-
- ByteCopy( &nPoints, pabyRec + 40 + 8, 4 );
- ByteCopy( &nParts, pabyRec + 36 + 8, 4 );
-
- nRecordSize = 52;
-
- /*
- * Write part start positions.
- */
- ByteCopy( psObject->panPartStart, pabyRec + 44 + 8,
- 4 * psObject->nParts );
- for( i = 0; i < psObject->nParts; i++ )
- {
- if( bBigEndian ) SwapWord( 4, pabyRec + 44 + 8 + 4*i );
- nRecordSize += 4;
- }
-
- /*
- * Write multipatch part types if needed.
- */
- if( psObject->nSHPType == ossim_SHPT_MULTIPATCH )
- {
- memcpy( pabyRec + nRecordSize, psObject->panPartType,
- 4*psObject->nParts );
- for( i = 0; i < psObject->nParts; i++ )
- {
- if( bBigEndian ) SwapWord( 4, pabyRec + nRecordSize );
- nRecordSize += 4;
- }
- }
-
- /*
- * Write the (x,y) vertex values.
- */
- for( i = 0; i < psObject->nVertices; i++ )
- {
- ByteCopy( psObject->padfX + i, pabyRec + nRecordSize, 8 );
- ByteCopy( psObject->padfY + i, pabyRec + nRecordSize + 8, 8 );
-
- if( bBigEndian )
- SwapWord( 8, pabyRec + nRecordSize );
-
- if( bBigEndian )
- SwapWord( 8, pabyRec + nRecordSize + 8 );
-
- nRecordSize += 2 * 8;
- }
-
- /*
- * Write the Z coordinates (if any).
- */
- if( psObject->nSHPType == ossim_SHPT_POLYGONZ
- || psObject->nSHPType == ossim_SHPT_ARCZ
- || psObject->nSHPType == ossim_SHPT_MULTIPATCH )
- {
- ByteCopy( &(psObject->dfZMin), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- ByteCopy( &(psObject->dfZMax), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- for( i = 0; i < psObject->nVertices; i++ )
- {
- ByteCopy( psObject->padfZ + i, pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
- }
- }
-
- /*
- * Write the M values, if any.
- */
- if( psObject->nSHPType == ossim_SHPT_POLYGONM
- || psObject->nSHPType == ossim_SHPT_ARCM
-#ifndef DISABLE_MULTIPATCH_MEASURE
- || psObject->nSHPType == ossim_SHPT_MULTIPATCH
-#endif
- || psObject->nSHPType == ossim_SHPT_POLYGONZ
- || psObject->nSHPType == ossim_SHPT_ARCZ )
- {
- ByteCopy( &(psObject->dfMMin), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- ByteCopy( &(psObject->dfMMax), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- for( i = 0; i < psObject->nVertices; i++ )
- {
- ByteCopy( psObject->padfM + i, pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
- }
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Extract vertices for a MultiPoint. */
-/* -------------------------------------------------------------------- */
- else if( psObject->nSHPType == ossim_SHPT_MULTIPOINT
- || psObject->nSHPType == ossim_SHPT_MULTIPOINTZ
- || psObject->nSHPType == ossim_SHPT_MULTIPOINTM )
- {
- int32 nPoints;
- int i;
-
- nPoints = psObject->nVertices;
-
- _SHPSetBounds( pabyRec + 12, psObject );
-
- if( bBigEndian ) SwapWord( 4, &nPoints );
- ByteCopy( &nPoints, pabyRec + 44, 4 );
-
- for( i = 0; i < psObject->nVertices; i++ )
- {
- ByteCopy( psObject->padfX + i, pabyRec + 48 + i*16, 8 );
- ByteCopy( psObject->padfY + i, pabyRec + 48 + i*16 + 8, 8 );
-
- if( bBigEndian ) SwapWord( 8, pabyRec + 48 + i*16 );
- if( bBigEndian ) SwapWord( 8, pabyRec + 48 + i*16 + 8 );
- }
-
- nRecordSize = 48 + 16 * psObject->nVertices;
-
- if( psObject->nSHPType == ossim_SHPT_MULTIPOINTZ )
- {
- ByteCopy( &(psObject->dfZMin), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- ByteCopy( &(psObject->dfZMax), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- for( i = 0; i < psObject->nVertices; i++ )
- {
- ByteCopy( psObject->padfZ + i, pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
- }
- }
-
- if( psObject->nSHPType == ossim_SHPT_MULTIPOINTZ
- || psObject->nSHPType == ossim_SHPT_MULTIPOINTM )
- {
- ByteCopy( &(psObject->dfMMin), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- ByteCopy( &(psObject->dfMMax), pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
-
- for( i = 0; i < psObject->nVertices; i++ )
- {
- ByteCopy( psObject->padfM + i, pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
- }
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Write point. */
-/* -------------------------------------------------------------------- */
- else if( psObject->nSHPType == ossim_SHPT_POINT
- || psObject->nSHPType == ossim_SHPT_POINTZ
- || psObject->nSHPType == ossim_SHPT_POINTM )
- {
- ByteCopy( psObject->padfX, pabyRec + 12, 8 );
- ByteCopy( psObject->padfY, pabyRec + 20, 8 );
-
- if( bBigEndian ) SwapWord( 8, pabyRec + 12 );
- if( bBigEndian ) SwapWord( 8, pabyRec + 20 );
-
- nRecordSize = 28;
-
- if( psObject->nSHPType == ossim_SHPT_POINTZ )
- {
- ByteCopy( psObject->padfZ, pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
- }
-
- if( psObject->nSHPType == ossim_SHPT_POINTZ
- || psObject->nSHPType == ossim_SHPT_POINTM )
- {
- ByteCopy( psObject->padfM, pabyRec + nRecordSize, 8 );
- if( bBigEndian ) SwapWord( 8, pabyRec + nRecordSize );
- nRecordSize += 8;
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Not much to do for null geometries. */
-/* -------------------------------------------------------------------- */
- else if( psObject->nSHPType == ossim_SHPT_NULL )
- {
- nRecordSize = 12;
- }
-
- else
- {
- /* unknown type */
- assert( FALSE );
- }
-
-/* -------------------------------------------------------------------- */
-/* Establish where we are going to put this record. If we are */
-/* rewriting and existing record, and it will fit, then put it */
-/* back where the original came from. Otherwise write at the end. */
-/* -------------------------------------------------------------------- */
- if( nShapeId == -1 || psSHP->panRecSize[nShapeId] < nRecordSize-8 )
- {
- if( nShapeId == -1 )
- nShapeId = psSHP->nRecords++;
-
- psSHP->panRecOffset[nShapeId] = nRecordOffset = psSHP->nFileSize;
- psSHP->panRecSize[nShapeId] = nRecordSize-8;
- psSHP->nFileSize += nRecordSize;
- }
- else
- {
- nRecordOffset = psSHP->panRecOffset[nShapeId];
- }
-
-/* -------------------------------------------------------------------- */
-/* Set the shape type, record number, and record size. */
-/* -------------------------------------------------------------------- */
- i32 = nShapeId+1; /* record # */
- if( !bBigEndian ) SwapWord( 4, &i32 );
- ByteCopy( &i32, pabyRec, 4 );
-
- i32 = (nRecordSize-8)/2; /* record size */
- if( !bBigEndian ) SwapWord( 4, &i32 );
- ByteCopy( &i32, pabyRec + 4, 4 );
-
- i32 = psObject->nSHPType; /* shape type */
- if( bBigEndian ) SwapWord( 4, &i32 );
- ByteCopy( &i32, pabyRec + 8, 4 );
-
-/* -------------------------------------------------------------------- */
-/* Write out record. */
-/* -------------------------------------------------------------------- */
- if( fseek( psSHP->fpSHP, nRecordOffset, 0 ) != 0
- || fwrite( pabyRec, nRecordSize, 1, psSHP->fpSHP ) < 1 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_FileIO,
- "Error in fseek() or fwrite() writing object to .shp file." );
-#endif
- free( pabyRec );
- return -1;
- }
-
- free( pabyRec );
-
-/* -------------------------------------------------------------------- */
-/* Expand file wide bounds based on this shape. */
-/* -------------------------------------------------------------------- */
- if( psSHP->adBoundsMin[0] == 0.0
- && psSHP->adBoundsMax[0] == 0.0
- && psSHP->adBoundsMin[1] == 0.0
- && psSHP->adBoundsMax[1] == 0.0 )
- {
- if( psObject->nSHPType == ossim_SHPT_NULL || psObject->nVertices == 0 )
- {
- psSHP->adBoundsMin[0] = psSHP->adBoundsMax[0] = 0.0;
- psSHP->adBoundsMin[1] = psSHP->adBoundsMax[1] = 0.0;
- psSHP->adBoundsMin[2] = psSHP->adBoundsMax[2] = 0.0;
- psSHP->adBoundsMin[3] = psSHP->adBoundsMax[3] = 0.0;
- }
- else
- {
- psSHP->adBoundsMin[0] = psSHP->adBoundsMax[0] = psObject->padfX[0];
- psSHP->adBoundsMin[1] = psSHP->adBoundsMax[1] = psObject->padfY[0];
- psSHP->adBoundsMin[2] = psSHP->adBoundsMax[2] = psObject->padfZ[0];
- psSHP->adBoundsMin[3] = psSHP->adBoundsMax[3] = psObject->padfM[0];
- }
- }
-
- for( i = 0; i < psObject->nVertices; i++ )
- {
- psSHP->adBoundsMin[0] = MIN(psSHP->adBoundsMin[0],psObject->padfX[i]);
- psSHP->adBoundsMin[1] = MIN(psSHP->adBoundsMin[1],psObject->padfY[i]);
- psSHP->adBoundsMin[2] = MIN(psSHP->adBoundsMin[2],psObject->padfZ[i]);
- psSHP->adBoundsMin[3] = MIN(psSHP->adBoundsMin[3],psObject->padfM[i]);
- psSHP->adBoundsMax[0] = MAX(psSHP->adBoundsMax[0],psObject->padfX[i]);
- psSHP->adBoundsMax[1] = MAX(psSHP->adBoundsMax[1],psObject->padfY[i]);
- psSHP->adBoundsMax[2] = MAX(psSHP->adBoundsMax[2],psObject->padfZ[i]);
- psSHP->adBoundsMax[3] = MAX(psSHP->adBoundsMax[3],psObject->padfM[i]);
- }
-
- return( nShapeId );
-}
-
-/************************************************************************/
-/* ossim_SHPReadObject() */
-/* */
-/* Read the vertices, parts, and other non-attribute information */
-/* for one shape. */
-/************************************************************************/
-
-ossim_SHPObject ossim_SHPAPI_CALL1(*)
- ossim_SHPReadObject( ossim_SHPHandle psSHP, int hEntity )
-
-{
- ossim_SHPObject *psShape;
-
-/* -------------------------------------------------------------------- */
-/* Validate the record/entity number. */
-/* -------------------------------------------------------------------- */
- if( hEntity < 0 || hEntity >= psSHP->nRecords )
- return( NULL );
-
-/* -------------------------------------------------------------------- */
-/* Ensure our record buffer is large enough. */
-/* -------------------------------------------------------------------- */
- if( psSHP->panRecSize[hEntity]+8 > psSHP->nBufSize )
- {
- psSHP->nBufSize = psSHP->panRecSize[hEntity]+8;
- psSHP->pabyRec = (uchar *) SfRealloc(psSHP->pabyRec,psSHP->nBufSize);
- }
-
-/* -------------------------------------------------------------------- */
-/* Read the record. */
-/* -------------------------------------------------------------------- */
- if( fseek( psSHP->fpSHP, psSHP->panRecOffset[hEntity], 0 ) != 0
- || fread( psSHP->pabyRec, psSHP->panRecSize[hEntity]+8, 1,
- psSHP->fpSHP ) != 1 )
- {
-#ifdef USE_CPL
- CPLError( CE_Failure, CPLE_FileIO,
- "Error in fseek() or fread() reading object from .shp file." );
-#endif
- return NULL;
- }
-
-/* -------------------------------------------------------------------- */
-/* Allocate and minimally initialize the object. */
-/* -------------------------------------------------------------------- */
- psShape = (ossim_SHPObject *) calloc(1,sizeof(ossim_SHPObject));
- psShape->nShapeId = hEntity;
-
- memcpy( &psShape->nSHPType, psSHP->pabyRec + 8, 4 );
- if( bBigEndian ) SwapWord( 4, &(psShape->nSHPType) );
-
-/* ==================================================================== */
-/* Extract vertices for a Polygon or Arc. */
-/* ==================================================================== */
- if( psShape->nSHPType == ossim_SHPT_POLYGON || psShape->nSHPType == ossim_SHPT_ARC
- || psShape->nSHPType == ossim_SHPT_POLYGONZ
- || psShape->nSHPType == ossim_SHPT_POLYGONM
- || psShape->nSHPType == ossim_SHPT_ARCZ
- || psShape->nSHPType == ossim_SHPT_ARCM
- || psShape->nSHPType == ossim_SHPT_MULTIPATCH )
- {
- int32 nPoints, nParts;
- int i, nOffset;
-
-/* -------------------------------------------------------------------- */
-/* Get the X/Y bounds. */
-/* -------------------------------------------------------------------- */
- memcpy( &(psShape->dfXMin), psSHP->pabyRec + 8 + 4, 8 );
- memcpy( &(psShape->dfYMin), psSHP->pabyRec + 8 + 12, 8 );
- memcpy( &(psShape->dfXMax), psSHP->pabyRec + 8 + 20, 8 );
- memcpy( &(psShape->dfYMax), psSHP->pabyRec + 8 + 28, 8 );
-
- if( bBigEndian ) SwapWord( 8, &(psShape->dfXMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfYMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfXMax) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfYMax) );
-
-/* -------------------------------------------------------------------- */
-/* Extract part/point count, and build vertex and part arrays */
-/* to proper size. */
-/* -------------------------------------------------------------------- */
- memcpy( &nPoints, psSHP->pabyRec + 40 + 8, 4 );
- memcpy( &nParts, psSHP->pabyRec + 36 + 8, 4 );
-
- if( bBigEndian ) SwapWord( 4, &nPoints );
- if( bBigEndian ) SwapWord( 4, &nParts );
-
- psShape->nVertices = nPoints;
- psShape->padfX = (double *) calloc(nPoints,sizeof(double));
- psShape->padfY = (double *) calloc(nPoints,sizeof(double));
- psShape->padfZ = (double *) calloc(nPoints,sizeof(double));
- psShape->padfM = (double *) calloc(nPoints,sizeof(double));
-
- psShape->nParts = nParts;
- psShape->panPartStart = (int *) calloc(nParts,sizeof(int));
- psShape->panPartType = (int *) calloc(nParts,sizeof(int));
-
- for( i = 0; i < nParts; i++ )
- psShape->panPartType[i] = ossim_SHPP_RING;
-
-/* -------------------------------------------------------------------- */
-/* Copy out the part array from the record. */
-/* -------------------------------------------------------------------- */
- memcpy( psShape->panPartStart, psSHP->pabyRec + 44 + 8, 4 * nParts );
- for( i = 0; i < nParts; i++ )
- {
- if( bBigEndian ) SwapWord( 4, psShape->panPartStart+i );
- }
-
- nOffset = 44 + 8 + 4*nParts;
-
-/* -------------------------------------------------------------------- */
-/* If this is a multipatch, we will also have parts types. */
-/* -------------------------------------------------------------------- */
- if( psShape->nSHPType == ossim_SHPT_MULTIPATCH )
- {
- memcpy( psShape->panPartType, psSHP->pabyRec + nOffset, 4*nParts );
- for( i = 0; i < nParts; i++ )
- {
- if( bBigEndian ) SwapWord( 4, psShape->panPartType+i );
- }
-
- nOffset += 4*nParts;
- }
-
-/* -------------------------------------------------------------------- */
-/* Copy out the vertices from the record. */
-/* -------------------------------------------------------------------- */
- for( i = 0; i < nPoints; i++ )
- {
- memcpy(psShape->padfX + i,
- psSHP->pabyRec + nOffset + i * 16,
- 8 );
-
- memcpy(psShape->padfY + i,
- psSHP->pabyRec + nOffset + i * 16 + 8,
- 8 );
-
- if( bBigEndian ) SwapWord( 8, psShape->padfX + i );
- if( bBigEndian ) SwapWord( 8, psShape->padfY + i );
- }
-
- nOffset += 16*nPoints;
-
-/* -------------------------------------------------------------------- */
-/* If we have a Z coordinate, collect that now. */
-/* -------------------------------------------------------------------- */
- if( psShape->nSHPType == ossim_SHPT_POLYGONZ
- || psShape->nSHPType == ossim_SHPT_ARCZ
- || psShape->nSHPType == ossim_SHPT_MULTIPATCH )
- {
- memcpy( &(psShape->dfZMin), psSHP->pabyRec + nOffset, 8 );
- memcpy( &(psShape->dfZMax), psSHP->pabyRec + nOffset + 8, 8 );
-
- if( bBigEndian ) SwapWord( 8, &(psShape->dfZMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfZMax) );
-
- for( i = 0; i < nPoints; i++ )
- {
- memcpy( psShape->padfZ + i,
- psSHP->pabyRec + nOffset + 16 + i*8, 8 );
- if( bBigEndian ) SwapWord( 8, psShape->padfZ + i );
- }
-
- nOffset += 16 + 8*nPoints;
- }
-
-/* -------------------------------------------------------------------- */
-/* If we have a M measure value, then read it now. We assume */
-/* that the measure can be present for any shape if the size is */
-/* big enough, but really it will only occur for the Z shapes */
-/* (options), and the M shapes. */
-/* -------------------------------------------------------------------- */
- if( psSHP->panRecSize[hEntity]+8 >= nOffset + 16 + 8*nPoints )
- {
- memcpy( &(psShape->dfMMin), psSHP->pabyRec + nOffset, 8 );
- memcpy( &(psShape->dfMMax), psSHP->pabyRec + nOffset + 8, 8 );
-
- if( bBigEndian ) SwapWord( 8, &(psShape->dfMMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfMMax) );
-
- for( i = 0; i < nPoints; i++ )
- {
- memcpy( psShape->padfM + i,
- psSHP->pabyRec + nOffset + 16 + i*8, 8 );
- if( bBigEndian ) SwapWord( 8, psShape->padfM + i );
- }
- }
-
- }
-
-/* ==================================================================== */
-/* Extract vertices for a MultiPoint. */
-/* ==================================================================== */
- else if( psShape->nSHPType == ossim_SHPT_MULTIPOINT
- || psShape->nSHPType == ossim_SHPT_MULTIPOINTM
- || psShape->nSHPType == ossim_SHPT_MULTIPOINTZ )
- {
- int32 nPoints;
- int i, nOffset;
-
- memcpy( &nPoints, psSHP->pabyRec + 44, 4 );
- if( bBigEndian ) SwapWord( 4, &nPoints );
-
- psShape->nVertices = nPoints;
- psShape->padfX = (double *) calloc(nPoints,sizeof(double));
- psShape->padfY = (double *) calloc(nPoints,sizeof(double));
- psShape->padfZ = (double *) calloc(nPoints,sizeof(double));
- psShape->padfM = (double *) calloc(nPoints,sizeof(double));
-
- for( i = 0; i < nPoints; i++ )
- {
- memcpy(psShape->padfX+i, psSHP->pabyRec + 48 + 16 * i, 8 );
- memcpy(psShape->padfY+i, psSHP->pabyRec + 48 + 16 * i + 8, 8 );
-
- if( bBigEndian ) SwapWord( 8, psShape->padfX + i );
- if( bBigEndian ) SwapWord( 8, psShape->padfY + i );
- }
-
- nOffset = 48 + 16*nPoints;
-
-/* -------------------------------------------------------------------- */
-/* Get the X/Y bounds. */
-/* -------------------------------------------------------------------- */
- memcpy( &(psShape->dfXMin), psSHP->pabyRec + 8 + 4, 8 );
- memcpy( &(psShape->dfYMin), psSHP->pabyRec + 8 + 12, 8 );
- memcpy( &(psShape->dfXMax), psSHP->pabyRec + 8 + 20, 8 );
- memcpy( &(psShape->dfYMax), psSHP->pabyRec + 8 + 28, 8 );
-
- if( bBigEndian ) SwapWord( 8, &(psShape->dfXMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfYMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfXMax) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfYMax) );
-
-/* -------------------------------------------------------------------- */
-/* If we have a Z coordinate, collect that now. */
-/* -------------------------------------------------------------------- */
- if( psShape->nSHPType == ossim_SHPT_MULTIPOINTZ )
- {
- memcpy( &(psShape->dfZMin), psSHP->pabyRec + nOffset, 8 );
- memcpy( &(psShape->dfZMax), psSHP->pabyRec + nOffset + 8, 8 );
-
- if( bBigEndian ) SwapWord( 8, &(psShape->dfZMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfZMax) );
-
- for( i = 0; i < nPoints; i++ )
- {
- memcpy( psShape->padfZ + i,
- psSHP->pabyRec + nOffset + 16 + i*8, 8 );
- if( bBigEndian ) SwapWord( 8, psShape->padfZ + i );
- }
-
- nOffset += 16 + 8*nPoints;
- }
-
-/* -------------------------------------------------------------------- */
-/* If we have a M measure value, then read it now. We assume */
-/* that the measure can be present for any shape if the size is */
-/* big enough, but really it will only occur for the Z shapes */
-/* (options), and the M shapes. */
-/* -------------------------------------------------------------------- */
- if( psSHP->panRecSize[hEntity]+8 >= nOffset + 16 + 8*nPoints )
- {
- memcpy( &(psShape->dfMMin), psSHP->pabyRec + nOffset, 8 );
- memcpy( &(psShape->dfMMax), psSHP->pabyRec + nOffset + 8, 8 );
-
- if( bBigEndian ) SwapWord( 8, &(psShape->dfMMin) );
- if( bBigEndian ) SwapWord( 8, &(psShape->dfMMax) );
-
- for( i = 0; i < nPoints; i++ )
- {
- memcpy( psShape->padfM + i,
- psSHP->pabyRec + nOffset + 16 + i*8, 8 );
- if( bBigEndian ) SwapWord( 8, psShape->padfM + i );
- }
- }
- }
-
-/* ==================================================================== */
-/* Extract vertices for a point. */
-/* ==================================================================== */
- else if( psShape->nSHPType == ossim_SHPT_POINT
- || psShape->nSHPType == ossim_SHPT_POINTM
- || psShape->nSHPType == ossim_SHPT_POINTZ )
- {
- int nOffset;
-
- psShape->nVertices = 1;
- psShape->padfX = (double *) calloc(1,sizeof(double));
- psShape->padfY = (double *) calloc(1,sizeof(double));
- psShape->padfZ = (double *) calloc(1,sizeof(double));
- psShape->padfM = (double *) calloc(1,sizeof(double));
-
- memcpy( psShape->padfX, psSHP->pabyRec + 12, 8 );
- memcpy( psShape->padfY, psSHP->pabyRec + 20, 8 );
-
- if( bBigEndian ) SwapWord( 8, psShape->padfX );
- if( bBigEndian ) SwapWord( 8, psShape->padfY );
-
- nOffset = 20 + 8;
-
-/* -------------------------------------------------------------------- */
-/* If we have a Z coordinate, collect that now. */
-/* -------------------------------------------------------------------- */
- if( psShape->nSHPType == ossim_SHPT_POINTZ )
- {
- memcpy( psShape->padfZ, psSHP->pabyRec + nOffset, 8 );
-
- if( bBigEndian ) SwapWord( 8, psShape->padfZ );
-
- nOffset += 8;
- }
-
-/* -------------------------------------------------------------------- */
-/* If we have a M measure value, then read it now. We assume */
-/* that the measure can be present for any shape if the size is */
-/* big enough, but really it will only occur for the Z shapes */
-/* (options), and the M shapes. */
-/* -------------------------------------------------------------------- */
- if( psSHP->panRecSize[hEntity]+8 >= nOffset + 8 )
- {
- memcpy( psShape->padfM, psSHP->pabyRec + nOffset, 8 );
-
- if( bBigEndian ) SwapWord( 8, psShape->padfM );
- }
-
-/* -------------------------------------------------------------------- */
-/* Since no extents are supplied in the record, we will apply */
-/* them from the single vertex. */
-/* -------------------------------------------------------------------- */
- psShape->dfXMin = psShape->dfXMax = psShape->padfX[0];
- psShape->dfYMin = psShape->dfYMax = psShape->padfY[0];
- psShape->dfZMin = psShape->dfZMax = psShape->padfZ[0];
- psShape->dfMMin = psShape->dfMMax = psShape->padfM[0];
- }
-
- return( psShape );
-}
-
-/************************************************************************/
-/* ossim_SHPTypeName() */
-/************************************************************************/
-
-const char ossim_SHPAPI_CALL1(*)
- ossim_SHPTypeName( int nSHPType )
-
-{
- switch( nSHPType )
- {
- case ossim_SHPT_NULL:
- return "NullShape";
-
- case ossim_SHPT_POINT:
- return "Point";
-
- case ossim_SHPT_ARC:
- return "Arc";
-
- case ossim_SHPT_POLYGON:
- return "Polygon";
-
- case ossim_SHPT_MULTIPOINT:
- return "MultiPoint";
-
- case ossim_SHPT_POINTZ:
- return "PointZ";
-
- case ossim_SHPT_ARCZ:
- return "ArcZ";
-
- case ossim_SHPT_POLYGONZ:
- return "PolygonZ";
-
- case ossim_SHPT_MULTIPOINTZ:
- return "MultiPointZ";
-
- case ossim_SHPT_POINTM:
- return "PointM";
-
- case ossim_SHPT_ARCM:
- return "ArcM";
-
- case ossim_SHPT_POLYGONM:
- return "PolygonM";
-
- case ossim_SHPT_MULTIPOINTM:
- return "MultiPointM";
-
- case ossim_SHPT_MULTIPATCH:
- return "MultiPatch";
-
- default:
- return "UnknownShapeType";
- }
-}
-
-/************************************************************************/
-/* ossim_SHPPartTypeName() */
-/************************************************************************/
-
-const char ossim_SHPAPI_CALL1(*)
- ossim_SHPPartTypeName( int nPartType )
-
-{
- switch( nPartType )
- {
- case ossim_SHPP_TRISTRIP:
- return "TriangleStrip";
-
- case ossim_SHPP_TRIFAN:
- return "TriangleFan";
-
- case ossim_SHPP_OUTERRING:
- return "OuterRing";
-
- case ossim_SHPP_INNERRING:
- return "InnerRing";
-
- case ossim_SHPP_FIRSTRING:
- return "FirstRing";
-
- case ossim_SHPP_RING:
- return "Ring";
-
- default:
- return "UnknownPartType";
- }
-}
-
-/************************************************************************/
-/* ossim_SHPDestroyObject() */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPDestroyObject( ossim_SHPObject * psShape )
-
-{
- if( psShape == NULL )
- return;
-
- if( psShape->padfX != NULL )
- free( psShape->padfX );
- if( psShape->padfY != NULL )
- free( psShape->padfY );
- if( psShape->padfZ != NULL )
- free( psShape->padfZ );
- if( psShape->padfM != NULL )
- free( psShape->padfM );
-
- if( psShape->panPartStart != NULL )
- free( psShape->panPartStart );
- if( psShape->panPartType != NULL )
- free( psShape->panPartType );
-
- free( psShape );
-}
-
-/************************************************************************/
-/* ossim_SHPRewindObject() */
-/* */
-/* Reset the winding of polygon objects to adhere to the */
-/* specification. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_SHPRewindObject( ossim_SHPHandle hSHP, ossim_SHPObject * psObject )
-
-{
- int iOpRing, bAltered = 0;
-
-/* -------------------------------------------------------------------- */
-/* Do nothing if this is not a polygon object. */
-/* -------------------------------------------------------------------- */
- if( psObject->nSHPType != ossim_SHPT_POLYGON
- && psObject->nSHPType != ossim_SHPT_POLYGONZ
- && psObject->nSHPType != ossim_SHPT_POLYGONM )
- return 0;
-
- if( psObject->nVertices == 0 || psObject->nParts == 0 )
- return 0;
-
-/* -------------------------------------------------------------------- */
-/* Process each of the rings. */
-/* -------------------------------------------------------------------- */
- for( iOpRing = 0; iOpRing < psObject->nParts; iOpRing++ )
- {
- int bInner, iVert, nVertCount, nVertStart, iCheckRing;
- double dfSum, dfTestX, dfTestY;
-
-/* -------------------------------------------------------------------- */
-/* Determine if this ring is an inner ring or an outer ring */
-/* relative to all the other rings. For now we assume the */
-/* first ring is outer and all others are inner, but eventually */
-/* we need to fix this to handle multiple island polygons and */
-/* unordered sets of rings. */
-/* -------------------------------------------------------------------- */
- dfTestX = psObject->padfX[psObject->panPartStart[iOpRing]];
- dfTestY = psObject->padfY[psObject->panPartStart[iOpRing]];
-
- bInner = FALSE;
- for( iCheckRing = 0; iCheckRing < psObject->nParts; iCheckRing++ )
- {
- int iEdge;
-
- if( iCheckRing == iOpRing )
- continue;
-
- nVertStart = psObject->panPartStart[iCheckRing];
-
- if( iCheckRing == psObject->nParts-1 )
- nVertCount = psObject->nVertices
- - psObject->panPartStart[iCheckRing];
- else
- nVertCount = psObject->panPartStart[iCheckRing+1]
- - psObject->panPartStart[iCheckRing];
-
- for( iEdge = 0; iEdge < nVertCount; iEdge++ )
- {
- int iNext;
-
- if( iEdge < nVertCount-1 )
- iNext = iEdge+1;
- else
- iNext = 0;
-
- if( (psObject->padfY[iEdge+nVertStart] < dfTestY
- && psObject->padfY[iNext+nVertStart] >= dfTestY)
- || (psObject->padfY[iNext+nVertStart] < dfTestY
- && psObject->padfY[iEdge+nVertStart] >= dfTestY) )
- {
- if( psObject->padfX[iEdge+nVertStart]
- + (dfTestY - psObject->padfY[iEdge+nVertStart])
- / (psObject->padfY[iNext+nVertStart]
- - psObject->padfY[iEdge+nVertStart])
- * (psObject->padfX[iNext+nVertStart]
- - psObject->padfX[iEdge+nVertStart]) < dfTestX )
- bInner = !bInner;
- }
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Determine the current order of this ring so we will know if */
-/* it has to be reversed. */
-/* -------------------------------------------------------------------- */
- nVertStart = psObject->panPartStart[iOpRing];
-
- if( iOpRing == psObject->nParts-1 )
- nVertCount = psObject->nVertices - psObject->panPartStart[iOpRing];
- else
- nVertCount = psObject->panPartStart[iOpRing+1]
- - psObject->panPartStart[iOpRing];
-
- dfSum = 0.0;
- for( iVert = nVertStart; iVert < nVertStart+nVertCount-1; iVert++ )
- {
- dfSum += psObject->padfX[iVert] * psObject->padfY[iVert+1]
- - psObject->padfY[iVert] * psObject->padfX[iVert+1];
- }
-
- dfSum += psObject->padfX[iVert] * psObject->padfY[nVertStart]
- - psObject->padfY[iVert] * psObject->padfX[nVertStart];
-
-/* -------------------------------------------------------------------- */
-/* Reverse if necessary. */
-/* -------------------------------------------------------------------- */
- if( (dfSum < 0.0 && bInner) || (dfSum > 0.0 && !bInner) )
- {
- int i;
-
- bAltered++;
- for( i = 0; i < nVertCount/2; i++ )
- {
- double dfSaved;
-
- /* Swap X */
- dfSaved = psObject->padfX[nVertStart+i];
- psObject->padfX[nVertStart+i] =
- psObject->padfX[nVertStart+nVertCount-i-1];
- psObject->padfX[nVertStart+nVertCount-i-1] = dfSaved;
-
- /* Swap Y */
- dfSaved = psObject->padfY[nVertStart+i];
- psObject->padfY[nVertStart+i] =
- psObject->padfY[nVertStart+nVertCount-i-1];
- psObject->padfY[nVertStart+nVertCount-i-1] = dfSaved;
-
- /* Swap Z */
- if( psObject->padfZ )
- {
- dfSaved = psObject->padfZ[nVertStart+i];
- psObject->padfZ[nVertStart+i] =
- psObject->padfZ[nVertStart+nVertCount-i-1];
- psObject->padfZ[nVertStart+nVertCount-i-1] = dfSaved;
- }
-
- /* Swap M */
- if( psObject->padfM )
- {
- dfSaved = psObject->padfM[nVertStart+i];
- psObject->padfM[nVertStart+i] =
- psObject->padfM[nVertStart+nVertCount-i-1];
- psObject->padfM[nVertStart+nVertCount-i-1] = dfSaved;
- }
- }
- }
- }
-
- return bAltered;
-}
diff --git a/Utilities/otbossim/src/ossim/vec/shptree.c b/Utilities/otbossim/src/ossim/vec/shptree.c
deleted file mode 100644
index 00f49cf92b..0000000000
--- a/Utilities/otbossim/src/ossim/vec/shptree.c
+++ /dev/null
@@ -1,986 +0,0 @@
-/******************************************************************************
- * $Id: shptree.c 9094 2006-06-13 19:12:40Z dburken $
- *
- * Project: Shapelib
- * Purpose: Implementation of quadtree building and searching functions.
- * 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 LICENSE.LGPL). 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$
- * Revision 1.6 2006/06/13 19:11:33 dburken
- * Global include path change.
- *
- * Revision 1.5 2005/10/08 12:57:26 gpotts
- * Added adiitional patches
- *
- * Revision 1.10 2005/01/03 22:30:13 fwarmerdam
- * added support for saved quadtrees
- *
- * Revision 1.9 2003/01/28 15:53:41 warmerda
- * Avoid build warnings.
- *
- * Revision 1.8 2002/05/07 13:07:45 warmerda
- * use qsort() - patch from Bernhard Herzog
- *
- * Revision 1.7 2002/01/15 14:36:07 warmerda
- * updated email address
- *
- * Revision 1.6 2001/05/23 13:36:52 warmerda
- * added use of ossim_SHPAPI_CALL
- *
- * Revision 1.5 1999/11/05 14:12:05 warmerda
- * updated license terms
- *
- * Revision 1.4 1999/06/02 18:24:21 warmerda
- * added trimming code
- *
- * Revision 1.3 1999/06/02 17:56:12 warmerda
- * added quad'' subnode support for trees
- *
- * Revision 1.2 1999/05/18 19:11:11 warmerda
- * Added example searching capability
- *
- * Revision 1.1 1999/05/18 17:49:20 warmerda
- * New
- *
- */
-
-#include <ossim/vec/shapefil.h>
-
-#include <math.h>
-#include <assert.h>
-#include <stdlib.h>
-#include <string.h>
-
-ossim_SHP_CVSID("$Id: shptree.c 9094 2006-06-13 19:12:40Z dburken $")
-
-#ifndef TRUE
-# define TRUE 1
-# define FALSE 0
-#endif
-
-static int bBigEndian = 0;
-
-
-/* -------------------------------------------------------------------- */
-/* If the following is 0.5, nodes will be split in half. If it */
-/* is 0.6 then each subnode will contain 60% of the parent */
-/* node, with 20% representing overlap. This can be help to */
-/* prevent small objects on a boundary from shifting too high */
-/* up the tree. */
-/* -------------------------------------------------------------------- */
-
-#define ossim_SHP_SPLIT_RATIO 0.55
-
-/************************************************************************/
-/* SfRealloc() */
-/* */
-/* A realloc cover function that will access a NULL pointer as */
-/* a valid input. */
-/************************************************************************/
-
-static void * SfRealloc( void * pMem, int nNewSize )
-
-{
- if( pMem == NULL )
- return( (void *) malloc(nNewSize) );
- else
- return( (void *) realloc(pMem,nNewSize) );
-}
-
-/************************************************************************/
-/* ossim_SHPTreeNodeInit() */
-/* */
-/* Initialize a tree node. */
-/************************************************************************/
-
-static ossim_SHPTreeNode *ossim_SHPTreeNodeCreate( double * padfBoundsMin,
- double * padfBoundsMax )
-
-{
- ossim_SHPTreeNode *psTreeNode;
-
- psTreeNode = (ossim_SHPTreeNode *) malloc(sizeof(ossim_SHPTreeNode));
-
- psTreeNode->nShapeCount = 0;
- psTreeNode->panShapeIds = NULL;
- psTreeNode->papsShapeObj = NULL;
-
- psTreeNode->nSubNodes = 0;
-
- if( padfBoundsMin != NULL )
- memcpy( psTreeNode->adfBoundsMin, padfBoundsMin, sizeof(double) * 4 );
-
- if( padfBoundsMax != NULL )
- memcpy( psTreeNode->adfBoundsMax, padfBoundsMax, sizeof(double) * 4 );
-
- return psTreeNode;
-}
-
-
-/************************************************************************/
-/* ossim_SHPCreateTree() */
-/************************************************************************/
-
-ossim_SHPTree ossim_SHPAPI_CALL1(*)
- ossim_SHPCreateTree( ossim_SHPHandle hSHP, int nDimension, int nMaxDepth,
- double *padfBoundsMin, double *padfBoundsMax )
-
-{
- ossim_SHPTree *psTree;
-
- if( padfBoundsMin == NULL && hSHP == NULL )
- return NULL;
-
-/* -------------------------------------------------------------------- */
-/* Allocate the tree object */
-/* -------------------------------------------------------------------- */
- psTree = (ossim_SHPTree *) malloc(sizeof(ossim_SHPTree));
-
- psTree->hSHP = hSHP;
- psTree->nMaxDepth = nMaxDepth;
- psTree->nDimension = nDimension;
- psTree->nTotalCount = 0;
-
-/* -------------------------------------------------------------------- */
-/* If no max depth was defined, try to select a reasonable one */
-/* that implies approximately 8 shapes per node. */
-/* -------------------------------------------------------------------- */
- if( psTree->nMaxDepth == 0 && hSHP != NULL )
- {
- int nMaxNodeCount = 1;
- int nShapeCount;
-
- ossim_SHPGetInfo( hSHP, &nShapeCount, NULL, NULL, NULL );
- while( nMaxNodeCount*4 < nShapeCount )
- {
- psTree->nMaxDepth += 1;
- nMaxNodeCount = nMaxNodeCount * 2;
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Allocate the root node. */
-/* -------------------------------------------------------------------- */
- psTree->psRoot = ossim_SHPTreeNodeCreate( padfBoundsMin, padfBoundsMax );
-
-/* -------------------------------------------------------------------- */
-/* Assign the bounds to the root node. If none are passed in, */
-/* use the bounds of the provided file otherwise the create */
-/* function will have already set the bounds. */
-/* -------------------------------------------------------------------- */
- if( padfBoundsMin == NULL )
- {
- ossim_SHPGetInfo( hSHP, NULL, NULL,
- psTree->psRoot->adfBoundsMin,
- psTree->psRoot->adfBoundsMax );
- }
-
-/* -------------------------------------------------------------------- */
-/* If we have a file, insert all it's shapes into the tree. */
-/* -------------------------------------------------------------------- */
- if( hSHP != NULL )
- {
- int iShape, nShapeCount;
-
- ossim_SHPGetInfo( hSHP, &nShapeCount, NULL, NULL, NULL );
-
- for( iShape = 0; iShape < nShapeCount; iShape++ )
- {
- ossim_SHPObject *psShape;
-
- psShape = ossim_SHPReadObject( hSHP, iShape );
- ossim_SHPTreeAddShapeId( psTree, psShape );
- ossim_SHPDestroyObject( psShape );
- }
- }
-
- return psTree;
-}
-
-/************************************************************************/
-/* ossim_SHPDestroyTreeNode() */
-/************************************************************************/
-
-static void ossim_SHPDestroyTreeNode( ossim_SHPTreeNode * psTreeNode )
-
-{
- int i;
-
- for( i = 0; i < psTreeNode->nSubNodes; i++ )
- {
- if( psTreeNode->apsSubNode[i] != NULL )
- ossim_SHPDestroyTreeNode( psTreeNode->apsSubNode[i] );
- }
-
- if( psTreeNode->panShapeIds != NULL )
- free( psTreeNode->panShapeIds );
-
- if( psTreeNode->papsShapeObj != NULL )
- {
- for( i = 0; i < psTreeNode->nShapeCount; i++ )
- {
- if( psTreeNode->papsShapeObj[i] != NULL )
- ossim_SHPDestroyObject( psTreeNode->papsShapeObj[i] );
- }
-
- free( psTreeNode->papsShapeObj );
- }
-
- free( psTreeNode );
-}
-
-/************************************************************************/
-/* ossim_SHPDestroyTree() */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPDestroyTree( ossim_SHPTree * psTree )
-
-{
- ossim_SHPDestroyTreeNode( psTree->psRoot );
- free( psTree );
-}
-
-/************************************************************************/
-/* ossim_SHPCheckBoundsOverlap() */
-/* */
-/* Do the given boxes overlap at all? */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_SHPCheckBoundsOverlap( double * padfBox1Min, double * padfBox1Max,
- double * padfBox2Min, double * padfBox2Max,
- int nDimension )
-
-{
- int iDim;
-
- for( iDim = 0; iDim < nDimension; iDim++ )
- {
- if( padfBox2Max[iDim] < padfBox1Min[iDim] )
- return FALSE;
-
- if( padfBox1Max[iDim] < padfBox2Min[iDim] )
- return FALSE;
- }
-
- return TRUE;
-}
-
-/************************************************************************/
-/* ossim_SHPCheckObjectContained() */
-/* */
-/* Does the given shape fit within the indicated extents? */
-/************************************************************************/
-
-static int ossim_SHPCheckObjectContained( ossim_SHPObject * psObject, int nDimension,
- double * padfBoundsMin, double * padfBoundsMax )
-
-{
- if( psObject->dfXMin < padfBoundsMin[0]
- || psObject->dfXMax > padfBoundsMax[0] )
- return FALSE;
-
- if( psObject->dfYMin < padfBoundsMin[1]
- || psObject->dfYMax > padfBoundsMax[1] )
- return FALSE;
-
- if( nDimension == 2 )
- return TRUE;
-
- if( psObject->dfZMin < padfBoundsMin[2]
- || psObject->dfZMax < padfBoundsMax[2] )
- return FALSE;
-
- if( nDimension == 3 )
- return TRUE;
-
- if( psObject->dfMMin < padfBoundsMin[3]
- || psObject->dfMMax < padfBoundsMax[3] )
- return FALSE;
-
- return TRUE;
-}
-
-/************************************************************************/
-/* ossim_SHPTreeSplitBounds() */
-/* */
-/* Split a region into two subregion evenly, cutting along the */
-/* longest dimension. */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPTreeSplitBounds( double *padfBoundsMinIn, double *padfBoundsMaxIn,
- double *padfBoundsMin1, double * padfBoundsMax1,
- double *padfBoundsMin2, double * padfBoundsMax2 )
-
-{
-/* -------------------------------------------------------------------- */
-/* The output bounds will be very similar to the input bounds, */
-/* so just copy over to start. */
-/* -------------------------------------------------------------------- */
- memcpy( padfBoundsMin1, padfBoundsMinIn, sizeof(double) * 4 );
- memcpy( padfBoundsMax1, padfBoundsMaxIn, sizeof(double) * 4 );
- memcpy( padfBoundsMin2, padfBoundsMinIn, sizeof(double) * 4 );
- memcpy( padfBoundsMax2, padfBoundsMaxIn, sizeof(double) * 4 );
-
-/* -------------------------------------------------------------------- */
-/* Split in X direction. */
-/* -------------------------------------------------------------------- */
- if( (padfBoundsMaxIn[0] - padfBoundsMinIn[0])
- > (padfBoundsMaxIn[1] - padfBoundsMinIn[1]) )
- {
- double dfRange = padfBoundsMaxIn[0] - padfBoundsMinIn[0];
-
- padfBoundsMax1[0] = padfBoundsMinIn[0] + dfRange * ossim_SHP_SPLIT_RATIO;
- padfBoundsMin2[0] = padfBoundsMaxIn[0] - dfRange * ossim_SHP_SPLIT_RATIO;
- }
-
-/* -------------------------------------------------------------------- */
-/* Otherwise split in Y direction. */
-/* -------------------------------------------------------------------- */
- else
- {
- double dfRange = padfBoundsMaxIn[1] - padfBoundsMinIn[1];
-
- padfBoundsMax1[1] = padfBoundsMinIn[1] + dfRange * ossim_SHP_SPLIT_RATIO;
- padfBoundsMin2[1] = padfBoundsMaxIn[1] - dfRange * ossim_SHP_SPLIT_RATIO;
- }
-}
-
-/************************************************************************/
-/* ossim_SHPTreeNodeAddShapeId() */
-/************************************************************************/
-
-static int
- ossim_SHPTreeNodeAddShapeId( ossim_SHPTreeNode * psTreeNode, ossim_SHPObject * psObject,
- int nMaxDepth, int nDimension )
-
-{
- int i;
-
-/* -------------------------------------------------------------------- */
-/* If there are subnodes, then consider wiether this object */
-/* will fit in them. */
-/* -------------------------------------------------------------------- */
- if( nMaxDepth > 1 && psTreeNode->nSubNodes > 0 )
- {
- for( i = 0; i < psTreeNode->nSubNodes; i++ )
- {
- if( ossim_SHPCheckObjectContained(psObject, nDimension,
- psTreeNode->apsSubNode[i]->adfBoundsMin,
- psTreeNode->apsSubNode[i]->adfBoundsMax))
- {
- return ossim_SHPTreeNodeAddShapeId( psTreeNode->apsSubNode[i],
- psObject, nMaxDepth-1,
- nDimension );
- }
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* Otherwise, consider creating four subnodes if could fit into */
-/* them, and adding to the appropriate subnode. */
-/* -------------------------------------------------------------------- */
-#if MAX_SUBNODE == 4
- else if( nMaxDepth > 1 && psTreeNode->nSubNodes == 0 )
- {
- double adfBoundsMinH1[4], adfBoundsMaxH1[4];
- double adfBoundsMinH2[4], adfBoundsMaxH2[4];
- double adfBoundsMin1[4], adfBoundsMax1[4];
- double adfBoundsMin2[4], adfBoundsMax2[4];
- double adfBoundsMin3[4], adfBoundsMax3[4];
- double adfBoundsMin4[4], adfBoundsMax4[4];
-
- ossim_SHPTreeSplitBounds( psTreeNode->adfBoundsMin,
- psTreeNode->adfBoundsMax,
- adfBoundsMinH1, adfBoundsMaxH1,
- adfBoundsMinH2, adfBoundsMaxH2 );
-
- ossim_SHPTreeSplitBounds( adfBoundsMinH1, adfBoundsMaxH1,
- adfBoundsMin1, adfBoundsMax1,
- adfBoundsMin2, adfBoundsMax2 );
-
- ossim_SHPTreeSplitBounds( adfBoundsMinH2, adfBoundsMaxH2,
- adfBoundsMin3, adfBoundsMax3,
- adfBoundsMin4, adfBoundsMax4 );
-
- if( ossim_SHPCheckObjectContained(psObject, nDimension,
- adfBoundsMin1, adfBoundsMax1)
- || ossim_SHPCheckObjectContained(psObject, nDimension,
- adfBoundsMin2, adfBoundsMax2)
- || ossim_SHPCheckObjectContained(psObject, nDimension,
- adfBoundsMin3, adfBoundsMax3)
- || ossim_SHPCheckObjectContained(psObject, nDimension,
- adfBoundsMin4, adfBoundsMax4) )
- {
- psTreeNode->nSubNodes = 4;
- psTreeNode->apsSubNode[0] = ossim_SHPTreeNodeCreate( adfBoundsMin1,
- adfBoundsMax1 );
- psTreeNode->apsSubNode[1] = ossim_SHPTreeNodeCreate( adfBoundsMin2,
- adfBoundsMax2 );
- psTreeNode->apsSubNode[2] = ossim_SHPTreeNodeCreate( adfBoundsMin3,
- adfBoundsMax3 );
- psTreeNode->apsSubNode[3] = ossim_SHPTreeNodeCreate( adfBoundsMin4,
- adfBoundsMax4 );
-
- /* recurse back on this node now that it has subnodes */
- return( ossim_SHPTreeNodeAddShapeId( psTreeNode, psObject,
- nMaxDepth, nDimension ) );
- }
- }
-#endif /* MAX_SUBNODE == 4 */
-
-/* -------------------------------------------------------------------- */
-/* Otherwise, consider creating two subnodes if could fit into */
-/* them, and adding to the appropriate subnode. */
-/* -------------------------------------------------------------------- */
-#if MAX_SUBNODE == 2
- else if( nMaxDepth > 1 && psTreeNode->nSubNodes == 0 )
- {
- double adfBoundsMin1[4], adfBoundsMax1[4];
- double adfBoundsMin2[4], adfBoundsMax2[4];
-
- ossim_SHPTreeSplitBounds( psTreeNode->adfBoundsMin, psTreeNode->adfBoundsMax,
- adfBoundsMin1, adfBoundsMax1,
- adfBoundsMin2, adfBoundsMax2 );
-
- if( ossim_SHPCheckObjectContained(psObject, nDimension,
- adfBoundsMin1, adfBoundsMax1))
- {
- psTreeNode->nSubNodes = 2;
- psTreeNode->apsSubNode[0] = ossim_SHPTreeNodeCreate( adfBoundsMin1,
- adfBoundsMax1 );
- psTreeNode->apsSubNode[1] = ossim_SHPTreeNodeCreate( adfBoundsMin2,
- adfBoundsMax2 );
-
- return( ossim_SHPTreeNodeAddShapeId( psTreeNode->apsSubNode[0], psObject,
- nMaxDepth - 1, nDimension ) );
- }
- else if( ossim_SHPCheckObjectContained(psObject, nDimension,
- adfBoundsMin2, adfBoundsMax2) )
- {
- psTreeNode->nSubNodes = 2;
- psTreeNode->apsSubNode[0] = ossim_SHPTreeNodeCreate( adfBoundsMin1,
- adfBoundsMax1 );
- psTreeNode->apsSubNode[1] = ossim_SHPTreeNodeCreate( adfBoundsMin2,
- adfBoundsMax2 );
-
- return( ossim_SHPTreeNodeAddShapeId( psTreeNode->apsSubNode[1], psObject,
- nMaxDepth - 1, nDimension ) );
- }
- }
-#endif /* MAX_SUBNODE == 2 */
-
-/* -------------------------------------------------------------------- */
-/* If none of that worked, just add it to this nodes list. */
-/* -------------------------------------------------------------------- */
- psTreeNode->nShapeCount++;
-
- psTreeNode->panShapeIds = (int *)
- SfRealloc( psTreeNode->panShapeIds,
- sizeof(int) * psTreeNode->nShapeCount );
- psTreeNode->panShapeIds[psTreeNode->nShapeCount-1] = psObject->nShapeId;
-
- if( psTreeNode->papsShapeObj != NULL )
- {
- psTreeNode->papsShapeObj = (ossim_SHPObject **)
- SfRealloc( psTreeNode->papsShapeObj,
- sizeof(void *) * psTreeNode->nShapeCount );
- psTreeNode->papsShapeObj[psTreeNode->nShapeCount-1] = NULL;
- }
-
- return TRUE;
-}
-
-/************************************************************************/
-/* ossim_SHPTreeAddShapeId() */
-/* */
-/* Add a shape to the tree, but don't keep a pointer to the */
-/* object data, just keep the shapeid. */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL
- ossim_SHPTreeAddShapeId( ossim_SHPTree * psTree, ossim_SHPObject * psObject )
-
-{
- psTree->nTotalCount++;
-
- return( ossim_SHPTreeNodeAddShapeId( psTree->psRoot, psObject,
- psTree->nMaxDepth, psTree->nDimension ) );
-}
-
-/************************************************************************/
-/* ossim_SHPTreeCollectShapesIds() */
-/* */
-/* Work function implementing ossim_SHPTreeFindLikelyShapes() on a */
-/* tree node by tree node basis. */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPTreeCollectShapeIds( ossim_SHPTree *hTree, ossim_SHPTreeNode * psTreeNode,
- double * padfBoundsMin, double * padfBoundsMax,
- int * pnShapeCount, int * pnMaxShapes,
- int ** ppanShapeList )
-
-{
- int i;
-
-/* -------------------------------------------------------------------- */
-/* Does this node overlap the area of interest at all? If not, */
-/* return without adding to the list at all. */
-/* -------------------------------------------------------------------- */
- if( !ossim_SHPCheckBoundsOverlap( psTreeNode->adfBoundsMin,
- psTreeNode->adfBoundsMax,
- padfBoundsMin,
- padfBoundsMax,
- hTree->nDimension ) )
- return;
-
-/* -------------------------------------------------------------------- */
-/* Grow the list to hold the shapes on this node. */
-/* -------------------------------------------------------------------- */
- if( *pnShapeCount + psTreeNode->nShapeCount > *pnMaxShapes )
- {
- *pnMaxShapes = (*pnShapeCount + psTreeNode->nShapeCount) * 2 + 20;
- *ppanShapeList = (int *)
- SfRealloc(*ppanShapeList,sizeof(int) * *pnMaxShapes);
- }
-
-/* -------------------------------------------------------------------- */
-/* Add the local nodes shapeids to the list. */
-/* -------------------------------------------------------------------- */
- for( i = 0; i < psTreeNode->nShapeCount; i++ )
- {
- (*ppanShapeList)[(*pnShapeCount)++] = psTreeNode->panShapeIds[i];
- }
-
-/* -------------------------------------------------------------------- */
-/* Recurse to subnodes if they exist. */
-/* -------------------------------------------------------------------- */
- for( i = 0; i < psTreeNode->nSubNodes; i++ )
- {
- if( psTreeNode->apsSubNode[i] != NULL )
- ossim_SHPTreeCollectShapeIds( hTree, psTreeNode->apsSubNode[i],
- padfBoundsMin, padfBoundsMax,
- pnShapeCount, pnMaxShapes,
- ppanShapeList );
- }
-}
-
-/************************************************************************/
-/* ossim_SHPTreeFindLikelyShapes() */
-/* */
-/* Find all shapes within tree nodes for which the tree node */
-/* bounding box overlaps the search box. The return value is */
-/* an array of shapeids terminated by a -1. The shapeids will */
-/* be in order, as hopefully this will result in faster (more */
-/* sequential) reading from the file. */
-/************************************************************************/
-
-/* helper for qsort */
-static int
-compare_ints( const void * a, const void * b)
-{
- return (*(int*)a) - (*(int*)b);
-}
-
-int ossim_SHPAPI_CALL1(*)
- ossim_SHPTreeFindLikelyShapes( ossim_SHPTree * hTree,
- double * padfBoundsMin, double * padfBoundsMax,
- int * pnShapeCount )
-
-{
- int *panShapeList=NULL, nMaxShapes = 0;
-
-/* -------------------------------------------------------------------- */
-/* Perform the search by recursive descent. */
-/* -------------------------------------------------------------------- */
- *pnShapeCount = 0;
-
- ossim_SHPTreeCollectShapeIds( hTree, hTree->psRoot,
- padfBoundsMin, padfBoundsMax,
- pnShapeCount, &nMaxShapes,
- &panShapeList );
-
-/* -------------------------------------------------------------------- */
-/* Sort the id array */
-/* -------------------------------------------------------------------- */
-
- qsort(panShapeList, *pnShapeCount, sizeof(int), compare_ints);
-
- return panShapeList;
-}
-
-/************************************************************************/
-/* ossim_SHPTreeNodeTrim() */
-/* */
-/* This is the recurve version of ossim_SHPTreeTrimExtraNodes() that */
-/* walks the tree cleaning it up. */
-/************************************************************************/
-
-static int ossim_SHPTreeNodeTrim( ossim_SHPTreeNode * psTreeNode )
-
-{
- int i;
-
-/* -------------------------------------------------------------------- */
-/* Trim subtrees, and free subnodes that come back empty. */
-/* -------------------------------------------------------------------- */
- for( i = 0; i < psTreeNode->nSubNodes; i++ )
- {
- if( ossim_SHPTreeNodeTrim( psTreeNode->apsSubNode[i] ) )
- {
- ossim_SHPDestroyTreeNode( psTreeNode->apsSubNode[i] );
-
- psTreeNode->apsSubNode[i] =
- psTreeNode->apsSubNode[psTreeNode->nSubNodes-1];
-
- psTreeNode->nSubNodes--;
-
- i--; /* process the new occupant of this subnode entry */
- }
- }
-
-/* -------------------------------------------------------------------- */
-/* We should be trimmed if we have no subnodes, and no shapes. */
-/* -------------------------------------------------------------------- */
- return( psTreeNode->nSubNodes == 0 && psTreeNode->nShapeCount == 0 );
-}
-
-/************************************************************************/
-/* ossim_SHPTreeTrimExtraNodes() */
-/* */
-/* Trim empty nodes from the tree. Note that we never trim an */
-/* empty root node. */
-/************************************************************************/
-
-void ossim_SHPAPI_CALL
- ossim_SHPTreeTrimExtraNodes( ossim_SHPTree * hTree )
-
-{
- ossim_SHPTreeNodeTrim( hTree->psRoot );
-}
-
-/************************************************************************/
-/* SwapWord() */
-/* */
-/* Swap a 2, 4 or 8 byte word. */
-/************************************************************************/
-
-static void SwapWord( int length, void * wordP )
-
-{
- int i;
- unsigned char temp;
-
- for( i=0; i < length/2; i++ )
- {
- temp = ((unsigned char *) wordP)[i];
- ((unsigned char *)wordP)[i] = ((unsigned char *) wordP)[length-i-1];
- ((unsigned char *) wordP)[length-i-1] = temp;
- }
-}
-
-/************************************************************************/
-/* ossim_SHPSearchDiskTreeNode() */
-/************************************************************************/
-
-static int
- ossim_SHPSearchDiskTreeNode( FILE *fp, double *padfBoundsMin, double *padfBoundsMax,
- int **ppanResultBuffer, int *pnBufferMax,
- int *pnResultCount, int bNeedSwap )
-
-{
- int i;
- int offset;
- int numshapes, numsubnodes;
- double adfNodeBoundsMin[2], adfNodeBoundsMax[2];
-
-/* -------------------------------------------------------------------- */
-/* Read and unswap first part of node info. */
-/* -------------------------------------------------------------------- */
- fread( &offset, 4, 1, fp );
- if ( bNeedSwap ) SwapWord ( 4, &offset );
-
- fread( adfNodeBoundsMin, sizeof(double), 2, fp );
- fread( adfNodeBoundsMax, sizeof(double), 2, fp );
- if ( bNeedSwap )
- {
- SwapWord( 8, adfNodeBoundsMin + 0 );
- SwapWord( 8, adfNodeBoundsMin + 1 );
- SwapWord( 8, adfNodeBoundsMax + 0 );
- SwapWord( 8, adfNodeBoundsMax + 1 );
- }
-
- fread( &numshapes, 4, 1, fp );
- if ( bNeedSwap ) SwapWord ( 4, &numshapes );
-
-/* -------------------------------------------------------------------- */
-/* If we don't overlap this node at all, we can just fseek() */
-/* pass this node info and all subnodes. */
-/* -------------------------------------------------------------------- */
- if( !ossim_SHPCheckBoundsOverlap( adfNodeBoundsMin, adfNodeBoundsMax,
- padfBoundsMin, padfBoundsMax, 2 ) )
- {
- offset += numshapes*sizeof(int) + sizeof(int);
- fseek(fp, offset, SEEK_CUR);
- return TRUE;
- }
-
-/* -------------------------------------------------------------------- */
-/* Add all the shapeids at this node to our list. */
-/* -------------------------------------------------------------------- */
- if(numshapes > 0)
- {
- if( *pnResultCount + numshapes > *pnBufferMax )
- {
- *pnBufferMax = (int) ((*pnResultCount + numshapes + 100) * 1.25);
- *ppanResultBuffer = (int *)
- SfRealloc( *ppanResultBuffer, *pnBufferMax * sizeof(int) );
- }
-
- fread( *ppanResultBuffer + *pnResultCount,
- sizeof(int), numshapes, fp );
-
- if (bNeedSwap )
- {
- for( i=0; i<numshapes; i++ )
- SwapWord( 4, *ppanResultBuffer + *pnResultCount + i );
- }
-
- *pnResultCount += numshapes;
- }
-
-/* -------------------------------------------------------------------- */
-/* Process the subnodes. */
-/* -------------------------------------------------------------------- */
- fread( &numsubnodes, 4, 1, fp );
- if ( bNeedSwap ) SwapWord ( 4, &numsubnodes );
-
- for(i=0; i<numsubnodes; i++)
- {
- if( !ossim_SHPSearchDiskTreeNode( fp, padfBoundsMin, padfBoundsMax,
- ppanResultBuffer, pnBufferMax,
- pnResultCount, bNeedSwap ) )
- return FALSE;
- }
-
- return TRUE;
-}
-
-/************************************************************************/
-/* ossim_SHPSearchDiskTree() */
-/************************************************************************/
-
-int ossim_SHPAPI_CALL1(*)
- ossim_SHPSearchDiskTree( FILE *fp,
- double *padfBoundsMin, double *padfBoundsMax,
- int *pnShapeCount )
-
-{
- int i, bNeedSwap, nBufferMax = 0;
- unsigned char abyBuf[16];
- int *panResultBuffer = NULL;
-
- *pnShapeCount = 0;
-
-/* -------------------------------------------------------------------- */
-/* Establish the byte order on this machine. */
-/* -------------------------------------------------------------------- */
- i = 1;
- if( *((unsigned char *) &i) == 1 )
- bBigEndian = FALSE;
- else
- bBigEndian = TRUE;
-
-/* -------------------------------------------------------------------- */
-/* Read the header. */
-/* -------------------------------------------------------------------- */
- fseek( fp, 0, SEEK_SET );
- fread( abyBuf, 16, 1, fp );
-
- if( memcmp( abyBuf, "SQT", 3 ) != 0 )
- return NULL;
-
- if( (abyBuf[3] == 2 && bBigEndian)
- || (abyBuf[3] == 1 && !bBigEndian) )
- bNeedSwap = FALSE;
- else
- bNeedSwap = TRUE;
-
-/* -------------------------------------------------------------------- */
-/* Search through root node and it's decendents. */
-/* -------------------------------------------------------------------- */
- if( !ossim_SHPSearchDiskTreeNode( fp, padfBoundsMin, padfBoundsMax,
- &panResultBuffer, &nBufferMax,
- pnShapeCount, bNeedSwap ) )
- {
- if( panResultBuffer != NULL )
- free( panResultBuffer );
- *pnShapeCount = 0;
- return NULL;
- }
-/* -------------------------------------------------------------------- */
-/* Sort the id array */
-/* -------------------------------------------------------------------- */
- qsort(panResultBuffer, *pnShapeCount, sizeof(int), compare_ints);
-
- return panResultBuffer;
-}
-
-/************************************************************************/
-/* ossim_SHPGetSubNodeOffset() */
-/* */
-/* Determine how big all the subnodes of this node (and their */
-/* children) will be. This will allow disk based searchers to */
-/* seek past them all efficiently. */
-/************************************************************************/
-
-static int ossim_SHPGetSubNodeOffset( ossim_SHPTreeNode *node)
-{
- int i;
- long offset=0;
-
- for(i=0; i<node->nSubNodes; i++ )
- {
- if(node->apsSubNode[i])
- {
- offset += 4*sizeof(double)
- + (node->apsSubNode[i]->nShapeCount+3)*sizeof(int);
- offset += ossim_SHPGetSubNodeOffset(node->apsSubNode[i]);
- }
- }
-
- return(offset);
-}
-
-/************************************************************************/
-/* ossim_SHPWriteTreeNode() */
-/************************************************************************/
-
-static void ossim_SHPWriteTreeNode( FILE *fp, ossim_SHPTreeNode *node)
-{
- int i,j;
- int offset;
- unsigned char *pabyRec = NULL;
-
- offset = ossim_SHPGetSubNodeOffset(node);
-
- pabyRec = (unsigned char *)
- malloc(sizeof(double) * 4
- + (3 * sizeof(int)) + (node->nShapeCount * sizeof(int)) );
-
- memcpy( pabyRec, &offset, 4);
-
- /* minx, miny, maxx, maxy */
- memcpy( pabyRec+ 4, node->adfBoundsMin+0, sizeof(double) );
- memcpy( pabyRec+12, node->adfBoundsMin+1, sizeof(double) );
- memcpy( pabyRec+20, node->adfBoundsMax+0, sizeof(double) );
- memcpy( pabyRec+28, node->adfBoundsMax+1, sizeof(double) );
-
- memcpy( pabyRec+36, &node->nShapeCount, 4);
- j = node->nShapeCount * sizeof(int);
- memcpy( pabyRec+40, node->panShapeIds, j);
- memcpy( pabyRec+j+40, &node->nSubNodes, 4);
-
- fwrite( pabyRec, 44+j, 1, fp );
- free (pabyRec);
-
- for(i=0; i<node->nSubNodes; i++ )
- {
- if(node->apsSubNode[i])
- ossim_SHPWriteTreeNode( fp, node->apsSubNode[i]);
- }
-}
-
-/************************************************************************/
-/* ossim_SHPWriteTree() */
-/************************************************************************/
-
-int ossim_SHPWriteTree(ossim_SHPTree *tree, const char *filename )
-{
- char signature[4] = "SQT";
- int i;
- char abyBuf[32];
- FILE *fp;
-
-/* -------------------------------------------------------------------- */
-/* Open the output file. */
-/* -------------------------------------------------------------------- */
- fp = fopen(filename, "wb");
- if( fp == NULL )
- {
- return FALSE;
- }
-
-/* -------------------------------------------------------------------- */
-/* Establish the byte order on this machine. */
-/* -------------------------------------------------------------------- */
- i = 1;
- if( *((unsigned char *) &i) == 1 )
- bBigEndian = FALSE;
- else
- bBigEndian = TRUE;
-
-/* -------------------------------------------------------------------- */
-/* Write the header. */
-/* -------------------------------------------------------------------- */
- memcpy( abyBuf+0, signature, 3 );
-
- if( bBigEndian )
- abyBuf[3] = 2; /* New MSB */
- else
- abyBuf[3] = 1; /* New LSB */
-
- abyBuf[4] = 1; /* version */
- abyBuf[5] = 0; /* next 3 reserved */
- abyBuf[6] = 0;
- abyBuf[7] = 0;
-
- fwrite( abyBuf, 8, 1, fp );
-
- fwrite( &(tree->nTotalCount), 4, 1, fp );
-
- /* write maxdepth */
-
- fwrite( &(tree->nMaxDepth), 4, 1, fp );
-
-/* -------------------------------------------------------------------- */
-/* Write all the nodes "in order". */
-/* -------------------------------------------------------------------- */
-
- ossim_SHPWriteTreeNode( fp, tree->psRoot );
-
- fclose( fp );
-
- return TRUE;
-}
--
GitLab