日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程语言 > python >内容正文

python

python shp地理坐标系转换平面坐标系-(涉及geopandas安装)

發布時間:2024/1/1 python 41 豆豆
生活随笔 收集整理的這篇文章主要介紹了 python shp地理坐标系转换平面坐标系-(涉及geopandas安装) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

首先需要正確安裝geopandas,步驟請參考:https://geoffboeing.com/2014/09/using-geopandas-windows/

話不多說,直接上代碼。例子寫的是4326轉2383,想要別的坐標系的就去代碼里提到的網站里找對應編碼對應的proj4字符串,替換進去。(長得類似于+proj=tmerc +lat_0=0.....一長串字符?)

我想了一下,之所以需要這樣寫,而不是直接用epsg=2383的原因是,pyproj里識別不出來2383是什么坐標系,所以只能直接輸入這個坐標系的proj4字符串。

import pandas as pd import numpy as np import geopandas as gpddef projTranslate(inputdir): ''' 功能:坐標系4326轉2383輸入:原始shp路徑輸出:轉換后的 shp,可以直接用于繪圖注意:geopandas依賴的pyproj新版本有些問題,網上有人用1.9.6版本的可以轉換成功,但是1.9.6版本的不支持python3.8以上。解決辦法:在to_crs的時候需要去 https://spatialreference.org/ 上找到坐標系對應的 proj4格式的字符串這樣才不會報錯。'''df=gpd.read_file(inputdir)df.set_crs(epsg=4326)df=df.to_crs("+proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=500000 +y_0=0 +a=6378140 +b=6356755.288157528 +units=m +no_defs ")return df

?

總結

以上是生活随笔為你收集整理的python shp地理坐标系转换平面坐标系-(涉及geopandas安装)的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。