[OSM-dev-fr] [SQL] calcul de surface

[Vincent Pottier]

On 06/10/2010 16:30, Vincent de Chateau-Thierry wrote:
> Bonjour,
>
>    
>> De : "Vincent Pottier"
>>
>> Pour calculer la surface couverte en landuse, j'utilise :
>>
>> SELECT count(*), SUM(ST_area( way)) AS area
>> FROM france_polygon
>> WHERE landuse LIKE '%'
>>
>> Ce qui donne :
>> 205217 44.5064532891517
>>
>> Mais je n'ai aucune idée de l'unité du résultat. Comment le traduire en
>> km² ?
>>      
> En effet, tu as obtenu des degrés carré, pas très pratique pour se faire une idée de la surface :-)
>
>    
>> De même :
>> SELECT osm_id, name, ROUND(ST_area( way))/1000000 AS area
>> FROM france_polygon
>> WHERE admin_level='2'
>>
>> me retourne systématiquement 0 pour les surfaces. Mais je peux trouver
>> la surface de la France métropolitaine sur wikipédia.
>>      
> Idem, la division par 1000000 des degrés carrés ramène vite à 0.
>
> Si tu es en PostGIS 1.5, tu peux pour obtenir "simplement" des mètres carrés, manipuler
> le type "Geography" à la place du type "Geometry". Plein d'explications ici :
> http://www.postgis.org/documentation/manual-1.5/ch04.html#PostGIS_Geography
>    
Je suis en 1.5. Ça marche avec EPSG : 3035

377100 km² de landuse. Je suis étonné de la faiblesse du résultat : ça 
fait un gros 50 % des 675 417 km² annoncés sur wikipédia pour la France 
métropolitaine (56%). Je m'attendais à plus.

On 06/10/2010 16:35, sly (sylvain letuffe) wrote:
>
>> De même :
>> SELECT osm_id, name, ROUND(ST_area( way))/1000000 AS area
>> FROM france_polygon
>> WHERE admin_level='2'
>>
>> me retourne systématiquement 0 pour les surfaces.
>>      
> T'es sûr que tu as quelque chose dans ta base ?
> select count(*) from france_polygon where admin_level='2' ?
>
>    
SELECT osm_id, name, ROUND(ST_area(st_transform(way,3035)))/1000000 AS area
FROM france_polygon
WHERE admin_level='2'

osm_id       name          area
-9407       "Andorra"    494.839945
-92863      "España"      12.813932
-1124039    "Monaco-territorial waters"    82.737834
-270009     "Guernsey"  3840.635954
-367988     "Jersey"    2404.757317
18887957    ""             0.892333
4543265     ""             0.004062

La base est le FranceLarge.osm de Pierren (de quelques jours...)