defra lidar viewer

Données lidar de Defra

Un article de blog décrivant ce travail: http://www.ianhopkinson.org.uk/2015/10/analysing-lidar-data-for-the-uk/

Usage

process.py {OS_grid_cell} {friendly_name} affichera une fenêtre affichant les données d'un carreau de carte lidar. OS_grid_cell est, par exemple, SJ46 pour Chester. Le {Friendly_name} sera plus tard utilisé dans l'affichage des menus et est facultatif.

Les données sont supposées être dans un fichier zip avec le nom du formulaire:

LIDAR-DSM-2M-{OS_grid_cell}.zip

Ces données doivent être téléchargées manuellement à partir de:

http://environment.data.gov.uk/ds/survey#/download

Jusqu'à présent, j'ai utilisé avec succès le terrain de résolution 2M et les fichiers de surface.

Courir ./process.py , et après une pause short (<1 minute), vous devriez voir une image en niveaux de gris de votre carreau de carte avec une hauteur codée sous forme d'ombre (noir = plus bas, blanc = plus élevé). Cette image est écrite à pleine résolution dans les images/ répertoires avec le nom {OS_grid_cell}.png .

Nominalement, il y a 100 sous-tariles à un ensemble, mais certains sont manquants et apparaissent sous forme de carrés noirs, certains emplacements ont Nodata, ceux-ci semblent également noirs dans l'image.

Cela suppose que vous avez requis les bibliothèques installées:

pip install -r requirements.txt

Il y a des tests minimaux qui peuvent être effectués en faisant nosetests , en supposant que vous avez installé nosetests

Il existe actuellement deux visualisations expérimentales HTML / JavaScript. Pour les exécuter, exécutez le serveur Web minimal dans Python 3:

python -m http.server 8888 &

Puis accédez à localhost:8888 pour une visualisation de superposition de cartes leaflet.js , les images disponibles et leurs boîtes de délimitation seront récupérées dans le fichier data_dict.json qui est mise à jour lorsque process.py est exécuté.

Naviguez vers localhost:8888/surface.html pour voir un rendu de surface 3D Whizzy, modifiez simplement le nom de fichier d'image approprié sur surface.html . (J'ai fait des cultures de 512x512 pixels des images de taille complète à l'aide de peinture .NET.) Tensionner la souris gauche (ou A) et le déplacement de la souris tourne la vue, le bouton de souris droit (ou D) plonge la vue, au milieu -Mououse-Button (ou S) zooms la vue.

FAIRE

• Pour les données de résolution de 25 cm et 50 cm, ne rendent pas une tuile entière de 10 kmx10 km - rendez les carreaux individuels ou ne faites pas une boîte de délimitation serrée?
• Essayant de traiter un ensemble de données 50cm, nous obtenons une erreur de mémoire à faible mémoire, puis une mémoire de mémoire sur l'iMShow, 20000x20000 pixel PNG créé avec succès cependant
• Essayant de traiter un ensemble de données de 25 cm, nous obtenons des erreurs hors un dans Yoffset (et l'image est trop grande pour Imshow)
• Résoudre le problème de l'inadéquation avec la visualisation de folidlet.js, pourrait-il être le problème:

https://help.openstreetmap.org/questions/2056/using-the-ordnance-survey-stanal-grid-with-openstreetmap

Ou il se peut que le système d'exploitation utilise le Mercator transverse et OSM utilise une projection de mercator sphérique

• Téléchargement automatique des données
• Sélectionnez dans leaflet.js et allez à la vue de surface, cela semble pratique:

https://github.com/heyman/leaflet-areaselect/

• Traitement Spinner pour surface.html Démarrage initial
• Embellissement avec bootstrap

Données

Les données sur le terrain sont une estimation du terrain sous-jacent, les données de surface sont la mesure réelle, y compris les bâtiments, les véhicules et la végétation.

Les données vivent ici:

http://environment.data.gov.uk/ds/survey#/download

http://environment.data.gov.uk/ds/survey#/download?grid=SJ36

Il y a des métadonnées JSON pour chaque grille:

http://www.geostore.com/environment-agency/rest/product/os_gb_10km/SJ46

L' id dans le JSON vous indique où vous pouvez télécharger le fichier zip, par exemple:

http://www.geostore.com/environment-agency/rest/product/download/6129

Références

Collecte pratique de liens

http://mapgubbins.tumblr.com/post/131424021480/open-data-release-of-aerial-lidar-data-for-for

Cet outil convertit la référence de la carte du système d'exploitation x, y en lat, GNL (et les montre sur une carte)

http://gridreferencefinder.com/

Code pour convertir OS X, Y en lat, LNG est par Hannah Fry, décrit ici:

http://www.hannahfry.co.uk/blog/2012/02/01/converting-british-stanal-grid-to-latitude-and-ongitude-ii

Obtenez l'origine des carreaux de la référence de la grille nationale de l'OS:

http://digimap.edina.ac.uk/webhelp/os/gazetteer_plus/grid_ref_conversion.htm

Nous pourrions faire une superposition d'image en utilisant Leafet.js:

http://leafletjs.com/reference.html#imageoverlay

https://www.mapbox.com/mapbox.js/example/v1.0.0/imageoverlay-georeferend/

Rendant une surface en trois.js

http://www.smartjava.org/content/threjs-rider-real-world-terrain-heightmap-using-open-data

Ajout d'un contrôle interactif de trackball à trois.js

http://stackoverflow.com/questions/18347256/trackballControls-in- trois-js

Le contrôleur d'Opacity est soulevé à partir d'une démo de mapbox:

https://www.mapbox.com/mapbox.js/example/v1.0.0/opacity/