码迷,mamicode.com
首页 > 其他好文 > 详细

GeoTools坐标转换(投影转换和仿射变换)

时间:2016-03-26 21:54:34      阅读:1535      评论:0      收藏:0      [点我收藏+]

标签:

GeoTools是在java下的gis开源软件,以下介绍坐标转换的两种方法:投影转换和仿射变换


  • 投影转换

这里以xian80经纬度坐标转xian80,3度分带 111中央经线平面坐标为例

转换函数如下:

1         Point pointXian80 = projectTransform(lon, lat, "EPSG:4610", "EPSG:2382");

 

 1     /**
 2      * 投影转换, lon=经度,lat=纬度,ESPG格式(例):EPSG:4610
 3      */
 4     public static Point projectTransform(double lon, double lat,
 5             String epsgSource, String epsgTarget) throws FactoryException,
 6             MismatchedDimensionException, TransformException {
 7         // 原始坐标点
 8         // PS:通常逻辑上理解经度应该是横坐标x,纬度是y,可是这里经度要填到y,纬度x,否则会报错
 9         Point sourcePoint = JtsHelper.createPoint(lat, lon);
10 
11         // 定义转换前和转换后的投影,可以用ESPG或者wkt
12         // "PROJCS[\"Xian_1980_3_Degree_GK_CM_111E\",GEOGCS[\"GCS_Xian_1980\",DATUM[\"D_Xian_1980\",SPHEROID[\"Xian_1980\",6378140.0,298.257]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Gauss_Kruger\"],PARAMETER[\"False_Easting\",500000.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",111.0],PARAMETER[\"Scale_Factor\",1.0],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]";
13         // CoordinateReferenceSystem mercatroCRS = CRS.parseWKT(strWKTMercator);
14         CoordinateReferenceSystem crsSource = CRS.decode(epsgSource);
15         CoordinateReferenceSystem crsTarget = CRS.decode(epsgTarget);
16         // 投影转换
17         MathTransform transform = CRS.findMathTransform(crsSource, crsTarget);
18         Point pointTarget = (Point) JTS.transform(sourcePoint, transform);
19 
20         return pointTarget;
21     }

关于定义坐标系的说明,GeoTools中定义坐标系有两种方法,一是坐标系的wkt,二是ESPG
坐标系的wkt可以从arcgis的prj文件中过去,用记事本打开prj文件,里面的内容就是该坐标系的wkt格式内容,例如:PROJCS[\"Xian_1980_3_Degree_GK_CM_111E\",GEOGCS[\"GCS_Xian_1980\",DATUM[\"D_Xian_1980\",SPHEROID[\"Xian_1980\",6378140.0,298.257]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Gauss_Kruger\"],PARAMETER[\"False_Easting\",500000.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",111.0],PARAMETER[\"Scale_Factor\",1.0],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]

获取到wkt后,可以通过 CoordinateReferenceSystem mercatroCRS = CRS.parseWKT(strWKTMercator); 获取到坐标系对象

另一种方法是ESPG,意思是每一个坐标系都有一个ESPG标准的号码,查询坐标系的ESPG可通过以下网站 http://spatialreference.org/

PS:当转换的两种坐标系的datum不同,会报错Bursa-Wolf parameters,暂时未解决

 


  • 仿射变换

代码如下:

首先是输入三组参考点,前三个是转换前,后三个是转换后

最后结果是pointReuslt

 1         // 参考点坐标
 2         Coordinate s1 = new Coordinate(429275.549, 2801455.153);
 3         Coordinate s2 = new Coordinate(428110.626, 2792148.620);
 4         Coordinate s3 = new Coordinate(428966.479, 2800016.622);
 5 
 6         Coordinate t1 = new Coordinate(4628.339, 801349.338);
 7         Coordinate t2 = new Coordinate(3515.906, 792036.308);
 8         Coordinate t3 = new Coordinate(4327.381, 799909.069);
 9         // 建立仿射变换对象
10         AffineTransformationBuilder afb = new AffineTransformationBuilder(s1,
11                 s2, s3, t1, t2, t3);
12         AffineTransformation atf = afb.getTransformation();
13 
14         Point pointReuslt = JtsHelper.createPoint(x,y);
15 
16         // 坐标转换
17         pointReuslt.apply(atf);

 

GeoTools坐标转换(投影转换和仿射变换)

标签:

原文地址:http://www.cnblogs.com/cannel/p/5323925.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!