Sunday, 15 March 2015

Coordinate system for accurately calculating areas of polygons that cross UTM Zones?


I have multiple polygons, all located offshore around the UK. I am trying to calculate the areas of these polygons in square kilometres, without splitting them by UTM zone.


Is it possible to do this accurately? I have been using World Mercator, but there are big distortions at this scale. If the polygons were limited to the east coast, say, I would use UTM31N. However, as they are all over the UK waters, I am unsure what to use.


Is there a coordinate system that covers this area? British National Grid does not extend far enough out to sea to be of use.



Answer



Any equal-area projection will do the job well. There are loads of equal-area projections that cover the entire earth (minus a point or two). Many of them are versions of a Cylindrical Equal-Area projection (such as the Gall-Peters).


You don't have to permanently reproject your polygons: create a temporary copy of the layer if you like, reproject it, calculate the areas, save the results as a table with just a pair of columns (feature id and area), and join the results to the original table.


(If you constantly project and unproject one master dataset, you risk accumulating small computational errors. One-off reprojections are often needed for analysis but shouldn't be used to modify the original data if you can help it.)


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...