[ferret_users] ferret isosurfaces in google earth with marching cubes algorithm

To: ferret users <ferret_users@xxxxxxxx>

Subject: [ferret_users] ferret isosurfaces in google earth with marching cubes algorithm

From: Patrick Brockmann <patrick.brockmann@xxxxxxxxxxxx>

Date: Mon, 4 Dec 2017 17:22:07 +0100 (CET)

Arc-authentication-results: i=3; mx.google.com; arc=pass (i=2); spf=softfail (google.com: domain of transitioning patrick.brockmann@xxxxxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=patrick.brockmann@xxxxxxxxxxxx

Arc-authentication-results: i=2; mx.google.com; arc=pass (i=1 spf=pass spfdomain=lsce.ipsl.fr); spf=softfail (google.com: domain of transitioning patrick.brockmann@xxxxxxxxxxxx does not designate 65.127.216.227 as permitted sender) smtp.mailfrom=patrick.brockmann@xxxxxxxxxxxx

Arc-authentication-results: i=1; mx.google.com; spf=pass (google.com: domain of patrick.brockmann@xxxxxxxxxxxx designates 194.254.241.32 as permitted sender) smtp.mailfrom=patrick.brockmann@xxxxxxxxxxxx

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=thread-topic:thread-index:mime-version:subject:to:from:date :arc-authentication-results:arc-message-signature :arc-authentication-results:arc-message-signature:message-id :arc-authentication-results; bh=7CgzG6omggEKGD96ytoKsmYo+TTvh9CK6LbOaqSnXpo=; b=XU2HL7ccJryB8Srn4ByktlW8kJmxarez0sBQxEMNWj5ZmuOcST4aVXIYj1YrovbSTn CT3RvZm07LRccGpjChTYVpnRFiHq1ueWqVLtF3GlA4TGhPtyM0VrNQJA3hfd9gBmoqjf QOtRibIFTYgJMbRtngg/n7ITHwucAsv9VL48YC5cONN0Bp7vqObGSoVLuI9EZp4rAQsD 4+4pFk3vL1c6SFgxxyxshLiAnMbIO5bgAr1EHP9FfWNo3nR5lJYyzrPBji9jl6XZKCdD MQGNOCwWHWsSO/VqVFwstM9N/pY7LmcklxRDYI8qndEKZ2G8wQJ8k/f+gc0pVKpotvO/ tTYg==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=thread-topic:thread-index:mime-version:subject:message-id:to:from :date:arc-authentication-results:arc-message-signature :arc-authentication-results; bh=7CgzG6omggEKGD96ytoKsmYo+TTvh9CK6LbOaqSnXpo=; b=sTef3qhyA8MJpgor8KRm6D/AtBArOIaA2NkzuvvKHwLAhkjPJyF32cVeCc84yCQJc/ PMxvWTFId3CB3LRttrxxStd1Y+GQ5XR/712YFIp7foZ+Yhr/rSDkOqYovF6nJ+kc5HZz kt5lGSklNHCJUlW1HqUByMuiJnmDcaG49bOUnjeG6toBhfy6KwnBGHKlZcHmR6TtBxYq +bTpp6tpikXQBbnlgjL4cZGCdPuZuQpR/1agzGvY6mHc5H5ol9AsN+Ivz5KIMfIdNXwS tEXCOU5S04KyM6qwZ6oF9DbJzBtdSd+ihlZ/M4ybpUQxI94P0qGa0Da8UftUIrZejniI kmWQ==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=thread-topic:thread-index:mime-version:subject:message-id:to:from :date:arc-authentication-results; bh=7CgzG6omggEKGD96ytoKsmYo+TTvh9CK6LbOaqSnXpo=; b=Bh4Yee1pMJT4bYgZLXxd49kdM8FPOd6yw+/KDA4bSO44soCry8DgeqOQV1L1ZwiHBs 9YR2ikLD2Q3AA53+FylbnWJ/F7OTh1DupBsBFJv4C2F5X+qwHUP0Hg7nit1J5LHdna8D r3O8yS4KV4cw7PH1r3cnTGMDEzAvmQfd2jLNMyCtJe7hcea23n9FqxcPDDMc4Yb0sZ3T OLum4Fk/U0Qw3ryQrb0adZ/5d/8wX2yX6ncZMUJuVXnp6WO/doCcaZmYOfiQy3RLNezK GmX6N+3d8tmB2hFnDk3Qn7Xq0twWiTKtZhjkhRJebdyfW5UWHIKGs4G1X/jR//myd1YV xMHg==

Arc-seal: i=3; a=rsa-sha256; t=1512404539; cv=pass; d=google.com; s=arc-20160816; b=iAgfk1HY26iJn7Lz7/07CeJHDryWqwrp/B+rn+V0THy1N5KmfnHyLettxvlGQHSnFY 7Ipsx+0DXW4tPYQfaAAxEQjqA1n4A6rW+aYYUyDgvYdDBOizwnmYDP7+JeCoZjXEy2j6 WCOKtAebvUpn1Mq+ur3PBfqllIsrRYAO8Py2fBP8+Cent1txao693zryMKpwCTSU4u0r 371FDQyLN482rrLnF/HpNJcCqRy/VCYliZRpJiUpdtxRTlE4hKgaw1Vf0bm0/2wUUlmG 4rj3ATmnS6d/PE0oVb85qvYjP8iIEqW+R8hscOolv5K1NunstYVL8/YiA1W6fcpBX/52 KUxw==

Arc-seal: i=2; a=rsa-sha256; t=1512404534; cv=pass; d=google.com; s=arc-20160816; b=Q4BunxvfACrGgEzLBZm3Z39blTK99UhxW/rsrBYZ9yjENYPbzFKbX3NgCNe91YKYTI QDAlBQ2/hU5p+SCkPQFCG2LZdRRYgEC/vAEBjvB2zzYkXYzx03olNr5DeVZ2j2AoisBN iMQAtLWWtahdziE1ImkILI5aCli0fFgHua+2DydCR0uu9vBEegybExvFIwQ2xm1xYPz2 6tpOmqmcYDdhKSLg+gIAd441oIrs1CmddhL4nsiORCVQLgOPzMpVv3/GhzSDMDQLk9Jf 0+IWfczlfKG8JNVxkYCjwsok/ZNj/Lf/WZqEuiEnwris5vXXA4PB55/V936rMn9NjkV4 EEJw==

Arc-seal: i=1; a=rsa-sha256; t=1512404528; cv=none; d=google.com; s=arc-20160816; b=xuK+nTy+LHJKRm24SvXfKqDr5j/86XGaPMfY/XJMXnD/OElxEmAgtjOFMQp24W/OiA j8+UMcAzdgQZZn7+KwsGBjXBviOujl6tyR5QvqnhGeb/P/P4SakXar8FK2kFQw8j4SQ+ moVl8oq++tMApplOLpc6gCeozDPVctnVYfpaeTd7nX43E8C9nWmHm+ktl4OUSyBWaWEK EdsWQcz9/y9F8ufS7O9ixfW/AuKzzNVFI3+fc0rxDmpcy+6sDr/xrynylyT5stECZkCG NwfZf6t1pyi4eMVbUlECe8CR7TZ7qbxbrAi6Y7XdozbHsZvU+nSm6mBitpShc0yHDwIZ dE/A==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa-gov.20150623.gappssmtp.com; s=20150623; h=message-id:date:from:to:subject:mime-version:thread-index :thread-topic:x-original-sender:x-original-authentication-results :precedence:mailing-list:list-id:list-post:list-help:list-archive; bh=7CgzG6omggEKGD96ytoKsmYo+TTvh9CK6LbOaqSnXpo=; b=C7spz6mAvoOh08uHFqrTHwanYNHgmd6H+ncc/Ov8QJSjXpu1YUmargDuOzZicv3a31 lzQyn+Vaim1TVLV0730V1dADhDxkn5NdFgukIeGINRh+KdQJ16bzo6gkzOlnHouiZo2u ggfwd/66H0QpZU4dISo0ztf4gZDyvnMiKaDjIS8pcoQODv3IxYubw5+I+5IXAxIlr3a9 U+0c66eO1+FL5x0gz/X5SAvjHK7QC5WoyZ8D6AxbzVTsWvmxgCu23KGMOnJidmqnkebb tK9zTkJRWU/DuKoZYmUNAEknI5ExSA12irtegF3Dsbxe83lCo0jxoye2bRzmLdK+iFWN efLw==

List-archive: <https://groups.google.com/a/noaa.gov/group/ferret_users/>

List-help: <https://support.google.com/a/noaa.gov/bin/topic.py?topic=25838>, <mailto:ferret_users+help@noaa.gov>

List-id: <ferret_users.noaa.gov>

List-post: <https://groups.google.com/a/noaa.gov/group/ferret_users/post>, <mailto:ferret_users@noaa.gov>

Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx

Sender: owner-ferret_users@xxxxxxxx

Thread-index: tdolCQl5jmroyOqXOicg9Zadkh5mQg==

Thread-topic: ferret isosurfaces in google earth with marching cubes algorithm

Hi ferreters,

I would be interested to share strategies or codes to plot from ferret (or pyferret) isosurfaces in google earth.

A colleague made a request on how to represent volcanic concentration isosurfaces nicelly (ie for outreach).

Also I have renewed my interest on Google Earth and have found a GooglePlot toolbox from matplot

that exactly provides what I would like to provide from ferret (https://publicwiki.deltares.nl/display/OET/Matlab+plotting+into+Google+Earth).

I have also found that the marching cube algorithm used to calculate isosurfaces

has been implemented from the scikit-image python module

(http://scikit-image.org/docs/dev/auto_examples/edges/plot_marching_cubes.html).

Now the kml export of polygons is more classic and should not be a difficult part.

I have found several related articles where this approach has been explored and implemented:

*Visualising volcanic gas plumes with virtual globes -2009

https://doi.org/10.1016/j.cageo.2009.02.005

* Isosurface Display of 3-D Scalar Fields from a Meteorological Model on Google Earth - 2013

http://www.dtic.mil/docs/citations/ADA582182

Any elements or share of existing initiatives on this topic will be welcome.

Regards

Patrick

--
Data Analysis and Visualization Engineer
LSCE/IPSL, CEA-CNRS-UVSQ laboratory
LSCE - Climate and Environment Sciences Laboratory
IPSL - Institut Pierre Simon Laplace
--