Question

我只是空间分析的先驱，我陷入了困境。

我有一个犯罪数据集，其中的点数以纬度和经度给出。我有另一个数据集（芝加哥的形状文件），我想使用形状文件中的多边形在地图上绘制所有纬度点。

问题是形状文件包含我不知道的不同格式的多边形信息。我从中检索形状文件 https://data.cityofchicago.org/Facilities-Geographic-Boundaries/Boundaries-Neighborhoods/9wp7-iasj

从上面的下载我使用了Neighborhoods_2012b.shp文件

犯罪数据的纬度经度：

POINT (-87.680162979 41.998718085)
POINT (-87.746717696 41.934629749)

Chicago Shapefile中的多边形形状:(所有都是正值）

POLYGON ((1182322.0429 1876674.730700001, 1182...
POLYGON ((1176452.803199999 1897600.927599996,...

我尝试将纬度和经度信息转换为不同的投影（墨卡托），例如（epsg：3857，epsg：3395），但这些投影给出了正值和负值

epsg:3857:
POINT (-9760511.095493518 5160787.421333898)
POINT (-9767919.932699846 5151192.321624438)

我甚至尝试将所有Lat-Long转换为UTM（使用python UTM库），这有望为我提供所有正面价值，但它仍然看起来不是正确的格式，因为情节的规模非常不同。 / p>

Using UTM python Library (utm.from_latlon)
POINT (4649857.621612935 443669.2483944244)
POINT (4642787.870839979 438095.1726599361)

我不知道如何处理这种情况，有没有办法知道在给定空间点的情况下使用什么类型的投影？

我很乐意提供任何帮助。

Answer 1

prj文件说：

PROJCS [＆＃34; NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet＆＃34;，GEOGCS [＆＃34; GCS_North_American_1983＆＃34;，DATUM [＆＃34; D_North_American_1983＆＃34;，SPHEROID [＆＃34; GRS_1980＆＃34;，6378137.0,298.257222101 ]，PRIMEM [＆＃34;格林威治＆＃34;，0.0]，单位[＆＃34;度＆＃34;，0.0174532925199433]，投影[＆＃34; Transverse_Mercator＆＃34]，PARAMETER [＆＃34; False_Easting＆＃34;，984250.0]，PARAMETER [＆＃34; False_Northing＆＃34;，0.0]，PARAMETER [＆＃34; Central_Meridian＆＃34;， - 88.33333333333333]，PARAMETER [＆＃34; Scale_Factor＆＃34;，0.999975] ，PARAMETER [＆＃34; Latitude_Of_Origin＆＃34;，36.66666666666666]，UNIT [＆＃34; Foot_US＆＃34;，0.3048006096012192]]

我用QGIS打开了图层，但没有直接使用prj文件。但是，使用prj文件的信息，您可以使用CRS选择器来检索它。搜索参数为：NAD83 Illinois East此处。根据prj文件的建议选择Feet中的那个。例如，EPSG = 6455是一个很好的。我想你现在有足够的信息可以继续...

如何确定给定空间点的投影或坐标参考系

1 个答案: