def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Hospitals/FeatureServer/0" return (esri2gpd.get(url).to_crs(epsg=EPSG).pipe( geocode, ZIPCodes.get()).pipe(geocode, Neighborhoods.get()).pipe( geocode, PUMAs.get()))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Land_Use/FeatureServer/0" return ( esri2gpd.get(url).to_crs(epsg=EPSG) # .rename(columns={"c_dig1": "Type"}) )
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Vacant_Indicators_Bldg/FeatureServer/0" gdf = esri2gpd.get(url) return (gdf.to_crs(epsg=EPSG).pipe(geocode, ZIPCodes.get()).pipe( geocode, Neighborhoods.get()).pipe(geocode, PUMAs.get()))
def get_zip_codes(): """ZIP Codes in Philadelphia.""" return (esri2gpd.get( "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Philadelphia_ZCTA_2018/FeatureServer/0", fields=["zip_code"], ).to_crs(epsg=EPSG).rename(columns={"zip_code": "zip"}))
def get_police_districts(): """Police Districts in Philadelphia.""" return (esri2gpd.get( "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Boundaries_District/FeatureServer/0", fields=["DIST_NUM"], ).to_crs(epsg=EPSG).rename(columns={"DIST_NUM": "police"}))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/ArcGIS/rest/services/CITY_LANDMARKS/FeatureServer/0" return (esri2gpd.get( url, where="NAME = 'City Hall' AND FEAT_TYPE = 'Municipal Building'"). to_crs(epsg=EPSG).assign(x=lambda df: df.geometry.x, y=lambda df: df.geometry.y))
def get_council_districts(): """Council Districts in Philadelphia.""" return (esri2gpd.get( "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Council_Districts_2016/FeatureServer/0/", fields=["DISTRICT"], ).rename(columns={ "DISTRICT": "council" }).to_crs(epsg=EPSG))
def get_school_catchments(): """Elementary school catchments in in Philadelphia.""" return (esri2gpd.get( "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Philadelphia_Elementary_School_Catchments_SY_2019_2020/FeatureServer/0", fields=["name"], ).rename(columns={ "name": "school" }).to_crs(epsg=EPSG))
def get_pa_house_districts(): """Elementary school catchments in in Philadelphia.""" return (esri2gpd.get( "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/PA_House_Districts/FeatureServer/0", fields=["district"], ).rename(columns={ "district": "house_district" }).to_crs(epsg=EPSG))
def get_neighborhoods(): """Neighborhoods in Philadelphia.""" return (esri2gpd.get( "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Philly_NTAs/FeatureServer/0", fields=["neighborhood"], ).rename(columns={ "neighborhood": "hood" }).to_crs(epsg=EPSG))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/PPR_Assets/FeatureServer/0" return (esri2gpd.get( url, fields=["OBJECTID", "ASSET_NAME", "SITE_NAME", "ADDRESS" ]).to_crs(epsg=EPSG).pipe(geocode, ZIPCodes.get()).pipe( geocode, Neighborhoods.get()).pipe(geocode, PUMAs.get()))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Census_Tracts_2010/FeatureServer/0" return (esri2gpd.get(url, fields=["GEOID10", "NAME10" ]).rename(columns={ "GEOID10": "census_tract_id", "NAME10": "census_tract_name" }).to_crs(epsg=EPSG))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/City_Facilities_pub/FeatureServer/0" return (esri2gpd.get( url, where="ASSET_SUBT1_DESC = 'Library Branch'", fields=["ASSET_NAME" ]).to_crs(epsg=EPSG).rename(columns={ "ASSET_NAME": "asset_name" }).assign(x=lambda df: df.geometry.x, y=lambda df: df.geometry.y))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Schools/FeatureServer/0" return (esri2gpd.get(url, fields=[ "LOCATION_ID" ]).to_crs(epsg=EPSG).rename(columns={ "LOCATION_ID": "ulcs_code" }).dropna(subset=["ulcs_code"]).assign( ulcs_code=lambda df: df.ulcs_code.astype(int), x=lambda df: df.geometry.x, y=lambda df: df.geometry.y, ))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Commercial_Corridors/FeatureServer/0" gdf = esri2gpd.get(url) return (gdf.to_crs(epsg=EPSG).pipe(geocode, ZIPCodes.get(), use_centroids=True).pipe( geocode, Neighborhoods.get(), use_centroids=True).pipe( geocode, PUMAs.get(), use_centroids=True))
def download(cls, **kwargs): # Get the year YEAR = kwargs.get("year", DEFAULT_YEAR) if YEAR >= 2010: URL = "https://tigerweb.geo.census.gov/arcgis/rest/services/Census2010/tigerWMS_Census2010/MapServer/0" WHERE = "STATE=42 AND PUMA LIKE '%032%'" # trim to philadelphia elif YEAR >= 2000: URL = "https://tigerweb.geo.census.gov/arcgis/rest/services/Census2010/tigerWMS_Census2000/MapServer/0" WHERE = "STATE=42 AND PUMA LIKE '%041%'" # trim to philadelphia else: raise ValueError("'year' must be greater than or equal to 2000") return ( esri2gpd.get(URL, where=WHERE, fields=["GEOID", "NAME"]) .rename(columns={"GEOID": "geo_id", "NAME": "geo_name"}) .sort_values("geo_id") .reset_index(drop=True) .to_crs(epsg=EPSG) )
def download(cls, **kwargs): fields = [ "ASSET_NAME", "ASSET_ADDR", "SITE_NAME", "Copy_of_Master_Site_List_9_24_2", "Copy_of_Master_Site_List_9_24_5", "Copy_of_Master_Site_List_9_24_7", "Copy_of_Master_Site_List_9_24_8", "Copy_of_Master_Site_List_9_24_9", "Copy_of_Master_Site_List_9_2_10", "Copy_of_Master_Site_List_9_2_11", "Copy_of_Master_Site_List_9_2_13", ] url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Rebuild_Sites/FeatureServer/0" return ( esri2gpd.get(url, fields=fields) .to_crs(epsg=EPSG) .rename( columns={ "Copy_of_Master_Site_List_9_24_2": "SITE_TYPE", "Copy_of_Master_Site_List_9_24_5": "PRELIMINARY_SITE_NEEDS", "Copy_of_Master_Site_List_9_24_7": "ROUGH_EST_COST", "Copy_of_Master_Site_List_9_24_8": "PROJECT_STATUS", "Copy_of_Master_Site_List_9_24_9": "PROJECT_TYPE", "Copy_of_Master_Site_List_9_2_10": "COUNCIL_PRIORITY", "Copy_of_Master_Site_List_9_2_11": "PPR_PRIORITY", "Copy_of_Master_Site_List_9_2_13": "HIGH_NEED", } ) .pipe(geocode, ZIPCodes.get()) .pipe(geocode, Neighborhoods.get()) .pipe(geocode, PUMAs.get()) )
def download(cls, **kwargs): """ Download the census tract boundaries """ # Get the year YEAR = kwargs.get("year", DEFAULT_YEAR) # trim to PA (42) and Philadelphia County (101) WHERE = "STATE=42 AND COUNTY=101" if YEAR >= 2010: URL = "https://tigerweb.geo.census.gov/arcgis/rest/services/Census2010/tigerWMS_Census2010/MapServer/10" elif YEAR >= 2000: URL = "https://tigerweb.geo.census.gov/arcgis/rest/services/Census2010/tigerWMS_Census2000/MapServer/6" else: raise ValueError("'year' must be greater than or equal to 2000") return ( esri2gpd.get(URL, where=WHERE, fields=["GEOID", "NAME"]) .rename(columns={"GEOID": "geo_id", "NAME": "geo_name"}) .sort_values("geo_id") .reset_index(drop=True) .to_crs(epsg=EPSG) )
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/ArcGIS/rest/services/Philly_Neighborhoods/FeatureServer/0" return (esri2gpd.get(url, fields=[ "MAPNAME" ]).to_crs(epsg=EPSG).rename(columns={"MAPNAME": "neighborhood"}))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Philadelphia_ZCTA_2018/FeatureServer/0" return esri2gpd.get(url, fields=["zip_code"]).to_crs(epsg=EPSG)
# -*- coding: utf-8 -*- """ Created on Wed Apr 1 16:56:41 2020 @author: Jason """ import esri2gpd fields = ['kode_prov,fungsi,fungsi_kws,sk_kawasan,lindung'] url = "https://gis.wwf.id/server/rest/services/National/Kawasan_Hutan_Indonesia_KLHK_2014/FeatureServer/0" gdf = esri2gpd.get(url, fields=fields, where="kode_prov='18'") gdf.head() gdf.to_file("D:/z_Temp/klhk/kh_2014/kh_prov13.shp")
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/ArcGIS/rest/services/Universities_Colleges/FeatureServer/0" return (esri2gpd.get(url, fields=["NAME"]).to_crs(epsg=EPSG).assign( x=lambda df: df.geometry.x, y=lambda df: df.geometry.y))
def test_limit(): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Philadelphia_ZCTA_2018/FeatureServer/0" gdf = esri2gpd.get(url, fields=["zip_code"], limit=1) assert len(gdf) == 1
def test_where(): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Philadelphia_ZCTA_2018/FeatureServer/0" gdf = esri2gpd.get(url, fields=["zip_code"], where="zip_code=19123") assert (gdf["zip_code"] == 19123).all()
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/City_Limits/FeatureServer/0" return esri2gpd.get(url).to_crs(epsg=EPSG)
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/ArcGIS/rest/services/PPR_Assets/FeatureServer/0" return (esri2gpd.get(url).to_crs(epsg=EPSG).assign( x=lambda df: df.geometry.x, y=lambda df: df.geometry.y))
def download(cls, **kwargs): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Boundaries_District/FeatureServer/0" return (esri2gpd.get(url, fields=[ "DIST_NUM" ]).to_crs(epsg=EPSG).rename(columns={"DIST_NUM": "police_district"}))
def test_fields(): url = "https://services.arcgis.com/fLeGjb7u4uXqeF9q/arcgis/rest/services/Philadelphia_ZCTA_2018/FeatureServer/0" gdf = esri2gpd.get(url, fields=["zip_code"]) assert all(col in gdf.columns for col in ["zip_code", "geometry"])