Conversion huge raster (50 000 000 cells) in to .xyz

[المشرف العام]

I do have DEM approx. 40x5km with cell size 2x2m, which gives approx. 50 000 000 points.

I need to export it to .xyz (1pt = 1 line in format x y z) file in order to use it in MIKE21.

It seemed like there is loads of tools that will do it (I do have access to Arcinfo 9.3.1 and 10), though I'm stuck for 2 days already. I tried some solutions and all have failed.

I did:1. SA/EXTRACTION/SAMPLE, receiving table, saving it as .dbf, and trying to import it in to Access 2007. I could not import .dbf produced in Arc into Access. Exporting table to .txt gives me always error.

1. Raster Calculator and Sample(). This worked great. Raster converted in 5 minutes producing 1.2GB file, but in order: ZXY not XYZ. MIKE do not allow to change order of columns, so I need to change order of columns. No clue how to do it. Tried Access 2007 again, but failed on import data.
   
2. Script Raster2XYZ ver. 1.1 from ESRI website (in VB). Can not install it properly now (gives an error), I used it before and it worked 6 months ago. However I remember that it worked really slow back then. I am not much of a programmer and can not find the bug.
   
3. Grass GIS and r.out.xyz have bug. Problem with double or float writing(can't remember now).
   
4. I have R script which does export .asc to .xyz. It take 2.5h to export 40 000 000 cells. Output file is 4.3GB. MIKE will not swallow it.
   

Any ideas how to export such raster to .xyz? Or how to reorder file from pt.2 to receive XYZ instead ZXY?

I always used Access to convert data (ie. .dbf to .txt(.xyz)), but now it seems like no option. Is there any other solution?

I do have HP workstation Z600 (Xeon 2.4Ghz, 8GB RAM, GF Quadro FX 3800).

Any sugestions?



أكثر...