开发者社区

文档建议反馈控制台

最新优惠活动

文章/答案/技术大牛

发布

在Python中将lon/lat转换为tmerc

在Python中将经纬度（lon/lat）转换为tmerc（Transverse Mercator）投影坐标可以使用pyproj库来实现。

pyproj是一个用于地理坐标转换的Python库，它提供了丰富的投影转换功能。下面是一个示例代码，演示如何使用pyproj将经纬度转换为tmerc投影坐标：

import pyproj

def lon_lat_to_tmerc(lon, lat):
    # 定义lon/lat的坐标系为WGS84
    wgs84 = pyproj.CRS('EPSG:4326')
    # 定义tmerc投影坐标系
    tmerc = pyproj.CRS('EPSG:32649')  # 32649为UTM Zone 49N的EPSG代码，可以根据实际情况选择合适的投影坐标系
    
    # 创建投影转换器
    transformer = pyproj.Transformer.from_crs(wgs84, tmerc, always_xy=True)
    
    # 进行坐标转换
    x, y = transformer.transform(lon, lat)
    
    return x, y

# 示例使用
lon = 116.3975
lat = 39.9086
x, y = lon_lat_to_tmerc(lon, lat)
print(f"lon/lat: ({lon}, {lat})")
print(f"tmerc: ({x}, {y})")

在上述示例代码中，我们首先定义了源坐标系（WGS84）和目标坐标系（tmerc）。然后，通过pyproj.Transformer.from_crs()方法创建了一个投影转换器。最后，使用转换器的transform()方法将经纬度转换为tmerc投影坐标。

需要注意的是，示例中使用的目标坐标系是UTM Zone 49N（EPSG:32649），你可以根据实际需求选择合适的投影坐标系。

关于腾讯云相关产品和产品介绍链接地址，由于要求不能提及具体的云计算品牌商，我无法提供相关链接。但你可以通过搜索引擎或腾讯云官方网站查找与地理信息相关的云服务产品。

页面内容是否对你有帮助？

有帮助

没帮助

相关·内容

python json转字符串_在python中将json转换为字符串

在python中将json转换为字符串时，请尝试使用str()和json.dumps()。

10.2K3 0

在 Python 中将数值变量转换为分类变量

这篇文章是今天发布的CTGAN的补充，我们可以使用pandas的cut函数将数据进行离散化、将连续变量进行分段汇总，这比写自定义函数要简单的多。

2.1K1 0

Basemap系列教程：使用shapefiles绘制地图

., resolution='i', projection='tmerc', lat_0 = 39.5, lon_0 = 1) map.drawmapboundary(fill_color...如果不是，可以使用 ogr2ogr 进行转换元素必须只有2维此例只有在元素是 polygons 或 polylines 时有效如图所示，结果仅是 polygons（polylines）的边界。...,urcrnrlat=43., resolution='i', projection='tmerc', lat_0 = 39.5, lon_0 = 1) map.drawmapboundary...,urcrnrlat=43., resolution='i', projection='tmerc', lat_0 = 39.5, lon_0 = 1) map.drawmapboundary...,urcrnrlat=43., resolution='i', projection='tmerc', lat_0 = 39.5, lon_0 = 1) map.drawmapboundary

4.7K2 0

Basemap大区域的UTM遥感图可视化（叠加矢量）

然后就要说说basemap了，写的真的不知道该怎么形容，突然羡慕起R里很多代码写的清清楚楚，心疼python3秒。。。basemap在画图上有些局限，从源码到应用，好多问题！就比如源码中： ?...其次是应用，太多东西要吐槽了，反正最想说的就是，我在ARCGIS里明明做的很好的东西，basemap就是不认！ ?...latst=np.linspace(x0,x1,ncols) lonst=np.linspace(y0,y1,nrows) lont, latt = np.meshgrid(latst, lonst) lon..._0t=(x0+x1)/2 lat_0t=(y0+y1)/2 map = Basemap(projection = 'tmerc',resolution='i',rsphere=(6378137.00,6356752.3142,298.257223563...),lat_0 = 33.86499240539949, lon_0 = 122.20569578770636,width=3915000,height=3600000)#我试了好多方法，最后是tmerc

8973 0

Basemap系列教程之基本函数

在地图上绘制一个点通常使用 plot 方法在地图上添加一个点： from mpl_toolkits.basemap import Basemap import matplotlib.pyplot as...plt map = Basemap(projection='ortho', lat_0=0, lon_0=0) map.drawmapboundary(fill_color...', lat_0=0, lon_0=3, llcrnrlon=1.819757266426611, llcrnrlat...', lat_0=0, lon_0=3, llcrnrlon=1.819757266426611, llcrnrlat...', lat_0=0, lon_0=3, llcrnrlon=1.819757266426611, llcrnrlat

2.9K1 0

Basemap工具函数(1)

在地理坐标系中使用会有一个问题，就是比例尺无法放到地图外 lon0 和 lat0 表示要计算比例尺的点 length 表示比例尺上呈现的千米数 barstyle 表示比例尺的类型。...,urcrnrlat=44., resolution='i', projection='tmerc', lat_0 = 39.5, lon_0 = -3.25) map.drawmapboundary...（python 2.7 ，basemap 1.0.7） ? gcpoints 沿着给定两个坐标点的线计算 n 个点。...gcpoints(lon1, lat1, lon2, lat2, npoints) lon1 和 lat1 是地理坐标系中的起始点 lon2 和 lat2 是地理坐标系中的终止点 npoints 表示要计算的点数...[注2] drawgreatcircle(lon1, lat1, lon2, lat2, del_s=100.0, **kwargs) lon1 和 lat1 是地理坐标系中的起始点 lon2 和 lat2

2.3K3 0

Basemap系列：管理投影

投影方式及其特征在Basemap对象创建时确定。这种方式和其他库（比如 GDAL）的做法具有很大不同。...Basemap支持的epsg代码在 /mpl_toolkits/basemap/data/epsg 文件中。...尽管epsg代码在文件中，但有时还是会出现一些问题。...,urcrnrlat=44., resolution='i', projection='tmerc', lat_0 = 39.5, lon_0 = -3.25) map.drawmapboundary...0 = 30., lon_0 = -27., llcrnrx=500000.

1.6K2 0

Basemap系列教程：绘图

_1=30., lat_2=60., lat_0=34.83158, lon_0=-98.)...,urcrnrlat=43., resolution='i', projection='tmerc', lat_0 = 39.5, lon_0 = 1) map.readshapefile...坐标被转换为地图单元第二个 logo 在全球图外部，因此没有地图坐标能对其坐标进行合适的转换 1）轴名为 axicon 的轴被创建，分别确定了轴的 x，y 位置及 width 和 height..._1=30., lat_2=60., lat_0=34.83158, lon_0=-98.)..._1=30., lat_2=60., lat_0=34.83158, lon_0=-98.)

4.3K1 0

Basemap系列教程：自定义colormap

import linspace from numpy import meshgrid from matplotlib.colors import LinearSegmentedColormap # 表示在[...LinearSegmentedColormap.from_list("my_colormap", ((0, 0, 0), (1, 1, 1)), N=6, gamma=1.0) map = Basemap(projection='tmerc...', lat_0=0, lon_0=3, llcrnrlon=1.819757266426611, llcrnrlat.../sample_files/colorfile.clr") map = Basemap(projection='tmerc', lat_0=0, lon_0=3,...plt.show() 使用 read_color_table 读取颜色文件，并返回等级划分和自定义的colormap 1）检查是否有注释行 2）原始值的范围是 [0 255]，转换为

1.1K4 0

Python | 绘制影像地图（Basemap）

在Basemap底层使用了GEOS库，用来将海岸线和边界特征剪切到所需的地图投影区域。随着python2.7的寿终正寝，Basemap的支持也结束了，后续会被Cartopy替代。...在创建Basemap类实例时，必须指定所需的地图投影。有两种方法提供矩形映射投影区域的四个角的每一个的纬度和经度值。提供地图投影区域中心的lat/lon值以及地图投影坐标中的区域的宽度和高度。..., lat = -104.237, 40.125 # Location of Boulder # 通过basemap对象，将地理坐标直接变换为投影坐标。...xpt,ypt = m(lon,lat) # 设定inverse参数为True，可以把投影坐标转换为地理坐标，再转回来。...,lat,mag分别代表经度、纬度、震级 for lon,lat,mag in zip(lons,lats,magnitudes): x,y = eq_map(lon,lat) msize

6.6K3 0

python实现地理位置类数据爬取与geohash应用初探

（经纬度）功能；逆向地理编码服务提供将坐标点（经纬度）转换为对应位置信息（如所在行政区划，周边地标点分布）功能。...""" 百度坐标系(BD-09)转火星坐标系(GCJ-02) 百度——>谷歌、高德 :param bd_lat:百度坐标纬度 :param bd_lon:百度坐标经度...= z * math.sin(theta) return [gg_lng, gg_lat] def wgs84_to_gcj02(lng, lat): """ WGS84转..., bd_lat): lon, lat = bd09_to_gcj02(bd_lon, bd_lat) return gcj02_to_wgs84(lon, lat) def wgs84..._to_bd09(lon, lat): lon, lat = wgs84_to_gcj02(lon, lat) return gcj02_to_bd09(lon, lat) def

3.6K4 0

Basemap系列教程：背景方法

,urcrnrlat=46., resolution='i', projection='cass', lat_0 = 39.5, lon_0 = 0.)...,urcrnrlat=50., resolution='i', projection='tmerc', lat_0 = 40., lon_0 = -80) map.drawmapboundary...,urcrnrlat=50., resolution='i', projection='tmerc', lat_0 = 40., lon_0 = -80) map.drawmapboundary...,projection='lcc',\ lat_1=45.,lat_2=55,lat_0=50,lon_0=-107.)...译注：如果你使用的 python2.x，在改变 tmpdir 变量的基础上，还需要更改 import Image 为 from PIL import Image wmsimage 使用 WMS 协议下载并绘制图像

3.4K2 1

Basemap系列教程：绘制子图及小地图

使用 matplotlib 中的 subplots 可以在同一个 figure 中绘制多个地图。...', lat_0 = 39.5, lon_0 = 3, ax=ax3) map3.drawmapboundary(fill_color='#9999FF') map3.fillcontinents(color...使用嵌入定位器可以在大地图中添加小地图，结果比在主地图中创建子图要好。嵌入定位器是一个非常酷的类，可以放大一个图的局部，并绘制在这个图上，从而展示某一块区域。注：比如用来在地图拐角显示南海地区。...'][:] lon = data.variables['lon'][:] lons, lone = lon[28], lon[57] lats, late = lat[35], lat[14] nx =..., llcrnrlat = lat1, urcrnrlon = lon2, urcrnrlat = lat2, projection = 'cyl', ax = ax) map.drawparallels

5K1 0

【Python】应用：pyproj地理计算库应用

1. pyproj介绍 pyproj 是 Python 中一个非常有用的库，用于执行地理坐标系统的转换和地图投影。...它是基于 proj 库的 Python 接口，因此继承了 proj 库的强大功能，但以 Python 的形式提供了更友好的接口。 2. 环境配置安装：pip install pyproj 3....= pyproj.Transformer.from_crs(lat_lon_crs, utm_crs, always_xy=True) # 定义经纬度坐标 lon, lat = 120.30550,..., lat = proj(easting, northing, inverse=True) return lat, lon # 示例 UTM 坐标：Zone 18N, Easting: 500000..., distance = geod.inv(lon1, lat1, lon2, lat2) print(f"Distance: {distance} meters") print(f"Initial

931 0

结合Geotools实现百度09，国测局02和经纬度的相互转换

4、从图中可以看出，gcj02和wgs84的区别不是很大在一些不是很精确地情况下可以认为是一样的，bd09的区别稍微大一点； # 实现思路由于坐标转换是单个点的，所以在处理一个shp的坐标转换的时候..., gg_lat}; } /** * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换 * 即谷歌、高德转百度 * @param gcj_lon...) + 0.006; return new double[]{bd_lon, bd_lat}; } /** * WGS84转GCj02 * @param...} } /** * GCJ02 转换为 WGS84 * @param gcj_lon * @param gcj_lat * @returns {*..., gcj_lat * 2 - mglat}; } } /** * 百度转换为wgs84 * @param bd_lon * @param

2.3K2 0

Python编程 | T-N波作用通量水平分量

