Skip to content

Instantly share code, notes, and snippets.

What would you like to do?
computo indice di boscosita tramite spatialite - gis
create table indice_boscosita as
select COMUNE,CastToMultiPolygon(st_union(st_intersection (clc.geom, com.geom))) as geom, count (*)as nro,
(sum (st_area(st_intersection (clc.geom, com.geom))) / st_area(com.geom))*100.00 as indice
from CLC_Molise clc, comuni_molise com
where liv_3 like 'Boschi%' and st_intersects (clc.geom, com.geom) = 1
group by 1
order by indice desc
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment