[0 0 0 0]
[0 0 0 0]
[0 0 0 0]
[0 0 0 0]

[0 0 0 0]
[0 0 0 0]
[6 8 0 0]
[5 3 0 0]

[0 0 0 0]
[0 0 0 0]
[6 8 6 6]
[5 3 4 5]

Code: Select all

  template <unsigned int Increment>
  void DTED<Increment>::read_dted_file ( boost::filesystem::path const& file,
                                        int position_index,
                                        boost::shared_ptr<WritePosition> position_obj )
  {
    LOG_FUNCTION
    log::record_debug ( "DTED File: %1%", file );
    log::record_debug ( "Position index: %1%", position_index );
    log::record_debug ( position_obj );

    // Open the target DTED data file
    GDALDataset* dataset = (GDALDataset*) GDALOpen ( file.string().c_str(), GA_ReadOnly );
    if ( dataset == NULL )
      {
        throw MissingFileError(file.string());
      }

    // Get Band for index position                                                                                                         
    // The GDALDataSet only contains one GDALRasterBand for DTED ata                 
    GDALRasterBand* band = dataset->GetRasterBand ( 1 );
    log::record_debug ( "Raster dimension: (%1%,%2%)", band->GetXSize(), band->GetYSize() );
   
    //Calculate the number of blocks of DTED data that must be read
    //from the GDALRasterBand
    int nXSize = 0;
    int nYSize = 0;
    band->GetBlockSize ( &nXSize, &nYSize );
    log::record_debug ( "Raster block dimension: (%1%,%2%)", nXSize, nYSize );
    int nXBlocks = ( band->GetXSize() + nXSize - 1 ) / nXSize;
    int nYBlocks = ( band->GetYSize() + nYSize - 1 ) / nYSize;
    log::record_debug ( "Raster block count: (%1%,%2%)", nXBlocks, nYBlocks );
    log::record_debug ( "Raster data type: %1%", band->GetRasterDataType() );

    // Create destination array
    short* dted_data = Eigen::ei_aligned_new<short>(nXSize * nYSize);

    // Read in blocks of data
    for ( int iYBlock = 0; iYBlock < nYBlocks; iYBlock++ )
      {
        for ( int iXBlock = 0; iXBlock < nXBlocks; iXBlock++ )
          {
            log::record_debug("Reading (%1%,%2%)", iXBlock, iYBlock );
            band->ReadBlock ( iXBlock, iYBlock, dted_data );
          }
      }

    // Line of code below is at DTED_T.inl:213
    // Point of Error: See description below
    Eigen::Map<Eigen::Matrix<int, Eigen::Dynamic, Eigen::Dynamic> > data_map ( dted_data );

    m_elevation_data.block ( position_obj->x, position_obj->Y, position_obj->rows, position_obj->columns ).swap ( data_map );
}


Code: Select all

Eigen::MatrixXi m_elevation_data;


Code: Select all

namespace dted {

  struct WritePosition
  {
    int X;
    int Y;
    int columns;
    int rows;
  };

} // namespace dted 

