Commande LOAD DEM

[camillemonchicourt]

Permet d'automatiser le chargement du MNT dans la BDD PostGIS à partir d'un fichier. Voir https://github.com/GeotrekCE/Geotrek-admin/blob/master/geotrek/altimetry/management/commands/loaddem.py

[camillemonchicourt]

Pour le moment, on télécharge la BDalti de l'IGN France metropole à 250m (75m trop lourd sur toute la France), on insert le MNT raster dans la BDD, puis on le vectorise pour pouvoir l'indexer :

GeoNature/install_db.sh

Lines 244 to 266 in 1fe9ee4

	if $install_default_dem
	then
	echo "Insert default French DEM (IGN 250m BD alti)"
	echo "" &>> log/install_db.log
	echo "" &>> log/install_db.log
	echo "--------------------" &>> log/install_db.log
	echo "Insert default French DEM (IGN 250m BD alti)" &>> log/install_db.log
	echo "--------------------" &>> log/install_db.log
	echo "" &>> log/install_db.log
	if [ ! -f '/tmp/BDALTIV2_2-0_250M_ASC_LAMB93-IGN69_FRANCE_2017-06-21.zip' ]
	then
	wget --cache=off http://geonature.fr/data/ign/BDALTIV2_2-0_250M_ASC_LAMB93-IGN69_FRANCE_2017-06-21.zip -P /tmp
	else
	echo "/tmp/BDALTIV2_2-0_250M_ASC_LAMB93-IGN69_FRANCE_2017-06-21.zip already exist"
	fi
	unzip /tmp/BDALTIV2_2-0_250M_ASC_LAMB93-IGN69_FRANCE_2017-06-21.zip -d /tmp
	#gdalwarp -t_srs EPSG:$srid_local /tmp/BDALTIV2_250M_FXX_0098_7150_MNT_LAMB93_IGN69.asc /tmp/dem.tif &>> log/install_db.log
	export PGPASSWORD=$user_pg_pass;raster2pgsql -s $srid_local -c -C -I -M -d -t 100x100 /tmp/BDALTIV2_250M_FXX_0098_7150_MNT_LAMB93_IGN69.asc ref_geo.dem|psql -h $db_host -U $user_pg -d $db_name &>> log/install_db.log
	echo "Vectorisation of dem raster. This may take a few minutes..."
	sudo -n -u postgres -s psql -d $db_name -c "INSERT INTO ref_geo.dem_vector (geom, val) SELECT (ST_DumpAsPolygons(rast)).* FROM ref_geo.dem;" &>> log/install_db.log
	echo "Create dem vector spatial index. This may take a few minutes..."
	sudo -n -u postgres -s psql -d $db_name -c "CREATE INDEX index_dem_vector_geom ON ref_geo.dem_vector USING gist (geom);" &>> log/install_db.log
	fi

La BDalti 75m ou 250m IGN est aussi disponible pour les DOM-TOM : http://professionnels.ign.fr/bdalti#tab-3

Sinon au niveau mondial, l'agence aérospatial japonaise propose depuis peu un MNT mondial à 30m :
http://www.eorc.jaxa.jp/ALOS/en/aw3d30/index.htm (TIF en 4326)

J'en ai téléchargé et mosaïqué quelques bouts pour tester sur la Guadeloupe : http://geonature.fr/data/jaxa/

[camillemonchicourt]

En attendant d'avoir une commande pour faire ça, la méthode manuelle a été décrite dans la doc :
http://docs.geonature.fr/admin-manual.html#referentiel-geographique

[jpm-cbna]

Il existe aussi les données de la NASA aussi précises à 30m: https://lpdaac.usgs.gov/products/srtmgl1v003/
(Cliquer sur le bouton "Access data", en haut de page, pour accéder aux sites de téléchargement des données)
Un autre site propose de télécharger les données via une carte: http://dwtkns.com/srtm30m/

Les données japonaises semblent plus récente, peut être sont elles plus précises.