一、基本知识了解
/** * 根据图幅计算经纬度 * 1:100万 A 经6度，纬4度; 1:50万 B 经3度，纬2度; 1:25万 C 经1度30分，纬1度; 1:10万 D 经30分 纬20分; * 1:5万 E 经15分 纬10分; 1：2.5万 F 经3分45秒 纬2分30秒;1：1万 G 经1分52.5秒 纬1分15秒; * 1：5000 H; 1:2000 J; 1:1000 K; 1:500 L; 1:20万 R; * 1:150万 S; 1:250万 T; 1:400万 U; 1:500万 V; 1:600万 W; 1:800万 X;1:1000万 Y * mappableUnit=J50D002002 */ public static MappableUnitHelper doComputerByMappableUnit(String mappableUnitStr){ // 分步计算 // y=h*y1-(h1-1)*y2,x=(l-31)*x1+(l1-1)*x2 // ,y代表纬度，y1代表比例尺100万地图的纬差，x代表经度，y1代表纬差，x1代表经差,y2代表纬差，x2代表经差 //下代表比例尺100万地图的经差。[]为商取整符号，如9.3则取9，h代表图号行号，l代表图号列号，h1代表 // h1代表图号的地形图在基础图幅内所位于的行号，l1代表图号的地形图在基础图幅内所位于的列号 if(StringUtils.isBlank(mappableUnitStr) ||mappableUnitStr.length() != 10 ){ logger.info("GeotoolsUtil.doComputerByMappableUnit.mappableUnit is not valid"); return null ; } //图号行号 String h = mappableUnitStr.substring(0,1) ; //图号列号 double l =Double.valueOf(mappableUnitStr.substring(1,3)) ; //比例尺代号 String mappableCode =mappableUnitStr.substring(3,4) ; //基础图幅内位于的行号 double h1 =Double.valueOf(mappableUnitStr.substring(4,7)) ; //基础图幅内位于的列号 double l1 =Double.valueOf(mappableUnitStr.substring(7,10)) ; //比例尺 Object bilichi = MappableUnitHelper.getMappableUnitMap().get(mappableCode) ; if(bilichi == null ){ return null ; } String xy = MappableUnitHelper.getxYDvalueMap().get(bilichi).toString() ; double x2 = Double.valueOf(xy.toString().split(",")[0]) ; double y2 = Double.valueOf(xy.toString().split(",")[1]) ; Object doubleLine = MappableUnitHelper.getMappableUnitLineMap().get(h) ; if(doubleLine == null){ return null ; } double hDouble = Double.valueOf(doubleLine.toString()) ; double ymax=(hDouble*14400-(h1-1)*y2)/3600 ; double xmin=((l-31)*21600+(l1-1)*x2)/3600; double ymin= ymax-y2/3600; double xmax=xmin+x2/3600; MappableUnitHelper mappableUnit = new MappableUnitHelper() ; mappableUnit.setXmin(xmin); mappableUnit.setXmax(xmax); mappableUnit.setYmin(ymin); mappableUnit.setYmax(ymax); return mappableUnit; } /** * 根据经纬度计算图幅 * 1:100万 A 经6度，纬4度; 1:50万 B 经3度，纬2度; 1:25万 C 经1度30分，纬1度; 1:10万 D 经30分 纬20分; * 1:5万 E 经15分 纬10分; 1：2.5万 F 经3分45秒 纬2分30秒;1：1万 G 经1分52.5秒 纬1分15秒; * 1：5000 H; 1:2000 J; 1:1000 K; 1:500 L; 1:20万 R; * 1:150万 S; 1:250万 T; 1:400万 U; 1:500万 V; 1:600万 W; 1:800万 X;1:1000万 Y * */ public static String changeToMappale(double x,double y,String scale){ // 分步计算 // h=[y/y1]+1;l=[x/x1]+31;h1=y1/y2-[(y/y1)/y2];l1=[(x/x1)/x2]+1; // ,y代表纬度，y1代表比例尺100万地图的纬差，x代表经度，y1代表纬差(1:100万),值为4度，x1代表经差(1:100万)，值为6度,y2代表纬差，x2代表经差 //下代表比例尺100万地图的经差。[]为商取整符号，如9.3则取9，（）为取余，h代表图号行号，l代表图号列号，h1代表 // h1代表图号的地形图在基础图幅内所位于的行号，l1代表图号的地形图在基础图幅内所位于的列号 //比例尺 String zeroLength="000" ; Object xy = MappableUnitHelper.getxYDvalueMap().get(scale); if(xy == null){ return null ; } double x2 = Double.valueOf(xy.toString().split(",")[0]) ; double y2 = Double.valueOf(xy.toString().split(",")[1]) ; String h = String.valueOf((int)Math.floor(y/4)+1) ; int l = (int)Math.floor(x/6)+31; String h1 = String.valueOf((int)(14400/y2-Math.floor(((y%4)*3600/y2)))); String l1 = String.valueOf((int)(Math.floor(x%6*3600/x2)+1)) ; String hkey = MappableUnitHelper.getKey(MappableUnitHelper.getMappableUnitLineMap(),h) ; String scaleCode = MappableUnitHelper.getKey(MappableUnitHelper.getMappableUnitMap(),scale) ; String hk1 = zeroLength.substring(0,zeroLength.length()-h1.length())+h1; String lk1 = zeroLength.substring(0,zeroLength.length()-l1.length())+l1; return hkey+l+scaleCode+hk1+lk1 ; }
package com.unis.dimensional.geotools.utils;
import java.util.HashMap;
import java.util.Map;
public class MappableUnitHelper {
private static Map<String,Object> mappableUnitMap = new HashMap<String,Object>() ;
private static Map<String,Object> xYDvalueMap = new HashMap<String,Object>() ;
private static Map<String,Object> mappableUnitLineMap = new HashMap<String,Object>() ;
private double xmin ;
private double xmax ;
private double ymin ;
private double ymax ;
public double getXmin() {
return xmin;
}
public void setXmin(double xmin) {
this.xmin = xmin;
}
public double getXmax() {
return xmax;
}
public void setXmax(double xmax) {
this.xmax = xmax;
}
public double getYmin() {
return ymin;
}
public void setYmin(double ymin) {
this.ymin = ymin;
}
public double getYmax() {
return ymax;
}
public void setYmax(double ymax) {
this.ymax = ymax;
}
static{
setMappableUnitMap() ;
setXYDvalueMap() ;
setMappableUnitLineMap();
}
private static void setMappableUnitMap(){
mappableUnitMap.put("A", "1:1000000") ;
mappableUnitMap.put("B", "1:500000") ;
mappableUnitMap.put("C", "1:250000") ;
mappableUnitMap.put("D", "1:100000") ;
mappableUnitMap.put("E", "1:50000") ;
mappableUnitMap.put("F", "1:25000") ;
mappableUnitMap.put("G", "1:10000") ;
mappableUnitMap.put("H", "1:5000") ;
}
private static void setMappableUnitLineMap(){
mappableUnitLineMap.put("A", "1") ;
mappableUnitLineMap.put("B", "2") ;
mappableUnitLineMap.put("C", "3") ;
mappableUnitLineMap.put("D", "4") ;
mappableUnitLineMap.put("E", "5") ;
mappableUnitLineMap.put("F", "6") ;
mappableUnitLineMap.put("G", "7") ;
mappableUnitLineMap.put("H", "8") ;
mappableUnitLineMap.put("I", "9") ;
mappableUnitLineMap.put("J", "10") ;
mappableUnitLineMap.put("K", "11") ;
mappableUnitLineMap.put("L", "12") ;
mappableUnitLineMap.put("M", "13") ;
mappableUnitLineMap.put("N", "13") ;
}
//插值，秒
private static void setXYDvalueMap(){
xYDvalueMap.put("1:1000000","21600,14400") ;
xYDvalueMap.put("1:500000", "10800,7200") ;
xYDvalueMap.put("1:250000","5400,3600") ;
xYDvalueMap.put("1:100000","1800,1200") ;
xYDvalueMap.put("1:50000","900,600") ;
xYDvalueMap.put("1:25000","450,300") ;
xYDvalueMap.put("1:10000","225,150") ;
xYDvalueMap.put("1:5000","112.5,75") ;
}
public static Map<String, Object> getMappableUnitMap() {
return mappableUnitMap;
}
public static Map<String, Object> getxYDvalueMap() {
return xYDvalueMap;
}
public static Map<String, Object> getMappableUnitLineMap() {
return mappableUnitLineMap;
}
public static String getKey(Map map, Object value) {
for (Object key : map.keySet()) {
if (map.get(key).equals(value)) {
return key.toString();
}
}
return null ;
}
}