Calculate ocean area using ETOPO

3 次查看（过去 30 天）
Ben
Ben 2022-7-9
回答： Anurag 2023-11-26
Hi,
I am wanting to calculate the area of a section of ocean (which also contains land) between two (lat,lon) points. I know I can import the ocean coverage and bathmetry in to Matlab from ETOPO, but unsure how to then get an area from this. Area would need to account for depth and bathymetry (i.e not just surface area). Is this even possible? If so I'd appreciate any pointers.
Thanks
  1 个评论
KSSV
KSSV 2022-7-9
How is your domain shape? If it is a perfect rectangle, simply use lengthxbreadth.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Anurag
Anurag 2023-11-26
Hi Ben,
I understand that you need to calculate the surface area of a section of the ocean between two (lat, lon) points also accounting for the depth and bathymetry. Considering both the surface area and the bathymetry, involves integrating over the three-dimensional space.
Follow the following steps to compute this:
[Z, R] = geotiffread('path_to_etopo_file.tif');
  • Define region of interest
latRange = [startLat, endLat];
lonRange = [startLon, endLon];
latIndices = find(R.LatitudeLimits(1) <= R.LatitudeLimits(2) & R.LatitudeLimits(1) <= latRange & latRange <= R.LatitudeLimits(2));
lonIndices = find(R.LongitudeLimits(1) <= R.LongitudeLimits(2) & R.LongitudeLimits(1) <= lonRange & lonRange <= R.LongitudeLimits(2));
subZ = Z(latIndices, lonIndices);
volume = trapz(trapz(subZ, 1), 2);
Please note that the above calculations have been done with the assumption that the region of interest is small enough to ignore the curvature of earth, else you may need to convert the latitudes and longitudes to Cartesian coordinates and accounting for the earth’s radius.
Hope this helps!
Regards,
Anurag

类别

在 Help Center 和 File Exchange 中查找有关 Oceanography and Hydrology 的更多信息

标签

产品


版本

R2022a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by