Python 匀称';geopandas中的s比例函数返回类似震级的点
我有一个.shp文件,读入geopandas数据帧。我将坐标参考系更改为2163,因为我正在制作一些矩形贴图,希望它们看起来有点正常Python 匀称';geopandas中的s比例函数返回类似震级的点,python,gis,scale,shapely,geopandas,Python,Gis,Scale,Shapely,Geopandas,我有一个.shp文件,读入geopandas数据帧。我将坐标参考系更改为2163,因为我正在制作一些矩形贴图,希望它们看起来有点正常 geo_df = geo.GeoDataFrame.from_file('path to shp files here') geo_df = geo_df.to_crs(epsg=2163) Dataframe如下所示: In [10]: geo_df.geometry Out[10]: 0 POLYGON ((1189879.121395004
geo_df = geo.GeoDataFrame.from_file('path to shp files here')
geo_df = geo_df.to_crs(epsg=2163)
Dataframe如下所示:
In [10]: geo_df.geometry
Out[10]:
0 POLYGON ((1189879.121395004 -1019103.847184072...
1 POLYGON ((1220434.289428635 -1303875.122589418...
2 POLYGON ((1210969.088787247 -1295221.772496042...
3 POLYGON ((1217371.162725744 -1300978.843646188...
现在我在计算一个成对的距离矩阵,但由于这些坐标是以米为单位给出的,所以数字相当大,我遇到了溢出问题。所以我试着这样做:
geo_df.geometry = geo_df.geometry.scale(xfact=1/10000, yfact=1/10000, zfact=1.0)
但是我的多边形的坐标并没有变小(只是移动了一点,这很好)。打印假定已重新缩放的数据帧会产生以下结果:
In [12]: geo_df.geometry.scale(xfact=1/10000, yfact=1/10000, zfact=1.0)
Out[12]:
0 POLYGON ((1197230.99656014 -1026296.564074459,...
1 POLYGON ((1221475.564597901 -1304490.752661498...
2 POLYGON ((1215982.083961291 -1290526.930856638...
3 POLYGON ((1223839.012413891 -1292585.80012823,...
我希望得到这样的东西
Out[12]:
0 POLYGON ((119.723099656014 -102.6296564074459,...
1 POLYGON ((122.1475564597901 -130.4490752661498...
2 POLYGON ((121.5982083961291 -129.0526930856638...
3 POLYGON ((122.3839012413891 -129.258580012823,...
也许我只是不明白应该做什么,尽管这看起来很简单。缩放是相对于中心点定义的。应用的缩放确实使多边形变小,但将中心点保持在多边形边界框的中心。
shapely.affinity.scale
中的origin
关键字参数默认为多边形中心;要将坐标值减少10000倍,可以设置原点=(0,0)
,如下所示:
geo_df.geometry = geo_df.geometry.scale(xfact=1/10000, yfact=1/10000, zfact=1.0, origin=(0, 0))