python将tif从wgs84转gcj02

# mark: 主要是dem从wgs84转成gcj02，步骤为wgs84的4326转成3857，之后投影进行偏移，再将偏移后的3857转成4326from osgeo import gdal, osr# 经纬度转投影
def tif4326To3857(input_file, output_file):options = gdal.WarpOptions(format='GTiff', srcSRS='EPSG:4326', dstSRS='EPSG:3857', width=1452)gdal.Warp(output_file, input_file, options=options)print('转换完成')# 投影转经纬度
def tif3857To4326(input_file, output_file):options = gdal.WarpOptions(format='GTiff', srcSRS='EPSG:3857', dstSRS='EPSG:4326', width=1452)gdal.Warp(output_file, input_file, options=options)print('转换完成')# wgs84投影偏移至gcj02
def tifWgs84ToGcj02(input_file, output_file):dataset = gdal.Open(input_file)im_width = dataset.RasterXSize #栅格矩阵的列数im_height = dataset.RasterYSize  # 栅格矩阵的列数im_array = dataset.GetRasterBand(1).ReadAsArray()im_projection = dataset.GetProjection()im_GeoTransform = dataset.GetGeoTransform()im_GeoTransform_new = (im_GeoTransform[0]+740.4591671, im_GeoTransform[1], im_GeoTransform[2], im_GeoTransform[3]+24.60901496, im_GeoTransform[4], im_GeoTransform[5])# 创建新的tif文件driver = gdal.GetDriverByName('GTiff')out_tif = driver.Create(output_file, im_width, im_height, 1, gdal.GDT_Float32)out_tif.SetGeoTransform(im_GeoTransform_new)out_tif.SetProjection(im_projection)  # 给新建图层赋予投影信息# # 获取地理坐标系统信息，用于选取需要的地理坐标系统# srs = osr.SpatialReference()# # 导入tif文件的坐标系# srs.ImportFromEPSG(3857)  # 定义输出的坐标系为"WGS 84"，AUTHORITY["EPSG","4326"]# out_tif.SetProjection(srs.ExportToWkt())  # 给新建图层赋予投影信息# 数据写出out_tif.GetRasterBand(1).WriteArray(im_array)  # 将数据写入内存，此时没有写入硬盘out_tif.FlushCache()  # 将数据写入硬盘out_tif = None  # 注意必须关闭tif文件print('偏移完成')if __name__ == '__main__':file1 = 'C:\\Users\\ZJZN\\Desktop\\tiff\\4326.tif'file2 = 'C:\\Users\\ZJZN\\Desktop\\tiff\\3857.tiff'file3 = 'C:\\Users\\ZJZN\\Desktop\\tiff\\3857-gcj02.tiff'file4 = 'C:\\Users\\ZJZN\\Desktop\\tiff\\4326-gcj02.tiff'tif4326To3857(file1, file2)tifWgs84ToGcj02(file2, file3)tif3857To4326(file3, file4)

tif4326To3857(file1, file2) 是将wgs84的经纬度转成投影
tifWgs84ToGcj02(file2, file3) 是将投影的进行偏移
tif3857To4326(file3, file4) 再将偏移后的投影转成经纬度

注意：im_GeoTransform_new 中的偏移量是各地根据实际计算得到的。

	gcj02-3857投影	wgs84-3857投影	差值
x	12587949.13	12587208.67	740.4591671
y	4337622.252	4337597.643	24.60901496

python将tif从wgs84转gcj02相关推荐

Unity经纬度相互转换（WGS-84、GCJ-02、BD-09）
/*** * * Title: MXFramework * 主题: GPS工具类 * Description: * 功能:1.各种坐标 ...
WGS84与GCJ02、BD09经纬度坐标转换介绍
在涉及到经纬度的地方,我们在编码的时候,经常需要进行转换.而且转换代码似乎非常复杂,真正理解,其实需要做一些工作,这里我将了解的经纬度坐标相关的知识做一个梳理. 首先了解经纬度坐标系统. 因为地球并不 ...
GPS坐标系转换（标准坐标系WGS84转GCJ-02火星坐标系）
GPS坐标系转换(标准坐标系WGS84转GCJ-02火星坐标系) 坐标系简介 WGS-84 标准坐标系 GCJ-02 - 国测局坐标 BD-09 - 百度坐标系干货前的说明 WGS84转GCJ-02 ...
python读取tif文件与png文件
python读取tif文件与png文件读取tif文件: import tifffile as tiff image = tiff.imread(image) 读取png文件: import nump ...
WGS84、GCJ-02、BD-09、图吧坐标简介及坐标转换实现（js版）
1.不同地图服务商的坐标系出于国家安全考虑不同的国家在规定地图服务提供商在提供地图服务的时候需要将地图数据进行加密.现在比较流行的坐标系有WGS84.GCJ-02.BD-09,且各个坐标之间存在偏差 ...
WGS84、GCJ-02、BD-09、图吧坐标简介及坐标转换实现（java版）
1.定位系统及如何定位现在全球有四个卫星定位系统:美国的全球定位系统GPS(Global Positioning System),俄罗斯的格洛纳斯GIONASS, 欧盟的伽利略系统,我国的北斗. 每 ...
Android开发：WGS-84、GCJ02坐标名词解释及坐标转换
WGS-84.GCJ02坐标名词解释及坐标转换目录一.什么是WSG-84坐标系? 二 .什么是GCJ-02坐标系? 三.GCJ-02坐标系使用场景四.函数说明五.工具类代码注:本文仅包含国测 ...
python读取tif图片时保留其16bit的编码格式
python读取tif图片时保留其16bit的编码格式 tif图片的编码格式一般是16bit的,在使用python-opencv读取tif文件时,为了保留其编码格式,我们需要用以下的方式: impor ...
转换WGS84、GCJ02、BD09坐标，解决百度高德坐标偏移
转换WGS84.GCJ02.BD09坐标,解决百度高德坐标偏移做过地图相关开发的同学肯定会遇到这样一个问题:同样的经纬度坐标,在百度地图和高德地图上位置不一样. 关于坐标系我们通常用经纬度来表示一 ...

python将tif从wgs84转gcj02

python将tif从wgs84转gcj02相关推荐

最新文章

热门文章