spatialreference.hpp

Go to the documentation of this file.

/******************************************************************************
 * $Id: lasspatialreference.hpp 1091 2009-03-10 20:17:35Z hobu $
 *
 * Project:  libLAS - http://liblas.org - A BSD library for LAS format data.
 * Purpose:  LAS Spatial Reference class 
 * Author:   Howard Butler, hobu.inc@gmail.com
 *
 ******************************************************************************
 * Copyright (c) 2009, Howard Butler
 *
 * All rights reserved.
 * 
 * Redistribution and use in source and binary forms, with or without 
 * modification, are permitted provided that the following 
 * conditions are met:
 * 
 *     * Redistributions of source code must retain the above copyright 
 *       notice, this list of conditions and the following disclaimer.
 *     * Redistributions in binary form must reproduce the above copyright 
 *       notice, this list of conditions and the following disclaimer in 
 *       the documentation and/or other materials provided 
 *       with the distribution.
 *     * Neither the name of the Martin Isenburg or Iowa Department 
 *       of Natural Resources nor the names of its contributors may be 
 *       used to endorse or promote products derived from this software 
 *       without specific prior written permission.
 * 
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 
 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 
 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS 
 * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE 
 * COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, 
 * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, 
 * BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS 
 * OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED 
 * AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, 
 * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT 
 * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY 
 * OF SUCH DAMAGE.
 ****************************************************************************/

#ifndef LIBLAS_LASSPATIALREFERENCE_HPP_INCLUDED
#define LIBLAS_LASSPATIALREFERENCE_HPP_INCLUDED

#include <liblas/detail/fwd.hpp>
#include <liblas/detail/private_utility.hpp>
#include <liblas/variablerecord.hpp>
#include <liblas/exception.hpp>
#include <liblas/capi/las_config.h>
#include <liblas/export.hpp>
#include <liblas/external/property_tree/ptree.hpp>

// std
#include <stdexcept> // std::out_of_range
#include <cstdlib> // std::size_t
#include <string>

// Fake out the compiler if we don't have libgeotiff includes already
#if !defined(__geotiff_h_) && !defined(LIBGEOTIFF_GEOTIFF_H_)
typedef struct GTIFS *GTIF;
#endif
#if !defined(__geo_simpletags_h_) && !defined(LIBGEOTIFF_GEO_SIMPLETAGS_H_)
typedef struct ST_TIFFS *ST_TIFF;
#endif

namespace liblas {

class LAS_DLL SpatialReference
{
public:
    enum WKTModeFlag
    {
        eHorizontalOnly = 1,
        eCompoundOK = 2
    };

    enum GeoVLRType
    {
        eGeoTIFF = 1,
        eOGRWKT = 2
    };

    SpatialReference();

    ~SpatialReference();

    SpatialReference(std::vector<VariableRecord> const& vlrs);

    SpatialReference(SpatialReference const& other);

    SpatialReference& operator=(SpatialReference const& rhs);

    bool operator==(const SpatialReference& other) const;
    bool operator!=(const SpatialReference& other) const { return !(*this == other); }
    
    const GTIF* GetGTIF();

    void SetGTIF(GTIF* pgtiff, ST_TIFF* ptiff);

    std::string GetWKT(WKTModeFlag mode_flag = eHorizontalOnly) const;
    std::string GetWKT(WKTModeFlag mode_flag, bool pretty) const;
    
    void SetWKT(std::string const& v);

    void SetVerticalCS(int32_t verticalCSType, 
                       std::string const& citation = std::string(0),
                       int32_t verticalDatum = -1,
                       int32_t verticalUnits = 9001);

    void SetFromUserInput(std::string const& v);
        
    std::string GetProj4() const;

    void SetProj4(std::string const& v);
    
    void SetVLRs(std::vector<VariableRecord> const& vlrs);
    
    void AddVLR(VariableRecord const& vlr);
    
    std::vector<VariableRecord> GetVLRs() const;

    void ClearVLRs( GeoVLRType eType );

    liblas::property_tree::ptree GetPTree() const;    
private:

    // FIXME: Define as shared_ptr<GTIF> with custom deleter to get rid of bloated mem management, unsafe anyway --mloskot
    GTIF*       m_gtiff;
    ST_TIFF*    m_tiff;

    std::string m_wkt;

    std::vector<VariableRecord> m_vlrs;
    bool IsGeoVLR(VariableRecord const& vlr) const;
    std::string GetGTIFFText() const;

    void ResetVLRs();
};

} // namespace liblas

LAS_DLL std::ostream& operator<<(std::ostream& ostr, const liblas::SpatialReference& srs);

LAS_C_START
#if defined(__geotiff_h_) || defined(LIBGEOTIFF_GEOTIFF_H_)
#if defined(GEO_NORMALIZE_H_INCLUDED)
char LAS_DLL * GTIFGetOGISDefn(GTIF*, GTIFDefn*);
#endif

int  LAS_DLL   GTIFSetFromOGISDefn(GTIF*, const char*);
void SetLinearUnitCitation(GTIF* psGTIF, char* pszLinearUOMName);

#if defined(_OGR_SRS_API_H_INCLUDED)
void SetGeogCSCitation(GTIF* psGTIF, OGRSpatialReference* poSRS, char* angUnitName, int nDatum, short nSpheroid);
#endif // defined _OGR_SRS_API_H_INCLUDED
#endif // defined __geotiff_h_ || defined LIBGEOTIFF_GEOTIFF_H_

LAS_C_END

#endif // LIBLAS_LASSPATIALREFERENCE_HPP_INCLUDED