天天看點

JAVA實作将GeoHash轉化為對應的經緯度坐标

轉載請注明出處:http://blog.csdn.net/xiaojimanman/article/details/50568428

http://www.llwjy.com/blogdetail/fc484929cff1efac413ec2524680c5d7.html

個人部落格站已經上線了,網址 www.llwjy.com ~歡迎各位吐槽~

-------------------------------------------------------------------------------------------------

      在部落格JAVA實作空間索引編碼(GeoHash)中介紹了什麼是GeoHash以及如何将坐标轉化為GeoHash,這篇部落格就介紹下,如何将GeoHash轉化為對應區域中心點的坐标。

Base32的對應關系

      在上次實作中介紹了數值和base32的對應關系,用數組儲存這種對應關系;現在我們需要找到一個資料結構,可以快速的查找base32字元和數值對應關系,這樣我們可以用map來儲存這種關系。

private static final char[] CHARS = {'0', '1', '2', '3', '4', '5', '6', '7', 
			'8', '9', 'b', 'c', 'd', 'e', 'f', 'g', 'h', 'j', 'k', 'm', 'n', 
			'p', 'q', 'r', 's', 't', 'u', 'v', 'w', 'x', 'y', 'z'};
private static HashMap<Character, Integer> CHARSMAP;

static {
	CHARSMAP = new HashMap<Character, Integer>();
	for (int i = 0; i < CHARS.length; i++) {
		CHARSMAP.put(CHARS[i], i);
	}
}
           

将數值轉化為二進制字元串

      當我們擷取每個字元對應的數字之後,需要将數字轉化為二進制字元串,我們可以直接調用Integer.toBinaryString方法來實作,由于我們要保證轉換之後的二進制要是5位,對于小于16的,我們需要在前部補0,為了實作這個功能,我們對數字統一加上32,在轉化成二進制字元串之後,再把第一位的1移除即可。

private String getBase32BinaryString(int i) {
	if (i < 0 || i > 31) {
		return null;
	}
	String str = Integer.toBinaryString(i + 32);
	return str.substring(1);
}
           

将GeoHash轉化為二進制字元串

      上一步是将數值轉化為五位的二進制字元串,下面我們就介紹下如何将GeoHash字元串轉化為對應的二進制字元串。通過上面介紹的Base32字元和數字的對應關系,我們可以知道每個字元對應的數字,然後将這個數字轉化為對應的五位二進制字元串,我們再把這些字元串拼接起來,就構成了一個完整的二進制字元串。

private String getGeoHashBinaryString(String geoHash) {
	if (geoHash == null || "".equals(geoHash)) {
		return null;
	}
	StringBuffer sb = new StringBuffer();
	for (int i = 0; i < geoHash.length(); i++) {
		char c = geoHash.charAt(i);
		if (CHARSMAP.containsKey(c)) {
			String cStr = getBase32BinaryString(CHARSMAP.get(c));
			if (cStr != null) {
				sb.append(cStr);
			}
		}
	}
	return sb.toString();
}
           

将GeoHash二進制字元串拆分成經緯度二進制字元串

      在上次部落格中介紹了GeoHash二進制字元串是由經緯度的二進制字元串整合而成的,其中奇數為是緯度,偶數為是經度,進而我們隻需要再把GeoHash字元串拆開即可。

StringBuffer lat = new StringBuffer();
StringBuffer lng = new StringBuffer();
for (int i = 0; i < geoHashBinaryStr.length(); i++) {
	if (i % 2 != 0) {
		lat.append(geoHashBinaryStr.charAt(i));
	} else {
		lng.append(geoHashBinaryStr.charAt(i));
	}
}
           

将經緯度二進制轉化為中點對應的經緯度

      上次部落格介紹如何将一個經緯度轉化為對應的二進制,就是看經緯度在對應的那個區域,這裡我們也來對其做你運算,如果二進制字元串第一位為1,則說明經緯度在大的區域,如果為0,則在小的區域,然後我們再對所在的區域進行分割,再來判斷二進制字元串的第二位是1還是0,以此類推,最後傳回最後一個區域中點的值。

private double getGeoHashMid(String binaryStr, double min, double max) {
	if (binaryStr == null || binaryStr.length() < 1) {
		return (min + max) / 2.0;
	}
	if (binaryStr.charAt(0) == '1') {
		return getGeoHashMid(binaryStr.substring(1), (min + max) / 2.0, max);
	} else {
		return getGeoHashMid(binaryStr.substring(1), min, (min + max) / 2.0);
	}
}
           

運作結果檢驗

JAVA實作将GeoHash轉化為對應的經緯度坐标

完整代碼

      對于DistanceUtil 和 LocationBean幾乎沒有做改變,這裡就不再介紹他們的源碼。

/**  
 *@Description: GeoHash實作經緯度的轉化
 */ 
package com.lulei.geo;  

import java.util.ArrayList;
import java.util.Arrays;
import java.util.HashMap;
import java.util.List;

import com.lulei.geo.bean.LocationBean;
import com.lulei.util.JsonUtil;
  
public class GeoHash {
	private LocationBean location;
	/**
	 * 1 2500km;2 630km;3 78km;4 30km
	 * 5 2.4km; 6 610m; 7 76m; 8 19m
	 */
	private int hashLength = 8; //經緯度轉化為geohash長度
	private int latLength = 20; //緯度轉化為二進制長度
	private int lngLength = 20; //經度轉化為二進制長度
	
	private double minLat;//每格緯度的機關大小
	private double minLng;//每個經度的機關大小
	private static final char[] CHARS = {'0', '1', '2', '3', '4', '5', '6', '7', 
				'8', '9', 'b', 'c', 'd', 'e', 'f', 'g', 'h', 'j', 'k', 'm', 'n', 
				'p', 'q', 'r', 's', 't', 'u', 'v', 'w', 'x', 'y', 'z'};
	private static HashMap<Character, Integer> CHARSMAP;
	
	static {
		CHARSMAP = new HashMap<Character, Integer>();
		for (int i = 0; i < CHARS.length; i++) {
			CHARSMAP.put(CHARS[i], i);
		}
	}
	
	public GeoHash(double lat, double lng) {
		location = new LocationBean(lat, lng);
		setMinLatLng();
	}
	
	public int gethashLength() {
		return hashLength;
	}
	
	/**
	 * @Author:lulei  
	 * @Description: 設定經緯度的最小機關
	 */
	private void setMinLatLng() {
		minLat = LocationBean.MAXLAT - LocationBean.MINLAT;
		for (int i = 0; i < latLength; i++) {
			minLat /= 2.0;
		}
		minLng = LocationBean.MAXLNG - LocationBean.MINLNG;
		for (int i = 0; i < lngLength; i++) {
			minLng /= 2.0;
		}
	}
	
	/**
	 * @return
	 * @Author:lulei  
	 * @Description: 求所在坐标點及周圍點組成的九個
	 */
	public List<String> getGeoHashBase32For9() {
		double leftLat = location.getLat() - minLat;
		double rightLat = location.getLat() + minLat;
		double upLng = location.getLng() - minLng;
		double downLng = location.getLng() + minLng;
		List<String> base32For9 = new ArrayList<String>();
		//左側從上到下 3個
		String leftUp = getGeoHashBase32(leftLat, upLng);
		if (!(leftUp == null || "".equals(leftUp))) {
			base32For9.add(leftUp);
		}
		String leftMid = getGeoHashBase32(leftLat, location.getLng());
		if (!(leftMid == null || "".equals(leftMid))) {
			base32For9.add(leftMid);
		}
		String leftDown = getGeoHashBase32(leftLat, downLng);
		if (!(leftDown == null || "".equals(leftDown))) {
			base32For9.add(leftDown);
		}
		//中間從上到下 3個
		String midUp = getGeoHashBase32(location.getLat(), upLng);
		if (!(midUp == null || "".equals(midUp))) {
			base32For9.add(midUp);
		}
		String midMid = getGeoHashBase32(location.getLat(), location.getLng());
		if (!(midMid == null || "".equals(midMid))) {
			base32For9.add(midMid);
		}
		String midDown = getGeoHashBase32(location.getLat(), downLng);
		if (!(midDown == null || "".equals(midDown))) {
			base32For9.add(midDown);
		}
		//右側從上到下 3個
		String rightUp = getGeoHashBase32(rightLat, upLng);
		if (!(rightUp == null || "".equals(rightUp))) {
			base32For9.add(rightUp);
		}
		String rightMid = getGeoHashBase32(rightLat, location.getLng());
		if (!(rightMid == null || "".equals(rightMid))) {
			base32For9.add(rightMid);
		}
		String rightDown = getGeoHashBase32(rightLat, downLng);
		if (!(rightDown == null || "".equals(rightDown))) {
			base32For9.add(rightDown);
		}
		return base32For9;
	}

	/**
	 * @param length
	 * @return
	 * @Author:lulei  
	 * @Description: 設定經緯度轉化為geohash長度
	 */
	public boolean sethashLength(int length) {
		if (length < 1) {
			return false;
		}
		hashLength = length;
		latLength = (length * 5) / 2;
		if (length % 2 == 0) {
			lngLength = latLength;
		} else {
			lngLength = latLength + 1;
		}
		setMinLatLng();
		return true;
	}
	
	/**
	 * @return
	 * @Author:lulei  
	 * @Description: 擷取經緯度的base32字元串
	 */
	public String getGeoHashBase32() {
		return getGeoHashBase32(location.getLat(), location.getLng());
	}
	
	/**
	 * @param lat
	 * @param lng
	 * @return
	 * @Author:lulei  
	 * @Description: 擷取經緯度的base32字元串
	 */
	private String getGeoHashBase32(double lat, double lng) {
		boolean[] bools = getGeoBinary(lat, lng);
		if (bools == null) {
			return null;
		}
		StringBuffer sb = new StringBuffer();
		for (int i = 0; i < bools.length; i = i + 5) {
			boolean[] base32 = new boolean[5];
			for (int j = 0; j < 5; j++) {
				base32[j] = bools[i + j];
			}
			char cha = getBase32Char(base32);
			if (' ' == cha) {
				return null;
			}
			sb.append(cha);
		}
		return sb.toString();
	}
	
	/**
	 * @param base32
	 * @return
	 * @Author:lulei  
	 * @Description: 将五位二進制轉化為base32
	 */
	private char getBase32Char(boolean[] base32) {
		if (base32 == null || base32.length != 5) {
			return ' ';
		}
		int num = 0;
		for (boolean bool : base32) {
			num <<= 1;
			if (bool) {
				num += 1;
			}
		}
		return CHARS[num % CHARS.length];
	}
	
	/**
	 * @param i
	 * @return
	 * @Author:lulei  
	 * @Description: 将數字轉化為二進制字元串
	 */
	private String getBase32BinaryString(int i) {
		if (i < 0 || i > 31) {
			return null;
		}
		String str = Integer.toBinaryString(i + 32);
		return str.substring(1);
	}
	
	/**
	 * @param geoHash
	 * @return
	 * @Author:lulei  
	 * @Description: 将geoHash轉化為二進制字元串
	 */
	private String getGeoHashBinaryString(String geoHash) {
		if (geoHash == null || "".equals(geoHash)) {
			return null;
		}
		StringBuffer sb = new StringBuffer();
		for (int i = 0; i < geoHash.length(); i++) {
			char c = geoHash.charAt(i);
			if (CHARSMAP.containsKey(c)) {
				String cStr = getBase32BinaryString(CHARSMAP.get(c));
				if (cStr != null) {
					sb.append(cStr);
				}
			}
		}
		return sb.toString();
	}
	
	/**
	 * @param geoHash
	 * @return
	 * @Author:lulei  
	 * @Description: 傳回geoHash 對應的坐标
	 */
	public LocationBean getLocation(String geoHash) {
		String geoHashBinaryStr = getGeoHashBinaryString(geoHash);
		if (geoHashBinaryStr == null) {
			return null;
		}
		StringBuffer lat = new StringBuffer();
		StringBuffer lng = new StringBuffer();
		for (int i = 0; i < geoHashBinaryStr.length(); i++) {
			if (i % 2 != 0) {
				lat.append(geoHashBinaryStr.charAt(i));
			} else {
				lng.append(geoHashBinaryStr.charAt(i));
			}
		}
		double latValue = getGeoHashMid(lat.toString(), LocationBean.MINLAT, LocationBean.MAXLAT);
		double lngValue = getGeoHashMid(lng.toString(), LocationBean.MINLNG, LocationBean.MAXLNG);
		LocationBean location = new LocationBean(latValue, lngValue);
		location.setGeoHash(geoHash);
		return location;
	}
	
	/**
	 * @param binaryStr
	 * @param min
	 * @param max
	 * @return
	 * @Author:lulei  
	 * @Description: 傳回二進制對應的中間值
	 */
	private double getGeoHashMid(String binaryStr, double min, double max) {
		if (binaryStr == null || binaryStr.length() < 1) {
			return (min + max) / 2.0;
		}
		if (binaryStr.charAt(0) == '1') {
			return getGeoHashMid(binaryStr.substring(1), (min + max) / 2.0, max);
		} else {
			return getGeoHashMid(binaryStr.substring(1), min, (min + max) / 2.0);
		}
	}
	
	/**
	 * @param lat
	 * @param lng
	 * @return
	 * @Author:lulei  
	 * @Description: 擷取坐标的geo二進制字元串
	 */
	private boolean[] getGeoBinary(double lat, double lng) {
		boolean[] latArray = getHashArray(lat, LocationBean.MINLAT, LocationBean.MAXLAT, latLength);
		boolean[] lngArray = getHashArray(lng, LocationBean.MINLNG, LocationBean.MAXLNG, lngLength);
		return merge(latArray, lngArray);
	}
	
	/**
	 * @param latArray
	 * @param lngArray
	 * @return
	 * @Author:lulei  
	 * @Description: 合并經緯度二進制
	 */
	private boolean[] merge(boolean[] latArray, boolean[] lngArray) {
		if (latArray == null || lngArray == null) {
			return null;
		}
		boolean[] result = new boolean[lngArray.length + latArray.length];
		Arrays.fill(result, false);
		for (int i = 0; i < lngArray.length; i++) {
			result[2 * i] = lngArray[i];
		}
		for (int i = 0; i < latArray.length; i++) {
			result[2 * i + 1] = latArray[i];
		}
		return result;
	}
	
	/**
	 * @param value
	 * @param min
	 * @param max
	 * @return
	 * @Author:lulei  
	 * @Description: 将數字轉化為geohash二進制字元串
	 */
	private boolean[] getHashArray(double value, double min, double max, int length) {
		if (value < min || value > max) {
			return null;
		}
		if (length < 1) {
			return null;
		}
		boolean[] result = new boolean[length];
		for (int i = 0; i < length; i++) {
			double mid = (min + max) / 2.0;
			if (value > mid) {
				result[i] = true;
				min = mid;
			} else {
				result[i] = false;
				max = mid;
			}
		}
		return result;
	}
	

	public static void main(String[] args) {
		// TODO Auto-generated method stub 
		GeoHash g = new GeoHash(40.221227, 116.24875);
		String geoHash = g.getGeoHashBase32();
		System.out.println(geoHash);
		LocationBean bean = g.getLocation(geoHash);
		System.out.println(JsonUtil.parseJson(bean));
		System.out.println(new GeoHash(bean.getLat(), bean.getLng()).getGeoHashBase32());
		System.out.println(DistanceUtil.getDistance(bean.getLat(), bean.getLng(), bean.getLat() - g.minLat, bean.getLng() - g.minLng));
	}

}
           

-------------------------------------------------------------------------------------------------

小福利

-------------------------------------------------------------------------------------------------

      個人在極客學院上《Lucene案例開發》課程已經上線了,歡迎大家吐槽~

第一課:Lucene概述

第二課:Lucene 常用功能介紹

第三課:網絡爬蟲

第四課:資料庫連接配接池

第五課:小說網站的采集

第六課:小說網站資料庫操作

第七課:小說網站分布式爬蟲的實作