In [ ]:

!pip install geopandas

# install geopandas

In [ ]:

import geopandas
from google.colab import drive
drive.mount('/content/drive')

#impor the libraries I need

In [ ]:

data = geopandas.read_file('/content/drive/MyDrive/geospatial training/toxic_release_pennsylvania.shp')

#read in data to variable 'data'

In [ ]:

data.head()

# View data

Out[ ]:

	YEAR	CITY	COUNTY	ST	LATITUDE	LONGITUDE	CHEMICAL	UNIT_OF_ME	TOTAL_RELE	geometry
0	2016	PHILADELPHIA	PHILADELPHIA	PA	40.005901	-75.072103	FORMIC ACID	Pounds	0.160	POINT (2718560.227 256380.179)
1	2016	PHILADELPHIA	PHILADELPHIA	PA	39.920120	-75.146410	ETHYLENE GLYCOL	Pounds	13353.480	POINT (2698674.606 224522.905)
2	2016	PHILADELPHIA	PHILADELPHIA	PA	40.023880	-75.220450	CERTAIN GLYCOL ETHERS	Pounds	104.135	POINT (2676833.394 261701.856)
3	2016	PHILADELPHIA	PHILADELPHIA	PA	39.913540	-75.198890	LEAD COMPOUNDS	Pounds	1730.280	POINT (2684030.004 221697.388)
4	2016	PHILADELPHIA	PHILADELPHIA	PA	39.913540	-75.198890	BENZENE	Pounds	39863.290	POINT (2684030.004 221697.388)

In [ ]:

data['CHEMICAL'].unique()

# The types of chemicals available

Out[ ]:

array(['FORMIC ACID', 'ETHYLENE GLYCOL', 'CERTAIN GLYCOL ETHERS',
       'LEAD COMPOUNDS', 'BENZENE', 'ETHYLBENZENE',
       '1,2,4-TRIMETHYLBENZENE',
       'SULFURIC ACID (1994 AND AFTER ACID AEROSOLS" ONLY)"',
       'ZINC COMPOUNDS', 'BARIUM COMPOUNDS', 'ANTIMONY COMPOUNDS',
       'MERCURY COMPOUNDS', 'SILVER', 'PHTHALIC ANHYDRIDE',
       'DICYCLOPENTADIENE', 'AMMONIA', 'STYRENE', 'NICKEL',
       'TETRABROMOBISPHENOL A', 'TOLUENE', 'N-BUTYL ALCOHOL', 'COPPER',
       'PHENOL', 'XYLENE (MIXED ISOMERS)', 'CUMENE', 'METHANOL',
       'ETHYLENE', 'HYDROGEN SULFIDE', 'NICKEL COMPOUNDS', 'PROPYLENE',
       'ACETOPHENONE', 'TRADE SECRET CHEMICAL', 'MANGANESE',
       'CUMENE HYDROPEROXIDE', 'LEAD', 'METHYL METHACRYLATE', 'CHROMIUM',
       'BENZO(G,H,I)PERYLENE', 'COPPER COMPOUNDS', 'N-HEXANE',
       'ACETALDEHYDE', 'MALEIC ANHYDRIDE', 'FORMALDEHYDE',
       '2,2-BIS(BROMOMETHYL)-1,3-PROPANEDIOL', 'MANGANESE COMPOUNDS',
       'HYDROGEN CYANIDE', 'DICHLOROMETHANE', 'PROPIONALDEHYDE',
       'COBALT COMPOUNDS', 'NAPHTHALENE', 'CYCLOHEXANE',
       'HYDROGEN FLUORIDE', 'POLYCYCLIC AROMATIC COMPOUNDS',
       'ZINC (FUME OR DUST)', 'TETRACHLOROETHYLENE', 'TRICHLOROETHYLENE',
       'BENZOYL PEROXIDE', 'MERCURY', 'NITRIC ACID',
       'N,N-DIMETHYLFORMAMIDE', 'N-METHYL-2-PYRROLIDONE', 'OXYFLUORFEN',
       'POLYCHLORINATED BIPHENYLS', 'METHYL ISOBUTYL KETONE',
       'ASBESTOS (FRIABLE)', 'NITRATE COMPOUNDS',
       'DI(2-ETHYLHEXYL) PHTHALATE', 'DIETHANOLAMINE',
       'METHYL TERT-BUTYL ETHER', 'CHLORODIFLUOROMETHANE',
       'DIBUTYL PHTHALATE', 'DIOXIN AND DIOXIN-LIKE COMPOUNDS', 'MIXTURE',
       'VANADIUM COMPOUNDS', '1,3-BUTADIENE', 'METHYL ETHYL KETONE',
       'CYANIDE COMPOUNDS', 'DECABROMODIPHENYL OXIDE', 'SODIUM NITRITE',
       '1,2-BUTYLENE OXIDE',
       'CHROMIUM COMPOUNDS(EXCEPT CHROMITE ORE MINED IN THE TRANSVAAL REGION)',
       'TERT-BUTYL ALCOHOL', '1,2-DICHLOROETHANE',
       'TOLUENE DIISOCYANATE (MIXED ISOMERS)', 'CHLORENDIC ACID',
       'BUTYL ACRYLATE', 'CHLORINE', 'SEC-BUTYL ALCOHOL',
       'ALUMINUM (FUME OR DUST)', 'POLYCHLORINATED ALKANES',
       'LITHIUM CARBONATE', 'PHOSPHORIC ACID', 'BERYLLIUM',
       'HYDROCHLORIC ACID (1995 AND AFTER ACID AEROSOLS" ONLY)"',
       'COBALT', 'CHLOROMETHANE', 'CADMIUM', '2-ETHOXYETHANOL',
       'O-XYLENE', 'FREON 113', '1,1,1-TRICHLOROETHANE',
       'ISOBUTYRALDEHYDE', 'ETHYLENE OXIDE', '1,2-DICHLOROPROPANE',
       'M-XYLENE', 'BENZYL CHLORIDE', 'ARSENIC COMPOUNDS',
       'PROPYLENE OXIDE', 'CARBON TETRACHLORIDE',
       'DICHLORODIFLUOROMETHANE', 'P-XYLENE', 'ACETONE',
       '1,4-DICHLOROBENZENE', 'CHLOROMETHYL METHYL ETHER', 'BROMOMETHANE',
       'TOLUENE-2,4-DIISOCYANATE', 'AMMONIUM SULFATE (SOLUTION)',
       'CHLOROPHENOLS', 'EPICHLOROHYDRIN', 'BIS(CHLOROMETHYL) ETHER',
       'CHLOROFORM', 'BIPHENYL',
       'ISOPROPYL ALCOHOL (MANUFACTURING,STRONG-ACID PROCESS ONLY,NO SUPPLIER)',
       'CHLOROBENZENE', 'DIETHYL PHTHALATE', 'MOLYBDENUM TRIOXIDE',
       '1,2,4-TRICHLOROBENZENE', 'ACRYLIC ACID', 'N-DIOCTYL PHTHALATE',
       'ALUMINUM OXIDE (FIBROUS FORMS)', 'PHOSPHORUS (YELLOW OR WHITE)',
       'SODIUM HYDROXIDE (SOLUTION)', 'ANTIMONY', 'BARIUM',
       'TOLUENE-2,6-DIISOCYANATE', 'SODIUM SULFATE (SOLUTION)',
       '1,2-DIBROMOETHANE', 'VINYL CHLORIDE'], dtype=object)

From assessment this data is only on Pennsylvania, I will drop non-relevant columns

In [ ]:

data = data.drop(['CITY','COUNTY','ST'], axis = 1)

In [ ]:

data. head()

Out[ ]:

	YEAR	LATITUDE	LONGITUDE	CHEMICAL	UNIT_OF_ME	TOTAL_RELE	geometry
0	2016	40.005901	-75.072103	FORMIC ACID	Pounds	0.160	POINT (2718560.227 256380.179)
1	2016	39.920120	-75.146410	ETHYLENE GLYCOL	Pounds	13353.480	POINT (2698674.606 224522.905)
2	2016	40.023880	-75.220450	CERTAIN GLYCOL ETHERS	Pounds	104.135	POINT (2676833.394 261701.856)
3	2016	39.913540	-75.198890	LEAD COMPOUNDS	Pounds	1730.280	POINT (2684030.004 221697.388)
4	2016	39.913540	-75.198890	BENZENE	Pounds	39863.290	POINT (2684030.004 221697.388)

In [ ]:

data['geometry'].plot()

Out[ ]:

<Axes: >

In [ ]: