Group Abstract Group Abstract

Message Boards Message Boards

How to estimate the area and volume of the dam reservoir from the satellite image?

Erfan Abdi
Posted 3 years ago

Hello all dear, I have the satellite image of Canton dam reservoir and I want to get its volume and area.

The actual specifications of the dam reservoir: volume=0.472km^3 & Surface area=32.01km^2

I did the following command to estimate the area, but the result was not acceptable and the unit is not clear. Please guide me to estimate the area and volume. Thank you

enter image description here

First, I manually obtained the shape below from the Mask Tool command.

enter image description here

imge =(image);

img = Binarize@imge~Erosion~1;
(m = MorphologicalComponents[img]) // Colorize

enter image description here

linecount = 
  3 /. (Binarize@imge // ColorNegate // Thinning // 
     ComponentMeasurements[#, "PerimeterCount"] &);
areacount = 
  2 /. ComponentMeasurements[MorphologicalComponents[img], "Count"];
areacount/linecount^2 // N

out put is : 0.111111
POSTED BY: Erfan Abdi
11 Replies
Sort By:

I realised that! Well at least the area is somewhat accurate. I'll try to improve the depth mechanism in sometime
POSTED BY: Navvye Anand
Erfan Abdi
Erfan Abdi
Posted 3 years ago

Thank you Navvye for share note book, but Volume of reservoir is 0.472km^3

POSTED BY: Erfan Abdi

Just multiply the area with the depth.

The entire code is uploaded as a nb
POSTED BY: Navvye Anand
Erfan Abdi
Erfan Abdi
Posted 3 years ago

Thanks for your efforts Navvye, is it possible to combined with (Gianluca command) and get the volume?
POSTED BY: Erfan Abdi
POSTED BY: Navvye Anand
Erfan Abdi
Erfan Abdi
Posted 3 years ago

Thanks, this is better. Can I get information from Mathematica about bathymetry to find depth and look for volume?
POSTED BY: Erfan Abdi

Great solution. I wonder if there's any way to recognise the dam automatically without manually using track-coordinates
POSTED BY: Navvye Anand
POSTED BY: Gianluca Gorni
Erfan Abdi
Erfan Abdi
Posted 3 years ago
POSTED BY: Erfan Abdi
POSTED BY: Navvye Anand

Your code amounts to this:

linecount = 3 /. {1 -> 587, 2 -> 76}
areacount = 2 /. {1 -> 11150, 2 -> 1, 3 -> 187734}

which does not make much sense to me.
POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard