Aquabiota Solutions AB mapcloud

## csv-to-shapefile-geopandas.py
import pandas as pd
from geopandas import GeoDataFrame
from shapely.geometry import Point
import fiona

df = pd.read_csv('data.csv')

geometry = [Point(xy) for xy in zip(df.x, df.y)]
crs = {'init': 'epsg:2263'} #http://www.spatialreference.org/ref/epsg/2263/
geo_df = GeoDataFrame(df, crs=crs, geometry=geometry)

## remove-empty-columns.csv

          
            foo
            bar
            baz

            
              a
              
              1

            
              b
              
              2

            
              c

            
              
              
              4

            
              e
              
              5

## compassbearing.py
# LICENSE: public domain

def calculate_initial_compass_bearing(pointA, pointB):
    """
    Calculates the bearing between two points.

    The formulae used is the following:
        θ = atan2(sin(Δlong).cos(lat2),
                  cos(lat1).sin(lat2) − sin(lat1).cos(lat2).cos(Δlong))
	import pandas as pd
	from geopandas import GeoDataFrame
	from shapely.geometry import Point
	import fiona

	df = pd.read_csv('data.csv')

	geometry = [Point(xy) for xy in zip(df.x, df.y)]
	crs = {'init': 'epsg:2263'} #http://www.spatialreference.org/ref/epsg/2263/
	geo_df = GeoDataFrame(df, crs=crs, geometry=geometry)
	# LICENSE: public domain

	def calculate_initial_compass_bearing(pointA, pointB):
	"""
	Calculates the bearing between two points.

	The formulae used is the following:
	θ = atan2(sin(Δlong).cos(lat2),
	cos(lat1).sin(lat2) − sin(lat1).cos(lat2).cos(Δlong))