Postgis - Invalid Geometries in Cutouts

0

Objective: To calculate the area of shp1 that intersects (st_intersection) with shp2 in continental proportion.

The shp1 and shp2 are present information for the Brazilian biome of the whole Cerrado, that is, they are relatively large.

shp1 contains 15,139 records and shp2 contains 45,322 records.

I tried some ways to not execute st_intersection knowing that it is a function that demands a lot of computational processing.

My first attempt was to make a cut through Qgis 2.18.19, but the invalid geometries were not considered and I need to count them, they represent almost 90% of the map that I want to get the area.

I made a st_buffer and also a st_makevalid for the invalid geometries and did not succeed.

In addition to: merge / dissolve layers, use v.clean of GRASS, st_union, st_force2d, st_buildarea, st_force2d algorithms.

Can you help me find a way to calculate the area of this intersection so that my machine will support? It's possible?

I have an IntelCore i5-7200U and 8gb memory.

The files have spatial indexes created, I tried to transform multipolygons into polygons, options like using st_intersects and st_within in the WHERE field and I can not do it. I really do not know what else to do.

The files used are:

shp1 - link ?

shp2 - link ?

Thank you.

    
asked by anonymous 05.07.2018 / 17:28

0 answers