相关疑难解决方法(0)

我有5个numpy形状nx,ny

lons.shape = (nx,ny)
lats.shape = (nx,ny)
reds.shape = (nx,ny)
greens.shape = (nx,ny)
blues.shape = (nx,ny)

红色,绿色和蓝色数组包含的值范围为0-255,lat/lon数组是纬度/经度像素坐标.

我的问题是如何将这些数据写入geotiff？

我最终想用底图绘制图像.

这是我到目前为止的代码,但是我得到一个巨大的GeoTIFF文件(~500MB),它出现空白(只是一个黑色图像).另请注意,nx,ny = 8120,5416.

from osgeo import gdal
from osgeo import osr
import numpy as np
import h5py
import os

os.environ['GDAL_DATA'] = "/Users/andyprata/Library/Enthought/Canopy_64bit/User/share/gdal"

# read in data
input_path = '/Users/andyprata/Desktop/modisRGB/'
with h5py.File(input_path+'red.h5', "r") as f:
    red = f['red'].value
    lon = f['lons'].value
    lat = f['lats'].value

with h5py.File(input_path+'green.h5', "r") as f:
    green = f['green'].value

with h5py.File(input_path+'blue.h5', "r") as f:
    blue = f['blue'].value

# convert …