tif_path,
'w',
driver='GTiff',
height=arr.shape[0],
width=arr.shape[1],
count=1,
dtype=arr.dtype,
crs=crs,
transform=transform,
nodata=-999.0
) as ds:
ds.write(arr, 1)
if __name__ == '__main__':
base_path = load_config()["base_path"]
try:
csv_path = sys.argv[1]
except:
# CSV input
csv_path = os.path.join(
base_path, 'results', 'proximity_results',
'population_gha_2019-07-01_proximity.csv')
try:
value_col = sys.argv[2]
except:
# Column
value_col = 'length_km'
"""Generate darthmouth events
"""
import os
import sys
import configparser
import pandas as pd
import geopandas as gpd
import subprocess
import numpy as np
from shapely.geometry import Point
import igraph as ig
import ghaa.config as config
import ghaa.analysis.spatial_and_network_functions as spnf
from tqdm import tqdm
BASE_PATH = config.load_config()['base_path']
tqdm.pandas()
def kind_of_access(x):
if x.nearest_hc_r_node == x.nearest_pop_r_node:
return 1
elif x.length_km > 0:
return 2
else:
return 0
def main():
crs = {'init': 'epsg:4326'} # Set the projection system
data_path = os.path.join(BASE_PATH,'data')
incoming_path = os.path.join(BASE_PATH,'incoming')
result_path = os.path.join(BASE_PATH,'results')