J'ai besoin de mosaïquer environ 550 Go d'imagerie tif ensemble et le logiciel que j'ai essayé échoue toujours. La zone a été divisée en zones afin que la plus petite ait environ 200 tuiles.
J'ai utilisé les dernières versions d'ERDAS (Imagine and Mapper), ArcINFO et Global Mapper sur un Intel Xeon E31245 3,30 gigahertz, DELL, 16 Go de RAM, Win 7 Professional 64 bits. Machine Mullti-core (4 au total), hyper-filetée (8 au total). Mon C a 700 Go de libre et D a 1,5 To.
Je cherche à utiliser Grass (jamais auparavant) mais i.image.mosaic ne semble gérer que 4 fichiers ... certains des miens ont 600 tuiles. D'autres options ou logiciels open source à essayer?
Désolé, nous devons ajouter que nous ne pouvons pas utiliser une mosaïque (ou l'équivalent dans d'autres logiciels) car nous devons créer des zones avec des zones sans données définies comme ecw afin qu'elles puissent être ouvertes dans n'importe quel logiciel SIG et combinées avec une résolution inférieure / ancienne données lorsque les nouvelles données n'existent pas de manière transparente.
Un exemple de l'apparence de certains fichiers mosaïqués dans différents logiciels. Global Mapper / ERDAS sont corrects mais ce n'est pas correct dans arcgis.
--- INFOS ANCIENNES ---
Désolé pour le dessin approximatif. Donc, avoir les zones colorées en 5 zones minimisera les zones sans données dans la plus grande AOI.
Dans arcgis, le code est le suivant (il est exécuté comme un modèle et non en python car je ne peux pas le faire prendre l'entrée tifList).
arcpy.MosaicToNewRaster_management(tifList+";" +mask,RootOutput,"Tile1.tif","PROJCS['GDA_1994_MGA_Zone_55',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',147.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]","16_BIT_UNSIGNED","0.5","3","MAXIMUM","#")
# Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script
# The following inputs are layers or table views: "test2"
arcpy.CopyRaster_management(OutputFile,RootOutput+"Tile1b.tif","#","256","256","NONE","NONE","16_BIT_UNSIGNED")
où tifList devrait être lu à partir d'un fichier csv mais cela n'a pas fonctionné en python donc j'exécute ce qui précède dans un modèle à la place ...
J'ai 1,5 To + d'espace libre sur mon lecteur mais le processus se bloque avec une erreur 9999.
100 carreaux seraient-ils traités? -devrions-nous envisager de diviser davantage les zones?