gnn · January 19, 2021 12:18
diff --git a/rasterize.py b/rasterize.py
 from functools import partial
 import json

 from geocube.api.core import make_geocube
 from geocube.rasterize import rasterize_image
 from rasterio.enums import MergeAlg
 from shapely.wkt import loads
 from shapely.geometry import box, mapping, MultiPoint
 import pandas as pd
 import geopandas as gpd


 df = pd.read_csv(
    "testdata.csv",
    index_col=False,
    usecols=[3, 8],
    converters={8: loads},
 )
 gdf = gpd.GeoDataFrame(
    df,
    crs="epsg:4326",
    geometry=df.geometry,
 )
 lcc = gdf.to_crs("epsg:8763")
 grid = make_geocube(
    vector_data=lcc,
    measurements=["speed"],
    resolution=(-12000, 12000),
    rasterize_function=partial(rasterize_image, merge_alg=MergeAlg.add),
    fill=0,
    align=(9000, 8000),
 )
diff --git a/testdata.csv b/testdata.csv
	from functools import partial
	import json

	from geocube.api.core import make_geocube
	from geocube.rasterize import rasterize_image
	from rasterio.enums import MergeAlg
	from shapely.wkt import loads
	from shapely.geometry import box, mapping, MultiPoint
	import pandas as pd
	import geopandas as gpd


	df = pd.read_csv(
	"testdata.csv",
	index_col=False,
	usecols=[3, 8],
	converters={8: loads},
	)
	gdf = gpd.GeoDataFrame(
	df,
	crs="epsg:4326",
	geometry=df.geometry,
	)
	lcc = gdf.to_crs("epsg:8763")
	grid = make_geocube(
	vector_data=lcc,
	measurements=["speed"],
	resolution=(-12000, 12000),
	rasterize_function=partial(rasterize_image, merge_alg=MergeAlg.add),
	fill=0,
	align=(9000, 8000),
	)