在遥感影像中，地理坐标的信息异常重要，但是对于numpy array数据来说，没有地理位置信息的，因此如果需要地理信息的情况下，如果再将遥感影像作为普通多波段RGB影像来进行处理就有点不合适了，那怎么合适呢？我就来了。。。。（欢呼）
使用gdal的第一步，就需要安装gdal，安装教程在：gdal安装教程
为一张无地理信息的遥感影像添加地理信息，一个函数就可以呀！
from osgeo import osr, gdaldef 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 = Nonesrc_ds = Nonesrc_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 = Noneassign_spatial_reference(image_path)可以自己构造一个地理坐标投影，然后按照你想要的地理位置加上相应的地理六参数即可，如果你对地理位置（经纬度坐标）都没有要求的话，只是为了增加地理坐标而加地理坐标，可以像我的代码一样，赋一个假的地理坐标即可。
希望能给这位朋友提供帮助。