平面投影坐标转经纬度坐标
平面投影坐标转经纬度坐标
本人在进行开发时，使用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进行检验正确。代码还不完整，不断更新中。
更新（2018年1月1日）：
当我把这个类移植到另一台电脑上并添加 gt-geotiff-17.0.jar后，恩，就添加了这一个jar包没报错，运行时报错：Exception in thread “main” java.lang.NoClassDefFoundError: javax/media/jai/I。个人觉得这是jar包不全所致，解决办法是先删除全部的import* 引用，然后一个一个的添加jar包，jar包源码在GitHub上有，有的方法名不止一个jar包，这个需特别注意。
本人原创，转载请附上本文连接。
欢迎技术交流。