geotif 添加坐标_利用Geotools来转换影像的坐标系[转]

使用Operations类下的resample方法(重采样)可以解决这个问题，它的方法定义如下：

Coverage org.geotools.coverage.processing.Operations.resample(Coverage source, CoordinateReferenceSystem crs) throws CoverageProcessingException

File file = new File("xxxx.tif");

if(file.exists()){

Reader br = new Reader();

GridCoverage2D old2D = br.getGridCoverage2D(file);

final CoordinateReferenceSystem WGS = CRS.decode("EPSG:3857");

final CoordinateReferenceSystem sourceCRS = old2D.getCoordinateReferenceSystem();

System.out.println(String.format("源坐标系为: %s", sourceCRS.getName()));

GridCoverage2D new2D = (GridCoverage2D) Operations.DEFAULT.resample(old2D, WGS);

System.err.println(String.format("目标坐标系为: %s", new2D.getCoordinateReferenceSystem().getName()))；

}

---------------------

作者：WilsonOnIsland

来源：CSDN

原文：https://blog.csdn.net/u013323965/article/details/79453383

1.计算指定像素坐标点的坡度

public float calcSlope(int cellX, int cellY, PlanarImage image) throws IOException {

DecimalFormat df = new DecimalFormat("#.0000");

final int[] dest = null;

int e = image.getTile(image.XToTileX(cellX), image.YToTileY(cellY)).getPixel(cellX, cellY, dest)[0];

int e1 = image.getTile(image.XToTileX(cellX - 1), image.YToTileY(cellY)).getPixel(cellX - 1, cellY, dest)[0];

int e2 = image.getTile(image.XToTileX(cellX), image.YToTileY(cellY - 1)).getPixel(cellX, cellY - 1, dest)[0];

int e3 = image.getTile(image.XToTileX(cellX + 1), image.YToTileY(cellY)).getPixel(cellX + 1, cellY, dest)[0];

int e4 = image.getTile(image.XToTileX(cellX), image.YToTileY(cellY + 1)).getPixel(cellX, cellY + 1, dest)[0];

int e5 = image.getTile(image.XToTileX(cellX - 1), image.YToTileY(cellY - 1)).getPixel(cellX - 1, cellY - 1,

dest)[0];

int e6 = image.getTile(image.XToTileX(cellX + 1), image.YToTileY(cellY - 1)).getPixel(cellX + 1, cellY - 1,

dest)[0];