{
   OGRRegisterAll();

   OGRDataSource       *poDS;

   poDS = OGRSFDriverRegistrar::Open("C:\\projects\\osgearth\\data\\shapefiles_dresden", FALSE);
   if (poDS == NULL)
   {
      printf("Open failed.\n");
      exit(1);
   }

   OGRLayer  *poLayer;

   poLayer = poDS->GetLayerByName("gis.osm_pois_v06");

   OGRFeature *poFeature;

   OGRPolygon* polygon = new OGRPolygon();
   OGRLinearRing* ring = new OGRLinearRing();

   ring->addPoint(51.089028, 13.617821);
   ring->addPoint(51.089675, 13.807336);
   ring->addPoint(51.025585, 13.807679);
   ring->addPoint(51.024937, 13.616105);

   polygon->addRing(ring);

   poLayer->SetSpatialFilter(polygon);

   poLayer->ResetReading();
   while ((poFeature = poLayer->GetNextFeature()) != NULL)
   {
      OGRFeatureDefn *poFDefn = poLayer->GetLayerDefn();

      OGRFeature::DestroyFeature(poFeature);
   }

   OGRDataSource::DestroyDataSource(poDS);
}