Suite

Extraire plusieurs grilles de HDF avec l'outil de sous-ensemble à l'aide d'ArcPy


Cela devrait être une chose assez simple mais qui me dérange pour le moment. J'ai ~5k fichiers HDF dans un dossier avec des noms :

MOD11A1.A2001001.h24v05.005.2006343042921 et ainsi de suite…

J'ai la suite qui les parcourt toutes pour extraire une grille de sous-ensemble de chaque fichier HDF. Pour une raison inconnue, le script python tombe en panne à chaque fois autour du numéro de fichier 3557. Je suis capable d'extraire manuellement plusieurs grilles autour de ce numéro à partir de fichiers HDF. Peut-être que j'atteins la limite de RAM. Mais même avec toutes les autres applications fermées, il présente le même comportement.

rasters = arcpy.ListRasters("*", "All") pour inLST dans les rasters : f = inLST[9:16] outLST = "D:/temp/subset/" + str(f) + '_lst' arcpy.ExtractSubDataset_management (inLST, outLST, "0")

Existe-t-il une alternative àarcpy.ListRasters, que puis-je utiliser ?


l'extraction de plus de 1k fichiers peut éventuellement avoir une mémoire saturée avec le nom du chemin de fichier"D:/temp/subset/" + str(f) + '_lst'. Cependant, j'ai pu l'exécuter en réduisant le chemin vers"D:/" + str(f) + '_lst'. Voici un extrait complet qui a fonctionné pour moi :

rasters = arcpy.ListRasters("*", "All") pour inLST dans les rasters : f = inLST[9:16] outLST = "D:/" + str(f) + '_lst' arcpy.ExtractSubDataset_management(inLST, outLST , "0")