Skip to content

Instantly share code, notes, and snippets.

View jachym's full-sized avatar

Jachym Cepicky jachym

View GitHub Profile
@jachym
jachym / jtsk.py
Created March 17, 2021 12:32
Výpočet rozlišení pro S-JTSK dlaždice 512px
try:
from tabulate import tabulate
except ModuleNotFoundError as e:
tabulate = None
print ("tabulate not installed, using `print`")
import json
dpi = 90
maxscaledenom = 7315200/2
@jachym
jachym / obce.py
Created October 29, 2020 08:09
Obce a koncovky
ice = 0
ov = 0
inn = 0
any = 0
with open("obce.csv") as obce:
for line in obce.readlines():
@jachym
jachym / basemaps-howto-cs.md
Last active July 9, 2020 09:28
MapServer Basemaps: instalace, publikace a vlastní styly

Cíle

  1. Publikovat mapovou službu nad daty OpenStreetMap, která bude publikovat tato data jako WMS v určitých stylech (Google, Default, BW, ...)
  2. Vytvořit vlastní styl, upravit, otestovat (např. v QGIS)

Linky

Příprava a stažení všeho potřebného

Gisquick (http://gisquick.org) je open source nástroj pro webovou prezentaci prostorových dat. Při jeho tvorbě náš slovensko-český tým zúročil dlouholeté zkušenosti při tvorbě webových mapových aplikací pro obce, zkušenosti s organizací vývoje otevřeného software a zkušenosti při nasazování těchto aplikací do produkčního prostředí. Gisquick umožňuje publikaci mapového projektu z prostředí QGIS na web – doslova na několik kliknutí. Protože na straně serveru je využit mapový server QGIS, jsou data na webu vykreslena přesně tak, jak byla nastylována v desktopovém prostředí. Součástí webové aplikace je pak vyhledávání prvků pomocí filtrů, měření, navigace, identifikace objektů a tisk mapy do PDF. Celé řešení ja navrženo tak, aby nedošlo k přílišnému zatížení serveru a prohlížeče – data jsou ukládána do vyrovnávací paměti, dlaždicována, atd. Celé prostředí je pak zabaleno do kontenjneru (Vagrant/Virtualbox) pro jeho rychlé nasazení v domácí infrastruktuře.

Celý projekt má tři části: zásuvný modul pro publikaci pr

@jachym
jachym / resources.json
Created March 14, 2018 12:24
melown resources klokan
[
{ "include": "examples.d/*.json" },
{
"comment" : "Klokan Satellite",
"credits" : [ ],
"definition" : {
"remoteUrl" : "https://maps.tilehosting.com/data/satellite/{loclod}/{locx}/{locy}.jpg?key=kSwrAdFeuIo6rD2Bm9dc"
},
"type" : "tms",
"driver" : "tms-raster-remote",
@jachym
jachym / zelenize2pokryti.py
Created September 16, 2017 23:08
Skript pro převod měření síly mobilního signálu na železnici z CSV do GeoJSON
#!/usr/bin/env python3
import csv
import os
import json
def getdata(value):
if value == 'x':
return None
else:
@jachym
jachym / dalnice2geojson.py
Created September 16, 2017 23:05
Skript pro převod pokrytí dálničí sítě mobilním signálem CSV do GeoJSON
#!/usr/bin/env python3
import csv
import os
import json
def get_features(filename):
with open(filename, newline='') as csvfile:
spamreader = csv.reader(csvfile, delimiter=';')
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
import rasterio as rio
from rasterio import features
from rasterio.rio.helpers import coords, write_features
import numpy as np
import json
from shapely.geometry import asShape
from shapely.geometry import mapping
def generate_jsons():
#!/usr/bin/env python3
"""
Script for updating data produced from LINZ
works on local 'linz' database
python3 update_temporal_table.py --csv ~/Data/linz/layer-804-changeset.csv --name nz-property-titles
nz-property-titles will be converted to schema and table name nz_property_titles.nz_property_titles
"""