GeoTools，Java：如何将 shapefile 数据从 DOR 转换为具有纬度和经度坐标的可用多边形答案

【问题标题】：GeoTools, Java: How to convert shapefile data from DOR to usable polygons with Latitude and Longitude coordinatesGeoTools，Java：如何将 shapefile 数据从 DOR 转换为具有纬度和经度坐标的可用多边形
【发布时间】：2018-01-25 06:27:28
【问题描述】：

我关心的 shapefile 可以在这里找到，尽管我确信这可以应用于任何使用基于状态平面坐标的坐标的 shapefile：

http://dor.wa.gov/content/FindTaxesAndRates/stshpdownloads.aspx

我目前正在使用 GeoTools，这是一个 Java 工具包，旨在读取此类 GIS shapefile 并显示有关它们的信息。我能够遍历文件并显示与每个税收管辖区相关的多边形，但这里有一个示例输出行：

Seattle|0.096|MULTIPOLYGON(((1157963.4622609913 862693.0463841558, 
1157982.575083658 862880.4607074857, 1158623.4100880772 
864826.6100303978, ... )))

根据我的研究，这些坐标似乎与任何典型的坐标系都不相符。我期望它们是从称为 State Plane Coordinates 的坐标系中派生出来的，但即使对于那些来说，目前似乎也不是正确的格式。我希望有人可能知道这些值的来源，以及是否有办法让 GeoTools 以更标准化的纬度/经度格式输出它们。

【问题讨论】：

如果我没记错的话，shapefile 会与它们一起存储投影。如果您可以从文件或投影的EPSG-code（投影的全局标识符）中获取投影，您就会自动知道如何处理坐标。
查看 geotools.org 上的 crs 教程，您将了解如何做到这一点

标签： java gis projection geotools spc

【解决方案1】：

我刚刚从该站点下载了最新的 shapefile 并将它们放入 ArcGIS 中，它显示以下内容：

Data Type:  Shapefile Feature Class 
Shapefile:  C:\Users\(...)\Downloads\Cities_17Q3\Cities_17Q3.shp
Geometry Type:  Polygon
Coordinates have Z values:  No 
Coordinates have measures:  No 

Projected Coordinate System:    NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet
Projection: Lambert_Conformal_Conic
False_Easting:  1640416,66666667
False_Northing: 0,00000000
Central_Meridian:   -120,50000000
Standard_Parallel_1:    45,83333333
Standard_Parallel_2:    47,33333333
Latitude_Of_Origin: 45,33333333
Linear Unit:    Foot_US

Geographic Coordinate System:   GCS_North_American_1983_HARN
Datum:  D_North_American_1983_HARN
Prime Meridian:     Greenwich
Angular Unit:   Degree

此外，西雅图的多边形似乎在视觉上排列：

这意味着您正在使用GCS_North_American_1983_HARN 坐标系和Lambert_Conformal_Conic 投影

【讨论】：