index.js 3.3 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117
  1. "use strict";
  2. Object.defineProperty(exports, "__esModule", { value: true });
  3. var meta_1 = require("@turf/meta");
  4. // Note: change RADIUS => earthRadius
  5. var RADIUS = 6378137;
  6. /**
  7. * Takes one or more features and returns their area in square meters.
  8. *
  9. * @name area
  10. * @param {GeoJSON} geojson input GeoJSON feature(s)
  11. * @returns {number} area in square meters
  12. * @example
  13. * var polygon = turf.polygon([[[125, -15], [113, -22], [154, -27], [144, -15], [125, -15]]]);
  14. *
  15. * var area = turf.area(polygon);
  16. *
  17. * //addToMap
  18. * var addToMap = [polygon]
  19. * polygon.properties.area = area
  20. */
  21. function area(geojson) {
  22. return meta_1.geomReduce(geojson, function (value, geom) {
  23. return value + calculateArea(geom);
  24. }, 0);
  25. }
  26. exports.default = area;
  27. /**
  28. * Calculate Area
  29. *
  30. * @private
  31. * @param {Geometry} geom GeoJSON Geometries
  32. * @returns {number} area
  33. */
  34. function calculateArea(geom) {
  35. var total = 0;
  36. var i;
  37. switch (geom.type) {
  38. case "Polygon":
  39. return polygonArea(geom.coordinates);
  40. case "MultiPolygon":
  41. for (i = 0; i < geom.coordinates.length; i++) {
  42. total += polygonArea(geom.coordinates[i]);
  43. }
  44. return total;
  45. case "Point":
  46. case "MultiPoint":
  47. case "LineString":
  48. case "MultiLineString":
  49. return 0;
  50. }
  51. return 0;
  52. }
  53. function polygonArea(coords) {
  54. var total = 0;
  55. if (coords && coords.length > 0) {
  56. total += Math.abs(ringArea(coords[0]));
  57. for (var i = 1; i < coords.length; i++) {
  58. total -= Math.abs(ringArea(coords[i]));
  59. }
  60. }
  61. return total;
  62. }
  63. /**
  64. * @private
  65. * Calculate the approximate area of the polygon were it projected onto the earth.
  66. * Note that this area will be positive if ring is oriented clockwise, otherwise it will be negative.
  67. *
  68. * Reference:
  69. * Robert. G. Chamberlain and William H. Duquette, "Some Algorithms for Polygons on a Sphere",
  70. * JPL Publication 07-03, Jet Propulsion
  71. * Laboratory, Pasadena, CA, June 2007 https://trs.jpl.nasa.gov/handle/2014/40409
  72. *
  73. * @param {Array<Array<number>>} coords Ring Coordinates
  74. * @returns {number} The approximate signed geodesic area of the polygon in square meters.
  75. */
  76. function ringArea(coords) {
  77. var p1;
  78. var p2;
  79. var p3;
  80. var lowerIndex;
  81. var middleIndex;
  82. var upperIndex;
  83. var i;
  84. var total = 0;
  85. var coordsLength = coords.length;
  86. if (coordsLength > 2) {
  87. for (i = 0; i < coordsLength; i++) {
  88. if (i === coordsLength - 2) {
  89. // i = N-2
  90. lowerIndex = coordsLength - 2;
  91. middleIndex = coordsLength - 1;
  92. upperIndex = 0;
  93. }
  94. else if (i === coordsLength - 1) {
  95. // i = N-1
  96. lowerIndex = coordsLength - 1;
  97. middleIndex = 0;
  98. upperIndex = 1;
  99. }
  100. else {
  101. // i = 0 to N-3
  102. lowerIndex = i;
  103. middleIndex = i + 1;
  104. upperIndex = i + 2;
  105. }
  106. p1 = coords[lowerIndex];
  107. p2 = coords[middleIndex];
  108. p3 = coords[upperIndex];
  109. total += (rad(p3[0]) - rad(p1[0])) * Math.sin(rad(p2[1]));
  110. }
  111. total = (total * RADIUS * RADIUS) / 2;
  112. }
  113. return total;
  114. }
  115. function rad(num) {
  116. return (num * Math.PI) / 180;
  117. }