JAVA使用gdal裁切

//裁切
    @Test
public void m2(){
    String srcFile="E:\\jinzita\\caiqie\\polygon11.tif";//原始数据
    String dstFile="E:\\jinzita\\caiqie\\clip2.tif";//裁切生成数据
    String shpFile="E:\\jinzita\\caiqie\\geojsonshp\\fhfh_3.shp";//shp文件
    String shpLayerName="fhfh_3";//shp文件名字
    gdal.AllRegister();
    // 读取要切的原图
    Dataset srcDs = gdal.Open(srcFile, gdalconstConstants.GA_ReadOnly);
    if (srcDs == null) {
        System.err.println("GDALOpen failed - " + gdal.GetLastErrorNo());
        System.err.println(gdal.GetLastErrorMsg());
        return;
    }
    // 构建warp参数
    Vector options = new Vector();
    options.add("-of");
    options.add("GTiff");
    options.add("-cutline");//按照矢量图像线裁剪
    options.add(shpFile);
    options.add("-cl");
    options.add(shpLayerName);
    options.add("-crop_to_cutline"); //裁剪后图层范围使用矢量图层大小
    options.add("-dstnodata");
    options.add("-10001.0"); //nulldata的值
    WarpOptions wo = new WarpOptions(options);
    Dataset warp = gdal.Warp(dstFile, new Dataset[]{srcDs}, wo,new buildOverViewCallBack());
    // 获取切图的原点坐标信息
    double[] warpTransform = warp.GetGeoTransform();
    System.out.println("warp原点坐标   = " + warpTransform[0] + "," + warpTransform[3]);
    System.out.println("warp像素坐标差 = " + warpTransform[1] + "," + warpTransform[5]);
    srcDs.delete();
}
    //进度回调
    static class buildOverViewCallBack extends ProgressCallback {
        @Override
        public int run(double dfComplete, String pszMessage) {
            if (StringUtils.isNotBlank(pszMessage)){
                System.out.println(pszMessage);
            }
            System.out.printf("%.2f%n",dfComplete*100);
            return 1;
        }
    }

原文地址：http://www.cnblogs.com/zcxxcvbn/p/16898603.html