テストデータ
- 1 TIF ファイル(159KB)
目標 3 倍:
raster2pgsqlQGISを使用してPostGISにラスターをロードし、視覚化します- IPython Notebook で PostGIS に接続し、ラスターを NumPy 配列に読み込みます
- 私の IPyhton Notebook では、Pandas を使用して、PostGIS に保存されている異なる時間ステップでラスターの 1 ピクセルの時系列を読み込みます
これまでのプロセスraster2pgsqlコマンドを使用して 1 つのラスター イメージを PostGIS Raster に取得し、DB マネージャーを使用して QGIS で視覚化することができました。
raster2pgsql -s 4326 -d -I -C -M -R -l 4 D:\Downloads\raster//modis.tif -F -t 100x100 public.ndvi | psql -U postgres -d rastertest -h localhost -p 5432
しかし、IPython Notebook 内からこのラスターにアクセス/クエリする方法は?
SQLALchemy と GeoAlchemy2 に関するこのプレゼンテーションを見つけました。また、PostGIS Raster もサポートしていることが言及されています。とても面白そうです！しかし、ドキュメントを使用すると、これをラスター データに適用する方法がわかりません
postgres=user次のコードを使用して PostGIS データベースに接続できるpassword=adminと思いますdatabase=rastertest。
from sqlalchemy import create_engine
engine = create_engine('postgresql://postgres:admin@localhost/rastertest', echo=True)
しかし、それから..どんなアドバイスも大歓迎です。