Commit fbd6b623 authored by Commandre Benjamin's avatar Commandre Benjamin
Browse files

Correction lat,long

1 merge request!4Develop
Showing with 12 additions and 10 deletions
+12 -10
...@@ -66,7 +66,7 @@ class Archive(): ...@@ -66,7 +66,7 @@ class Archive():
self.liste_archive = [] self.liste_archive = []
self.url = '' self.url = ''
def vector_projection_with_epsg(self, chemin_vecteur, output_path, epsg=32631): def vector_projection_with_epsg(self, chemin_vecteur, output_path, epsg=32630):
driver = ogr.GetDriverByName('ESRI Shapefile') driver = ogr.GetDriverByName('ESRI Shapefile')
...@@ -155,7 +155,7 @@ class Archive(): ...@@ -155,7 +155,7 @@ class Archive():
area_coord_corner = str(extent_[2]) + ',' + str(extent_[0]) + ',' + str(extent_[3]) + ',' + str(extent_[1]) area_coord_corner = str(extent_[2]) + ',' + str(extent_[0]) + ',' + str(extent_[3]) + ',' + str(extent_[1])
return "{},{},{},{}".format(extent_[2], extent_[0], extent_[3], extent_[1]) return "{},{},{},{}".format(extent_[0], extent_[2], extent_[1], extent_[3])
def listing_by_tile(self): def listing_by_tile(self):
...@@ -172,8 +172,10 @@ class Archive(): ...@@ -172,8 +172,10 @@ class Archive():
def listing_by_coord(self): def listing_by_coord(self):
self.requete['box'] = self.coord_box_dd() self.requete['box'] = self.coord_box_dd()
self.url = "{0}/{1}/api/collections/{2}/search.json?{3}".format(self.serveur, self.resto, self.capteur, urlencode(self.requete)) self.url = "{0}/{1}/api/collections/{2}/search.json?{3}".format(self.serveur, self.resto, self.capteur, urlencode(self.requete))
print(self.url)
self.get_liste_images() self.get_liste_images()
...@@ -287,7 +289,7 @@ class Archive(): ...@@ -287,7 +289,7 @@ class Archive():
sauvegarde = configparser.ConfigParser() sauvegarde = configparser.ConfigParser()
# Pour toutes les images disponible # Pour toutes les images disponible
for l in self.liste_archive : for l in self.liste_archive :
# Récupération de la date quand l'image a été prise # Récupération de la date quand l'image a été prise
...@@ -299,7 +301,7 @@ class Archive(): ...@@ -299,7 +301,7 @@ class Archive():
# Lecture du fichier de sauvegarde # Lecture du fichier de sauvegarde
sauvegarde.read("{}/sauvegarde.ini".format(self.dossier_sortie)) sauvegarde.read("{}/sauvegarde.ini".format(self.dossier_sortie))
# Si la date existe dans le fichier de sauvegarde ... # Si la date existe dans le fichier de sauvegarde ...
if str(date) in sauvegarde.keys() : if str(date) in sauvegarde.keys() :
# ... mais pas l'image, on la rajoute dans la liste des images à traiter. # ... mais pas l'image, on la rajoute dans la liste des images à traiter.
...@@ -373,7 +375,7 @@ class Archive(): ...@@ -373,7 +375,7 @@ class Archive():
del liste_content del liste_content
# Mis à jour du fichier de sauvegarde # Mis à jour du fichier de sauvegarde
# Lecture du fichier de sauvegarde # Lecture du fichier de sauvegarde
sauvegarde.read("{}/sauvegarde.ini".format(self.dossier_sortie)) sauvegarde.read("{}/sauvegarde.ini".format(self.dossier_sortie))
...@@ -486,7 +488,7 @@ class Archive(): ...@@ -486,7 +488,7 @@ class Archive():
tuiles_image = [] tuiles_image = []
tuiles_nuage = [] tuiles_nuage = []
# Options Gdal pour la fusion des bandes # Options Gdal pour la fusion des bandes
options_vrt = gdal.ParseCommandLine('-resolution highest -srcnodata -10000 -vrtnodata NaN -separate') options_vrt = gdal.ParseCommandLine('-resolution highest -srcnodata -10000 -vrtnodata NaN -separate')
self.logger.info("Extraction des images") self.logger.info("Extraction des images")
...@@ -516,8 +518,8 @@ class Archive(): ...@@ -516,8 +518,8 @@ class Archive():
# liste_bandes.append(gdal.Open(mmap_name)) # liste_bandes.append(gdal.Open(mmap_name))
liste_bandes.append(mmap_name) liste_bandes.append(mmap_name)
# On fusionne les différentes bandes en une seule image # On fusionne les différentes bandes en une seule image
# on découpe l'image selon l'emprise # on découpe l'image selon l'emprise
# et on conserve l'image obtenue # et on conserve l'image obtenue
vrt = gdal.BuildVRT("", liste_bandes, options=options_vrt) vrt = gdal.BuildVRT("", liste_bandes, options=options_vrt)
...@@ -568,7 +570,7 @@ class Archive(): ...@@ -568,7 +570,7 @@ class Archive():
self.logger.debug("Dossier image : {0}".format(dossier)) self.logger.debug("Dossier image : {0}".format(dossier))
# On effectue une mosaïque des images qu'on découpe selon l'emprise. # On effectue une mosaïque des images qu'on découpe selon l'emprise.
Outils.clip(tuiles_image, self.emprise, form="GTiff", dst="{0}/{1}.tif".format(dossier,date)) Outils.clip(tuiles_image, self.emprise, form="GTiff", dst="{0}/{1}.tif".format(dossier,date))
del tuiles_nuage del tuiles_nuage
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment