The problem at hand includes 3 shapefiles. 2 Land Use classified satellite images and 1 regular polygon shapefile that covers the area of a National Park in my country. The software used in this analysis is QGIS 2.14.7. All shapefiles are using the SIRGAS 2000 CRS (EPSG: 4674)

My assignment was rather simple, I had to calculate the area for each class on the Land Use inside the park. The mold's area (the shapefile with the boundaries of the national park) is of 870.000 hectares.

I made sure both satellite images didn't overlap and merged them together. After that, I clipped the merged shapefile using the boundaries shapefiles as the mold. The output shapefile was created with the exact shape of the mold. Apparently everything was ok.

At last, I use the statistics for categories tool (located on the toolbox) to get a small table with the sums of each class. When I sum all the sums, the area goes to approximately 938.000 ha.

Any clues on why this might have happened?

I tried the topology checker plugin. I found 20 errors close to the limits of the satellite images, but they are small and wouldn't make up for almost 60.000 ha. Is it possible that just by the presence of errors, QGIS goes cuckoo?

I know I could dissolve the satellite shapes, but I'm not on a very powerful computer and the processing takes too long. Would it make a big difference anyway?

They're all on the same CRS, including the project.

  • You should immediately validate your polygons. Intersection code generally cheats and assumes certain topological relationships that are not true if either of the two polygons is ordered incorrectly. – Vince May 8 '17 at 23:27
  • I'm sorry, I'm not entirely sure I got that. I mean, I didn't use the intersection tool. Are you saying that one of those processes I listed before (I'm guessing Clip) utilises the same code of Intersection? – Eric Lino May 19 '17 at 20:02
  • Yes, a clip is a topological operation which would be impacted by incorrectly wound polygons. Just a single 4-vertex polygon could be responsible for an area error exceeding the envelope of the dataset (this is the essence of undefined behavior). – Vince May 20 '17 at 0:27

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.