?? gdalinfo.java
字號:
/****************************************************************************** * $Id: gdalinfo.java 9112 2006-02-02 21:04:35Z collinsb $ * * Name: gdalinfo.java * Project: GDAL SWIG Interface * Purpose: Java port of gdalinfo application * Author: Benjamin Collins, The MITRE Corporation * * * $Log$ * Revision 1.1 2006/02/02 21:04:35 collinsb * Expand support for SWIG Java bindings * * **/import java.util.Enumeration;import java.util.Hashtable;import java.util.Vector;import org.gdal.gdal.Band;import org.gdal.gdal.ColorTable;import org.gdal.gdal.Dataset;import org.gdal.gdal.Driver;import org.gdal.gdal.GCP;import org.gdal.gdal.gdal;import org.gdal.gdalconst.gdalconstConstants;import org.gdal.osr.CoordinateTransformation;import org.gdal.osr.SpatialReference;public class gdalinfo { /************************************************************************/ /* Usage() */ /************************************************************************/ public static void Usage() { System.out .println("Usage: gdalinfo [--help-general] [-mm] [-nogcp] [-nomd] " + "datasetname"); System.exit(1); } /************************************************************************/ /* main() */ /************************************************************************/ public static void main(String[] args) { { Dataset hDataset; Band hBand; int i, iBand; double[] adfGeoTransform = new double[6]; Driver hDriver; Hashtable papszMetadata; boolean bComputeMinMax = false, bSample = false; boolean bShowGCPs = true, bShowMetadata = true; boolean bStats = false; String pszFilename = null; gdal.AllRegister(); if (args.length < 1) { Usage(); System.exit(0); } /* -------------------------------------------------------------------- */ /* Parse arguments. */ /* -------------------------------------------------------------------- */ for (i = 0; i < args.length; i++) { if (args[i].equals("-mm")) bComputeMinMax = true; else if (args[i].equals("-stats")) bStats = true; else if (args[i].equals("-sample")) bSample = true; else if (args[i].equals("-nogcp")) bShowGCPs = false; else if (args[i].equals("-nomd")) bShowMetadata = false; else if (args[i].startsWith("-")) Usage(); else if (pszFilename == null) pszFilename = args[i]; else Usage(); } if (pszFilename == null) Usage(); /* -------------------------------------------------------------------- */ /* Open dataset. */ /* -------------------------------------------------------------------- */ hDataset = gdal.Open(pszFilename, gdalconstConstants.GA_ReadOnly); if (hDataset == null) { System.err .println("GDALOpen failed - " + gdal.GetLastErrorNo()); System.err.println(gdal.GetLastErrorMsg()); //gdal.DumpOpenDatasets( stderr ); //gdal.DestroyDriverManager(); //gdal.DumpSharedList( null ); System.exit(1); } /* -------------------------------------------------------------------- */ /* Report general info. */ /* -------------------------------------------------------------------- */ hDriver = hDataset.GetDriver(); System.out.println("Driver: " + hDriver.getShortName() + "/" + hDriver.getLongName()); System.out.println("Size is " + hDataset.getRasterXSize() + ", " + hDataset.getRasterYSize()); /* -------------------------------------------------------------------- */ /* Report projection. */ /* -------------------------------------------------------------------- */ if (hDataset.GetProjectionRef() != null) { SpatialReference hSRS; String pszProjection; pszProjection = hDataset.GetProjectionRef(); hSRS = new SpatialReference(pszProjection); if (hSRS != null) { String[] pszPrettyWkt = new String[1]; hSRS.ExportToPrettyWkt(pszPrettyWkt, 0); System.out.println("Coordinate System is:"); System.out.println(pszPrettyWkt[0]); //gdal.CPLFree( pszPrettyWkt ); } else System.out.println("Coordinate System is `" + hDataset.GetProjectionRef() + "'"); hSRS.delete(); } /* -------------------------------------------------------------------- */ /* Report Geotransform. */ /* -------------------------------------------------------------------- */ hDataset.GetGeoTransform(adfGeoTransform); { if (adfGeoTransform[2] == 0.0 && adfGeoTransform[4] == 0.0) { System.out.println("Origin = (" + adfGeoTransform[0] + "," + adfGeoTransform[3] + ")"); System.out.println("Pixel Size = (" + adfGeoTransform[1] + "," + adfGeoTransform[5] + ")"); } else System.out.println("GeoTransform ="); System.out.println(" " + adfGeoTransform[0] + ", " + adfGeoTransform[1] + ", " + adfGeoTransform[2]); System.out.println(" " + adfGeoTransform[3] + ", " + adfGeoTransform[4] + ", " + adfGeoTransform[5]); } /* -------------------------------------------------------------------- */ /* Report GCPs. */ /* -------------------------------------------------------------------- */ if (bShowGCPs && hDataset.GetGCPCount() > 0) { System.out.println("GCP Projection = " + hDataset.GetGCPProjection()); int count = 0; Vector GCPs = new Vector(); hDataset.GetGCPs(GCPs); Enumeration e = GCPs.elements(); while (e.hasMoreElements()) { GCP gcp = (GCP) e.nextElement(); System.out.println("GCP[" + (count++) + "]: Id=" + gcp.getId() + ", Info=" + gcp.getInfo()); System.out.println(" (" + gcp.getGCPPixel() + "," + gcp.getGCPLine() + ") (" + gcp.getGCPX() + "," + gcp.getGCPY() + "," + gcp.getGCPZ() + ")"); } } /* -------------------------------------------------------------------- */ /* Report metadata. */ /* -------------------------------------------------------------------- */ papszMetadata = hDataset.GetMetadata_Dict(""); if (bShowMetadata && papszMetadata.size() > 0) { Enumeration keys = papszMetadata.keys(); System.out.println("Metadata:"); while (keys.hasMoreElements()) { String key = (String) keys.nextElement(); System.out.println(" " + key + "=" + papszMetadata.get(key)); } } /* -------------------------------------------------------------------- */ /* Report subdatasets. */ /* -------------------------------------------------------------------- */ papszMetadata = hDataset.GetMetadata_Dict("SUBDATASETS"); if (papszMetadata.size() > 0) { System.out.println("Subdatasets:"); Enumeration keys = papszMetadata.keys(); while (keys.hasMoreElements()) { String key = (String) keys.nextElement(); System.out.println(" " + key + "=" + papszMetadata.get(key)); } } /* -------------------------------------------------------------------- */ /* Report corners. */ /* -------------------------------------------------------------------- */ System.out.println("Corner Coordinates:\n"); GDALInfoReportCorner(hDataset, "Upper Left", 0.0, 0.0); GDALInfoReportCorner(hDataset, "Lower Left", 0.0, hDataset .getRasterYSize()); GDALInfoReportCorner(hDataset, "Upper Right", hDataset .getRasterXSize(), 0.0); GDALInfoReportCorner(hDataset, "Lower Right", hDataset .getRasterXSize(), hDataset.getRasterYSize()); GDALInfoReportCorner(hDataset, "Center", hDataset.getRasterXSize() / 2.0, hDataset.getRasterYSize() / 2.0); /* ==================================================================== */ /* Loop over bands. */ /* ==================================================================== */ for (iBand = 0; iBand < hDataset.getRasterCount(); iBand++) { Double[] pass1 = new Double[1], pass2 = new Double[1]; double[] adfCMinMax = new double[2]; ColorTable hTable; hBand = hDataset.GetRasterBand(iBand + 1); /*if( bSample ) { float[] afSample = new float[10000]; int nCount; nCount = hBand.GetRandomRasterSample( 10000, afSample ); System.out.println( "Got " + nCount + " samples." ); }*/ System.out.println("Band " + (iBand+1) + " Type=" + gdal.GetDataTypeName(hBand.getDataType()) + ", ColorInterp=" + gdal.GetColorInterpretationName(hBand .GetRasterColorInterpretation())); String hBandDesc = hBand.GetDescription(); if (hBandDesc != null && hBandDesc.length() > 0) System.out.println(" Description = " + hBandDesc); hBand.GetMinimum(pass1); hBand.GetMaximum(pass2); if(pass1[0] != null || pass2[0] != null || bComputeMinMax) { System.out.println(" Min=" + pass1[0] + " Max=" + pass2[0]); } if (bComputeMinMax) { hBand.ComputeRasterMinMax(adfCMinMax, 0); System.out.println(" Computed Min/Max=" + adfCMinMax[0] + "," + adfCMinMax[1]); } /*eErr = hBand.GetRasterStatistics( hBand, false, bStats, dfMin, dfMax, dfMean, dfStdDev ); if( eErr == CE_None ) { System.out.println( " Minimum=%.3f, Maximum=%.3f, Mean=%.3f, StdDev=%.3f\n", dfMin, dfMax, dfMean, dfStdDev ); }*/ hBand.GetNoDataValue(pass1); if(pass1[0] != null) { System.out.println(" NoData Value=" + pass1[0]); } if (hBand.GetOverviewCount() > 0) { int iOverview; System.out.println(" Overviews: "); for (iOverview = 0; iOverview < hBand.GetOverviewCount(); iOverview++) { Band hOverview; if (iOverview != 0) System.out.print(", "); hOverview = hBand.GetOverview(iOverview); System.out.print(hOverview.getXSize() + "x" + hOverview.getYSize()); } System.out.println(""); } /*if( GDALHasArbitraryOverviews( hBand ) ) { System.out.println( " Overviews: arbitrary\n" ); }*/ /*if( strlen(GDALGetRasterUnitType(hBand)) > 0 ) { System.out.println( " Unit Type: %s\n", GDALGetRasterUnitType(hBand) ); }*/ /*if( GDALGetRasterCategoryNames(hBand) != null ) { String[][] papszCategories = GDALGetRasterCategoryNames(hBand); int i; System.out.println( " Categories:\n" ); for( i = 0; papszCategories[i] != null; i++ ) System.out.println( " %3d: %s\n", i, papszCategories[i] ); }*/ hBand.GetScale(pass1); if(pass1[0] != null) { System.out.print(" Offset: " + pass1[0]); } hBand.GetOffset(pass1); if(pass1[0] != null) { System.out.println(", Scale:" + pass1[0]); } papszMetadata = hBand.GetMetadata_Dict(""); if( bShowMetadata && papszMetadata.size() > 0 ) { Enumeration keys = papszMetadata.keys(); System.out.println("Metadata:"); while (keys.hasMoreElements()) { String key = (String) keys.nextElement(); System.out.println(" " + key + "=" + papszMetadata.get(key)); } } if (hBand.GetRasterColorInterpretation() == gdalconstConstants.GCI_PaletteIndex && (hTable = hBand.GetRasterColorTable()) != null) { int count; System.out.println(" Color Table (" + gdal.GetPaletteInterpretationName(hTable .GetPaletteInterpretation()) + " with " + hTable.GetCount() + " entries)"); for (count = 0; count < hTable.GetCount(); count++) { System.out.println(" " + count + ": " + hTable.GetColorEntry(count)); } } } hDataset.delete(); //CSLDestroy( argv ); //GDALDumpOpenDatasets( stderr ); //gdal.DestroyDriverManager(); //CPLDumpSharedList( null ); //CPLCleanupTLS(); System.exit(0); } } /************************************************************************/ /* GDALInfoReportCorner() */ /************************************************************************/ static boolean GDALInfoReportCorner(Dataset hDataset, String corner_name, double x, double y) { double dfGeoX, dfGeoY; String pszProjection; double[] adfGeoTransform = new double[6]; CoordinateTransformation hTransform = null; System.out.print(corner_name + " "); /* -------------------------------------------------------------------- */ /* Transform the point into georeferenced coordinates. */ /* -------------------------------------------------------------------- */ hDataset.GetGeoTransform(adfGeoTransform); { pszProjection = hDataset.GetProjectionRef(); dfGeoX = adfGeoTransform[0] + adfGeoTransform[1] * x + adfGeoTransform[2] * y; dfGeoY = adfGeoTransform[3] + adfGeoTransform[4] * x + adfGeoTransform[5] * y; } if (adfGeoTransform[0] == 0 && adfGeoTransform[1] == 0 && adfGeoTransform[2] == 0 && adfGeoTransform[3] == 0 && adfGeoTransform[4] == 0 && adfGeoTransform[5] == 0) { System.out.println("(" + x + "," + y + ")"); return false; } /* -------------------------------------------------------------------- */ /* Report the georeferenced coordinates. */ /* -------------------------------------------------------------------- */ System.out.print("(" + dfGeoX + "," + dfGeoY + ") "); /* -------------------------------------------------------------------- */ /* Setup transformation to lat/long. */ /* -------------------------------------------------------------------- */ if (pszProjection != null && pszProjection.length() > 0) { SpatialReference hProj, hLatLong = null; hProj = new SpatialReference(pszProjection); if (hProj != null) hLatLong = hProj.CloneGeogCS(); if (hLatLong != null) { //CPLPushErrorHandler( gdalconstConstants.CPLQuietErrorHandler ); hTransform = new CoordinateTransformation(hProj, hLatLong); //CPLPopErrorHandler(); hLatLong.delete(); } if (hProj != null) hProj.delete(); } /* -------------------------------------------------------------------- */ /* Transform to latlong and report. */ /* -------------------------------------------------------------------- */ if (hTransform != null) { double[] transPoint = new double[3]; hTransform.TransformPoint(transPoint, dfGeoX, dfGeoY, 0); System.out.print("(" + gdal.DecToDMS(transPoint[0], "Long", 2)); System.out .print("," + gdal.DecToDMS(transPoint[1], "Lat", 2) + ")"); } if (hTransform != null) hTransform.delete(); System.out.println(""); return true; }}?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -