पॉलीगनों की एक geopandas GeoDataFrameश्रृंखला को देखते हुए , मैं अपनी सूची में प्रत्येक सुविधा के किमी वर्ग में क्षेत्र प्राप्त करना चाहूंगा।
यह एक बहुत ही आम समस्या है, और अतीत में सामान्य रूप से सुझाए गए समाधान का उपयोग करना shapelyऔर pyprojसीधे (जैसे यहां और यहां ) किया गया है।
क्या शुद्ध में ऐसा करने का कोई तरीका है geopandas?
जवाबों:
यदि GeoDataFrame के crs को जाना जाता है (EPSG: 4326 यूनिट = डिग्री, यहाँ), तो आपको अपनी स्क्रिप्ट में Shapely, और pyproj की आवश्यकता नहीं है क्योंकि GeoPandas उनका उपयोग करता है)।
import geopandas as gpd
test = gpd.read_file("test_wgs84.shp")
print test.crs
test.head(2)
अब अपने GeoDataFrame को कॉपी करें और प्रोजेक्शन सिस्टम में प्रोजेक्शन को बदलें (EPSG: 3857, unit = m जैसा कि ResMar के उत्तर में है)
tost = test.copy()
tost= tost.to_crs({'init': 'epsg:3857'})
print tost.crs
tost.head(2)
अब वर्ग किलोमीटर में क्षेत्र
tost["area"] = tost['geometry'].area/ 10**6
tost.head(2)
लेकिन मर्केटर प्रोजेक्शन में सतहें सही नहीं हैं, इसलिए मीटर में अन्य प्रोजेक्शन के साथ।
tost= tost.to_crs({'init': 'epsg:32633'})
tost["area"] = tost['geometry'].area/ 10**6
tost.head(2)
.to_crsसमारोह के लिए पारित हो जाता है pyprojवैसे भी। समान क्षेत्र प्रक्षेपण का एक अच्छा उदाहरण: proj4.org/projections/cea.html जिसे निम्नानुसार पारित किया जा सकता है:.to_crs({'proj':'cea'})
{'proj':'cea'}निकटतम क्षेत्र का अनुमान है।
मुझे विश्वास है हाँ। निम्नलिखित को काम करना चाहिए:
gdf['geometry'].to_crs({'init': 'epsg:3395'})\
.map(lambda p: p.area / 10**6)यह ज्यामिति को एक समान क्षेत्र के प्रक्षेपण में परिवर्तित करता है, क्षेत्र को लाता है shapely(मीटर ^ 2 में लौटता है), और नक्शे पर एक किमी ^ 2 (यह अंतिम चरण वैकल्पिक है)।
epsg:3395सीआरएस को फिट करने के लिए इस उत्तर को संशोधित किया है । धन्यवाद।
epsg:3857, लेकिन आपका कोड हैepsg:3395, दोनों में से कौन सही है?