FOSS4G Asia 2023でオンライン発表したGEOS-WASMの話題など

この記事は、 [FOSS4G Advent Calendar 2023](https://qiita.com/advent-calendar/2023/foss4g) の4日目の内容になります。(が、遅れてしまって申し訳ありません...:bow::sweat_drops:)

# はじめに

先週の2023/11/29(水)にFOSS4G Asia 2023で、今年の夏頃から取り組んでいるGEOS-WASMの話をオンラインで発表してきましたので、スライドの概要と補足説明、イベントの感想などを記載しようと思います。

<script async class="docswell-embed" src="https://www.docswell.com/assets/libs/docswell-embed/docswell-embed.min.js" data-src="https://www.docswell.com/slide/Z1JEL3/embed" data-aspect="0.5625"></script><div class="docswell-link"><a href="https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118">[FOSS4G Asia 2023] GEOS runs on web browsers - WebAssembly power for geospatial analysis by @sanak</a></div>

# スライドの概要

## アジェンダ
[![geos-wasm/3](https://img.docswell.com/page/LJLZ6XX17R.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p3)

冒頭の以下の項目が、プレゼンで最低限紹介したかった箇所、
* What is GEOS ?
* What is GEOS-WASM ?
* GEOS-WASM using library/app

続く以下の項目が、背景、経緯と特徴、
* Background
* How the project started
* Characteristics of GEOS-WASM

最後の以下の項目が、最近の話題とWebAssemblyの将来の可能性、謝辞となります。
* Recent related topics
* Future possibilities of WebAssembly
* Ackowledgements

## GEOSって何？
[![geos-wasm/4](https://img.docswell.com/page/47MV5LL57W.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p4)
FOSS4G界隈の方はご存知だと思いますが、初見の方もいるかもしれなかったので、[GEOS](https://libgeos.org/)の簡単な説明をしました。
* Webサイト: https://libgeos.org/
* 計算幾何学のためのC/C++ライブラリ
* OGCの[シンプルフィーチャ](https://en.wikipedia.org/wiki/Simple_Features)のジオメトリモデルと空間関数を実装
* [PostGIS](https://postgis.net/)、[QGIS](https://qgis.org/)、[GDAL](https://gdal.org/)、[Shapely](https://shapely.readthedocs.io/en/stable/project.html)などのバックエンドで利用されている

## GEOS-WASMって何？
[![geos-wasm/5](https://img.docswell.com/page/P7R23KKZJ9.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p5)
このプレゼンのメインテーマ、[GEOS-WASM](https://chrispahm.github.io/geos-wasm/) の紹介になります。
* Webサイト: https://chrispahm.github.io/geos-wasm/
* [GEOS](https://libgeos.org/)ライブラリのWebAssemblyビルドで、ブラウザやNode.js環境で利用可能
* GEOSの低水準のC-APIを、抽象化なしでそのままJavaScriptに公開
* 高レベルのライブラリのビルディングブロックや、特定のユースケースでの利用を意図

Webサイトで左側メニューの **Examples** から [Buffer](https://chrispahm.github.io/geos-wasm/examples/buffer.html) を選択すると、実際に動作するサンプルが見れます。

<!-- [![geos-wasm/6](https://img.docswell.com/page/PJXD4LL1EX.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p6) -->

インストールは簡単で、Node.jsなら
```bash
npm install geos-wasm
```
ブラウザなら
```html
<script type="module">
  import initGeosJs from 'https://cdn.skypack.dev/geos-wasm';
</script>
```

ただし、利用の際には、初期化されたGEOSオブジェクトを取得後、GEOSの[C-API](https://libgeos.org/usage/c_api/)と同様な呼び出しが必要になります。
```js
import initGeosJs from 'geos-wasm'
// initGeosJsはGEOSオブジェクトへのPromiseを返す
const geos = await initGeosJs()

// GEOSオブジェクトを利用してGEOS関数を呼び出す
// 例: ポリゴンの面積を取得する
// WKTリーダを作成
const reader = geos.GEOSWKTReader_create()
const wkt = 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))'
// WKT文字列を読み込んでGEOSジオメトリのポインタを取得
const geomPtr = geos.GEOSWKTReader_read(reader, wkt)
// 面積が書き込まれるポインター領域を作成
const areaPtr = geos.Module._malloc(8)
// ジオメトリの面積を計算
geos.GEOSArea(geomPtr, areaPtr)
// 面積をポインタからJavaScriptの数値に変換
const area = geos.Module.getValue(areaPtr, 'double')

console.log(area) // area = 1

// WKTリーダ、ジオメトリ、ポインタを解放
geos.GEOSWKTReader_destroy(reader)
geos.GEOSGeom_destroy(geomPtr)
geos.GEOSFree(areaPtr)
```

## GEOS-WASMを使ったライブラリ・アプリ

[![geos-wasm/7](https://img.docswell.com/page/3JKXZLLMED.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p7)
GEOS-WASMを使ったライブラリ・アプリの1つ目は [@neocarto](https://github.com/neocarto) さん開発の `geotoolbox` になります。
* GitHub: [neocarto/geotoolbox](https://github.com/neocarto/geotoolbox)
* デモサイト: https://observablehq.com/@neocartocnrs/hello-geotoolbox
* 様々な動作可能なサンプルが含まれている
* GEOS-WASMはv2.0からいくつかの関数(`Buffer`, `Clip`, `Densify` と `Union`)で利用されている

[![geos-wasm/8](https://img.docswell.com/page/LE35D332J5.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p8)
GEOS-WASMを使ったライブラリ・アプリの2つ目は拙作の `GEOS WASM Tester` になります。
* GitHub: [sanak/geos-wasm-tester](https://github.com/sanak/geos-wasm-tester)
* デモサイト: https://sanak.github.io/geos-wasm-tester/
* デスクトップ環境で動作する [JTS Test Builder](https://www.eclipse.org/community/eclipse_newsletter/2018/december/jts.php) を参考に
* GEOS-WASMはgit submoduleとして利用

どちらもデモサイトから動作を確認できます。

## 背景 - 2次元の地理空間解析ライブラリ

[![geos-wasm/10](https://img.docswell.com/page/V7PXNMVZE8.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p10)
背景の1つ目として、2次元の地理空間解析ライブラリを開発開始年順とバインディング、利用アプリ・ライブラリ・DBと合わせて記載しています。
GEOSはJTS Topology Suiteのバインディングですが、最も利用されているライブラリと言えると思います。
また、2010年頃までは、Java・C++などの枯れた言語による実装が主でしたが、2010年以降は、JavaScriptやRustなどの新しい言語での実装・バインディングが進んでいるように思います。
太字の箇所は以降のスライドで名前が出てきます。

## 背景 - WebAssembly

[![geos-wasm/11](https://img.docswell.com/page/2JVLK9ZM7Q.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p11)
WebAssemblyの概要と経緯について簡単に記載しています。
全般に、一つ前のスライドの2次元の地理空間解析ライブラリが出揃った後(2015年以降)にWebAssemblyが登場したイメージです。

[![geos-wasm/12](https://img.docswell.com/page/5EGXQZ4XJL.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p12)
FOSS4GイベントでのWebAssembly関連の発表をまとめています。
最初のPirminさんの[発表](https://blog.sourcepole.ch/assets/2019/wasm-foss4g19.pdf)と、最後のMichaelさんの[発表](https://boardend.github.io/foss4g-2022-webassembly4g/1)(太字箇所)は、スライド資料へのリンクになってますので、WebAssemblyに興味のある方は是非一読頂けると良いと思います。

[![geos-wasm/13](https://img.docswell.com/page/4JQ8KLG57P.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p13)
Geo系でWebAssemblyに対応したオープンソースをリストアップしています。
既にEmscripten (C/C++)、Rustを利用した様々なプロジェクトがありますが、Emscripten (C/C++)ベースのプロジェクトは既存ライブラリをWebAssemblyにポーティングしたもの、Rustベースのプロジェクトは新しい実装になると思います。
太字の箇所は以降のスライドで名前が出てきます。

## プロジェクトが始まった経緯

[![geos-wasm/14](https://img.docswell.com/page/K742PD5VJ1.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p14)
プロジェクトのアイデア自体は[@tmcw](https://github.com/tmcw)さんのブレインダンプ(https://github.com/tmcw/geos-wasm) から始まりました。
* TurfはWeb開発で最も快適
* ただし、ロバストネス(堅牢性)には、JTSからのポートであるJSTSのような、他の安定した実装が必要
* GEOSと同じ機能・バグが得られることが望ましい
* GEOS/WASMのアイデアを提起

[![geos-wasm/15](https://img.docswell.com/page/LJ1M9ZW4EG.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p15)
[@chrispahm](https://github.com/chrispahm)さんがリポジトリを見つけ、[概念実証(PoC)バージョン](https://github.com/chrispahm/geos-wasm/tree/8c2679536bd33205377b6359a2e7bdc03d69db14)をGEOSとEmscriptenを用いて実装し、Issue([geos-wasm - Is it worth the effort? · Issue #2 · tmcw/geos-wasm](https://github.com/tmcw/geos-wasm/issues/2)) を起票
* 完全に動作する `buffer` デモと極めて詳細なドキュメント
* コンパイルプロセスは [bugra9/gdal3.js](https://github.com/bugra9/gdal3.js) を参考
* シリアライゼーションパフォーマンスの問題について記載

[![geos-wasm/16](https://img.docswell.com/page/GJW3Q96ZJ2.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p16)
[@kylebarron](https://github.com/kylebarron)さんがブログ([Thoughts on GEOS in WebAssembly](https://kylebarron.dev/blog/geos-wasm))を投稿し、以下に言及
* Apache Parquetバインディングとの類似性
* WebAssemblyのパフォーマンスの概念の詳細
* シリアライゼーションパフォーマンスの問題の詳細
* GeoRustの将来の可能性
  * 後に [geoarrow/geoarrow-rs](https://github.com/geoarrow/geoarrow-rs) プロジェクトを開始し、
  * Observableで記事([Prototyping GeoRust + GeoArrow in WebAssembly / Kyle Barron | Observable](https://observablehq.com/@kylebarron/prototyping-georust-geoarrow-in-webassembly))を投稿

[![geos-wasm/17](https://img.docswell.com/page/4EZVW9YLJ3.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p17)
私はX(Twitter)経由で @kylebarron　さんのブログ記事を見つけ、このプロジェクトに興奮を覚えました。なぜなら、
* 過去、Windows+IE (ActiveX)環境向けに [GEOS COM Wrapper](https://github.com/sanak/geos-com) の開発経験があったが、動作可能な環境を広げたいと思っていた。
* 最近になって [node-geos](https://github.com/kashif/node-geos) へのコントリビューション(貢献)を試みたが、Node (v8)のC++APIが私には難しかったのと、NodeはサポートできてもWebブラウザはサポートできないので諦めていた。
* 上記から、このプロジェクトへのコントリビューションを開始

[![geos-wasm/18](https://img.docswell.com/page/Y765GKDMEV.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p18)
最初の課題はシリアライゼーションパフォーマンス問題の解消でした。
[Preferred data transfer between JS and WASM context → shouldn't it be WKB? · Issue #4 · chrispahm/geos-wasm](https://github.com/chrispahm/geos-wasm/issues/4)
* JavaScriptとGEOS-WASM間でジオメトリを渡す時に、テキストベースのフォーマット(WKTや文字列化されたGeoJSON)だと遅い問題
* WKBフォーマットを [@syncpoint/wkx](https://www.npmjs.com/package/@syncpoint/wkx) パッケージと合わせて利用することで、パフォーマンスを改善
* (後に、PR(プルリクエスト)[#19](https://github.com/chrispahm/geos-wasm/pull/19)で別のヘルパーメソッドも追加されています。)

[![geos-wasm/19](https://img.docswell.com/page/G75YVP3Q74.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p19)
2つ目の課題はAPI設計でした。
[API preferences? · Issue #5 · chrispahm/geos-wasm](https://github.com/chrispahm/geos-wasm/issues/5)
* 元々は、Turf.jsと互換性のある `buffer` のような高水準の関数と組み合わせる設計となっていた。
* 個人的には、ShapelyのようなAPIが好ましいと思っていたが、JavaScriptのガベージコレクションでのメモリ解放が難しそうだった...
* 議論の結果、低水準のGEOS C-API関数のみ提供する形に。

## GEOS-WASMの特徴

[![geos-wasm/20](https://img.docswell.com/page/9J2G36ZWJR.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p20)
* GEOSのバインディングの中では、EmscriptenによりC/C++コードをコンパイルする点がユニーク
  * GEOS C/C++ライブラリへのダイナミックリンクを利用する、他のShapelyなどのGEOSバインディングとは異なる
  * JTSのバインディング(JSTS、dart_jts、NTSなど)は、Javaのコードを各言語にポートしているので、GEOS-WASMはそちらに少し似ている。
  * GEOSのパッチバージョンアップ時は、WebAssemblyへのフルコンパイルが必要
  * ライセンスはGEOSと同様、LGPL 2.1
  * フルコンパイルのため、npmパッケージサイズはTurf.jsより大きく、JSTSと同じレベル
    * [geos-wasm](https://www.npmjs.com/package/geos-wasm): 2.98 MB
    * [jsts](https://www.npmjs.com/package/jsts): 3.42 MB
    * [@turf/turf](https://www.npmjs.com/package/@turf/turf): 658 kB

## 最近の関連トピック

[![geos-wasm/21](https://img.docswell.com/page/DEYV19R9EM.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p21)
WebAssembly関連では
* [The State of WebAssembly 2023](https://blog.scottlogic.com/2023/10/18/the-state-of-webassembly-2023.html)
* [WebAssembly Garbage Collection (WasmGC) enabled by default in Chrome](https://developer.chrome.com/en/blog/wasmgc/)
  * 関連したFlutterのニュース: [Support for WebAssembly (Wasm) | Flutter](https://docs.flutter.dev/platform-integration/web/wasm)

FOSS4GでのWebAssembly関連では
* [qgis/qgis-js](https://github.com/qgis/qgis-js) ([Demo](https://qgis.github.io/qgis-js/))
  * QGIS(のコア機能)がEmscriptenビルドでWebブラウザ上で動作するというもの
* SotM EU 2023 ワークショップ: [Ungap the Map: designing moonshot cycling networks](https://stateofthemap.eu/program/ungap-the-map-designing-moonshot-cycling-networks) by Dustin Carlino
  * カスタムのコスト関数を定義可能なクライアント側の経路探索 ([リンク](https://twitter.com/CarlinoDustin/status/1719393159560524040))

## WebAssemblyの将来の可能性

[![geos-wasm/22](https://img.docswell.com/page/VJNM2LDD78.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p22)
* オープンな3D建物データとクライアン側経路探索による日陰を通る経路探索など
  * Rustの既存のリソース ([dabreegster/route_snapper](https://github.com/dabreegster/route_snapper)) を利用
  * もしくは、 [pgRouting](https://github.com/pgRouting) のコアの C/C++ 箇所をEmscriptenでWebAssemblyにエクスポートするなど？
* 言語間の相互運用性の向上
  * Michael SchmukiさんのFOSS4G 2022の[WebAssembly4G: Where we are, and where we're heading](https://boardend.github.io/foss4g-2022-webassembly4g/20?clicks=1)のスライドの [p.20](https://boardend.github.io/foss4g-2022-webassembly4g/20?clicks=1) [p.37](https://boardend.github.io/foss4g-2022-webassembly4g/37?clicks=2) で言及されていた [Component Model](https://component-model.bytecodealliance.org/)

## 謝辞

[![geos-wasm/23](https://img.docswell.com/page/YE9N14G8E3.jpg)](https://www.docswell.com/s/sanak/Z1JEL3-2023-11-30-100118#p23)
最後に以下の方々に感謝を述べます。
* GEOS-WASMの著者とコントリビュータ(貢献者)へ
  * 彼らのアウトプット・努力がなければ、プロジェクトは始まらなかったので
* GEOSコミュニティへ
  * 20年を超える堅牢なジオメトリ演算機能の提供・維持を感謝して
* そしてスライド内で言及した他の全てのプロジェクトと方々に

# イベントの感想

## 事前準備

Call for Presentationが8月末に終わって、イベントまで時間は十分あったはずなのに、スライドを書き始めたのが発表の1週間前で、かつ、アジェンダ部分がなかなか整理できずに焦りばかりつのってしまう状況でした。
会社で希望者が利用できる英語のレッスン(基本は週1で30分)を再開し、イベントまで集中的に予約して、アジェンダの部分もアドバイス頂いたことで、何とかスライド作成・発表までこぎつけました... :sweat_drops:

## 当日発表

オンライントークは日本からは私だけで、Zoomで司会をされていた方は中国科学院大学のSongさんという方でした。
前後がコミュニティ・学術系の発表で、少し気後れ・緊張しましたが、所々英語発音をかみながらも、何とか無事済ませられました。

20分の発表で18分ほどで終わってしまい、少し時間が余ったので司会のSongさんから、「Extensionのように動作するのか？」との質問があったんですが、想定環境が良く分からずでうまく回答できず、しどろもどろになってしまいました... :sweat:

他にも、スライド資料でリンクを多用していたこともあって、ZoomのQAでリンクが分からないという質問がありましたので、発表終了後に一旦Googleスライドのリンクを共有しましたが、今後発表の際は、スライド資料自体のURLをQRコードなどで提供するようにします...。

ZoomのQAでは、プロジェクト参加の動機に関する質問もあり、「GEOSがWebブラウザ上で動作すること自体が面白いと思ったから」のように回答したと思います。

## オンライン視聴

1日目(11/29)はプレゼン発表の後、力尽きてしまったので、2日目(11/30)の午後からYouTubeでライブ配信のメインホールの動画を視聴していました。
https://www.youtube.com/watch?v=WCV90Fpt4Mg&t=14440s

3日目(12/1)のクロージングセレモニーも印象的でしたので、頭出ししたリンクを上げておきます。
https://www.youtube.com/watch?v=V1EKvxzhSbI&t=21029s

# おわりに

以上、少し冗長・単調になったかもしれませんが、GEOS-WASMのスライドの概要とイベントの感想を記載しました。

個人的な反省点としては、普段からアウトプットをしていないと、いざという時に言葉が出てこないことがあるように思いましたので、今後は色々とアウトプットを心がけていきたいと思います。 :bow::sweat_drops:

最後に、今年も残りわずかとなりましたが、今後とも、よろしくお願いいたします！