From 2584136f75a82969c46c6c8a5bf196dc9714acb6 Mon Sep 17 00:00:00 2001
From: Otmane Lahlou <otmane.lahlou@c-s.fr>
Date: Fri, 20 Aug 2010 11:44:03 +0200
Subject: [PATCH] ADD : add the MapServer products writer classes
---
Code/IO/otbMapFileProductWriter.h | 186 ++++++++++
Code/IO/otbMapFileProductWriter.txx | 538 ++++++++++++++++++++++++++++
2 files changed, 724 insertions(+)
create mode 100644 Code/IO/otbMapFileProductWriter.h
create mode 100644 Code/IO/otbMapFileProductWriter.txx
diff --git a/Code/IO/otbMapFileProductWriter.h b/Code/IO/otbMapFileProductWriter.h
new file mode 100644
index 0000000000..1f541d82b5
--- /dev/null
+++ b/Code/IO/otbMapFileProductWriter.h
@@ -0,0 +1,186 @@
+/*=========================================================================
+
+ Program: ORFEO Toolbox
+ Language: C++
+ Date: $Date$
+ Version: $Revision$
+
+
+ Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
+ See OTBCopyright.txt for details.
+
+
+ This software is distributed WITHOUT ANY WARRANTY; without even
+ the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+ PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+#ifndef __otbMapFileProductWriter_h
+#define __otbMapFileProductWriter_h
+
+#include "itkProcessObject.h"
+#include "itkObjectFactory.h"
+
+// Image Tiling
+#include "otbMultiChannelExtractROI.h"
+#include "otbImageFileWriter.h"
+#include "otbVectorRescaleIntensityImageFilter.h"
+#include <ogr_spatialref.h>
+#include "otbGenericRSTransform.h"
+#include "otbStreamingShrinkImageFilter.h"
+
+// FileMap necessary includes
+#include "otbVectorData.h"
+#include "otbVectorDataFileWriter.h"
+
+
+namespace otb
+{
+
+/** \class MapFileProductWriter
+ * \brief This class writes Map Product file format (MapFile and MapFiles)
+ *
+ * \ingroup IO
+ *
+ */
+template <class TInputImage>
+class ITK_EXPORT MapFileProductWriter : public itk::ProcessObject
+{
+public:
+ /** Standard class typedefs. */
+ typedef MapFileProductWriter Self;
+ typedef itk::ProcessObject Superclass;
+ typedef itk::SmartPointer<Self> Pointer;
+ typedef itk::SmartPointer<const Self> ConstPointer;
+
+ /** Method for creation through the object factory. */
+ itkNewMacro(Self);
+
+ /** Run-time type information (and related methods). */
+ itkTypeMacro(MapFileProductWriter, itk::ProcessObject);
+
+ /** Some convenient typedefs. */
+ typedef TInputImage InputImageType;
+ typedef typename InputImageType::InternalPixelType InternalPixelType;
+ typedef typename InputImageType::SizeType SizeType;
+ typedef typename InputImageType::IndexType IndexType;
+ typedef unsigned char OutputPixelType;
+ typedef typename InputImageType::Pointer InputImagePointer;
+ typedef typename InputImageType::RegionType InputImageRegionType;
+ typedef typename InputImageType::PixelType InputImagePixelType;
+
+ typedef VectorData<double,2> VectorDataType;
+ typedef typename VectorDataType::DataNodeType DataNodeType;
+ typedef typename DataNodeType::PolygonType PolygonType;
+ typedef typename PolygonType::VertexType VertexType;
+
+ typedef VectorDataFileWriter<VectorDataType> VectorDataFileWriterType;
+
+
+ /// Multi channels Extract ROI filter
+ typedef MultiChannelExtractROI<InternalPixelType, OutputPixelType> VectorImageExtractROIFilterType;
+
+ // Writer
+ typedef ImageFileWriter< VectorImage<OutputPixelType> > VectorWriterType;
+
+ // Resampler
+ typedef StreamingShrinkImageFilter<InputImageType,InputImageType > StreamingShrinkImageFilterType;
+
+ // Intensity Rescale
+ typedef VectorRescaleIntensityImageFilter<InputImageType,
+ InputImageType> VectorRescaleIntensityImageFilterType;
+
+ // Transformer
+ typedef GenericRSTransform<> TransformType;
+ typedef TransformType::InputPointType InputPointType;
+ typedef TransformType::OutputPointType OutputPointType;
+
+ /** Dimension of input image. */
+ itkStaticConstMacro(InputImageDimension, unsigned int,
+ InputImageType::ImageDimension);
+
+ virtual void SetInput( const InputImageType *image);
+ virtual void SetInput( unsigned int, const TInputImage * image);
+ const InputImageType * GetInput(void);
+ const InputImageType * GetInput(unsigned int idx);
+
+ // Accessors macros
+ itkSetStringMacro(FileName);
+ itkSetMacro(TileSize,unsigned int);
+ itkGetMacro(TileSize,unsigned int);
+
+ /** Update Method */
+ virtual void Update()
+ {
+ this->Write();
+ }
+
+protected:
+ MapFileProductWriter();
+ virtual ~MapFileProductWriter();
+ void PrintSelf(std::ostream& os, itk::Indent indent) const;
+
+ /**Method for Tiling the input image*/
+ virtual void Tiling();
+
+ /** Method for writing */
+ virtual void Write();
+
+private:
+ MapFileProductWriter(const Self &); //purposely not implemented
+ void operator =(const Self&); //purposely not implemented
+
+ /** Mehtod to initialize the variables*/
+ virtual void Initialize();
+
+ /**MapFile IndexTule build*/
+ virtual void AddBBoxToIndexTile(OutputPointType lowerLeftCorner,
+ OutputPointType lowerRightCorner,
+ OutputPointType upperRightCorner,
+ OutputPointType upperLeftCorner,
+ unsigned int x, unsigned int y);
+
+ /** Method To Generate the mapFile*/
+ virtual void GenerateMapFile();
+
+ std::string m_Path;
+ bool m_UseExtendMode;
+ InputImagePointer m_VectorImage;
+ InputImagePointer m_ResampleVectorImage;
+
+ // Extract ROI
+ typename VectorImageExtractROIFilterType::Pointer m_VectorImageExtractROIFilter;
+
+ // Writer
+ typename VectorWriterType::Pointer m_VectorWriter;
+
+ // Resampler
+ typename StreamingShrinkImageFilterType::Pointer m_StreamingShrinkImageFilter;
+
+ // Rescale intensity
+ typename VectorRescaleIntensityImageFilterType::Pointer m_VectorRescaleIntensityImageFilter;
+
+ // Transformer
+ typename TransformType::Pointer m_Transform;
+
+ // VectorData Pointer Type for tileindex shapefile
+ typename VectorDataType::Pointer m_VectorDataIndexTile;
+ typename DataNodeType::Pointer m_Polygon;
+ typename DataNodeType::Pointer m_Folder;
+
+ // Tile size
+ int m_TileSize;
+ int m_CurrentDepth;
+
+ // File and path name
+ std::string m_FileName;
+ std::string m_IndexShapeFileName;
+};
+
+} // end namespace otb
+
+#ifndef OTB_MANUAL_INSTANTIATION
+#include "otbMapFileProductWriter.txx"
+#endif
+
+#endif
diff --git a/Code/IO/otbMapFileProductWriter.txx b/Code/IO/otbMapFileProductWriter.txx
new file mode 100644
index 0000000000..8bbf915245
--- /dev/null
+++ b/Code/IO/otbMapFileProductWriter.txx
@@ -0,0 +1,538 @@
+/*=========================================================================
+
+ Program: ORFEO Toolbox
+ Language: C++
+ Date: $Date$
+ Version: $Revision$
+
+
+ Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
+ See OTBCopyright.txt for details.
+
+
+ This software is distributed WITHOUT ANY WARRANTY; without even
+ the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+ PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+#ifndef __otbMapFileProductWriter_txx
+#define __otbMapFileProductWriter_txx
+
+#include "otbMapFileProductWriter.h"
+#include "itksys/SystemTools.hxx"
+
+namespace otb
+{
+template <class TInputImage>
+MapFileProductWriter<TInputImage>
+::MapFileProductWriter(): m_TileSize(256),m_UseExtendMode(true),m_Polygon(),m_Folder()
+{
+ // Modify superclass default values, can be overridden by subclasses
+ this->SetNumberOfRequiredInputs(1);
+}
+
+template <class TInputImage>
+MapFileProductWriter<TInputImage>
+::~MapFileProductWriter()
+{
+}
+
+/**
+ *
+ */
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::SetInput(const InputImageType *input)
+{
+ // Process object is not const-correct so the const_cast is required here
+ this->ProcessObject::SetNthInput(0,
+ const_cast< InputImageType * >( input ) );
+}
+
+
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::SetInput( unsigned int index, const TInputImage * image )
+{
+ // Process object is not const-correct so the const_cast is required here
+ this->ProcessObject::SetNthInput(index,
+ const_cast< TInputImage *>( image ) );
+}
+
+/**
+ *
+ */
+template <class TInputImage>
+const typename MapFileProductWriter<TInputImage>::InputImageType *
+MapFileProductWriter<TInputImage>
+::GetInput(void)
+{
+ if (this->GetNumberOfInputs() < 1)
+ {
+ return 0;
+ }
+
+ return static_cast<const TInputImage * >
+ (this->ProcessObject::GetInput(0) );
+}
+
+/**
+ *
+ */
+template <class TInputImage>
+const typename MapFileProductWriter<TInputImage>::InputImageType *
+MapFileProductWriter<TInputImage>
+::GetInput(unsigned int idx)
+{
+ return static_cast< const TInputImage * >
+ (this->ProcessObject::GetInput(idx));
+}
+
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::Write()
+{
+ // Get the input Image
+ m_VectorImage = const_cast<TInputImage *>(this->GetInput());
+
+ // Initialize the filename, the vectordatas
+ this->Initialize();
+
+ // Do the tiling
+ this->Tiling();
+
+ // Generate the mapFile
+ this->GenerateMapFile();
+
+ // Write the IndexTile
+ typename VectorDataFileWriterType::Pointer writer = VectorDataFileWriterType::New();
+ writer->SetFileName(m_IndexShapeFileName);
+ writer->SetInput(m_VectorDataIndexTile);
+ writer->Update();
+}
+
+
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::Initialize()
+{
+ // check that the right extension is given : expected .map
+ if (itksys::SystemTools::GetFilenameLastExtension(m_FileName) != ".map")
+ {
+ itkExceptionMacro(<<itksys::SystemTools::GetFilenameLastExtension(m_Path)<<" is a wrong Extension FileName : Expected .map");
+ }
+
+ // the path to the directory where the mapfile will be written
+ m_Path = itksys::SystemTools::GetFilenamePath(m_FileName);
+
+ // the index shapefile filename
+ std::ostringstream tempIndexShapeName;
+ tempIndexShapeName << m_Path<<"/";
+ tempIndexShapeName << itksys::SystemTools::GetFilenameWithoutExtension(m_FileName);
+ tempIndexShapeName << "_index.shp";
+ m_IndexShapeFileName = tempIndexShapeName.str();
+
+ // Intitialize the vectordata to build the indexTile
+ m_VectorDataIndexTile = VectorDataType::New();
+ m_VectorDataIndexTile->SetProjectionRef(m_VectorImage->GetProjectionRef());
+ DataNodeType::Pointer root = m_VectorDataIndexTile->GetDataTree()->GetRoot()->Get();
+ DataNodeType::Pointer document = DataNodeType::New();
+ m_Folder = DataNodeType::New();
+
+ document->SetNodeType(otb::DOCUMENT);
+ m_Folder->SetNodeType(otb::FOLDER);
+
+ document->SetNodeId("DOCUMENT");
+ m_Folder->SetNodeId("FOLDER");
+
+ m_VectorDataIndexTile->GetDataTree()->Add(document, root);
+ m_VectorDataIndexTile->GetDataTree()->Add(m_Folder, document);
+}
+
+
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::Tiling()
+{
+ unsigned int numberOfChannel = m_VectorImage->GetNumberOfComponentsPerPixel();
+ unsigned int curIdx = 0;
+
+ /** Image statistics*/
+ typename InputImageType::PixelType inMin(numberOfChannel), inMax(numberOfChannel), outMin(numberOfChannel), outMax(
+ numberOfChannel);
+ outMin.Fill(0);
+ outMax.Fill(255);
+
+ // Build wgs ref to compute long/lat
+ OGRSpatialReference oSRS;
+ oSRS.SetWellKnownGeogCS("WGS84");
+ char * wgsRef = NULL;
+ oSRS.exportToWkt(&wgsRef);
+
+ // Update image base information
+ m_VectorImage->UpdateOutputInformation();
+
+ // Get the image size
+ SizeType size;
+ size = m_VectorImage->GetLargestPossibleRegion().GetSize();
+
+ int sizeX = size[0];
+ int sizeY = size[1];
+
+ // Compute max depth
+ int maxDepth =
+ static_cast<int>(max(vcl_ceil(vcl_log(static_cast<float>(sizeX) / static_cast<float>(m_TileSize)) / vcl_log(2.0)),
+ vcl_ceil(vcl_log(static_cast<float>(sizeY) / static_cast<float>(m_TileSize)) / vcl_log(2.0))));
+
+// // Compute nbTile : Keep this for progress Bar ?
+// int nbTile = 0;
+
+// for (int i = 0; i <= maxDepth; i++)
+// {
+// int ratio = static_cast<int>(vcl_pow(2., (maxDepth - i)));
+// nbTile += (((sizeX / ratio) / m_TileSize) + 1) * (((sizeY / ratio) / m_TileSize) + 1);
+// }
+
+ // Extract size & index
+ SizeType extractSize;
+ IndexType extractIndex;
+
+ std::cout << "maxDepth " << maxDepth<< std::endl;
+
+ for (int depth = 0; depth <= maxDepth; depth++)
+ {
+
+ // update the attribute value Current Depth
+ m_CurrentDepth = depth;
+
+ // Resample image to the max Depth
+ int sampleRatioValue = static_cast<int>(vcl_pow(2., (maxDepth - depth)));
+
+ if (sampleRatioValue > 1)
+ {
+ m_StreamingShrinkImageFilter = StreamingShrinkImageFilterType::New();
+
+ m_StreamingShrinkImageFilter->SetShrinkFactor(sampleRatioValue);
+ m_StreamingShrinkImageFilter->SetInput(m_VectorImage);
+
+ m_VectorRescaleIntensityImageFilter = VectorRescaleIntensityImageFilterType::New();
+ m_VectorRescaleIntensityImageFilter->SetInput(m_StreamingShrinkImageFilter->GetOutput()); //m_ResampleVectorImage);
+ m_VectorRescaleIntensityImageFilter->SetOutputMinimum(outMin);
+ m_VectorRescaleIntensityImageFilter->SetOutputMaximum(outMax);
+
+ if (depth == 0)
+ {
+ m_VectorRescaleIntensityImageFilter->Update();
+ inMin = m_VectorRescaleIntensityImageFilter->GetInputMinimum();
+ inMax = m_VectorRescaleIntensityImageFilter->GetInputMaximum();
+ }
+ else
+ {
+ m_VectorRescaleIntensityImageFilter->SetInputMinimum(inMin);
+ m_VectorRescaleIntensityImageFilter->SetInputMaximum(inMax);
+ m_VectorRescaleIntensityImageFilter->SetAutomaticInputMinMaxComputation(false);
+ }
+
+ // New resample vector image
+ m_ResampleVectorImage = m_VectorRescaleIntensityImageFilter->GetOutput();
+ }
+ else
+ {
+ m_VectorRescaleIntensityImageFilter = VectorRescaleIntensityImageFilterType::New();
+ m_VectorRescaleIntensityImageFilter->SetInput(m_VectorImage);
+ m_VectorRescaleIntensityImageFilter->SetOutputMinimum(outMin);
+ m_VectorRescaleIntensityImageFilter->SetOutputMaximum(outMax);
+
+ m_VectorRescaleIntensityImageFilter->SetInputMinimum(inMin);
+ m_VectorRescaleIntensityImageFilter->SetInputMaximum(inMax);
+ m_VectorRescaleIntensityImageFilter->SetAutomaticInputMinMaxComputation(false);
+
+ m_ResampleVectorImage = m_VectorRescaleIntensityImageFilter->GetOutput();
+ }
+
+ m_ResampleVectorImage->UpdateOutputInformation();
+
+ // Get the image size
+ size = m_ResampleVectorImage->GetLargestPossibleRegion().GetSize();
+
+ sizeX = size[0];
+ sizeY = size[1];
+
+ int x = 0;
+ int y = 0;
+
+ // Tiling resample image
+ for (int tx = 0; tx < sizeX; tx += m_TileSize)
+ {
+ for (int ty = 0; ty < sizeY; ty += m_TileSize)
+ {
+ if ((tx + m_TileSize) >= sizeX)
+ {
+ extractIndex[0] = tx;
+ extractSize[0] = sizeX - tx;
+ }
+ else
+ {
+ extractIndex[0] = tx;
+ extractSize[0] = m_TileSize;
+ }
+
+ if ((ty + m_TileSize) >= sizeY)
+ {
+ extractIndex[1] = ty;
+ extractSize[1] = sizeY - ty;
+ }
+ else
+ {
+ extractIndex[1] = ty;
+ extractSize[1] = m_TileSize;
+ }
+
+ // Create directory where to store generated tiles
+ std::ostringstream path;
+ path << m_Path<<"/tiles";
+
+ ossimFilename cachingDir(path.str());
+ cachingDir.createDirectory();
+
+ // Generate Tile filename
+ std::ostringstream ossFileName;
+ ossFileName << m_Path<<"/";
+ ossFileName << "tiles/tile_";
+ ossFileName << m_CurrentDepth<<"_";
+ ossFileName << x<<"_";
+ ossFileName << y;
+ ossFileName << ".tif";
+
+ // Extract ROI
+ m_VectorImageExtractROIFilter = VectorImageExtractROIFilterType::New();
+
+ // Set extract roi parameters
+ m_VectorImageExtractROIFilter->SetStartX(extractIndex[0]);
+ m_VectorImageExtractROIFilter->SetStartY(extractIndex[1]);
+ m_VectorImageExtractROIFilter->SetSizeX(extractSize[0]);
+ m_VectorImageExtractROIFilter->SetSizeY(extractSize[1]);
+
+ // Set Channel to extract
+ m_VectorImageExtractROIFilter->SetChannel(1);
+ m_VectorImageExtractROIFilter->SetChannel(2);
+ m_VectorImageExtractROIFilter->SetChannel(3);
+
+ // Set extract roi input
+ m_VectorImageExtractROIFilter->SetInput(m_ResampleVectorImage);
+
+ // Configure writer
+ m_VectorWriter = VectorWriterType::New();
+ m_VectorWriter->SetFileName(ossFileName.str().c_str());
+ m_VectorWriter->SetInput(m_VectorImageExtractROIFilter->GetOutput());
+ m_VectorWriter->WriteGeomFileOn();
+ m_VectorWriter->Update();
+
+ /** TODO : Generate KML for this tile */
+ // Search Lat/Lon box
+ m_Transform = TransformType::New();
+ m_Transform->SetInputKeywordList(m_ResampleVectorImage->GetImageKeywordlist());
+ m_Transform->SetInputProjectionRef(m_VectorImage->GetProjectionRef());
+ m_Transform->SetOutputProjectionRef(wgsRef);
+ m_Transform->InstanciateTransform();
+
+ InputPointType inputPoint;
+ OutputPointType outputPoint;
+ IndexType indexTile;
+ SizeType sizeTile, demiSizeTile;
+
+ sizeTile = extractSize;
+
+ /** GX LAT LON **/
+ // Compute lower left corner
+ indexTile[0] = extractIndex[0];
+ indexTile[1] = extractIndex[1] + sizeTile[1];
+ m_ResampleVectorImage->TransformIndexToPhysicalPoint(indexTile, inputPoint);
+ outputPoint = m_Transform->TransformPoint(inputPoint);
+ OutputPointType lowerLeftCorner = outputPoint;
+
+ // Compute lower right corner
+ indexTile[0] = extractIndex[0] + sizeTile[0];
+ indexTile[1] = extractIndex[1] + sizeTile[1];
+ m_ResampleVectorImage->TransformIndexToPhysicalPoint(indexTile, inputPoint);
+ outputPoint = m_Transform->TransformPoint(inputPoint);
+ OutputPointType lowerRightCorner = outputPoint;
+
+ // Compute upper right corner
+ indexTile[0] = extractIndex[0] + sizeTile[0];
+ indexTile[1] = extractIndex[1];
+ m_ResampleVectorImage->TransformIndexToPhysicalPoint(indexTile, inputPoint);
+ outputPoint = m_Transform->TransformPoint(inputPoint);
+ OutputPointType upperRightCorner = outputPoint;
+
+ // Compute upper left corner
+ indexTile[0] = extractIndex[0];
+ indexTile[1] = extractIndex[1];
+ m_ResampleVectorImage->TransformIndexToPhysicalPoint(indexTile, inputPoint);
+ outputPoint = m_Transform->TransformPoint(inputPoint);
+ OutputPointType upperLeftCorner = outputPoint;
+
+ // Build The indexTile
+ this->AddBBoxToIndexTile(lowerLeftCorner,
+ lowerRightCorner,
+ upperRightCorner,
+ upperLeftCorner,x,y);
+
+ /** END GX LAT LON */
+ y++;
+ }
+ x++;
+ y = 0;
+ }
+ }
+}
+
+/**
+*/
+
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::AddBBoxToIndexTile(OutputPointType lowerLeftCorner,
+ OutputPointType lowerRightCorner,
+ OutputPointType upperRightCorner,
+ OutputPointType upperLeftCorner,
+ unsigned int x, unsigned int y)
+{
+ // From PointType to VertexType
+ VertexType pLL,pLR,pUR,pUL;
+
+ pLL[0]=lowerLeftCorner[0];
+ pLL[1]=lowerLeftCorner[1];
+
+ pLR[0]=lowerRightCorner[0];
+ pLR[1]=lowerRightCorner[1];
+
+ pUR[0]=upperRightCorner[0];
+ pUR[1]=upperRightCorner[1];
+
+ pUL[0]=upperLeftCorner[0];
+ pUL[1]=upperLeftCorner[1];
+
+ // Form a polygon with the vertices
+ PolygonType::Pointer poly = PolygonType::New();
+ poly->AddVertex(pLL);
+ poly->AddVertex(pLR);
+ poly->AddVertex(pUR);
+ poly->AddVertex(pUL);
+
+ // Add the polygon to the datanode
+ m_Polygon = DataNodeType::New();
+ m_Polygon->SetNodeId("FEATURE_POLYGON");
+ m_Polygon->SetNodeType(otb::FEATURE_POLYGON);
+ m_Polygon->SetPolygonExteriorRing(poly);
+
+ itk::OStringStream oss;
+ oss << "tiles/tile_";
+ oss <<m_CurrentDepth <<"_"<< x <<"_"<< y <<".tif";
+
+ // Add the field "LOCATION" used in mapserver clients
+ // to get the path to tiles
+ m_Polygon->SetFieldAsString("LOCATION",oss.str());
+
+ // Add the to vectordata
+ m_VectorDataIndexTile->GetDataTree()->Add(m_Polygon,m_Folder);
+}
+
+/**
+ */
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::GenerateMapFile()
+{
+ std::ofstream file(m_FileName.c_str());
+ file << fixed << setprecision(6);
+
+ file <<"MAP" << std::endl;
+ file <<" NAME Level0" << std::endl;
+ file <<" # Map image size" << std::endl;
+ file <<" SIZE "<< m_TileSize <<" "<< m_TileSize << std::endl;
+ file <<" UNITS dd" << std::endl<<std::endl;
+ file <<" EXTENT -180 -90 180 90" << std::endl;
+ file <<" PROJECTION" << std::endl;
+ file <<" \"init=epsg:4326\"" << std::endl;
+ file <<" END" << std::endl;
+
+ file <<" # Background color for the map canvas -- change as desired" << std::endl;
+ file <<" IMAGECOLOR 192 192 192" << std::endl;
+ file <<" IMAGEQUALITY 95" << std::endl;
+ file <<" IMAGETYPE PNG" << std::endl;
+ file <<" OUTPUTFORMAT" << std::endl;
+ file <<" NAME PNG" << std::endl;
+ file <<" DRIVER 'GD/PNG'" << std::endl;
+ file <<" MIMETYPE 'image/png'" << std::endl;
+ file <<" IMAGEMODE RGB" << std::endl;
+ file <<" FORMATOPTION INTERLACE=OFF" << std::endl;
+ file <<" EXTENSION 'png'" << std::endl;
+ file <<" END" << std::endl;
+
+
+ file <<" # Web interface definition. Only the template parameter" << std::endl;
+ file <<" # is required to display a map. See MapServer documentation" << std::endl;
+ file <<" WEB" << std::endl;
+ file <<" # Set IMAGEPATH to the path where MapServer should" << std::endl;
+ file <<" # write its output." << std::endl;
+ //file <<" #IMAGEPATH \'D:\OSGeo4W_bis2/tmp/ms_tmp/\'" << std::endl;
+
+ file <<" # Set IMAGEURL to the url that points to IMAGEPATH" << std::endl;
+ file <<" # as defined in your web server configuration" << std::endl;
+ file <<" #IMAGEURL '/ms_tmp/'" << std::endl;
+
+ file <<" # WMS server settings" << std::endl;
+ file <<" METADATA" << std::endl;
+ file <<" 'wms_title' 'Level0'" << std::endl;
+ // XXXX Set the path to this file
+ file <<" \'wms_onlineresource\' \'http://127.0.0.1/cgi-bin/mapserv.exe?map=D:\\PARTAGE\\testsmapserver2\\niveau1_raster.map&\'" << std::endl;
+ file <<" \'wms_srs\' \'EPSG:4326\'" << std::endl;
+ file <<" END" << std::endl;
+ file <<" END" << std::endl;
+
+ // Get the name of the layer
+ std::ostringstream tempIndexShapeName;
+ tempIndexShapeName << itksys::SystemTools::GetFilenameWithoutExtension(m_FileName);
+
+ file <<" LAYER" << std::endl;
+ file <<" NAME '"<<tempIndexShapeName.str()<<"'" << std::endl;
+ file <<" #GROUP 'earthsat'" << std::endl;
+ file <<" TYPE RASTER" << std::endl;
+ file <<" TILEITEM 'LOCATION'" << std::endl;
+ file <<" TILEINDEX \'"<<m_IndexShapeFileName<<"\'" << std::endl;
+ file <<" METADATA" << std::endl;
+ file <<" 'wms_title' 'earthsat'" << std::endl;
+ file <<" 'wms_name' 'earthsat'" << std::endl;
+ file <<" END" << std::endl;
+ file <<" PROCESSING \"RESAMPLE=AVERAGE\"" << std::endl;
+ file <<" STATUS OFF" << std::endl;
+ file <<" TRANSPARENCY 100" << std::endl;
+ file <<" PROJECTION" << std::endl;
+ file <<" \"init=epsg:4326\"" << std::endl;
+ file <<" END" << std::endl;
+ file <<" MINSCALE 250000" << std::endl;
+ file <<" END" << std::endl;
+ file <<"END" << std::endl;
+
+ file.close();
+}
+
+template <class TInputImage>
+void
+MapFileProductWriter<TInputImage>
+::PrintSelf(std::ostream& os, itk::Indent indent) const
+{
+ // Call the superclass implementation
+ Superclass::PrintSelf(os, indent);
+}
+
+}
+#endif
--
GitLab