Java各坐标系之间的转换（⾼斯、WGS84经纬度、Web墨卡托、⽡⽚坐标）本⽂整理了⼀些地理坐标系之间的转换（Java代码）
pom依赖
<dependency>
<groupId>com.vividsolutions</groupId>
<artifactId>jts</artifactId>
<version>1.13</version>
</dependency>
<dependency>
<groupId>org.osgeo</groupId>
<artifactId>proj4j</artifactId>
<version>0.1.0</version>
</dependency>
坐标转换⼯具类
package;
import Pixel;
import Tile;
import Coordinate;
import Envelope;
import BasicCoordinateTransform;
import CRSFactory;
import CoordinateReferenceSystem;
import ProjCoordinate;
/
*************************************
*Class Name:GeoTransform
*Description:<；坐标转换⼯具类>
*@since 1.0.0
*************************************/
public class GeoTransform {
/**
* ⾚道半径
*/
private final static double EarthRadius =6378137.0;
/**
* 地球周长
*/
private final static double EarthPerimeter =2* Math.PI * EarthRadius;
/**
* ⽡⽚⼤⼩，默认256
*/
private final static int tileSize =256;
/**
* 初始像素分辨率.
*/
private final static double initialResolution = EarthPerimeter / tileSize;
/
**
* 坐标原点
*/
private final static Coordinate origin =new Coordinate(-EarthPerimeter /2.0, EarthPerimeter /2.0);
private final static BasicCoordinateTransform transform1;
private final static BasicCoordinateTransform transform2;
private final static CRSFactory crsFactory =new CRSFactory();
private final static CoordinateReferenceSystem WGS84CRS = ateFromName("EPSG:4326");
private final static CoordinateReferenceSystem WebMercatorCRS = ateFromName("EPSG:3857");
static{
transform1 =new BasicCoordinateTransform(WGS84CRS, WebMercatorCRS);
transform1 =new BasicCoordinateTransform(WGS84CRS, WebMercatorCRS);
transform2 =new BasicCoordinateTransform(WebMercatorCRS, WGS84CRS);
}
/**
* 缩放级别换算地图分辨率
*
* @param zoom 级别
*/
public double zoomToResolution(int zoom){
return initialResolution / Math.pow(2, zoom);
}
/**
* 经纬度转墨卡托
*
* @param pt 经纬度坐标
* @return 墨卡托坐标
*/
public Coordinate geographic2Mercator(Coordinate pt){
synchronized(transform1){
ProjCoordinate pt1 =new ProjCoordinate(pt.x, pt.y);
ProjCoordinate pt2 =new ProjCoordinate();
return new Coordinate(pt2.x, pt2.y);
}
}
/**
* 墨卡托转经纬度
*
* @param pt 墨卡托坐标
* @return 经纬度坐标
*/
public Coordinate mercator2Geographic(Coordinate pt){
synchronized(transform2){
ProjCoordinate pt1 =new ProjCoordinate(pt.x, pt.y);
ProjCoordinate pt2 =new ProjCoordinate();
return new Coordinate(pt2.x, pt2.y);
}
}
/**
* ⾼斯转经纬度
*
* @param pt ⾼斯坐标
* @param d 度带号(3度带)
* @return 经纬度坐标
*/
public Coordinate gk2Geographic(Coordinate pt,int d){
synchronized(crsFactory){
CoordinateReferenceSystem GKCRS = ateFromParameters("WGS84", String.format("+proj=tmerc +lat_0=0 +lon_0=%d +k=1 +x_0=5 00000 +y_0=0 +ellps=WGS84 +units=m +no_defs", d *3));
BasicCoordinateTransform transform =new BasicCoordinateTransform(GKCRS, WGS84CRS);
ProjCoordinate pt1 =new ProjCoordinate(pt.x, pt.y);
ProjCoordinate pt2 =new ProjCoordinate();
return new Coordinate(pt2.x, pt2.y);
}
}
/**
* 经纬度转⾼斯
*
reference group* @param pt 经纬度坐标
* @return ⾼斯坐标
*/
*/
public Coordinate geographic2GK(Coordinate pt){
synchronized(crsFactory){
int d =(int) Math.floor((pt.y +1.5)/3);
CoordinateReferenceSystem GKCRS = ateFromParameters("WGS84", String.format("+proj=tmerc +lat_0=0 +lon_0=%d +k=1 +x_0=5 00000 +y_0=0 +ellps=WGS84 +units=m +no_defs", d *3));
BasicCoordinateTransform transform =new BasicCoordinateTransform(WGS84CRS, GKCRS);
ProjCoordinate pt1 =new ProjCoordinate(pt.x, pt.y);
ProjCoordinate pt2 =new ProjCoordinate();
return new Coordinate(pt2.x, pt2.y);
}
}
/**
* ⾼斯转web墨卡托
*
* @param pt ⾼斯坐标
* @param d 度带好(3度带)
* @return 墨卡托坐标
*/
public Coordinate gk2Mercator(Coordinate pt,int d){
synchronized(crsFactory){
CoordinateReferenceSystem GKCRS = ateFromParameters("WGS84", String.format("+proj=tmerc +lat_0=0 +lon_0=%d +k=1 +x_0=5 00000 +y_0=0 +ellps=WGS84 +units=m +no_defs", d *3));
BasicCoordinateTransform transform =new BasicCoordinateTransform(GKCRS, WebMercatorCRS);
ProjCoordinate pt1 =new ProjCoordinate(pt.x, pt.y);
ProjCoordinate pt2 =new ProjCoordinate();
return new Coordinate(pt2.x, pt2.y);
}
}
/**
* 墨卡托转像素
*
* @param pt 墨卡托坐标
* @param zoom 缩放级别
* @return 像素坐标
*/
public Pixel mercator2Pixel(Coordinate pt,int zoom){
double res =zoomToResolution(zoom);
Double px =(pt.x - origin.x)/ res;
Double py =-(pt.y - origin.y)/ res;
//System.out.println(px+","+py);
//fixme 精度向下取整
return new Pixel((long) Math.floor(px),(long) Math.floor(py));
}
/
**
* 像素转墨卡托
*
* @param pixel 像素坐标
* @param zoom 缩放级别
* @return 墨卡托坐标
*/
public Coordinate pixel2Mercator(Pixel pixel,int zoom){
double res =zoomToResolution(zoom);
double x = X()* res + origin.x;
double y = origin.y - Y()* res;
return new Coordinate(x, y);
}
/**
* 像素坐标所在⽡⽚
*
* @param pixel 像素坐标
* @return ⽡⽚坐标
*/
*/
public Tile pixelAtTile(Pixel pixel){
long tileX = X()/ tileSize;
long tileY = Y()/ tileSize;
return new Tile(tileX, tileY);
}
/**
* 像素转⽡⽚内像素
*
* @param pixel 像素坐标
* @param tile ⽡⽚坐标
* @return ⽡⽚内像素坐标
*/
public Pixel pixel2Tile(Pixel pixel, Tile tile){
long pX = X()- X()* tileSize;
long pY = Y()- Y()* tileSize;
return new Pixel(pX, pY);
}
/**
* ⽡⽚内像素转像素
*
* @param p ⽡⽚内像素坐标
* @param tile ⽡⽚坐标
* @return 像素坐标
*/
public Pixel tile2Pixel(Pixel p, Tile tile){
long pixelX = p.getX()+ X()* tileSize;
long pixelY = p.getY()+ Y()* tileSize;
return new Pixel(pixelX, pixelY);
}
/**
* 墨卡托转⽡⽚内像素
*
* @param pt 墨卡托坐标
* @param tile ⽡⽚坐标
* @param zoom 缩放级别
* @return ⽡⽚内像素坐标
*/
public Pixel mercator2Tile(Coordinate pt, Tile tile,int zoom){ Pixel p =mercator2Pixel(pt, zoom);
Pixel pixel =pixel2Tile(p, tile);
return pixel;
}
/**
* 经纬度转像素
*
* @param pt 经纬度坐标
* @param zoom 缩放级别
* @return 像素坐标
*/
public Pixel geographic2Pixel(Coordinate pt,int zoom){
Coordinate mpt =geographic2Mercator(pt);
Pixel pixel =mercator2Pixel(mpt, zoom);
return pixel;
}
/**
* 经纬度转⽡⽚内像素
*
* @param pt 经纬度坐标
* @param tile ⽡⽚坐标
* @param zoom 缩放级别
* @return ⽡⽚内像素坐标
*/
public Pixel geographic2Tile(Coordinate pt, Tile tile,int zoom){ Pixel pixel =geographic2Pixel(pt, zoom);
Pixel p =this.pixel2Tile(pixel, tile);
Pixel p =this.pixel2Tile(pixel, tile);
return p;
}
/**
* 像素转经纬度
*
* @param pixel 像素坐标
* @param zoom 缩放级别
* @return 经纬度坐标
*/
public Coordinate pixel2Geographic(Pixel pixel,int zoom){ Coordinate mpt =pixel2Mercator(pixel, zoom);
Coordinate lonlat =ator2Geographic(mpt);
return lonlat;
}
/**
* Tile坐标转换为所在的Tile的矩形
*
* @param tile ⽡⽚坐标
* @return 矩形
*/
public Envelope tile2Envelope(Tile tile){
long px = X()* tileSize;
long py = Y()* tileSize;
Pixel pixel1 =new Pixel(px, py +256);//左下
Pixel pixel2 =new Pixel(px +256, py);//右上
Coordinate sw =pixel2Geographic(pixel1, Z()); Coordinate ne =pixel2Geographic(pixel2, Z());
return new Envelope(sw, ne);
}
/**
* ⽡⽚坐标转换为QuadKey四叉树键值
*
* @param tile ⽡⽚坐标
* @return String QuadKey四叉树键值
*/
public String tile2QuadKey(Tile tile){
long tileX = X();
long tileY = Y();
StringBuilder quadKey =new StringBuilder();
for(int i = Z(); i >0; i--){
char digit ='0';
int mask =1<<(i -1);
if((tileX & mask)!=0){
digit++;
}
if((tileY & mask)!=0){
digit++;
digit++;
}
quadKey.append(digit);
}
String();
}
/**
* QuadKey四叉树键值转换为⽡⽚坐标
*
* @param quadKey QuadKey四叉树键值
* @return ⽡⽚坐标
*/
public Tile quadKey2Tile(String quadKey){
long tileX =0;
long tileY =0;
int levelOfDetail = quadKey.length();
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论