I’m performing a carbon analysis for a province in Indonesia. I’m not sure whether my output is correct. Can anybody help me? I attach here the carbon pools excel file + the LULC shapefile I used + the output’s text document. I also attach a screen of my outputs in QGIS.
Furthermore, I have two QGIS-related questions:
for the InVEST analysis do I have to set a specific CRS?
How can I merge a LULC raster file with a vector file (roads)
Yes, you need to have your spatial inputs in the same projected coordinate system, where the units are in meters (see the User Guide Data Needs section for more information.) The model uses the area of each cell, along with the carbon pool values (which are Mg C/hectare) to calculate carbon storage. I am guessing that your land use map has a geographic coordinate system, where the units are in degrees, which is causing those very small output values.
There are probably several ways of merging rasters with vectors, but the way I just did the same thing is to turn my vector file into a raster, with the same cell size as the LULC, then use Raster Calculator to combine the two. In ArcGIS I use the Con tool for this, I’m not sure what the equivalent is in QGIS.
Since the log file doesn’t show anything wrong, and your carbon pool table looks fine (thanks for posting those!), it’s not obvious what’s happening. One thing is that I believe the underlying libraries don’t support Web Mercator, and I know yours isn’t exactly that, but I’d try reprojecting to a different coordinate system (UTM is a safe choice, and a lot of others work too) and see if that changes anything.
Thanks, Marco. I just ran InVEST 3.7 Carbon with the data you sent, and got the kind of values you’d expect (0 - 1476 MgC/pixel.) Everything looks good in terms of coordinate system and units in your inputs. Just in case, I looked at the result in both ArcGIS and QGIS and it looked fine in both, I am not seeing the e-07 values that you are.
So I’m not sure what’s happening on your end. Have you played with the raster symbology for tot_c_cur, in case it’s just a symbology issue? If you use the Identify tool and click around the raster, do the pixels have reasonable values? What’s the min/max in the raster properties?
Sorry, but I’m not sure what else to do to help since I can’t replicate the issue.
@swolny I also had a problem with my values. I assigned values for my land cover map with number from (1-25) and I assigned the carbon storage values to the land cover numbers. However, the result was all over the place. For example water was supposed to be zero but I got the values of open wetlands, the settlement values became zero. I tried to check my excel file but it was correct, and I did not get any error when I ran the model. I also checked my land cover map and it was assigned correctly