受“甲方”委托，我写了一个计算T-N波作用通量水平分量的Python脚本。...位势与位势高度位势高度的单位是gpm位势米或者m米，gpm在metpy的单位系统中与m是等价的。位势的单位是m²/s²，数值上约是位势高度的9.8倍。其实就是位势高度乘以重力加速度等于位势。...# 位势高度转位势 Φ_day = mpcalc.height_to_geopotential(hgt_day) Φ_mon = mpcalc.height_to_geopotential(hgt_mon...({'lon': lon_rad, 'lat': lat_rad}) v_climatic = v_climatic.assign_coords({'lon': lon_rad, 'lat': lat_rad...= lon_deg.values lon[lon>180] -= 360 fx = fx.assign_coords({'lon': lon, 'lat': lat_deg}).sortby(['lon

5.4K5 1

绘图技巧 | 我用python画出了北上广深的地铁路线动态图

幸好我在网上搜到了GCJ-02转BD-09的公式，并用python实现此公式： #需要的两个常量先设置好 pi = 3.1415926535897932384 #π r_pi = pi * 3000.0.../180.0 def gcj02_bd09(lon_gcj02,lat_gcj02): b = math.sqrt(lon_gcj02 * lon_gcj02 + lat_gcj02 * lat_gcj02...(lon_gcj02 * r_pi) lon_bd09 = b * math.cos(o) + 0.0065 lat_bd09 = b * math.sin(o) + 0.006...return [lon_bd09,lat_bd09] 这样我们就写好了一个python将GCJ-02坐标系转成BD-09的函数，调用这个函数，就可以将高德地图获取的坐标点集合统统转换成百度坐标系。...上图中的访问应用（AK），保存好，这在后续的可视化中将要用到。

1.6K3 0

WebGis技术汇总

var dLat = this.transformLat(lon - 105.0, lat - 35.0); var dLon = this.transformLon(lon - 105.0...d = this.delta(wgsLat, wgsLon); return { 'lat': wgsLat + d.lat, 'lon': wgsLon + d.lon };..., lon) { if (lon 137.8347) return true; if (lat < 0.8293...; return [bd_lng, bd_lat] }; /** * WGS-84 转 GCJ-02 * @param lng * @param lat *...用 fromLonLat() 方法把4326的坐标转换为3857的坐标。

7973 0

python练习题（3）——地球数据的计算

而Python作为一门功能强大且易学的编程语言，为我们提供了处理地球数据的优秀工具和技术。在本文中，我们将探索一道关于地球数据的Python程序练习题，通过实践和应用，发现Python编程的魅力。...编程练习有助于我们更深入地理解Python的各种功能和概念，同时帮助我们发现和纠正在实践中可能犯的错误。只有在实际应用中，我们才能真正体会到编程的乐趣和挑战。...本文的Python程序练习题本文的Python程序练习题将涉及地球数据的计算，帮助我们更好地理解Python在处理地球数据方面的应用。题目如下：地球数据的计算。...首先，我们需要将经纬度转换为弧度，然后使用球面三角函数来计算两点之间的距离。Python的math库提供了对数学计算的支持，包含了我们所需的三角函数。..., lon1, lat2, lon2): earth_radius = 6371.0 # 地球平均半径，单位：公里 # 将经纬度转换为弧度 lat1_radians = degrees_to_radians

5531 0

matlab、python矩阵导入ArcGIS（绘制空间图）

对于一些nc数据或者遥感影像处理时，虽然一些第三方软件可以出图，但我们往往需要借助python或者matlab软件进行数据处理，但最后保存下来数据如何导入arcgis进行分析呢？...格网的空间分辨率，以上面为例就是分辨率为0.0833333 NODATA_value代表没有数据的值，通常为-999等，还是看别人当初怎么定义的注意：NODATA_value不能是nan，如果是nan值，建议转换为...= data.variables['lat'][:].data # 读取纬度lon = data.variables['lon'][:].data # 读取经度time = data.variables...lines.insert(5, str6)with open('poyang_SRI.txt', 'w') as f: f.writelines(lines)print("over")另外一种方式：参考博客在arcgis...中ASCII转栅格（ASCII to Raster）后面根据需要进行后续操作，比如我的需要重分类一下才能看出区别。

1931 0

点击加载更多

扫码

添加站长进交流群

领取专属 10元无门槛券

手把手带您无忧上云

扫码加入开发者社群

相关资讯

热门标签

活动推荐

运营活动

活动名称

广告关闭