[R-sig-Geo] bottom area of a lake calculation

[Thomas Adams]

Hello Sophie,

Yes, what you show was my understanding. I'm sure what you want to do can
be done in R, but I'm not sure what packages to use. For doing such a
thing, I use GRASS GIS. Regardless, I think you should re-project your data
from latitude-longitude to Albers Equal Area or Lambert Conic Conformal,
etc., that is area preserving. Eventually, you need to do this to get an
area anyhow.

An additional step would be to reproject, then spatially interpolate the
point locations you have in your spreadsheet into a grid to Albers Equal
Area or Lambert Conic Conformal, etc.

I could take a stab at this and show you...

Regards,
Tom

On Wed, Aug 8, 2018 at 10:39 AM, <sophie.leblanc using neuf.fr> wrote:

> Dear Tom,
>
> thanks for your reply.
> Just to be sure I have explain correctly.
> With marmap I can get the green area but I also want to calculate the red
> one.
> I have attached the kind of data I have.
>
> ‌
> I do not have any GIS software.
> I wanted to convert my spherical coordinate to cartesian ones with
>
> x = radius * cos(latitude) * cos(longitude)
> y = radius * cos(latitude) * sin(longitude)
> z = radius * sin(latitude)
> radius= earth_radius-depth : not fully sure...
> Do you think it is correct?
>
> The slope is not the same everywhere around the lake.
> Is it possible to calculate area with your first explanations?
> what R package should I use?
>
> thanks,
> Sophie Leblanc
>
>
>
>
>
>
>
> *http://sophieleblanc6.wix.com/photographie
> <http://sophieleblanc6.wix.com/photographie>*
> De : "Thomas Adams"
> A : sophie.leblanc using neuf.fr,"r-sig-Geo using r-project.org"
> Envoyé: mercredi 8 août 2018 16:07
> Objet : Re: [R-sig-Geo] bottom area of a lake calculation
>
>
> Sophie,
>
> My apologies, I was thinking in terms of using GRASS GIS... but the idea
> is the same
>
> Tom
>
> On Wed, Aug 8, 2018 at 10:06 AM, Thomas Adams <tea3rd using gmail.com> wrote:
>>
>> Sophie,
>>
>> I think a reasonable approximation might be obtained:
>>
>> (1) if you were to calculate the slope at each grid point using
>> r.slope.aspect
>> (2) use r.mapcalc to calculate A/cos(slope), where A is the pixel
>> projected area (e.g., 100mx100m = 10000 m^2) -- lat-long coordinates will
>> not work.
>>
>> Consequently, if the slope were approximately 0, cos(slope) approaches 1
>> and A/cos(slope) = A; as slope approaches 90 deg, cos(slope) becomes very
>> small, so A becomes very large. This works where the slope direction is
>> orthogonal to a pixel edge, so a correction to account for other slope
>> directions may be needed. But, with smaller grids, the error may be
>> negligible for you...
>>
>> Tom
>>
>>
>> On Wed, Aug 8, 2018 at 4:55 AM, <sophie.leblanc using neuf.fr> wrote:
>>>
>>> <span style="font-family:arial,helvetica,sans-serif;
>>> font-size:12px">‌</span>Dear all,<br>
>>> <br>
>>> I would like to calculate the bottom area of a lake. I have already used
>>> marmap to calculate projected area, i.e. surface area.<br>
>>> But i have not found anything to calculate a non projected area. I have
>>> the coordinates of points (longitude, latitude and depth) I can convert
>>> them into xyz. But then I don't know how to calculate the surface I have
>>> found polygonal area calculation but I have a 3D object, a polyheron and I
>>> think I can not use polygonal area.<br>
>>> Do you have any idea? do you know an R package that could help?<br>
>>> <br>
>>> thanks,<br>
>>> Sophie Leblanc<br>
>>> <br>
>>> <!-- +Signature --><font size="2"><em><a href="http://sophieleblanc6.wi
>>> x.com/photographie">http://sophieleblanc6.wix.com/photographie</a></em></font>
>>> <!-- -Signature -->
>>>         [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo using r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>>
>>
>>
>>
>>
>>
>
>
>
>
>
>



-- 
Thomas E Adams, III
1724 Sage Lane
Blacksburg, VA 24060
tea3rd using gmail.com (personal)
tea using terrapredictions.org (work)

1 (513) 739-9512 (cell)

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20180808/9115b88e/attachment.html>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: pasteImage1533738760024.png
Type: image/png
Size: 14034 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20180808/9115b88e/attachment.png>