Q&A
すみません、pytohnです。
shpを読み込んで、shp内のポリゴンを結合し、再び出力したいです。
以下のようなエラーが出ます。
AttributeError: 'MultiPolygon' object has no attribute 'to_file'
恐らく、shp内のポリゴンを結合したら、マルチポリゴンになり、to_fileでは出力できないという意味と思っています。
マルチポリゴンでも出力するにはどうすれば良いでしょうか。
コード
gdf_res = gpd.read_file(dir_input1/input_file, encoding='cp932')
gdf_res=gdf_res["geometry"].unary_union
gdf_res.to_file(dir_output/output_file,encoding='cp932')
print('Output: ',dir_output/output_file)
回答1件
0
ベストアンサー
python
1gdf_res1 = gpd.read_file(input_file, encoding='cp932') 2gdf_res2=gdf_res1["geometry"].unary_union 3gdf_res3=gpd.GeoDataFrame(geometry=[gdf_res2]) 4gdf_res3.to_file(output_file,encoding='cp932')
で行けるかと思います。
まず、gdf_res1はGeoDataFrameです。
gdf_res2はマルチポリゴンです。
to_fileはGeoDataFrameのメソッドなので、マルチポリゴンに対してでは処理できません。
マルチポリゴンをジオメトリの要素と持つGeoDataFrameを作って
それをto_fileする必要があります。そのGeoDataFrameがgdf_res3です。
投稿2022/01/17 14:22
総合スコア1497
回答ありがとうございます。
出力されて、結合もされていたのですが、gpd.read_file(input_file, encoding='cp932')が持っていた座標系が変わってしまいました。（元々は平面直角座標7系）
以下のように変更してみたのですが、合っていますでしょうか。
gdf_res1 = gpd.read_file(input_file, encoding='cp932')
gdf_res2=gdf_res1["geometry"].unary_union
gdf_res3=gpd.GeoDataFrame(geometry=[gdf_res2])
gdf_res3.crs = from_epsg(2449)
gdf_res3.to_file(output_file,encoding='cp932')
よろしくお願いします。
これもインポートします。
from fiona.crs import from_epsg
座標系が変わったと言うより無くなったのだと思います。
質問者さんの提示のコードでも出来るかも知れませんが、座標系をコピーする方がよいと思います。
gdf_res3.crs = gdf_res1.crs
という感じで。
ありがとうございます。
上手くできました。
ご協力ありがとうございました。
あなたの回答
tips
プレビュー
