Automating Geo-processing tasks with GeoPandas and Python Scripting (Thiessen Weights Computation)

[@manishbhusal9299]

Thank you very much from Nepal🇳🇵

[@mteixeira87]

Brilliant, thank you for sharing this knowledge with all of us :D

[@geodeltalabs]

You are welcome

[@gustavoalvarezbea5855]

Great video, thank you very much! But wouldn't it be necessary to search for the area in the new 'district_areas' dataframe because it was already in 'intersect_poly' starting from the intersection?

[@RajanPudasainiz]

Love from Nepal

[@ShaibalAhmed]

This video is very well done and explained very well. Thank you! Subscribed! Can you please do a video where we can read in subsets of a large shapefile and save it to disk with predefined names? This would help me immensely since I am doing this manually with QGIS.

[@geodeltalabs]

Thanks for your comments and for sharing your interests. I will try to make a video to cover what you have requested :)

[@saii1990]

Very nice explanation.

[@geodeltalabs]

Thanks :)

[@VikasKumar-bh8kx]

hi, I have a CSV file filled with building footprints, which is quite heavy and can not be opened in the application directly, I want to clip my area of interest from that .csv file, how can I do that in any application? csv file size is 17gb

[@haniklakhe3434]

Interesting

[@mdmasroor3494]

Thank you so much. can you please tell me why you have divided area ten to the power six ?

[@geodeltalabs]

To convert the area into Sq. Km.

[@fulladehane6015]

Thanks

[@tubuobaggio]

hello how can i download the tutorial files? i have been trying to download the files to no avail thanks very much

[@vatsalyaanand8559]

Just a doubt isn't Thiessen Weights area of the Thiessen polygon divided by the whole area You divided Thiessen polygon with just the area of the district

[@geodeltalabs]

Hi, you may divide by the whole area if you would want to find out the contributions from the stations to the entire area, considered as one single unit (in this case, the whole country). But if you want to find out the weights by each district/ province/ catchment, then you have to divide it by the respective boundary area for each.

[@vatsalyaanand8559]

Thanks for clarification

[@lumafe1975]

I think you can simplify your code without creating a new dataframe in pandas........................

[@lumafe1975]

# Extracting the District and Area columns from the geodataframe to a new pandas dataframe # district_areas = pandas.DataFrame(districts[['District', 'Area']]) # district_areas.set_index('District', inplace = True) # Adding a new column to the intersect_poly attributes table intersect_poly['weights'] = 0 # Assign the thiessen weights to the new column for index_a, row_a in intersect_poly.iterrows(): # intersect_poly.loc[index_a, 'weights'] = row_a['Shape_Area']/(district_areas.loc[row_a['District']]['Area']) intersect_poly.loc[index_a, 'weights'] = row_a['Shape_Area']/row_a['Area']