python GDAL 影像处理系列之给遥感影像加地理信息（二）

  在遥感影像中，地理坐标的信息异常重要，但是对于numpy array数据来说，没有地理位置信息的，因此如果需要地理信息的情况下，如果再将遥感影像作为普通多波段RGB影像来进行处理就有点不合适了，那怎么合适呢？我就来了。。。。（欢呼）

  使用gdal的第一步，就需要安装gdal，安装教程在：gdal安装教程

为一张无地理信息的遥感影像添加地理信息，一个函数就可以呀！

from osgeo import osr, gdal


def assign_spatial_reference_byfile(src_path, dst_path):
    src_ds = gdal.Open(src_path, gdal.GA_ReadOnly)
    sr = osr.SpatialReference()
    sr.ImportFromWkt(src_ds.GetProjectionRef())
    geoTransform = src_ds.GetGeoTransform()
    dst_ds = gdal.Open(dst_path, gdal.GA_Update)
    dst_ds.SetProjection(sr.ExportToWkt())
    dst_ds.SetGeoTransform(geoTransform)
    dst_ds = None
    src_ds = None
    
src_path = 'F:/yi_2009.tif' 
dst_path = 'C:/Users/Dell/Desktop/gg_building.tif'  
assign_spatial_reference_byfile(src_path, dst_path)

     调用方法非常简单，src_path为带地理坐标的影像，dst_path为无地理信息的遥感影像，当然转完的影像dst_path与src_path的地理坐标是一致的。

---

很开心大家能提到非常宝贵的意见，收到大家的评论，欢迎讨论留言！！

有个朋友说如果没有参考的地理坐标，怎么加地理坐标？如下代码段所示：

def assign_spatial_reference(filepath):
    
    sr = osr.SpatialReference()
    sr.ImportFromEPSG(4326)                    #地理坐标投影
    
    ds = gdal.Open(filepath, gdal.GA_Update)
    ds.SetProjection(sr.ExportToWkt())
    ds.SetGeoTransform([0,1,0,0,0,1])          #地理6参数
    
    ds = None

assign_spatial_reference(image_path)

可以自己构造一个地理坐标投影，然后按照你想要的地理位置加上相应的地理六参数即可，如果你对地理位置（经纬度坐标）都没有要求的话，只是为了增加地理坐标而加地理坐标，可以像我的代码一样，赋一个假的地理坐标即可。

希望能给这位朋友提供帮助。

整理不易，欢迎点赞、评论、收藏。