[ 50%] Building CXX object CMakeFiles/dted.dir/dted/Interface.cpp.o
In file included from /home/storri/Documents/Code/dted/dted/DTED_T.hpp:89,
from /home/storri/Documents/Code/dted/dted/DTED_Level_One.hpp:4,
from /home/storri/Documents/Code/dted/dted/DTED_Factory.hpp:5,
from /home/storri/Documents/Code/dted/dted/Interface.cpp:3:
/home/storri/Documents/Code/dted/dted/DTED_T.inl: In member function ‘void dted::DTED<Increment>::read_dted_file(const boost::filesystem2::path&, int, boost::shared_ptr<dted::WritePosition>) [with unsigned int Increment = 1200u]’:
/home/storri/Documents/Code/dted/dted/DTED_T.inl:96: instantiated from ‘void dted::DTED<Increment>::init_impl(std::string, double, double, unsigned int, unsigned int) [with unsigned int Increment = 1200u]’
/home/storri/Documents/Code/dted/dted/DTED_Level_One.hpp:25: instantiated from here
/home/storri/Documents/Code/dted/dted/DTED_T.inl:213: error: no matching function for call to ‘Eigen::Map<Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, 1>::Map(short int*&)’
/usr/local/include/eigen2/Eigen/src/Core/Map.h:78: note: candidates are: Eigen::Map<MatrixType, PacketAccess>::Map(const typename Eigen::ei_traits<Eigen::Map<MatrixType, _PacketAccess> >::Scalar*, int, int) [with MatrixType = Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, int PacketAccess = 1]
/usr/local/include/eigen2/Eigen/src/Core/Map.h:76: note: Eigen::Map<MatrixType, PacketAccess>::Map(const typename Eigen::ei_traits<Eigen::Map<MatrixType, _PacketAccess> >::Scalar*, int) [with MatrixType = Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, int PacketAccess = 1]
/usr/local/include/eigen2/Eigen/src/Core/Map.h:74: note: Eigen::Map<MatrixType, PacketAccess>::Map(const typename Eigen::ei_traits<Eigen::Map<MatrixType, _PacketAccess> >::Scalar*) [with MatrixType = Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, int PacketAccess = 1]
/usr/local/include/eigen2/Eigen/src/Core/Map.h:61: note: Eigen::Map<Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, 1>::Map(const Eigen::Map<Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, 1>&)
/home/storri/Documents/Code/dted/dted/DTED_T.inl: In member function ‘void dted::DTED<Increment>::read_dted_file(const boost::filesystem2::path&, int, boost::shared_ptr<dted::WritePosition>) [with unsigned int Increment = 3600u]’:
/home/storri/Documents/Code/dted/dted/DTED_T.inl:96: instantiated from ‘void dted::DTED<Increment>::init_impl(std::string, double, double, unsigned int, unsigned int) [with unsigned int Increment = 3600u]’
/home/storri/Documents/Code/dted/dted/DTED_Level_Two.hpp:25: instantiated from here
/home/storri/Documents/Code/dted/dted/DTED_T.inl:213: error: no matching function for call to ‘Eigen::Map<Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, 1>::Map(short int*&)’
/usr/local/include/eigen2/Eigen/src/Core/Map.h:78: note: candidates are: Eigen::Map<MatrixType, PacketAccess>::Map(const typename Eigen::ei_traits<Eigen::Map<MatrixType, _PacketAccess> >::Scalar*, int, int) [with MatrixType = Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, int PacketAccess = 1]
/usr/local/include/eigen2/Eigen/src/Core/Map.h:76: note: Eigen::Map<MatrixType, PacketAccess>::Map(const typename Eigen::ei_traits<Eigen::Map<MatrixType, _PacketAccess> >::Scalar*, int) [with MatrixType = Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, int PacketAccess = 1]
/usr/local/include/eigen2/Eigen/src/Core/Map.h:74: note: Eigen::Map<MatrixType, PacketAccess>::Map(const typename Eigen::ei_traits<Eigen::Map<MatrixType, _PacketAccess> >::Scalar*) [with MatrixType = Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, int PacketAccess = 1]
/usr/local/include/eigen2/Eigen/src/Core/Map.h:61: note: Eigen::Map<Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, 1>::Map(const Eigen::Map<Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, 1>&)
make[2]: *** [CMakeFiles/dted.dir/dted/Interface.cpp.o] Error 1
make[1]: *** [CMakeFiles/dted.dir/all] Error 2
make: *** [all] Error 2

storri wrote:/home/storri/Documents/Code/dted/dted/DTED_T.inl:213: error: no matching function for call to ‘Eigen::Map<Eigen::Matrix<int, 10000, 10000, 2, 10000, 10000>, 1>::Map(short int*&)’

bjacob wrote:In Eigen, the type suffix 'i' means int, not short int. Eigen 2 has limited support for other integer types out of the box (but it's very easy to add support for more types). Eigen 3, on the other hand, supports every integer type out of the box, so you can without problem do:

typedef Matrix<short int, Dynamic, Dynamic> MatrixXsi;

and then use that type with Map<MatrixXsi>.

Code: Select all

In file included from /usr/local/include/eigen2/Eigen/Core:130,
                 from dted_stl_map_test.cpp:2:
/usr/local/include/eigen2/Eigen/src/Core/Matrix.h: In instantiation of ‘Eigen::ei_traits<Eigen::Matrix<short int, 10000, 10000, 2, 10000, 10000> >’:
/usr/local/include/eigen2/Eigen/src/Core/MatrixBase.h:62:   instantiated from ‘Eigen::MatrixBase<Eigen::Matrix<short int, 10000, 10000, 2, 10000, 10000> >’
/usr/local/include/eigen2/Eigen/src/Core/Matrix.h:129:   instantiated from ‘Eigen::Matrix<short int, 10000, 10000, 2, 10000, 10000>’
dted_stl_map_test.cpp:218:   instantiated from here
/usr/local/include/eigen2/Eigen/src/Core/Matrix.h:116: error: incomplete type ‘Eigen::NumTraits<short int>’ used in nested name specifier


storri wrote:
I used the example you gave above with the Eigen 2 code. I now get an error with the NumTraits which is what I think you said Eigen 3 takes care of for me. Which version of Eigen was I suppose to use with the typedef for MatrixXsi you gave me?