30 lines
740 B
Python
Executable File
30 lines
740 B
Python
Executable File
#!/usr/bin/python
|
|
from math import pi
|
|
from pyproj import CRS
|
|
import numpy as np
|
|
import geopandas as gpd
|
|
from shapely.geometry import LineString, MultiLineString
|
|
|
|
INTERVAL = 0.1
|
|
TAIL_LEN = 10
|
|
SINE_LEN = 7
|
|
|
|
TAILS = np.zeros(int(TAIL_LEN / INTERVAL))
|
|
|
|
|
|
def main():
|
|
sin_range = np.arange(-pi/4, SINE_LEN, INTERVAL)
|
|
amplitude = (np.sin(sin_range * pi / 2) + 1)*2
|
|
y = np.concatenate([TAILS, amplitude, TAILS])
|
|
x = np.arange(-TAIL_LEN - pi/4, SINE_LEN + TAIL_LEN, INTERVAL)
|
|
lines = LineString(zip(x, y))
|
|
geom = MultiLineString([lines])
|
|
df = gpd.GeoDataFrame(crs=CRS(3346))
|
|
df['geometry'] = None
|
|
df.loc[0, 'geometry'] = geom
|
|
df.to_file("sinewave.gpkg", driver='GPKG')
|
|
|
|
|
|
if __name__ == '__main__':
|
|
main()
|