2020年12月17日木曜日

5次メッシュを3次メッシュに集約

AVS30(5次メッシュ)を 3次メッシュ内で平均化したいなあと思案。

全国だと EXCEL では行数が足りないので、正攻法の ArcGISで処理。

時間がかかるので Python でも追っかけ処理。
shpを読んだ後の主な処理はこれだけ。

  1. #AVS30=0(水域)を削除
  2. gdf=gdf[gdf['5th_CODE'] != 0]
  3. #3次メッシュコード作成
  4. gdf['3rd_CODE']=gdf['5th_CODE'].astype(np.int64)//100
  5. #3次メッシュにおける平均を求める
  6. df=gdf.groupby('3rd_CODE').mean()

空間演算が必要ないので、Python の方が早く終わりました。
便利。

0 件のコメント:

コメントを投稿