I'm attempting to calculate acreage of the various "prime farmland" classifications in the GSSURGO dataset for the Omaha Reservation in Nebraska using QGIS but when I sum the total it's about 86 times more acreage than the total area of the reservation. The steps I followed to get to this calculation was:
- Load the GSSURGO MUPOLYGON layer into QGIS and reproject to EPSG 5070.
- Join with MAPUNIT Table (that contains prime farmland classifications for each soil unit in MUPOLYGON) using shared MUSYM field and export as joined shapefile.
- Clip joined GSSURGO data to include data within exterior boundaries of Omaha Reservation.
- Add new Decimal String field and Calculate acres for each polygon within the clipped layer using the Calculate Geometry polygon.
- Generate Summary Statistics using Vector>Analysis Tools>Basic Statistic for Fields.
The result it gives me is 17,0120,077 acres. However, the area of Omaha Reservation (calculated using boundary file used in clip operation in step 3) is only 198,395 acres. I've double checked I was calculating acres instead of square meters or square feet).
Why are these totals so wildly off?