经纬度坐标转平面投影坐标

本人在进行开发时，使用geotools的GeoTiffReader类读取卫星影像的顶点坐标，得到的是平面坐标，单位是米，打开arcGIS，加载影像，右下角会有坐标，在视图->数据框属性->常规->单位，可以更改坐标显示单位，一开始默认单位是米，我改成十进制度，如图：。
现想将平面坐标转化为经纬度，在网上看到一片博客：经纬度坐标和投影坐标的转换。参考该代码写类自己的，我的目的正好与博客写的相反。然后我就参考geotools-17的用户文档（
geotools/geotools-17.0-userguide/geotools-17.0/library/referencing/crs.html?highlight=transform）
现摘取部分文档（）：

程序主要靠这俩方法：CRS.findMathTransform()方法定义转换的坐标系，JTS.transform()进行坐标的转换。

package Test02;import org.geotools.coverage.grid.GridCoverage2D;import org.geotools.gce.geotiff.GeoTiffReader;import org.geotools.geometry.jts.JTS;import org.geotools.referencing.CRS;import org.geotools.referencing.crs.DefaultGeographicCRS;import org.opengis.geometry.MismatchedDimensionException;import org.opengis.referencing.FactoryException;import org.opengis.referencing.crs.CoordinateReferenceSystem;import org.opengis.referencing.operation.MathTransform;import org.opengis.referencing.operation.TransformException;import com.vividsolutions.jts.geom.Coordinate;import com.vividsolutions.jts.geom.GeometryFactory;import com.vividsolutions.jts.geom.Point;public class CoordinateConversion02 {    /**     * 平面坐标转经纬度     */    public static double[] convert(double x, double y,CoordinateReferenceSystem crs)             throws FactoryException, MismatchedDimensionException, TransformException {    Coordinate sourceCoord = new Coordinate(x, y);    GeometryFactory geoFactory = new GeometryFactory();    Point sourcePoint = geoFactory.createPoint(sourceCoord);    MathTransform transform = CRS.findMathTransform(crs, DefaultGeographicCRS.WGS84, false);    Point targetPoint = (Point) JTS.transform(sourcePoint, transform);    double[] targetCoord = {targetPoint.getX(), targetPoint.getY()};    return targetCoord;    }    /**     * 测试     */    public static void main( String[] args ) throws Exception    {        double x = -132202.486;        double y = 1657014.274;        //第一种获取方式，imageFile为卫星影像的文件地址，目的是获取影像的投影//        GeoTiffReader reader = new GeoTiffReader(imageFile);//        GridCoverage2D coverage = reader.read(null);//        //获取投影//        CoordinateReferenceSystem crs = coverage.getCoordinateReferenceSystem2D();        //这是第二种获取方式（参考博文），        final String wkt="PROJCS[\"unnamed\","                                 +"GEOGCS[\"WGS 84\","                                       + "DATUM[\"World Geodetic System 1984\", "                                         +"SPHEROID[\"WGS 84\", 6378137.0, 298.257223563, AUTHORITY[\"EPSG\",\"7030\"]],"                                         + "AUTHORITY[\"EPSG\",\"6326\"]], "                                     +  "PRIMEM[\"Greenwich\", 0.0, AUTHORITY[\"EPSG\",\"8901\"]],"                                      + " UNIT[\"degree\", 0.017453292519943295], "                                      + "AXIS[\"Geodetic longitude\", EAST], "                                      +" AXIS[\"Geodetic latitude\", NORTH], "                                      + "AUTHORITY[\"EPSG\",\"4326\"]], "                                   + " PROJECTION[\"Albers_Conic_Equal_Area\"], "                                   +  "PARAMETER[\"central_meridian\", 110.0], "                                   + " PARAMETER[\"latitude_of_origin\", 12.0], "                                    +" PARAMETER[\"standard_parallel_1\", 25.0], "                                   +  "PARAMETER[\"false_easting\", 0.0], "                                   +  "PARAMETER[\"false_northing\", 0.0],"                                     + "PARAMETER[\"standard_parallel_2\", 47.0],"                                     + "UNIT[\"m\", 1.0], "                                    + "AXIS[\"Easting\", EAST], "                                   +  "AXIS[\"Northing\", NORTH]]";        CoordinateReferenceSystem crs = CRS.parseWKT(wkt);        double[] coordinate = convert(x, y,crs);        System.out.println("X: " + coordinate[0] + ", Y: " + coordinate[1]);    }}";

好了，一切完成，用arcGIS进行检验正确。代码还不完整，不断更新中。本人原创，转载请附上本文连接，并通知本人邮箱。邮箱：[e-mail]
欢迎技术交流。