シェイプファイルからnetCDF4データを抽出しますか？

Question

python 3.6では、グローバルな降水値を含むnetCDF4ファイルをインポートしました。コロラド川流域の形状を含むシェイプファイルもインポートしました。私の目標は、シェイプファイル内の降水データのみを読み込み/抽出することです。私は複数の例を探しましたが、誰も本当に助けてくれませんでした。ここで

は、これまでの私のコードです：

from netCDF4 import Dataset 
import numpy as np 
import geopandas as gpd 

nc = Dataset('filename.nc') 
long = nc.variables['lon'][:] 
lati = nc.variables['lat'][:] 
rainfall = nc.variables['precip'][:] 

shapefile=gpd.read_file('filename.shp') 

は、上記のコード上のエラーメッセージはありません。

Answer 1

ああ、家の水文学者！ ;）

これまでのところ、コードをあまり使っていないので、ファイルをメモリに読み込んだだけです。私は（唯一のGRIBファイル付き）同じ分析を実行しようとしていたとき

、私は、まさにこのような目的のための偉大なPythonライブラリを見つけRasterStatsと呼ばれます。

これは、ndarrayラスタオブジェクトとGDALでサポートされているほとんどのラスタファイルタイプ（netCDFでもかまいません）の処理をサポートしており、必要なものを正確に生成します。より多くの場合

、非常に巧妙manualを見て、あなたがどこかで動けなくなるなら、私に知らせて！