Python GDAL学习笔记(二)

最新推荐文章于 2024-04-17 23:11:06 发布

原创

最新推荐文章于 2024-04-17 23:11:06 发布 · 7.7k 阅读

本文详细介绍了GIS中自定义空间坐标系的方法，包括使用EPSG、WKT及PROJ字符串初始化坐标系统，遥感影像的几何校正过程，空间参考系统转换至Albers投影坐标，以及GDAL Warp函数的参数详解。

一、自定义空间坐标系

使用EPSG或EPSGA编号进行初始化

EPSG：欧洲石油勘探组织（European Petroleum Survey Group，EPSG)，EPSG编译并传播了EPSG大地参数集、广泛使用的地球椭球体数据库、大地基准、地理和投影坐标系、测量单位等。该数据集是坐标参考系统和坐标转换定义的集合，这些定义可以应用在全球、区域、国家或局部地区。目前，这个数据集被广泛接受并使用。EPSG通过一个Web发布平台http://www.epsg.org/对该数据集进行分发。在该数据库中存储的椭球体，投影坐标系等不同组合都对应着不同的ID号，这个号在EPSG中被称为EPSGcode，它代表特定的椭球体、单位、地理坐标系或投影坐标系等信息。

from osgeo import gdal, osr

srs = osr.SpatialReference()
srs.ImportFromEPSG(4509) #WGS 84，其编码为4326
print(srs)

PROJCS["CGCS2000 / Gauss-Kruger CM 117E",
    GEOGCS["China Geodetic Coordinate System 2000",
        DATUM["China_2000",
            SPHEROID["CGCS2000",6378137,298.257222101,
                AUTHORITY["EPSG","1024"]],
            AUTHORITY["EPSG","1043"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4490"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",117],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Northing",NORTH],
    AXIS["Easting",EAST],
    AUTHORITY["EPSG","4509"]]

使用WKT字符串进行初始化

首先获取一幅已经进行投影的proj的Wkt信息。

from osgeo import gdal, osr

filename = r'F:\GDAL learning\Landsat8.tif'
dataset = gdal.Open(filename, gdal.GA_ReadOnly)
proj = dataset.GetProjection()
print(proj)

PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["A