Le portail francophone de la géomatique

SELECT id_pts, the_geom, moy(valeur) as moy, max(valeur) as max FROM mesures WHERE substance LIKE 'C14' AND date_mesure > 'jan-2009' AND date_mesure < 'dec-2009'

   
    //---------------------------------------------------//
    // Create a feature source with point data.
    // (The Sheboygan sample data does not contain such data,
    // so we'll create it.)
    
    // Create a feature class definition for the new feature source
    $classDefinition = new MgClassDefinition();
    $classDefinition->SetName("Points");
    $classDefinition->SetDescription("Points");
    $classDefinition->SetDefaultGeometryPropertyName("GEOM");
    
    // Create an identify property
    $identityProperty = new MgDataPropertyDefinition("KEY");
    $identityProperty->SetDataType(MgPropertyType::Int32);
    $identityProperty->SetAutoGeneration(true);
    $identityProperty->SetReadOnly(true);    
    // Add the identity property to the class definition
    $classDefinition->GetIdentityProperties()->Add($identityProperty);
    $classDefinition->GetProperties()->Add($identityProperty);
    
    // Create a name property
    $nameProperty = new MgDataPropertyDefinition("NAME");
    $nameProperty->SetDataType(MgPropertyType::String);    
    // Add the name property to the class definition
    $classDefinition->GetProperties()->Add($nameProperty);
    
        // Create a measure property
    $measureProperty = new MgDataPropertyDefinition("MEASURE");
    $measureProperty->SetDataType(MgPropertyType::Double);    
    // Add the name property to the class definition
    $classDefinition->GetProperties()->Add($measureProperty);
    
    // Create a geometry property
    $geometryProperty = new MgGeometricPropertyDefinition("GEOM");
    $geometryProperty->SetGeometryTypes(MgFeatureGeometricType::Point);    
    // Add the geometry property to the class definition
    $classDefinition->GetProperties()->Add($geometryProperty);
    
    // Create a feature schema
    $featureSchema = new MgFeatureSchema("PointSchema", "Point schema");
    // Add the feature schema to the class definition
    $featureSchema->GetClasses()->Add($classDefinition);             
    
    // Create the feature source
    $featureSourceName = 'Library://Analyse_RSE/Points.FeatureSource'; 
    $resourceIdentifier = new MgResourceIdentifier($featureSourceName);
    $srs = $map->GetMapSRS();
    $sdfParams = new MgCreateSdfParams("ArbitraryXY", $srs, $featureSchema);
    $featureService->CreateFeatureSource($resourceIdentifier, $sdfParams);
    
    // We need to add some data to the sdf before using it.  The spatial context
    // reader must have an extent.
    $batchPropertyCollection = new MgBatchPropertyCollection();
    $wktReaderWriter = new MgWktReaderWriter();
    $agfReaderWriter = new MgAgfReaderWriter();
    $geometryFactory = new MgGeometryFactory();

    // Make four points     
    $propertyCollection = MakePoint("Point A", 5.2, "POINT(838005 6355540)");
    $batchPropertyCollection->Add($propertyCollection);
    unset($propertyCollection);    
    
    $propertyCollection = MakePoint( "Point B", 9.6, "POINT(838350 6355590)");
    $batchPropertyCollection->Add($propertyCollection);
    unset($propertyCollection);    
            
        
    // Add the batch property collection to the feature source
    $cmd = new MgInsertFeatures("Points", $batchPropertyCollection); 
    $featureCommandCollection = new MgFeatureCommandCollection();
    $featureCommandCollection->Add($cmd);
    
    // Execute the "add" commands
    $featureService->UpdateFeatures($resourceIdentifier, $featureCommandCollection, false);    
    
    // ... 
    //---------------------------------------------------//
    // Create a new layer - à partir de là c'est comme si le code n'était pas pris en compte
    
    $factory = new LayerDefinitionFactory();
    // Create a mark symbol
    $resourceId = "Library://Analyse_RSE/BasicSymbols.SymbolLibrary";
    $symbolName = "cube";
    $width = "24";  // unit = points
    $height = "24"; // unit = points
    $color = "FFFF0000";
    $markSymbol = $factory->CreateMarkSymbol($resourceId, 
      $symbolName, $width, $height, $color);   
    
    // Create a text symbol
    $text = "ID";
    $fontHeight="12";
    $foregroundColor = "FF000000";
    $textSymbol = $factory->CreateTextSymbol($text, 
      $fontHeight, $foregroundColor);

    // Create a point rule.
    $legendLabel = "trees";
    $filter = "";
    $pointRule = $factory->CreatePointRule($legendLabel, 
      $filter, $textSymbol, $markSymbol);
      
    // Create a point type style.
    $pointTypeStyle = $factory->
      CreatePointTypeStyle($pointRule);
      
    // Create a scale range.
    $minScale = "0";
    $maxScale = "1000000000000";
    $pointScaleRange = $factory->CreateScaleRange(
      $minScale, $maxScale, $pointTypeStyle);      

    // Create the layer definiton.
   $featureName = "PointSchema:Points";
    $geometry = "GEOM";
    $layerDefinition = $factory->
        CreateLayerDefinition($featureClass, $featureName,
        $geometry, $pointScaleRange);
        
    //---------------------------------------------------//
    // ...
    // Add the layer to the map - je n'ai pas ajout de cette couche 'Points' dans mon map..
    
    $newLayer = add_layer_definition_to_map($layerDefinition, "Points", "Points", $sessionId, $resourceService, $map);
    add_layer_to_group($newLayer, "Nouveau", "Nouveau", $map);

    // --------------------------------------------------//
    // Turn on the visibility of this layer.
    // (If the layer does not already exist in the map, it will be visible by default when it is added.
    // But if the user has already run this script, he or she may have set the layer to be invisible.)
    $layerCollection = $map->GetLayers();
    if ($layerCollection->Contains("Points"))
    {
        $pointsLayer =$layerCollection->GetItem("Points");
        $pointsLayer->SetVisible(true);
    }
    
    //---------------------------------------------------//
    //  Save the map back to the session repository
    $sessionIdName = "Session:$sessionId//$mapName.Map";
    $sessionResourceID = new MgResourceIdentifier($sessionIdName);
    $sessionResourceID->Validate();
    $map->Save($resourceService, $sessionResourceID);

    //---------------------------------------------------//
}
catch (MgException $e)
{
    echo "<script language=\"javascript\" type=\"text/javascript\"> \n";
    echo "    alert(\" " . $e->GetExceptionMessage() . " \"); \n";
    echo "</script> \n";
}

///////////////////////////////////////////////////////////////////////////////////
function MakePoint($name, $measure, $the_geom)
{    
    $propertyCollection = new MgPropertyCollection();
    $nameProperty = new MgStringProperty("NAME", $name);
    $propertyCollection->Add($nameProperty);
    
    $measureProperty = new MgDoubleProperty("MEASURE", $measure);
    $propertyCollection->Add($measureProperty);

    $wktReaderWriter = new MgWktReaderWriter();
    $agfReaderWriter = new MgAgfReaderWriter();    
    
    $geometry = $wktReaderWriter->Read($the_geom); 
    $geometryByteReader = $agfReaderWriter->Write($geometry); 
    $geometryProperty = new MgGeometryProperty("GEOM", $geometryByteReader);
    $propertyCollection->Add($geometryProperty);
    
    return $propertyCollection;
}
?>