Skip to content
Snippets Groups Projects
Commit 2ebc577d authored by Emmanuel Christophe's avatar Emmanuel Christophe
Browse files

ENH: remove the projection conversion responsibility from the reader and writer

parent aabd3111
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Code/IO/otbKMLVectorDataIO.txx
+ 5
− 76
View file @ 2ebc577d
......@@ -540,46 +540,17 @@ namespace otb
{
// Retrieve data required for georeferencing
typename VectorDataType::PointType origin;
typename VectorDataType::SpacingType spacing;
origin.Fill(0.0);
spacing.Fill(1.0);
std::string projectionRefWkt = data->GetProjectionRef();
bool projectionInformationAvailable = !projectionRefWkt.empty();
ossimProjection * projection = NULL;
//Currently only handling map projections. Should be able to add the sensor models
//quite easily
if (projectionInformationAvailable)
if( projectionRefWkt.compare("GEOGCS[\"GCS_WGS_1984\", DATUM[\"WGS_1984\", SPHEROID[\"WGS_1984\",6378137,298.257223563]], PRIMEM[\"Greenwich\",0], UNIT[\"Degree\",0.017453292519943295]]") )
{
//TODO a mechanism similar to the SensorModel would be needed here for the MapProjection
// to instanciate directly from the keyword list staying at the OTB level.
ossimKeywordlist kwl;
ossimOgcWktTranslator wktTranslator;
projectionInformationAvailable = wktTranslator.toOssimKwl(projectionRefWkt, kwl);
projection = ossimMapProjectionFactory::instance()->createProjection(kwl);
if (!projection)
{
projectionInformationAvailable = false;
}
origin = data->GetOrigin();
spacing = data->GetSpacing();
otbMsgDevMacro(<< "Projection information : " << projectionRefWkt);
otbMsgDevMacro(<< " - Origin : " << origin);
otbMsgDevMacro(<< " - Spacing : " << spacing);
itkWarningMacro(<<"Vector data should be reprojected in geographic coordinates"
<< " before saving to KML. Please use otbVectorDataProjectionFilter. "
<< " The projection information is currently: " << projectionRefWkt
<< "\n We assume that you know what you're doing and proceed to save the KML file.");
}
if (!projectionInformationAvailable)
{
otbMsgDevMacro(<< "Projection information unavailable => spacing set to 1 and origin to 0");
}
//Create the factory
KmlFactory* factory = KmlFactory::GetFactory();
if(factory == NULL)
......@@ -643,13 +614,6 @@ namespace otb
// Create <coordinates>
CoordinatesPtr coordinates = factory->CreateCoordinates();
PointType pointCoord = it.Get()->GetPoint();
if (projectionInformationAvailable)
{
ossimDpt cartoPoint(pointCoord[0] * spacing[0] + origin[0], pointCoord[1] * spacing[1] + origin[1]);
ossimGpt geoPoint = projection->inverse(cartoPoint);
pointCoord[0] = static_cast<typename PointType::ValueType>(geoPoint.lond());
pointCoord[1] = static_cast<typename PointType::ValueType>(geoPoint.latd());
}
if(DataNodeType::Dimension>2)
{
coordinates->add_latlngalt(pointCoord[1],pointCoord[0],pointCoord[2]);
......@@ -698,13 +662,6 @@ namespace otb
while(vIt != vertexList->End())
{
VertexType pointCoord = vIt.Value();
if (projectionInformationAvailable)
{
ossimDpt cartoPoint(pointCoord[0] * spacing[0] + origin[0], pointCoord[1] * spacing[1] + origin[1]);
ossimGpt geoPoint = projection->inverse(cartoPoint);
pointCoord[0] = geoPoint.lond();
pointCoord[1] = geoPoint.latd();
}
if(DataNodeType::Dimension>2)
{
coordinates->add_latlngalt(pointCoord[1],pointCoord[0],pointCoord[2]);
......@@ -761,13 +718,6 @@ namespace otb
while(vIt != vertexList->End())
{
VertexType pointCoord = vIt.Value();
if (projectionInformationAvailable)
{
ossimDpt cartoPoint(pointCoord[0] * spacing[0] + origin[0], pointCoord[1] * spacing[1] + origin[1]);
ossimGpt geoPoint = projection->inverse(cartoPoint);
pointCoord[0] = geoPoint.lond();
pointCoord[1] = geoPoint.latd();
}
if(DataNodeType::Dimension>2)
{
coordinates->add_latlngalt(pointCoord[1],pointCoord[0],pointCoord[2]+1);//Drawing polygon 1m above ground to avoid z-buffer issues
......@@ -784,13 +734,6 @@ namespace otb
//Adding the first point again to close the polygon
vIt = vertexList->Begin();
VertexType pointCoord = vIt.Value();
if (projectionInformationAvailable)
{
ossimDpt cartoPoint(pointCoord[0] * spacing[0] + origin[0], pointCoord[1] * spacing[1] + origin[1]);
ossimGpt geoPoint = projection->inverse(cartoPoint);
pointCoord[0] = geoPoint.lond();
pointCoord[1] = geoPoint.latd();
}
if(DataNodeType::Dimension>2)
{
coordinates->add_latlngalt(pointCoord[1],pointCoord[0],pointCoord[2]+1);//Drawing polygon 1m above ground to avoid z-buffer issues
......@@ -816,13 +759,6 @@ namespace otb
{
vIt = vertexList->Begin();
VertexType pointCoord = vIt.Value();
if (projectionInformationAvailable)
{
ossimDpt cartoPoint(pointCoord[0] * spacing[0] + origin[0], pointCoord[1] * spacing[1] + origin[1]);
ossimGpt geoPoint = projection->inverse(cartoPoint);
pointCoord[0] = geoPoint.lond();
pointCoord[1] = geoPoint.latd();
}
if(DataNodeType::Dimension>2)
{
coordinates->add_latlngalt(pointCoord[1],pointCoord[0],pointCoord[2]+1);//Drawing polygon 1m above ground to avoid z-buffer issues
......@@ -839,13 +775,6 @@ namespace otb
//Adding the first point again to close the polygon
vIt = vertexList->Begin();
VertexType pointCoord = vIt.Value();
if (projectionInformationAvailable)
{
ossimDpt cartoPoint(pointCoord[0] * spacing[0] + origin[0], pointCoord[1] * spacing[1] + origin[1]);
ossimGpt geoPoint = projection->inverse(cartoPoint);
pointCoord[0] = geoPoint.lond();
pointCoord[1] = geoPoint.latd();
}
if(DataNodeType::Dimension>2)
{
coordinates->add_latlngalt(pointCoord[1],pointCoord[0],pointCoord[2]+1);//Drawing polygon 1m above ground to avoid z-buffer issues
......
Code/IO/otbSHPVectorDataIO.txx
+ 26
− 46
View file @ 2ebc577d
......@@ -137,17 +137,11 @@ namespace otb
if (projectionInformationAvailable)
{
otbMsgDevMacro(<< "Projection information : " << projectionRefWkt);
otbMsgDevMacro(<< " - Origin : " << this->m_Origin);
otbMsgDevMacro(<< " - Spacing : " << this->m_Spacing);
}
else
{
this->m_Origin.Fill(0.0);
this->m_Spacing.Fill(1.0);
otbMsgDevMacro(<< "Projection information unavailable => spacing set to 1 and origin to 0");
otbMsgDevMacro(<< "Projection information unavailable");
}
data->SetOrigin(this->m_Origin);
data->SetSpacing(this->m_Spacing);
// For each layer
for(int layerIndex = 0;layerIndex<m_DataSource->GetLayerCount();++layerIndex)
......@@ -443,17 +437,15 @@ namespace otb
{
itkGenericExceptionMacro(<<"Failed to convert OGRGeometry to OGRPoint");
}
SpacingType spacing = this->m_Spacing;
OriginType origin = this->m_Origin;
PointType otbPoint;
otbPoint.Fill(0);
otbPoint[0] = static_cast<typename DataNodeType::PrecisionType>((ogrPoint->getX()-origin[0])/spacing[0]);
otbPoint[1] = static_cast<typename DataNodeType::PrecisionType>((ogrPoint->getY()-origin[1])/spacing[1]);
otbPoint[0] = static_cast<typename DataNodeType::PrecisionType>(ogrPoint->getX());
otbPoint[1] = static_cast<typename DataNodeType::PrecisionType>(ogrPoint->getY());
if(DataNodeType::Dimension > 2)
{
otbPoint[2]=static_cast<typename DataNodeType::PrecisionType>((ogrPoint->getZ()-origin[2])/spacing[2]);
otbPoint[2]=static_cast<typename DataNodeType::PrecisionType>(ogrPoint->getZ());
}
DataNodePointerType node = DataNodeType::New();
......@@ -473,8 +465,6 @@ namespace otb
{
itkGenericExceptionMacro(<<"Failed to convert OGRGeometry to OGRLine");
}
SpacingType spacing = this->m_Spacing;
OriginType origin = this->m_Origin;
LinePointerType line = LineType::New();
......@@ -487,12 +477,12 @@ namespace otb
typename LineType::VertexType vertex;
vertex[0] = (ogrTmpPoint->getX()-origin[0])/spacing[0];
vertex[1] = (ogrTmpPoint->getY()-origin[1])/spacing[1];
vertex[0] = ogrTmpPoint->getX();
vertex[1] = ogrTmpPoint->getY();
if(DataNodeType::Dimension > 2)
{
vertex[2]= (ogrTmpPoint->getZ()-origin[2])/spacing[2];
vertex[2]= ogrTmpPoint->getZ();
}
line->AddVertex(vertex);
......@@ -517,8 +507,6 @@ namespace otb
{
itkGenericExceptionMacro(<<"Failed to convert OGRGeometry to OGRPolygon");
}
SpacingType spacing = this->m_Spacing;
OriginType origin = this->m_Origin;
OGRPoint * ogrTmpPoint = new OGRPoint();
......@@ -530,12 +518,12 @@ namespace otb
{
ogrRing->getPoint(pIndex,ogrTmpPoint);
typename PolygonType::VertexType vertex;
vertex[0] = (ogrTmpPoint->getX()-origin[0])/spacing[0];
vertex[1] = (ogrTmpPoint->getY()-origin[1])/spacing[1];
vertex[0] = ogrTmpPoint->getX();
vertex[1] = ogrTmpPoint->getY();
if(DataNodeType::Dimension > 2)
{
vertex[2]= (ogrTmpPoint->getZ()-origin[2])/spacing[2];
vertex[2]= ogrTmpPoint->getZ();
}
extRing->AddVertex(vertex);
......@@ -552,11 +540,11 @@ namespace otb
ogrRing->getPoint(pIndex,ogrTmpPoint);
typename PolygonType::VertexType vertex;
vertex[0] = (ogrTmpPoint->getX()-origin[0])/spacing[0];
vertex[1] = (ogrTmpPoint->getY()-origin[1])/spacing[1];
vertex[0] = ogrTmpPoint->getX();
vertex[1] = ogrTmpPoint->getY();
if(DataNodeType::Dimension > 2)
{
vertex[2]= (ogrTmpPoint->getZ()-origin[2])/spacing[2];
vertex[2]= ogrTmpPoint->getZ();
}
ring->AddVertex(vertex);
}
......@@ -626,25 +614,17 @@ namespace otb
}
// Retrieve data required for georeferencing
OriginType origin;
SpacingType spacing;
origin.Fill(0.0);
spacing.Fill(1.0);
std::string projectionRefWkt = data->GetProjectionRef();
bool projectionInformationAvailable = !projectionRefWkt.empty();
if (projectionInformationAvailable)
{
origin = data->GetOrigin();
spacing = data->GetSpacing();
otbMsgDevMacro(<< "Projection information : " << projectionRefWkt);
otbMsgDevMacro(<< " - Origin : " << origin);
otbMsgDevMacro(<< " - Spacing : " << spacing);
}
else
{
otbMsgDevMacro(<< "Projection information unavailable => spacing set to 1 and origin to 0");
otbMsgDevMacro(<< "Projection information unavailable");
}
// Retrieving root node
......@@ -721,13 +701,13 @@ namespace otb
case FEATURE_POINT:
{
OGRPoint ogrPoint;
ogrPoint.setX(it.Get()->GetPoint()[0] * spacing[0] + origin[0]);
ogrPoint.setY(it.Get()->GetPoint()[1] * spacing[1] + origin[1]);
ogrPoint.setX(it.Get()->GetPoint()[0]);
ogrPoint.setY(it.Get()->GetPoint()[1]);
if(DataNodeType::Dimension>2)
{
ogrPoint.setZ(it.Get()->GetPoint()[2] * spacing[2] + origin[2]);
ogrPoint.setZ(it.Get()->GetPoint()[2]);
}
if(ogrCollection == NULL)
......@@ -752,11 +732,11 @@ namespace otb
while(vIt != vertexList->End())
{
OGRPoint ogrPoint;
ogrPoint.setX(vIt.Value()[0] * spacing[0] + origin[0]);
ogrPoint.setY(vIt.Value()[1] * spacing[1] + origin[1]);
ogrPoint.setX(vIt.Value()[0]);
ogrPoint.setY(vIt.Value()[1]);
if(DataNodeType::Dimension>2)
{
ogrPoint.setZ(vIt.Value()[2] * spacing[2] + origin[2]);
ogrPoint.setZ(vIt.Value()[2]);
}
ogrLine.addPoint(&ogrPoint);
++vIt;
......@@ -785,11 +765,11 @@ namespace otb
while(vIt != vertexList->End())
{
OGRPoint ogrPoint;
ogrPoint.setX(vIt.Value()[0] * spacing[0] + origin[0]);
ogrPoint.setY(vIt.Value()[1] * spacing[1] + origin[1]);
ogrPoint.setX(vIt.Value()[0]);
ogrPoint.setY(vIt.Value()[1]);
if(DataNodeType::Dimension>2)
{
ogrPoint.setZ(vIt.Value()[2] * spacing[2] + origin[2]);
ogrPoint.setZ(vIt.Value()[2]);
}
ogrExternalRing->addPoint(&ogrPoint);
......@@ -809,11 +789,11 @@ namespace otb
while(vIt != vertexList->End())
{
OGRPoint ogrPoint;
ogrPoint.setX(vIt.Value()[0] * spacing[0] + origin[0]);
ogrPoint.setY(vIt.Value()[1] * spacing[1] + origin[1]);
ogrPoint.setX(vIt.Value()[0]);
ogrPoint.setY(vIt.Value()[1]);
if(DataNodeType::Dimension>2)
{
ogrPoint.setZ(vIt.Value()[2] * spacing[2] + origin[2]);
ogrPoint.setZ(vIt.Value()[2]);
}
ogrInternalRing->addPoint(&ogrPoint);
++vIt;
......
Code/IO/otbVectorDataFileReader.h
+ 0
− 28
View file @ 2ebc577d
......@@ -106,27 +106,6 @@ public :
itkSetStringMacro(FileName);
itkGetStringMacro(FileName);
/** Set the origin for the conversion when reading data.
* \sa GetOrigin() */
itkSetMacro(Origin, PointType);
virtual void SetOrigin( const double origin[VDimension] );
virtual void SetOrigin( const float origin[VDimension] );
itkGetConstReferenceMacro(Origin, PointType);
/** Set the spacing (size of a pixel) for the conversion when reading dat.
* \sa GetSpacing() */
virtual void SetSpacing (const SpacingType & spacing);
virtual void SetSpacing (const double spacing[VDimension]);
virtual void SetSpacing (const float spacing[VDimension]);
itkGetConstReferenceMacro(Spacing, SpacingType);
/** Set/Get the target projection for the vector data */
itkSetStringMacro(TargetProjection);
itkGetStringMacro(TargetProjection);
/** Set/Get the VectorDataIO helper class. Often this is created via the object
* factory mechanism that determines whether a particular VectorDataIO can
* read a certain file. This method provides a way to get the VectorDataIO
......@@ -156,13 +135,6 @@ protected:
std::string m_FileName; // The file to be read
/** Use to fill up information when data are read
* and eventually to do a conversion from geo
* coordinates to image coordinates */
SpacingType m_Spacing;
PointType m_Origin;
std::string m_TargetProjection;
private:
VectorDataFileReader(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
......
Code/IO/otbVectorDataFileReader.txx
+ 0
− 69
View file @ 2ebc577d
......@@ -38,8 +38,6 @@ VectorDataFileReader<TOutputVectorData>
m_VectorDataIO = 0;
m_FileName = "";
m_UserSpecifiedVectorDataIO = false;
m_Spacing.Fill(1);
m_Origin.Fill(0);
}
/**
* Destructor
......@@ -50,69 +48,6 @@ VectorDataFileReader<TOutputVectorData>
{
}
//----------------------------------------------------------------------------
template <class TOutputVectorData>
void
VectorDataFileReader<TOutputVectorData>
::SetSpacing(const SpacingType & spacing )
{
itkDebugMacro("setting Spacing to " << spacing);
if( this->m_Spacing != spacing )
{
this->m_Spacing = spacing;
this->Modified();
}
}
//----------------------------------------------------------------------------
template <class TOutputVectorData>
void
VectorDataFileReader<TOutputVectorData>
::SetSpacing(const double spacing[VDimension] )
{
SpacingType s(spacing);
this->SetSpacing(s);
}
//----------------------------------------------------------------------------
template <class TOutputVectorData>
void
VectorDataFileReader<TOutputVectorData>
::SetSpacing(const float spacing[VDimension] )
{
itk::Vector<float, VDimension> sf(spacing);
SpacingType s;
s.CastFrom( sf );
this->SetSpacing(s);
}
//----------------------------------------------------------------------------
template <class TOutputVectorData>
void
VectorDataFileReader<TOutputVectorData>
::SetOrigin(const double origin[VDimension] )
{
PointType p(origin);
this->SetOrigin( p );
}
//----------------------------------------------------------------------------
template <class TOutputVectorData>
void
VectorDataFileReader<TOutputVectorData>
::SetOrigin(const float origin[VDimension] )
{
itk::Point<float, VDimension> of(origin);
PointType p;
p.CastFrom( of );
this->SetOrigin( p );
}
/** Test whether the given filename exist and it is readable,
this is intended to be called before attempting to use
VectorDataIO classes for actually reading the file. If the file
......@@ -236,10 +171,6 @@ VectorDataFileReader<TOutputVectorData>
}
m_VectorDataIO->SetFileName(m_FileName.c_str());
// m_VectorDataIO->ReadVectorDataInformation();
m_VectorDataIO->SetOrigin(m_Origin);
m_VectorDataIO->SetSpacing(m_Spacing);
m_VectorDataIO->SetTargetProjection(m_TargetProjection);
//Copy MetaDataDictionary from instantiated reader to output VectorData.
output->SetMetaDataDictionary(m_VectorDataIO->GetMetaDataDictionary());
......
Code/IO/otbVectorDataIOBase.h
+ 0
− 27
View file @ 2ebc577d
......@@ -76,27 +76,6 @@ public:
itkGetStringMacro(FileName);
/** Set the origin for the conversion when reading data.
* \sa GetOrigin() */
itkSetMacro(Origin, PointType);
virtual void SetOrigin( const double origin[VDimension] );
virtual void SetOrigin( const float origin[VDimension] );
itkGetConstReferenceMacro(Origin, PointType);
/** Set the spacing (size of a pixel) for the conversion when reading dat.
* \sa GetSpacing() */
virtual void SetSpacing (const SpacingType & spacing);
virtual void SetSpacing (const double spacing[VDimension]);
virtual void SetSpacing (const float spacing[VDimension]);
itkGetConstReferenceMacro(Spacing, SpacingType);
/** Set/Get the target projection for the vector data */
itkSetStringMacro(TargetProjection);
itkGetStringMacro(TargetProjection);
/** Enums used to specify byte order; whether Big Endian or Little Endian.
* Some subclasses use this, some ignore it. */
typedef enum {BigEndian,LittleEndian,OrderNotApplicable} ByteOrder;
......@@ -192,12 +171,6 @@ protected:
/** Return the object to an initialized state, ready to be used */
virtual void Reset(const bool freeDynamic = true);
/** Use to fill up information when data are read
* and eventually to do a conversion from geo
* coordinates to image coordinates */
SpacingType m_Spacing;
PointType m_Origin;
std::string m_TargetProjection;
private:
VectorDataIOBase(const Self&); //purposely not implemented
......
Code/IO/otbVectorDataIOBase.txx
+ 0
− 65
View file @ 2ebc577d
......@@ -31,74 +31,9 @@ VectorDataIOBase<TData>
m_ByteOrder(OrderNotApplicable)
{
Reset(false);
m_Spacing.Fill(1);
m_Origin.Fill(0);
}
//----------------------------------------------------------------------------
template <class TData>
void
VectorDataIOBase<TData>
::SetSpacing(const SpacingType & spacing )
{
itkDebugMacro("setting Spacing to " << spacing);
if( this->m_Spacing != spacing )
{
this->m_Spacing = spacing;
this->Modified();
}
}
//----------------------------------------------------------------------------
template <class TData>
void
VectorDataIOBase<TData>
::SetSpacing(const double spacing[VDimension] )
{
SpacingType s(spacing);
this->SetSpacing(s);
}
//----------------------------------------------------------------------------
template <class TData>
void
VectorDataIOBase<TData>
::SetSpacing(const float spacing[VDimension] )
{
itk::Vector<float, VDimension> sf(spacing);
SpacingType s;
s.CastFrom( sf );
this->SetSpacing(s);
}
//----------------------------------------------------------------------------
template <class TData>
void
VectorDataIOBase<TData>
::SetOrigin(const double origin[VDimension] )
{
PointType p(origin);
this->SetOrigin( p );
}
//----------------------------------------------------------------------------
template <class TData>
void
VectorDataIOBase<TData>
::SetOrigin(const float origin[VDimension] )
{
itk::Point<float, VDimension> of(origin);
PointType p;
p.CastFrom( of );
this->SetOrigin( p );
}
template <class TData>
void
VectorDataIOBase<TData>
......
Code/Projections/otbVectorDataProjectionFilter.txx
+ 24
− 13
View file @ 2ebc577d
......@@ -332,20 +332,7 @@ namespace otb
itk::ExposeMetaData<std::string>(inputDict, MetaDataKey::ProjectionRefKey, m_InputProjectionRef );
}
//If the projection information for the output is provided, propagate it
OutputVectorDataPointer output = this->GetOutput();
itk::MetaDataDictionary & outputDict = output->GetMetaDataDictionary();
if (m_OutputKeywordList.GetSize() != 0)
{
ossimKeywordlist kwl;
m_OutputKeywordList.convertToOSSIMKeywordlist (kwl);
itk::EncapsulateMetaData<ossimKeywordlist>(outputDict, MetaDataKey::OSSIMKeywordlistKey, kwl );
}
if (m_InputProjectionRef.empty())
{
itk::EncapsulateMetaData<std::string>(outputDict, MetaDataKey::ProjectionRefKey, m_OutputProjectionRef );
}
otbMsgDevMacro(<< "Information to instanciate transform: ");
......@@ -358,6 +345,8 @@ namespace otb
otbMsgDevMacro(<< " - Output Origin: " << m_OutputOrigin);
otbMsgDevMacro(<< " - Output Spacing: " << m_OutputSpacing);
bool firstTransformGiveGeo = true;
//*****************************
//Set the input transformation
//*****************************
......@@ -391,6 +380,7 @@ namespace otb
if(m_InputTransform.IsNull())//default if we didn't manage to instantiate it before
{
m_InputTransform = itk::IdentityTransform< double, 2 >::New();
bool firstTransformGiveGeo = false;
otbMsgDevMacro(<< "Input projection set to identity")
}
......@@ -427,9 +417,30 @@ namespace otb
if(m_OutputTransform.IsNull())//default if we didn't manage to instantiate it before
{
m_OutputTransform = itk::IdentityTransform< double, 2 >::New();
if (firstTransformGiveGeo)
{
m_OutputProjectionRef = "GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.017453292519943295]]";
}
otbMsgDevMacro(<< "Output projection set to identity")
}
//If the projection information for the output is provided, propagate it
OutputVectorDataPointer output = this->GetOutput();
itk::MetaDataDictionary & outputDict = output->GetMetaDataDictionary();
if (m_OutputKeywordList.GetSize() != 0)
{
ossimKeywordlist kwl;
m_OutputKeywordList.convertToOSSIMKeywordlist (kwl);
itk::EncapsulateMetaData<ossimKeywordlist>(outputDict, MetaDataKey::OSSIMKeywordlistKey, kwl );
}
if (m_InputProjectionRef.empty())
{
itk::EncapsulateMetaData<std::string>(outputDict, MetaDataKey::ProjectionRefKey, m_OutputProjectionRef );
}
m_Transform->SetFirstTransform(m_InputTransform);
m_Transform->SetSecondTransform(m_OutputTransform);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment