Application: gvSIG desktop

/**********************************************************************

 * $Id: Writetif.java 804 2004-12-30 13:16:13Z igbrotru $

 *

 * Name:     writetif.java

 * Project:  gvSIG - Generalitat Valenciana SIG

 * Purpose:  Ejemplo para la prueba del interfaz java de la libreria 

 *                         gdal de Frank Warmerdam

 * Author:   Nacho Brodin, brodin_ign@gva.es

 *

 **********************************************************************/

import es.gva.cit.jgdal.*;

public class Writetif{

        public static void main(String[] args){

          Gdal migdal=new Gdal();

          Gdal dset_dstno=null;

          GeoTransform gt=null;

          int ancho=0,alto=0;

          GdalRasterBand mirasterband=null;

          GdalRasterBand rband=null;

          int nxsize=0;

          int nysize=0;

          int rastercount=0;

          float total=0;

          try{

            if(args.length==3){

                    migdal.open(args[0],Gdal.GA_ReadOnly);

                    ancho=Integer.parseInt(args[1]);

                    alto=Integer.parseInt(args[2]);

            }else {

                    System.out.println("Params: filename Xsize Ysize"); 

                    System.exit(1);

            }

          }catch(Exception ge){

                        ge.printStackTrace();

                        //...

          }

          try{

                  gt=migdal.getGeoTransform();

                  System.out.println("Origin = ("+gt.adfgeotransform[0]+","+gt.adfgeotransform[3]+")");

                  System.out.println("Pixel Size = ("+gt.adfgeotransform[1]+","+gt.adfgeotransform[5]+")");

          }catch(GdalException e){

                          System.out.println("I can't obtain the array geoTransform for this image");

                        e.printStackTrace();

                        //...........

          }

          try{

                  rastercount=migdal.getRasterCount();

            nxsize = migdal.getRasterXSize();

            nysize = migdal.getRasterYSize();

            System.out.println("N BANDAS="+rastercount+" SIZEX="+nxsize+" SIZEY="+nysize);

            System.out.println("WINDOWX="+ancho+" WINDOWY="+alto);

            System.out.println(migdal.getProjectionRef());

            int ngcp=migdal.getGCPCount();

            System.out.println("NGCP="+ngcp);

          }catch(GdalException ge){

                                                          ge.printStackTrace();

                                                          //...

                                                     }

          try{

                          //Obtenemos el driver y creamos el dataset del destino

                          GdalDriver drv=Gdal.getDriverByName("GTiff");

                          String[] params={"TILED=YES","PHOTOMETRIC=RGB","TFW=WORLDFILE"};

                          dset_dstno=drv.create("salidatif.tif",ancho,alto,rastercount,Gdal.GDT_Byte,null);

                          //A?adimos la informaci?n de proyecci?n (esta es de ejemplo. Hay que formatear la de la im?gen de entrada)

                          OGRSpatialReference oSRS=new OGRSpatialReference();

                          oSRS.setUTM(11,true);

                          oSRS.setWellKnownGeogCS("NAD27");

                          System.out.println("Nueva Proyecci?n ==> "+oSRS.exportToWkt());

                          dset_dstno.setProjection(oSRS.exportToWkt());

                          //A?adimos el vector geoTransform de la imag?n fuente a la de destino

                          dset_dstno.setGeoTransform(gt);

                          byte[] buf=null;

                          for(int iBand=0;iBand<rastercount;iBand++){

                                  System.out.println("Lectura/Escritura banda..."+iBand);

                                  mirasterband=migdal.getRasterBand(iBand+1);

                                  rband=dset_dstno.getRasterBand(iBand+1);

                                  buf = mirasterband.readRaster(0, 0, ancho, alto, ancho, alto, Gdal.GDT_Byte);

                                  long inicio=System.currentTimeMillis();                          

                                  rband.writeRaster(0,0, ancho, alto, buf,Gdal.GDT_Byte);

                                  long fin=System.currentTimeMillis();

                                  total+=(fin-inicio);                 

                          }

                          System.out.println("TIEMPO DE ESCRITURA ="+(total/1000));                          

                          dset_dstno.close();

                          migdal.close();

                          oSRS=null;

                  }catch(Exception e){

                          e.printStackTrace( );

                          //...

                  }

        }  

}