polygeohasher is a python package to implement polygon to geohash and vice versa with optimisation of geohash levels as per the user requirement, with error rate being controlled by the user.
Use the package manager pip to install polygeohasher.
pip3 install polygeohasherPoetry: Add package to lock file
poetry add polygeohasherfrom polygeohasher import polygeohasher
import geopandas as gpd
# read geojson(geometry) file
gdf = gpd.read_file("your geospatial file format") # read your geometry file here
# initialize polygeohasher
pgh = polygeohasher.Polygeohasher(gdf)
# declare geohash levels
INPUT_GEOHASH_LEVEL = 6
MINIMUM_GEOHASH_LEVEL = 5
MAXIMUM_GEOHASH_LEVEL = 7
# create a dataframe with list of geohashes for each geometry
initial_df = pgh.create_geohash_list(INPUT_GEOHASH_LEVEL,inner=False)
# get a dataframe with optimized list of geohashes
final_df = pgh.geohash_optimizer(initial_df, MINIMUM_GEOHASH_LEVEL, MAXIMUM_GEOHASH_LEVEL, INPUT_GEOHASH_LEVEL)
# prints optimization summary
pgh.optimization_summary(initial_df, final_df)
# convert geohash to geometry
geo_df = pgh.geohashes_to_geometry(final_df, "geohash_column_name")
# write file in desired spatial file format
geo_df.to_file("your write path.format",driver = "GeoJSON") Following is the optimization summary:
--------------------------------------------------
OPTIMIZATION SUMMARY
--------------------------------------------------
Total Counts of Initial Geohashes : 2597
Total Counts of Final Geohashes : 837
Percent of optimization : 67.77 %
--------------------------------------------------
Study are consist of division of City of Bengaluru in India.
Primary Output of geohashes without any optimisation.
Final Output of geohashes with optimization of number of geohashes at different levels to cover an area.
Pull requests are welcome. For major changes, please open an issue first to discuss what you would like to change.
Please make sure to update tests as appropriate.