How to calculate the area under a curve in km^2?
1 次查看(过去 30 天)
显示 更早的评论
Hello everyone!
I am a geologist working on this area of the Moon characterized by this fault-related landforms (see image), for some future calculations I need to know the so-called "excess area", so I need to quantify in m^2 or km^2 the area under the red line (this is a topographic profile across one of the above-mentioned structures) and above the black line whiich is my reference level. Since I am working on a GIS-based software with a known coordinate system I can measure all the distances I want to have a reference for calculation.
I read in a paper that "cross-sectional area of ridges was calculated from each profile by using a trapezoidal numerical integration" but I really don't know where to start...
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1390584/image.png)
Thank you in advance,
Jack
3 个评论
Torsten
2023-5-22
Then
trapz(x,y-ybase)
is the area where x and y are as above and ybase is the y-coordinate of the black line.
The unit of the area will be [unit_x]*[unit_y].
采纳的回答
Nathan Hardenberg
2023-5-22
编辑:Nathan Hardenberg
2023-5-22
Check out the trapz-function (https://de.mathworks.com/help/matlab/ref/trapz.html). This is exactly what you need.
You have to cut of the data from your starting point to your endpoint. And your "black line" has to be 0 level. because the trapz-function is just integrating numerically. If your data is not euqally spaced, you can also provide this in the function (see documentation)
% creating radom data
y = zeros(1,20);
y(5:14) = rand(1,10)+5
plot(y); % visualizing
xline(4, '--');
xline(15, '--');
y_cut = y(4:15); % cutting off data
area = trapz(y_cut)
0 个评论
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Geology 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!