Le portail francophone de la géomatique

// layer: le layer à copier, 
//destDs: la DataSource qui devra contenir le nouveau layer, 
//destSrs: la SpatialReference de ma nouvelle DataSource.
private static void InsertLayer(Layer layer, DataSource destDs, SpatialReference destSrs)
{
   // Création du nouveau Layer
   Layer destLayer = destDs.CreateLayer(layer.GetName(), destSrs, layer.GetGeomType(), null);
   if (destLayer == null)
      Trace.WriteLine("error CreateLayer");

   Feature tmpFeature = layer.GetNextFeature();
   if (tmpFeature == null)
      Trace.WriteLine("error no Feature");

   // Création des champs existants dans le nouveau Layer
   for (int i = 0; i < tmpFeature.GetFieldCount(); i++)
      if (destLayer.CreateField(tmpFeature.GetFieldDefnRef(i), 0) != Ogr.OGRERR_NONE)
         Trace.WriteLine("error");

   // Copie de chaque feature du layer à copier
   while (tmpFeature != null)
   {
      Feature newFeature = new Feature(destLayer.GetLayerDefn());
      newFeature.SetFrom(tmpFeature, 0);
      newFeature.SetFID(tmpFeature.GetFID());

      // On change la géométrie de la nouvelle feature (en rapport avec la SpatialRefercence)
      Geometry destGeometry = newFeature.GetGeometryRef();
      if (destGeometry != null)
         if (destGeometry.TransformTo(destSrs) != Ogr.OGRERR_NONE)
            Trace.WriteLine("error TransformTo");

      if (destLayer.SetFeature(newFeature) != Ogr.OGRERR_NONE)
         Trace.WriteLine("error SetFeature");
   }

   tmpFeature = layer.GetNextFeature();
}