利用gis通過shp文件將經緯度解析為城市名

1. android端通過gps解析出所在的城市名，可以使用Geocoder來完成，但是Geocoder函數和android api的版本有關。2.1以上不能用
2. 通過google map方式如：http://maps.googleapis.com/maps/api/geocode/xml?latlng=40.714224,-73.961452&sensor=true可以完成，也是比較簡單的一種形式，但是公司網關做了限制，不能訪問。
3. 沒有辦法，只能通過gis解析城市名。  
4. 先來了解一下什么是gis：即地理信息系統（Geographic Information System),下面分點來介紹下所用到的開源的java GIS系統。

• 第一個用到的是Geotools，他是一個開源的Java GIS工具包,可利用它來開發符合標准的地理信息系統。Geotools提供了OGC(Open Geospatial Consortium)規范的一個實現來作為他們的開發。
• 第二個用到的是JTS，他是加拿大的 Vivid Solutions 做的一套開放源碼的 Java API。它提供了一套空間數據操作的核心算法,為在兼容OGC標准的空間對象模型中進行基礎的幾何操作提供2D空間謂詞API。提供一個資料鏈接：http://www.vividsolutions.com/jts/jtshome.htm 
• 什么是shp文件：Shapefile文件是美國環境系統研究所（ESRI）所研制的GIS文件系統格式文件，是工業標准的矢量數據文件。 Shapefile將空間特征表中的非拓撲幾何對象和屬性信息存儲在數據集中，特征表中的幾何對象存為以坐標點集表示的圖形文件—SHP文件，Shapefile文件並不含拓撲（Topological）數據結構。一個Shape文件包括三個文件：一個主文件(*.shp)，一個索引文件(*.shx)，和一個dBASE(*.dbf)表。主文件是一個直接存取，變長度記錄的文件，其中每個記錄描述構成一個地理特征（Feature）的所有vertices坐標值。在索引文件中，每條記錄包含對應主文件記錄距離主文件頭開始的偏移量，dBASE表包含SHP文件中每一個Feature的特征屬性，表中幾何記錄和屬性數據之間的一一對應關系是基於記錄數目的ID。在dBASE文件中的屬性記錄必須和主文件中的記錄順序是相同的。圖形數據和屬性數據通過索引號建立一一對應的關系。相應的資料參考鏈接:http://blog.csdn.net/zapzqc/article/details/3123149

1. 進入正題：通過shp文件 ,並返回一個集合。

    // 讀取shp文件,返回所有文件信息
    public static FeatureIterator<SimpleFeature> readSHP()
    {
     ShapefileDataStore shpDataStore = null;
     long time1 = System.currentTimeMillis();
     try
     {
      // shp文件路徑
      shpDataStore = new ShapefileDataStore(new File("E:\\geo\\geo\\region.shp").toURI().toURL());
      shpDataStore.setStringCharset(Charset.forName("GBK"));
      // 文件名稱
      String typeName = shpDataStore.getTypeNames()[0];
      FeatureSource<SimpleFeatureType, SimpleFeature> featureSource = null;
      featureSource = (FeatureSource<SimpleFeatureType, SimpleFeature>) shpDataStore.getFeatureSource(typeName);
      FeatureCollection<SimpleFeatureType, SimpleFeature> result = featureSource.getFeatures();
      itertor = result.features();
     }
     catch (Exception e)
     {
      e.printStackTrace();
     }
     System.out.println(System.currentTimeMillis() - time1 + "讀取shp文件時間");
     return itertor;
    }

    

2. 傳入經緯度信息解析出城市名。

   // 根據個gps坐標獲取對應城市
    public static String readInfo(String poingX, String pointY)
    {
     StringBuffer mBuffer = null;
     String mCityName = "";
     try
     {
      if (null == itertor)
      {
       itertor = readSHP();
      }
      mBuffer = new StringBuffer();
      // 迭代iterator
      while (itertor.hasNext())
      {
       SimpleFeature feature = itertor.next();
       Collection<Property> p = feature.getProperties();
       Iterator<Property> it = p.iterator();
       while (it.hasNext())
       {
        String value = null;
        Property pro = it.next();
        value = pro.getValue().toString();
        mBuffer.append(value + "\r\n");
       }
       count++;
       mList.add(mBuffer.toString());
       mBuffer.delete(0, mBuffer.length());
      }
      itertor.close();
      String data = null;
      String[] arr = null;
      // 如果緩存mGpsBuffer為空，則從新到所有集合中循環。
      if (null != mGpsBuffer)
      {
       WKTReader reader = new WKTReader(geometryFactory);
       String point1 = "POINT(" + poingX + " " + pointY + ")";
       Point point = (Point) reader.read(point1);
       MultiPolygon multiPolygon = (MultiPolygon) reader.read(mGpsBuffer[0]);
       // 傳過來的gps經緯度坐標，為一個點
       Geometry geometry1 = geometryFactory.createGeometry(point);
       // 區域內所有點圍城的一個多邊形
       Geometry geometry2 = geometryFactory.createGeometry(multiPolygon);
       // 如果該點包含在多邊形里面則解析出城市名
       if (geometry2.contains(geometry1))
       {
        mCityName = mGpsBuffer[1];
        System.out.println("cityname1 = " + mGpsBuffer[1]);
        return mCityName;
       }
      }
      //循環所有點集合，判斷點是否在其中的多邊形里面。
      for (int i = 0; i < mList.size(); i++)
      {
       data = mList.get(i);
       arr = data.split("\r\n");
       try
       {
        mGpsBuffer = new String[arr.length];
        WKTReader reader = new WKTReader(geometryFactory);
        String point1 = "POINT(" + poingX + " " + pointY + ")";
        Point point = (Point) reader.read(point1);
        MultiPolygon multiPolygon = (MultiPolygon) reader.read(arr[0]);
        Geometry geometry1 = geometryFactory.createGeometry(point);
        Geometry geometry2 = geometryFactory.createGeometry(multiPolygon);
        if (geometry2.contains(geometry1))
        {
         mCityName = arr[1];
         //將多邊形的信息存入緩存中。
         mGpsBuffer = arr;
         System.out.println("cityname = " + arr[1]);
         break;
        }
       }
       catch (Exception e)
       {
        e.printStackTrace();
       }
      }
     }
     catch (Exception e)
     {
      e.printStackTrace();
      // TODO: handle exception
     }
     return mCityName;
    }

    

3. 至此解析工作已完成。進過測試，解析速度還是可以的，如果再優化效率還會提升。