cquest

# génération auto par gdal du fichier XML de description du service WMTS
gdal_translate "WMTS:https://wxs.ign.fr/an7nvfzojv5wa96dsga5nk8w/geoportail/wmts?Service=WMTS&Request=GetCapabilities&Version=1.0.0,layer=TRANSPORTS.DRONES.RESTRICTIONS" drones.xml -of WMTS

# téléchargement des tuiles et sortie en geotiff
gdal_translate drones.xml drones.tif -co COMPRESS=DEFLATE -co TILED=YES

cquest

                   events_what                    |    nb    |        dernier_cree        |      dernier_modifie       
--------------------------------------------------+----------+----------------------------+----------------------------
 public_transport.delay                           |   620959 | 2017-07-30 16:20:24.990843 | 2017-07-30 16:20:26.337225
 weather.warning.thunderstorm                     |    22168 | 2017-07-30 16:18:01.224928 | 2017-07-30 16:18:01.224928
 traffic.jam                                      |    91352 | 2017-07-30 16:17:06.586786 | 2017-07-30 16:21:06.260598
 traffic                                          |   243016 | 2017-07-30 16:17:04.457172 | 2017-07-30 16:21:06.744111
 traffic.accident                                 |    91085 | 2017-07-30 16:17:03.87314  | 2017-07-30 16:21:10.437644
 public_transport.delay.major                     |   105467 | 2017-07-30 16:15:14.104144 | 2017-07-30 16:20:26.681857
 traffic.nogaz                                    |    43494 | 2017-07

cquest

unzip erdas-ecw-sdk-5.3.0-linux.zip
chmod +x ERDAS_ECWJP2_SDK-5.3.0.bin
./ERDAS_ECWJP2_SDK-5.3.0.bin

sudo cp -r ~/hexagon/ERDAS-ECW_JPEG_2000_SDK-5.3.0/Desktop_Read-Only /usr/local/hexagon
ldconfig /usr/local/hexagon

cd gdal
CXXFLAGS="-D_GLIBCXX_USE_CXX11_ABI=0" ./configure --with-ecw=/usr/local/hexagon --without-jpeg12
make clean

cquest

wget http://download.geofabrik.de/europe/france/lorraine-latest.osm.pbf
wget "http://polygons.openstreetmap.fr/get_poly.py?id=51856&params=0" -O 54.poly
osmconvert lorraine-latest.osm.pbf -B=54.poly -o=54.pbf
osmconvert 54.pbf --out-osm | grep 'user=".*"' -o | sort | uniq -c | sort -n > stats.txt

cquest

# script permettant d'automatiser le téléchargement de la couche hydro de la BD Topo (licence Ouverte)
# les liens sur le site web IGN ne sont pas simplement téléchargeables (protection sur referer ou autre), d'où ce script
# ce script est sous licence WTFPL

MILLESIME=2017-03  # millesime à télécharger (figure dans les URL des fichiers)

# récup de la liste des URL sur la page http://professionnels.ign.fr/bdtopo-hydrographie
FILES=$(curl http://professionnels.ign.fr/bdtopo-hydrographie -s | sed 's/"/\n"\n/g' | grep "telechargement/inspire.*$MILLESIME")

for url in $FILES; do

cquest

cquest

cquest

psql -Atc "select insee from insee_cog_2015;" | \
	parallel -j 8 \
		wget -q -nc "https://www.geoportail-urbanisme.gouv.fr/document/info/?partition=DU_{}" -O {}.json

: > gpu.sjson; for f in *.json; do insee=$(echo $f |sed 's/.json//'); sed "s/\"success\"/\"insee\":\"$insee\", \"success\"/" $f >> gpu.sjson; done

jq -sc '.' gpu.sjson > gpu.json

cquest

select * from (
select case
  when id is null
  then format('<error class=\"33\" subclass=\"1\"><location lat=\"%s\" lon=\"%s\" /><text lang=\"fr\" value=\"%s (%s)\" /></error>',
    lat,lon,voie_cadastre,fantoir)
  when id_noname is not null and id_noname not like '%,%' and (l_geom-l_ways_noname<100) and ((l_noname > 0.5 and l2_noname<100) or (l_noname > 0.75)) and upper(voie_cadastre)!=voie_cadastre
  then format('<error class=\"32\" subclass=\"1\"><location lat=\"%s\" lon=\"%s\" /><text lang=\"fr\" value=\"%s (%s)\" /><way id=\"%s\"></way><fixes><fix><way id=\"%s\"><tag action=\"create\" k=\"name\" v=\"%s\" /></way></fix></fixes></error>',
    lat,lon,voie_cadastre,fantoir,id_noname,id_noname,voie_cadastre)
  when id is not null and id not like '%,%' and (l_geom-l_ways<100) and ((l > 0.5 and l2 < 100) or (l>0.75)) and names is null and upper(voie_cadastre)!=voie_cadastre
    then format('<error class=\"32\" subclass=\"1\"><location lat=\"%s\" lon=\"%s\" /><text lang=\"fr\" value=\"%s (%s)\" /><way id=\"%s\

cquest