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"])
