Discussion :

[Alaopi]

Bonjour,

Je suis un étudiant en classe prépa scientifique et dans le cadre de mon TIPE (projet scientifique mené sur un an et qui sera évalué aux concours) sur l'optimisation et le dimensionnement d'une centrale hydroélectrique en Martinique, j'ai été amené à soliciter votre aide.
Afin de créer un programme python d'optimisation du tracé des conduites d'eau forcées pour ma centrale hydroélectrique hypothétique, j'ai besoin des données sur l'altitude de cette zone. J'ai déjà consulté Géoportail où j'ai pu utiliser leur carte topographique pour repérer les dénivelés et même observer les variations d'altitude grâce à des profils altimétriques. Cependant, vu que pour mon programme, je compte modéliser la "zone de tracage" par un tableau de valeur où chaque indice correspondrait à une altitude (donc je subdiviserais le plan en "petits carrés"), il me serait profitable d'avoir accès directement aux relevés d'altitude avec lesquels ont été créés ces cartes. Je ne sais pas comment trouver un tel fichier, ou même s'il existe. J'ai cherché sur le site d'IGN pour obtenir plus de renseignements mais tout ceci m'a l'air très complexe...

Merci beaucoup de votre aide,

Alex Dejoie

[azarzelli]

Bonjour,

Ce que vous recherchez est un modèle numérique de terrain (MNT) : un raster (grille) qui à chaque case géoréférencée (dont on connaît les coordonnées) associe une altitude. https://fr.wikipedia.org/wiki/Mod%C3...que_de_terrain

L'IGN propose des MNT via ses produits RGE ALTI 5m (plus précis) et BD ALTI

Le MNT (RGE ALTI) de la Martinique est disponible en téléchargement par ftp via ftp://RGE_ALTI_ext:Thae5eerohsei8ve@..._2015-10-21.7z

Dans cette archive, dans le Dossier RGEALTI / 1_DONNES_LIVRAISON / RGEALTI_MNT_5M, vous retrouverez des tuiles couvrant l'emprise de la Martinique au format ASC, donc les spécifications se situent ici : https://desktop.arcgis.com/fr/arcmap...ter-format.htm

Ces données sont fournies dans le système de coordonnées correspondant à la zone 20N de l'UTM (EPSG:32620 https://epsg.io/32620).

Ces fichiers sont exploitables par un SIG (par exemple QGIS), ou bien directement en lecture par un script, car ils sont au format texte.

Bon courage pour votre TIPE,

Cordialement,


EDIT:
Ayant récemment travaillé sur le sujet, j'ai écrit une fonction en Python pour l'ouverture de ce genre de fichier. Je vous laisse reprendre le code à votre guise :

Code python :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
import numpy as np
 
def openAsc(filename):
  """
  Ouvre un fichier de MNT .asc et renvoie un dictionnaire correspondant à toutes les informations contenues dans le fichier
 
  Parameters
  ----------
  filename : str
    nom du ficher .asc à ouvrir
  Returns
  -------
  dict
    description du fichier.
    Clés du dictionnaire :
      ncols: nombre de colonnes, int
      nrows: nombre de lignes, int
      xllcorner: coordonnée X du coin inférieur gauche, float 
      yllcorner: coordonnée Y du coin inférieur gauche, float
      cellsize: longueur du côté d'une case en m, float
      NODATA_value: valeur d'une case sans donnée, float
      data: données du fichier MNT, np.array
  """
  result = {}
  try:
    with open(filename, 'r') as file:
      lines = file.readlines()
  except IOError:
    raise IOError("Le fichier MNT n'a pas été correctement ouvert")
 
  for i in range(6):
    try:
      line = lines[i]
    except IndexError:
      raise IOError("Le fichier .asc est mal formé : moins de 6 lignes d'en-tête")
 
    key, value = line.split()
    if i <= 1 :
      result[key] = int(value)
    else:
      result[key] = float(value)
 
  dem_list = []
  for line in lines[6:]:
    dem_list.append(np.array(line.split(), dtype=float))
 
  result['data'] = np.array(dem_list)
 
  return result

[azarzelli]

Bonjour,

Ma réponse vous convient-elle ? Dans l'absence de réponse à ce sujet, je le marquerai comme résolu et le fermerai.

Cordialement,

[Alaopi]

Bonjour,

Oui, votre réponse m'a été d'une grande aide. Veuillez m'excuser de ne pas avoir répondu plus tôt. Merci encore une fois.

Cordialement,