[Qgis-user] Getting landuse information for several circular polygons
Micha Silver
micha at arava.co.il
Mon Jul 16 12:13:06 PDT 2012
On 16/07/2012 20:22, Manuel Spínola wrote:
> Dear list members,
>
> I have a landuse layer and another layer with 20 circular polygons
> (buffers) and I want to obtain the amount of each landuse category
> within each circular polygon. I tried "clip" from Geoprocessing Tools
> but I dont get the information individualized for each circular
> polygon. Is it possible to get the landuse within each circular
> polygon in one step?
>
Not exactly "one step" but...
If you don't mind to drop the vectors into spatialite, then a query can
update the buffers layer with area of each landuse category as follows:
Suppose "buffers" has a label column, and there are landuse categories
lu1, lu2, lu3, .... then:
SELECT b.label AS Buffer,
ST_Area(ST_Intersection(b.Geometry, l.Geometry)) AS LandUse1
FROM buffers AS b, landuse AS l
WHERE l.category='lu1' AND
ST_Intersects(b.Geometry, l.Geometry)
GROUP BY b.buff_id;
If you add columns to the buffers table to hold total area for each
landuse category::
ALTER TABLE buffers ADD COLUMN lu1_area float;
ALTER TABLE buffers ADD COLUMN lu2_area float;
....
then you'd update each column like so:
UPDATE buffers SET lu1_area=(
SELECT ST_Area(ST_Intersection(b.Geometry, l.Geometry)) AS LandUse1
FROM buffers AS b, landuse AS l
WHERE l.category='lu1' AND
ST_Intersects(b.Geometry, l.Geometry) AND
b.label=buffers.label
);
HTH,
Micha
> Best,
>
> Manuel
>
> --
> *Manuel Spínola, Ph.D.*
> Instituto Internacional en Conservación y Manejo de Vida Silvestre
> Universidad Nacional
> Apartado 1350-3000
> Heredia
> COSTA RICA
> mspinola at una.ac.cr <mailto:mspinola at una.ac.cr>
> mspinola10 at gmail.com <mailto:mspinola10 at gmail.com>
> Teléfono: (506) 2277-3598
> Fax: (506) 2237-7036
> Personal website: Lobito de río
> <https://sites.google.com/site/lobitoderio/>
> Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>
>
> This mail was received via Mail-SeCure System.
>
>
> _______________________________________________
> Qgis-user mailing list
> Qgis-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/qgis-user
>
> This mail was received via Mail-SeCure System.
>
>
More information about the Qgis-user
mailing list