Ortho4XP

[atanorr]

voir post #251

Sorry, je suis passé à côté , Vraiment Grand merci pour votre aide c'est nickel

Quand je dit que j'ai la sorcière ... pffff
Etape 2
--------

-> Chargement des données altimétriques
Traceback (most recent call last):
File "./Ortho4XP_alpha2.py", line 2067, in <module>
alt_dem=load_altitude_matrix(lat,lon)
File "./Ortho4XP_alpha2.py", line 520, in load_altitude_matrix
return numpy.where(alt1==-32767,alt2,alt1)
ValueError: shape mismatch: objects cannot be broadcast to a single shape

@Odevillée
As tu le même problème sur la tuile 51 2 ?

[jasum]

Bravo OscarPilote et merci de nous faire partager vos compétences.

j'ai testé:
la tuile +48-6 zoom 18 BI avec la config par défaut aucun problème,
la tuile +45+6 17 BI pour essayer le mailleur,et là je ne reconnais pas la vue de Chamonix sur l'Aiguille Verte( les Drus semblent 'rabotés' de même que les Aiguilles de Chamonix), j'ai essayé en variant les options de config sans résultat (parfois pire) par contre les altiports sont presque parfaits comparés aux meshs d'AlpilotX qui présentent des aberrations importantes au niveau des altiports et aérodromes de montagnes.

Quelles options de maillage pour optimiser le mesh me conseilleriez vous;

PS:j'ai un warning à l'etape2:

Etape 2
--------

-> Chargement des données altimétriques
-> Calcul de la matrice de courbure du terrain
./Ortho4XP:499: RuntimeWarning: invalid value encountered in sqrt
hess_max_eig=abs(a+b+numpy.sign(a+b)*numpy.sqrt((a-b)**2+4*c**2))/2.0

[atanorr]

Bravo OscarPilote et merci de nous faire partager vos compétences.

j'ai testé:
la tuile +48-6 zoom 18 BI avec la config par défaut aucun problème,
la tuile +45+6 17 BI pour essayer le mailleur,et là je ne reconnais pas la vue de Chamonix sur l'Aiguille Verte( les Drus semblent 'rabotés' de même que les Aiguilles de Chamonix), j'ai essayé en variant les options de config sans résultat (parfois pire) par contre les altiports sont presque parfaits comparés aux meshs d'AlpilotX qui présentent des aberrations importantes au niveau des altiports et aérodromes de montagnes.

Quelles options de maillage pour optimiser le mesh me conseilleriez vous;

Odevillée ou Oscar ,
mais je tenterais ceci:
curvature_tol = 10

hmin = 5

hmax = 6000 , puisque hmax = 1000 si j'ai bien compris , le mailleur cesse de raffiner

Essaie en zoom 14 , juste le hmax pour commencer..

[jasum]

Merci atanor ,
j'essaie de suite.

[Oscar Pilote]

la tuile +45+6 17 BI pour essayer le mailleur,et là je ne reconnais pas la vue de Chamonix sur l'Aiguille Verte( les Drus semblent 'rabotés' de même que les Aiguilles de Chamonix)

C'est ma tuile "cible" depuis le début, celle que j'essaierai de faire tendre vers ce dont je peux rêver.
Qu'entends-tu par rabotés ? Dans les essais que j'ai déjà faits on distingue quand même bien les drus
et les aiguilles de Chamonix. Bien sûr, les fichiers d'altitude perdent en précision dans les zones très accidentées,
c'est inhérent à la manière dont ils sont obtenus. Il y aurait le RGE Alti de l'ign qui devrait faire bien mieux, mais
je n'y ai pas accès, je cherche depuis un bon moment une solution alternative. On ne distingue pas la différence
entre le grand et le petit dru par exemple, et l'arête ecclésiastique ne laisse pas apparaître tous ses membres,
mais c'est déjà pas mal (j'essaierai de faire une vidéo, il faudra que je refasse la tuile, c'est un peu le chantier
dans mes répertoires en ce moment...).
Je te conseille vivement le ZL19 en FR juste sur la zone en question, tu auras des zones toutes blanches de
l'autre côté des Jorasses car FR s'arrête à la frontière, mais le survol de la mer de glace et de l'envers des
aiguilles est bluffant.

Pour les paramètres :
hmax pas très important (1000 si tu veux)
hmin à 10 ou 20
max_curv à 10 par exemple (regarde le nombre de triangles que cela engendre à l'étape du maillage, plus max_curv
est petit plus tu auras de triangles. Pour comparaison les hdv3 on environ 2.5 millions de triangles, les uhd 6 millions
si ma mémoire est bonne. A nombre de triangles égal tu devrais avoir un meilleur rendu avec le mailleur d'ortho4xp
( = le mailleur de Jonathan Shewchuk pour l'essentiel) qu'avec celui de Laminar.

PS) Gentoo, cela ne rigole pas ;-)

[Oscar Pilote]

-> Chargement des données altimétriques
Traceback (most recent call last):
File "./Ortho4XP_alpha2.py", line 2067, in <module>
alt_dem=load_altitude_matrix(lat,lon)
File "./Ortho4XP_alpha2.py", line 520, in load_altitude_matrix
return numpy.where(alt1==-32767,alt2,alt1)
ValueError: shape mismatch: objects cannot be broadcast to a single shape

On en a parlé hier, le SRTM de la 51 2 (comme elle ne contient que de l'eau essentiellement)
a été taillé de moitié par ses auteurs...
J'ai fait la modif nécessaire pour pallier à cela dans le code et uploadé sur
sourceforge, si tu ne veux pas retélécharger il suffit de remplacer :

return numpy.where(alt1==-32767,alt2,alt1)

par

try:
    return numpy.where(alt1==-32767,alt2,alt1)
except:
    return alt2

là où cette ligne apparaît (ligne 520), attention à respecter l'indentation.

[Odevillée]

Ben moi je reste calé à l'étape 2 avec n'importe queele tuile.

Etape 2
--------

-> Chargement des données altimétriques
Traceback (most recent call last):
File "./Ortho4XP_alpha2.py", line 2067, in <module>
alt_dem=load_altitude_matrix(lat,lon)
File "./Ortho4XP_alpha2.py", line 520, in load_altitude_matrix
return numpy.where(alt1==-32767,alt2,alt1)
ValueError: shape mismatch: objects cannot be broadcast to a single shape

[jasum]

Merci pour ta réponse

Ce n'est pas catastrophique, oui on distingue les 2 drus mais le col est vraiment peu marqué.(manque la forme en pain de sucre du grand dru)

j'essayerai en ZL19 FR ça va être lourd mais j'ai de la place sur mes disques.

[atanorr]

return numpy.where(alt1==-32767,alt2,alt1)

par

try:
    return numpy.where(alt1==-32767,alt2,alt1)
except:
    return alt2

là où cette ligne apparaît (ligne 520), attention à respecter l'indentation.

Vraiment merci oscar pour ta rapidité à résoudre les soucis..
le site de source etant down pour le moment ( sorcière au dessus de ma tête comme d hab )
j'ai rentré la ligne

comme ceci:
def load_altitude_matrix(lat,lon):
print('-> Chargement des données altimétriques')
filename1=downloaded_dem_filename(lat,lon,'SRTM')
filename2=downloaded_dem_filename(lat,lon,'GDEM')
try:
ds1=gdal.Open(filename1)
ds2=gdal.Open(filename2)
except:
print("Le ou les fichiers DEM au format .tif sont manquants.")
alt1=numpy.array(ds1.GetRasterBand(1).ReadAsArray( ))
alt2=numpy.array(ds2.GetRasterBand(1).ReadAsArray( ))
try:
return numpy.where(alt1==-32767,alt2,alt1)
except:
return alt2

mais sorcière... ==>
yves@ubuntu:~/ortho4xp$ ./Ortho4XP_alpha2.py build_tile 51 2
File "./Ortho4XP_alpha2.py", line 521
return numpy.where(alt1==-32767,alt2,alt1)
^
SyntaxError: 'return' outside function , mal écris probablement..?

[Odevillée]

Ben moi je reste calé à l'étape 2 avec n'importe queele tuile.

Etape 2
--------

-> Chargement des données altimétriques
Traceback (most recent call last):
File "./Ortho4XP_alpha2.py", line 2067, in <module>
alt_dem=load_altitude_matrix(lat,lon)
File "./Ortho4XP_alpha2.py", line 520, in load_altitude_matrix
return numpy.where(alt1==-32767,alt2,alt1)
ValueError: shape mismatch: objects cannot be broadcast to a single shape

correction effectuée, test en cours ça mouline.

[atanorr]

correction effectuée, test en cours ça mouline.

ça mouline également , après un bon verre de rosé :D

Mais pas longtemps ahahaha :D

Etape 2
--------

-> Chargement des données altimétriques
-> Calcul de la matrice de courbure du terrain
-> Lancement du mailleur proprement dit :
./Utils/Triangle4XP: 1: ./Utils/Triangle4XP: Syntax error: "(" unexpected

Etape 3
--------
-> Calcul des altitudes des sommets du maillage par interpolation
Traceback (most recent call last):
File "./Ortho4XP_alpha2.py", line 2080, in <module>
vertices=build_3D_vertex_array(lat,lon,alt_dem)
File "./Ortho4XP_alpha2.py", line 540, in build_3D_vertex_array
f_node = open(node_filename,'r')
FileNotFoundError: [Errno 2] No such file or directory: '+51+002.1.node'

je vais reprendre un verre...

[Odevillée]

ça mouline également , après un bon verre de rosé :D

Mais pas longtemps ahahaha :D

Etape 2
--------

-> Chargement des données altimétriques
-> Calcul de la matrice de courbure du terrain
-> Lancement du mailleur proprement dit :
./Utils/Triangle4XP: 1: ./Utils/Triangle4XP: Syntax error: "(" unexpected

Etape 3
--------
-> Calcul des altitudes des sommets du maillage par interpolation
Traceback (most recent call last):
File "./Ortho4XP_alpha2.py", line 2080, in <module>
vertices=build_3D_vertex_array(lat,lon,alt_dem)
File "./Ortho4XP_alpha2.py", line 540, in build_3D_vertex_array
f_node = open(node_filename,'r')
FileNotFoundError: [Errno 2] No such file or directory: '+51+002.1.node'

je vais reprendre un verre...

Pareil pour moi tuile +51+002 OK, par contre tuile +50+004 erreur.

Etape 3
--------
-> Calcul des altitudes des sommets du maillage par interpolation
Traceback (most recent call last):
File "./Ortho4XP_alpha2.py", line 2080, in <module>
vertices=build_3D_vertex_array(lat,lon,alt_dem)
File "./Ortho4XP_alpha2.py", line 540, in build_3D_vertex_array
f_node = open(node_filename,'r')
FileNotFoundError: [Errno 2] No such file or directory: '+50+004.1.node'

On vas finir bourré si ça continue .

[atanorr]

c'est étrange.. pourquoi chez toi la tuile +51 +002 passe et pas chez moi...
je n'ai rien changé juste comme le montre oscar afin de créer le mask , le zoom en 14

On va attendre le maître de la lumière la dessus :D
ça mérite un verre tout ça..

[Odevillée]

c'est étrange.. pourquoi chez toi la tuile +51 +002 passe et pas chez moi...
je n'ai rien changé juste comme le montre oscar afin de créer le mask , le zoom en 14

On va attendre le maître de la lumière la dessus :D
ça mérite un verre tout ça..

Le l'ai fait sans masque pour essayer, direct en z17, paramètres water_overlay=True.

[atanorr]

Le l'ai fait sans masque pour essayer, direct en z17, paramètres water_overlay=True.

et avec water_overlay = false en zoom14?
le but de la manœuvre est de travailler avec le masque qui devrais permettre une découpe parfaite ,