Suite

Impossible de déformer les fichiers HDF5

Impossible de déformer les fichiers HDF5


J'ai de nombreux fichiers HDF5 (.nc) dont j'ai besoin pour traiter le warp par lots à l'aide de l'un des utilitaires gdal gdalwarp. Lorsque j'ai essayé de déformer les fichiers, une erreur s'est produite :

CONTRIBUTION:

gdalwarp -geoloc -te 109.975 3.475 135.025 25.025 HDF5:"@file"://geophysical_data/chlor_a %out_path%@fname.tif"

RÉSULTAT:

ERREUR 1 : Impossible de calculer une transformation basée sur GEOLOC_ARRAY entre pixel/ligne et coordonnées géoréférencées pour HDF5:A2015045060500.L2_LAC_OC.nc://geophysical _data/chlor_a.

Mise à jour1 :

Juste pour que ce soit clair, voulez-vous dire dans lat.vrt, lon.vrt et chlor.vrt, je devrais supprimer les identifiants GCP et la clé MDI et insérer cette section :

 oc-long.vrt 1 oc-lat.vrt 1 0 0 1 1

entre cette section ?

 ###### section métadonnées ici #######   HDF5:A2015194044000.L2_LAC.SeAHABS.nc://geophysical_data/chlor_a 1   

Après quelques tests, je pense que le geoloc ne fonctionne pas correctement. J'ai donc utilisé la méthode alternative en utilisant des fichiers vrt créés manuellement :

  1. Créer un fichier nommélon.vrt:
 GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTORITÉ["EPSG","7030"]],TOWGS84[0,0,0,0,0,0, 0],AUTORITÉ["EPSG","6326"]],PRIMEM["Greenwich",0,AUTORITÉ["EPSG","8901"]],UNIT["degré",0.0174532925199433,AUTORITÉ["EPSG"," 9108"]],AUTORITÉ["EPSG","4326"]]   HDF5:A2015045060000.L2_LAC_OC.nc://navigation_data/longitude 1   
  1. Idem pour les latitudes enlat.vrt:
 GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTORITÉ["EPSG","7030"]],TOWGS84[0,0,0,0,0,0, 0],AUTORITÉ["EPSG","6326"]],PRIMEM["Greenwich",0,AUTORITÉ["EPSG","8901"]],UNIT["degré",0.0174532925199433,AUTORITÉ["EPSG"," 9108"]],AUTORITÉ["EPSG","4326"]]   HDF5:A2015045060000.L2_LAC_OC.nc://navigation_data/latitude 1   
  1. et pour les donnéeschlor.vrt:
  lon.vrt 1 lat.vrt 1 0 0 1 1   HDF5:A2015045060000.L2_LAC_OC.nc://geophysical_data/chlor_a 1   
  1. Faire le gauchissement avec :

gdalwarp -geoloc -t_srs EPSG:4326 chlor.vrt chlor-out.tif

et le résultat s'adapte aux rivages autour de Bornéo :


Alternativement à la création manuelle des vrts, vous pouvez les créer avec GDAL :

gdal_translate -de VRT HDF5:A2015045060000.L2_LAC_OC.nc://navigation_data/longitude lon.vrt gdal_translate -de VRT HDF5:A2015045060000.L2_LAC_OC.nc://navigation_data/latitude lat.vrt.vrt gdal2_translate_OC201504 HDF60000.L2_LAC_OC.nc://navigation_data/latitude lat.vrt.L.L. nc://geophysical_data/chlor_a chlor.vrt

Avec un bon éditeur de texte, supprimez les listes GCP de toutes et insérez seul dans le chlor.vrt cette section à la place :

 lon.vrt 1 lat.vrt 1 0 0 1 1

Puis cours

gdalwarp -geoloc -t_srs EPSG:4326 -écraser chlor.vrt chlor-vrt.tif

pour obtenir la même image que ci-dessus.


Une autre solution, fonctionnant avec des points GCP modifiés manuellement, peut être trouvée dans ma réponse pour Utiliser GDALwarp pour reprojeter le fichier netCDF ?


Voir la vidéo: An Introduction to HDF5