AVS30(5次メッシュ)を 3次メッシュ内で平均化したいなあと思案。
全国だと EXCEL では行数が足りないので、正攻法の ArcGISで処理。
時間がかかるので Python でも追っかけ処理。
shpを読んだ後の主な処理はこれだけ。
#AVS30=0(水域)を削除 gdf=gdf[gdf['5th_CODE'] != 0] #3次メッシュコード作成 gdf['3rd_CODE']=gdf['5th_CODE'].astype(np.int64)//100 #3次メッシュにおける平均を求める df=gdf.groupby('3rd_CODE').mean()
空間演算が必要ないので、Python の方が早く終わりました。
便利。
0 件のコメント:
コメントを投稿