willemvanopstal

process;region;aandeel
original;0;0.0017
original;1;0.0051
original;2;0.0034
original;3;0.0187
original;4;0.0192
original;5;0.0186
original;6;0.0448
original;7;0.0673
original;8;0.2201

willemvanopstal

# 2018-02-14 | 15:41
# Author: Willem van Opstal
# License: public

dummyDict = {1:[2,3,4],2:[],3:[5,8],4:[6,10],5:[],6:[7,9],7:[],8:[],9:[],10:[],11:[],12:[11]}
childrenList = []
idOfInterest = 1

def complexStuff(entitity):
	#This is just a dummy-function, you can neglect it if you want.

willemvanopstal

<?xml version="1.0"?>
<gpx version="1.1" creator="willem" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.topografix.com/GPX/1/1" xsi:schemaLocation="http://www.topografix.com/GPX/1/1 http://www.topografix.com/GPX/1/1/gpx.xsd">
<metadata><bounds minlat="51.326305272762710" minlon="3.292901827520408" maxlat="52.146054421000002" maxlon="4.727864003228013"/></metadata>                    

<wpt lat="51.37579504006839" lon="4.214571758064346">
  <name>-3.0m</name><sym>empty</sym>
</wpt>
<wpt lat="51.382056010186766" lon="4.125308333848052">
  <name>-2.95m</name><sym>empty</sym>
</wpt>

willemvanopstal

<?xml version="1.0"?>
<gpx version="1.1" creator="willem" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ogr="http://osgeo.org/gdal" xmlns="http://www.topografix.com/GPX/1/1" xsi:schemaLocation="http://www.topografix.com/GPX/1/1 http://www.topografix.com/GPX/1/1/gpx.xsd">
<metadata><bounds minlat="52.758011117000002" minlon="4.461771476959187" maxlat="53.860092322218186" maxlon="7.215335698745105"/></metadata>                    
<wpt lat="53.239276498061059" lon="7.202685131430939">
  <name>-2.5m</name>
  <sym>empty</sym>
</wpt>
<wpt lat="53.256001686818472" lon="7.175547327184482">
  <name>-2.45m</name>
  <sym>empty</sym>

willemvanopstal

#!/usr/bin/python
# Initialisierung
import RPi.GPIO as GPIO
import time
import os
GPIO.setmode(GPIO.BCM)
 
# Hier den entsprechden GPIO-PIN auswaehlen
GPIO_TPIN = 21
 

willemvanopstal

import fiona
import shapely.geometry as geometry
from shapely.ops import cascaded_union, polygonize
from scipy.spatial import Delaunay
import numpy as np
import math
        

def alpha_shape(points, alpha):
    if len(points) < 4:

willemvanopstal

999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 999.000 

willemvanopstal

import numpy as np
from scipy.interpolate import griddata
import matplotlib.pyplot as plt
import numpy.ma as ma
from numpy.random import uniform, seed

points = np.loadtxt('npcsv.csv',delimiter=";")
xi = np.linspace(589000,603500,1000)
yi = np.linspace(6965000,6980000,1000)

willemvanopstal

import gdal

gridopt = gdal.GridOptions(format='GTiff',algorithm='linear') # this results in one grey blob. whatever radius I try.
#gridopt = gdal.GridOptions(format='GTiff',algorithm='invdist:power=7') # this works perfectly

output = gdal.Grid('outputlinear.tif','newpoints.vrt',options=gridopt)

willemvanopstal

230817
250817