python - 如何确定给定空间点的投影或坐标引用系

Question

我只是 Spatial Analysis 的新手，而且卡在了某个点上。

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

问题是形状文件包含我不知道的不同格式的多边形信息。我从
检索形状文件 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,...

我尝试将纬度和经度信息转换为不同的投影 (Mercator)，例如 (epsg:3857, epsg:3395)，但这些投影给了我正值和负值

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

我什至尝试将所有 Lat-Long 转换为 UTM(使用 python UTM 库)，希望这能给我所有积极的值(value)，但它似乎仍然不是正确的格式，因为绘图的比例非常不同。

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

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

我很乐意提供任何帮助。

Answer 1

prj 文件说:

PROJCS["NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",984250.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-88.33333333333333],PARAMETER["Scale_Factor",0.999975],PARAMETER["Latitude_Of_Origin",36.66666666666666],UNIT["Foot_US",0.3048006096012192]]

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