Annonce

Réduire
Aucune annonce.

LIDAR SHOM Reunion Ortho4XP130 : ERROR: in reading .tif

Réduire
X
 
  • Filtre
  • Heure
  • Afficher
Tout nettoyer
nouveaux messages

  • LIDAR SHOM Reunion Ortho4XP130 : ERROR: in reading .tif

    Salut,

    J'ai récupéré les fichiers .asc de la Réunion : http://diffusion.shom.fr/pro/amenage...union2016.html

    J'ai créé la mosaic virtuelle via la command

    D:\LIDAR\Reunion\MNT5m>gdalbuildvrt renuion.vrt *.asc (2817 fichiers)
    0...10...20...30...40...50...60...70...80...90...1 00 - done.
    Puis j'ai généré le fichier .tif sachant que le format d'origine n'est pas du Lambert93 comme pour Ouessant https://www.x-plane.fr/showthread.ph...l=1#post537078

    D:\LIDAR\Reunion\MNT5m>gdalwarp -s_srs IGNF:RGR92UTM40S -t_srs EPSG:4326 reunion.vrt reunion.tif
    0...10...20...30...40...50...60...70...80...90...1 00 - done.
    L'image fait 631*Mo


    Mais sur Ortho4XP

    Step 1
    Loading elevation data and smoothing it over airports.
    ERROR: in reading D:/LIDAR/Reunion/MNT5m/reunion.tif -> replaced with zero altitude.
    Ou Step 3
    -> Construction of the masks
    ERROR: in reading D:/LIDAR/Reunion/MNT5m/reunion.tif -> replaced with zero altitude.
    * Min altitude: 0.0 , Max altitude: 0.0 , Mean: 0.0
    Pourquoi Ortho4XP130 n'arrive pas à lire ce fichier ? Trop gros pour 681Mo ? Mauvaise conversion (exemple de .asc : LITTO3D_REU_0376_7648_MNT5_20120427_UTM40S_RGR92_R EUN89 )?
    PPL en cours sur LFPN avec l'Aéro-club AIR FRANCE

    Machine : Amstrad CPC464 sous Vulkan
    https://ailgorbot.wordpress.com/

  • #2
    Envoyé par ailgorbot Voir le message
    Salut,

    J'ai récupéré les fichiers .asc de la Réunion : http://diffusion.shom.fr/pro/amenage...union2016.html

    J'ai créé la mosaic virtuelle via la command



    Puis j'ai généré le fichier .tif sachant que le format d'origine n'est pas du Lambert93 comme pour Ouessant https://www.x-plane.fr/showthread.ph...l=1#post537078



    L'image fait 631*Mo


    Mais sur Ortho4XP

    Step 1


    Ou Step 3
    -> Construction of the masks


    Pourquoi Ortho4XP130 n'arrive pas à lire ce fichier ? Trop gros pour 681Mo ? Mauvaise conversion (exemple de .asc : LITTO3D_REU_0376_7648_MNT5_20120427_UTM40S_RGR92_R EUN89 )?
    J'ai essayé sur 4 Lidar au niveau de La Possession et je n'ai pas eu de problème.
    Vérifie ton tif: gdalinfo xxx.tif
    Gentoo Linux x86_64 RYZEN 7900X RTX2070 X-Plane12 & Hackintosh Somona (iGPU intelHD630) IntelCore i7 CPU 7700 @ 3.60GHz X-Plane 11

    Commentaire


    • #3
      J'ai lancé la commande. Le fichier .tif me semble correcte, même si je n'y connaît rien

      D:\LIDAR\Reunion\MNT5m>gdalinfo reunion.tif
      Driver: GTiff/GeoTIFF
      Files: reunion.tif
      Size is 13903, 11897
      Coordinate System is:
      GEOGCS["WGS 84",
      DATUM["WGS_1984",
      SPHEROID["WGS 84",6378137,298.257223563,
      AUTHORITY["EPSG","7030"]],
      AUTHORITY["EPSG","6326"]],
      PRIMEM["Greenwich",0],
      UNIT["degree",0.0174532925199433],
      AUTHORITY["EPSG","4326"]]
      Origin = (55.196091909454360,-20.844639423434963)
      Pixel Size = (0.000046784769716,-0.000046784769716)
      Metadata:
      AREA_OR_POINT=Area
      Image Structure Metadata:
      INTERLEAVE=BAND
      Corner Coordinates:
      Upper Left ( 55.1960919, -20.8446394) ( 55d11'45.93"E, 20d50'40.70"S)
      Lower Left ( 55.1960919, -21.4012378) ( 55d11'45.93"E, 21d24' 4.46"S)
      Upper Right ( 55.8465406, -20.8446394) ( 55d50'47.55"E, 20d50'40.70"S)
      Lower Right ( 55.8465406, -21.4012378) ( 55d50'47.55"E, 21d24' 4.46"S)
      Center ( 55.5213162, -21.1229386) ( 55d31'16.74"E, 21d 7'22.58"S)
      Band 1 Block=13903x1 Type=Float32, ColorInterp=Gray
      NoData Value=-9999.900390625

      D:\LIDAR\Reunion\MNT5m>
      via cet URL tous les .asc de la Réunion : https://www.dropbox.com/s/ktxcdw3k4tv9ugz/MNT5m.7z?dl=0
      PPL en cours sur LFPN avec l'Aéro-club AIR FRANCE

      Machine : Amstrad CPC464 sous Vulkan
      https://ailgorbot.wordpress.com/

      Commentaire


      • #4
        Bonsoir aligorbot

        J'ai construit le tif en WGS84 avec les fichiers que tu as fournis, ce tif passe sans probléme chez moi.


        Regarde si il passe chez toi.
        Dernière modification par jasum, 23 octobre 2018, 14h13.
        Gentoo Linux x86_64 RYZEN 7900X RTX2070 X-Plane12 & Hackintosh Somona (iGPU intelHD630) IntelCore i7 CPU 7700 @ 3.60GHz X-Plane 11

        Commentaire


        • #5
          On obtient bien le même fichier .tif

          Ton .tif : LaReunion_WGS84.tif = 631*Mo (661*711*534 octets)
          Mon .tif :reunion.tif = 631*Mo (661*711*534 octets)

          Et aussi la même erreur.
          Loading elevation data and smoothing it over airports.
          ERROR: in reading D:/LIDAR/Reunion/MNT5m/LaReunion_WGS84.tif -> replaced with zero altitude.
          Tu as pu créer les tuiles sans erreur à la lecture du .tif ?
          PPL en cours sur LFPN avec l'Aéro-club AIR FRANCE

          Machine : Amstrad CPC464 sous Vulkan
          https://ailgorbot.wordpress.com/

          Commentaire


          • #6
            Envoyé par ailgorbot Voir le message
            Tu as pu créer les tuiles sans erreur à la lecture du .tif ?
            Oui, sans aucun problème!!!
            Gentoo Linux x86_64 RYZEN 7900X RTX2070 X-Plane12 & Hackintosh Somona (iGPU intelHD630) IntelCore i7 CPU 7700 @ 3.60GHz X-Plane 11

            Commentaire


            • #7
              OK je vais réinstaller ortho4xp 1.3. Il y avoir un louper.
              PPL en cours sur LFPN avec l'Aéro-club AIR FRANCE

              Machine : Amstrad CPC464 sous Vulkan
              https://ailgorbot.wordpress.com/

              Commentaire


              • #8
                Finalement même avec la nouvelle version de ortho4xp fraîchement télécharger, il échoue à la lecture du fichier

                Step 1 : Building vector data for tile -22+055 :
                --------

                -> Dealing with airports
                * Recycling OSM data from .\OSM_data\-30+050\-22+055\-22+055_airports.osm.bz2
                FMEP Pierrefonds Airport 2 runways, lat= -21.32, lon= 55.42
                **** **** 1 runway , lat= -21.03, lon= 55.26
                Loading elevation data and smoothing it over airports.
                ERROR: in reading D:/LIDAR/Reunion/MNT5m/LaReunion_WGS84.tif -> replaced with zero altitude.
                * Min altitude: 0.0 , Max altitude: 0.0 , Mean: 0.0

                Auto-patched 3 runways and 12 pieces of taxiway.
                Flattened 22 helipads.
                Number of edges at this point: 801
                -> Dealing with roads
                * Recycling OSM data from .\OSM_data\-30+050\-22+055\-22+055_big_roads.osm.bz2
                * Checking which large roads need levelling.
                * Buffering banked road network as multipolygon.
                Encoding it.
                * Encoding the remaining primary road network as linestrings.
                Number of edges at this point: 9293
                Process interrupted
                __________________________________________________ _______________________________________________

                Pourtant je pense être bon au niveau des librairies.

                $ pip install pyproj
                Requirement already satisfied: pyproj in c:\users\ailgorbot\appdata\local\programs\python\p ython36-32\lib\site-packages (1.9.5.1)

                $ pip install Numpy
                Requirement already satisfied: Numpy in c:\users\ailgorbot\appdata\local\programs\python\p ython36-32\lib\site-packages (1.15.3)

                $ pip install Gdal
                Requirement already satisfied: Gdal in c:\users\ailgorbot\appdata\local\programs\python\p ython36-32\lib\site-packages (2.3.2)

                $ pip install Shapely
                Requirement already satisfied: Shapely in c:\users\ailgorbot\appdata\local\programs\python\p ython36-32\lib\site-packages (1.6.4.post1)

                $ pip install Rtree
                Requirement already satisfied: Rtree in c:\users\ailgorbot\appdata\local\programs\python\p ython36-32\lib\site-packages (0.8.3)
                Requirement already satisfied: setuptools in c:\users\ailgorbot\appdata\local\programs\python\p ython36-32\lib\site-packages (from Rtree) (39.0.1)

                $ pip install Pillow
                Requirement already satisfied: Pillow in c:\users\ailgorbot\appdata\local\programs\python\p ython36-32\lib\site-packages (5.3.0)

                $ python3 --version
                Python 3.6.0
                PPL en cours sur LFPN avec l'Aéro-club AIR FRANCE

                Machine : Amstrad CPC464 sous Vulkan
                https://ailgorbot.wordpress.com/

                Commentaire


                • #9
                  Bonjour aligorbot

                  C'est normal, tu fais une erreur de manipulation:

                  A l'étape1, Ortho4XP a besoin d'DEM de 1°x1° d'arc terrestre, ce qui n'est pas le cas de ce tif.
                  La manip correcte est de lancer Ortho4XP avec le DEM par default (il est recherché automatiquement)
                  L'étape2 terminée, tu fais un iterate= 1 en indiquant ton tif et tu relance l'étape2, ce tif vient "rafiner" le mesh dejà construit ( n'oublie pas cleaning_level= 0 dans ta configuration)
                  Gentoo Linux x86_64 RYZEN 7900X RTX2070 X-Plane12 & Hackintosh Somona (iGPU intelHD630) IntelCore i7 CPU 7700 @ 3.60GHz X-Plane 11

                  Commentaire


                  • #10
                    Bonjour tous les deux,

                    Il n'est pas nécessaire que le DEM suive les bords de la tuile, ni même qu'il la recouvre totalement (dans ce cas c'est étendu par zéro en dehors, ce qui n'est pas gênant pour une île) . Pour des données où toute la mer (comme dans ce tif) est mise en no_data il faut mettre fill_nodata à False, et le programme transforme alors no_data en 0. Mais je pense que l'erreur d'aligorbot est ailleurs (installation de gdal je dirais), pourrait-il essayer avec la version .exe pour vérifier. Sinon pour comprendre ce qui bloque lance python dans une fenêtre de commande, puis "import gdal", puis "ds=gdal.Open('Reunion.tif')", et enfin "arr=ds.GetRasterBand(1).ReadAsArray()". Il devrait te dire l'erreur qu'il rencontre.

                    Par ailleurs les données du SHOM contiennent parfois des absurdités en quelques points (ici il y a deux pics de deux pixels de large), je vous ferai parvenir les données corrigées.
                    Dernière modification par Oscar Pilote, 24 octobre 2018, 10h28.
                    Linux Debian sid - Intel i7 4Ghz - 16Gb RAM - Nvidia GTX 970
                    Ortho4XP v130 : Dépôt github (mises à jour au fil de l'eau (avant la période de sécheresse)), ou Version clés en main sous Windows.

                    Commentaire


                    • #11
                      Merci Oscar pour ces précisions

                      A propos des aberrations, une fois celles-ci repérées,corriges-tu le mesh par un patch osm ou bien édites-tu le MNT dans QGIS (par exemple) et corriges les "entités" concernées.
                      Dernière modification par jasum, 24 octobre 2018, 14h13.
                      Gentoo Linux x86_64 RYZEN 7900X RTX2070 X-Plane12 & Hackintosh Somona (iGPU intelHD630) IntelCore i7 CPU 7700 @ 3.60GHz X-Plane 11

                      Commentaire


                      • #12
                        YES!!!, le lancement via le .exe corrige la situation.

                        Pour info
                        Python 3.6.5 (v3.6.5:f59c0932b4, Mar 28 2018, 16:07:46) [MSC v.1900 32 bit (Intel)] on win32
                        Type "help", "copyright", "credits" or "license" for more information.
                        >>> import gdal
                        >>> ds=gdal.Open('D:\LIDAR\Reunion\MNT5m\Reunion.tif')
                        >>> arr=ds.GetRasterBand(1).ReadAsArray()
                        >>> print(arr)
                        [[-9999.9 -9999.9 -9999.9 ... -9999.9 -9999.9 -9999.9]
                        [-9999.9 -9999.9 -9999.9 ... -9999.9 -9999.9 -9999.9]
                        [-9999.9 -9999.9 -9999.9 ... -9999.9 -9999.9 -9999.9]
                        ...
                        [-9999.9 -9999.9 -9999.9 ... -9999.9 -9999.9 -9999.9]
                        [-9999.9 -9999.9 -9999.9 ... -9999.9 -9999.9 -9999.9]
                        [-9999.9 -9999.9 -9999.9 ... -9999.9 -9999.9 -9999.9]]
                        >>>
                        Merci à vous deux.

                        J'ai vu ici et là quelques imperfections. A regarder de plus près
                        PPL en cours sur LFPN avec l'Aéro-club AIR FRANCE

                        Machine : Amstrad CPC464 sous Vulkan
                        https://ailgorbot.wordpress.com/

                        Commentaire

                        Chargement...
                        X