সান গিলিজের রাস্টেরিও ব্যবহার করুন । এটা সহজেই সঙ্গে মিলিত হতে পারে Fiona, (পড়া এবং লেখা shapefiles) এবং সুষম একই লেখকের।
স্ক্রিপ্টে রাস্টেরিও_পলিজোনাইজ.পি
শুরু হয়
import rasterio
from rasterio.features import shapes
mask = None
with rasterio.drivers():
with rasterio.open('a_raster') as src:
image = src.read(1) # first band
results = (
{'properties': {'raster_val': v}, 'geometry': s}
for i, (s, v)
in enumerate(
shapes(image, mask=mask, transform=src.affine)))
ফলাফলটি জিওজেএসন বৈশিষ্ট্যগুলির একটি জেনারেটর
geoms = list(results)
# first feature
print geoms[0]
{'geometry': {'type': 'Polygon', 'coordinates': [[(202086.577, 90534.3504440678), (202086.577, 90498.96207), (202121.96537406777, 90498.96207), (202121.96537406777, 90534.3504440678), (202086.577, 90534.3504440678)]]}, 'properties': {'raster_val': 170.52000427246094}}
যে আপনি সুদৃশ্য জ্যামিতিতে রূপান্তর করতে পারেন
from shapely.geometry import shape
print shape(geoms[0]['geometry'])
POLYGON ((202086.577 90534.35044406779, 202086.577 90498.96206999999, 202121.9653740678 90498.96206999999, 202121.9653740678 90534.35044406779, 202086.577 90534.35044406779))
জিওপ্যান্ডাস ডেটাফ্রেম তৈরি করুন এবং স্থানিক যোগদান, প্লট করা, জিওজসন, ইএসআরআই শেফফাইল ইত্যাদি হিসাবে সংরক্ষণের কার্যকারিতা সহজেই সক্ষম করুন
geoms = list(results)
import geopandas as gp
gpd_polygonized_raster = gp.GeoDataFrame.from_features(geoms)