tabulate area line 242

[Saadi4469]

Hello sir,

In line2 42 you are tabulating the area with respect to the DEM and the grid, my area units are coming off. I think that is because both my grid and DEM are in WGS84. I have attached a screenshot for further reference. Can you please guide me on this Dr. Rathore?

[Saadi4469]

My DEM is of 30 meter resolution and my GRID has a resolution of 0.0625 degree (approximately 6 km). Both of them are in WGS84 coordinates.

[LSRathore]

Hi, can you manually check if the sum of all the bands' areas equal to the total area of a grid?

[Saadi4469]

For calculating the area, I will have to change the projection of the grid to UTM right? And if that is so then for the tabulate area step, both the DEM and grid have to be in projected coordinate system. Please correct me if I am wrong!

[LSRathore]

Yes, you have to change the projection of the fishnet file to UTM to calculate the area of each pixel, you can do that in Arcgis and check if the output of line242 matches with what you generate. The script requires the shapefiles and rasters to be in GCS. I am not sure if it would work if inputs are in UTM projection. I do not have access to my data so I can not cross-check it, sorry.

[Saadi4469]

Hi Sir,
If I upload my data (a very small area) and share the link here, can you please have a look at it?

Thank you,

[LSRathore]

I can do it but not today, maybe tomorrow. Share the link with me and I’ll have a look.

[Saadi4469]

Thank you so much sir!
VIC_GRID_fishnet.zip
The DEM can be downloaded from the following link:
https://app.box.com/s/8qg5e8qif9r8hqr7han440zdw61mlxfn
Please let me know if you data access issues!

[Saadi4469]

Hi Sir, I figured out the solution but now I am getting an error in line 271. I will post a new issue for that.
Thank you for feedback!

[LSRathore]

That's good, I downloaded the data but could not run because my computer has an issue with arcpy.