конвертировать TIFF в NetCDF - PullRequest
Новая
       21

конвертировать TIFF в NetCDF

0 голосов
/ 02 сентября 2018

Я пытаюсь преобразовать TIFF в NetCDF файл. Ошибка говорит об ошибке индекса: 

import numpy as np
from netCDF4 import Dataset
import rasterio


with rasterio.drivers():
    src=rasterio.open(r"ia.tiff","r")
    dst_transform=src.transform
    dst_width=src.width
    dst_height=src.height
    print (dst_transform)
    xmin = dst_transform[0]
    xmax = dst_transform[0] + dst_transform[1]*dst_width
    print (xmax)
    min = dst_transform[3] + dst_transform[5]*dst_height
    print(ymin)
    ymax = dst_transform[3]

    dst_width=dst_width+1
    dst_height=dst_height+1
    outf=Dataset(r'ia.nc','w',format='NETCDF4_CLASSIC')
    lats=np.linspace(ymin,ymax,dst_width)
    lons=np.linspace(xmin,xmax,dst_height)


    lat=outf.createDimension('lon',len(lats))
    lon=outf.createDimension('lat',len(lons))
    longitude=outf.createVariable('longitude',np.float64,('lon',))
    latitude=outf.createVariable('latitude',np.float64,('lat',))
    SHIA=outf.createVariable('SHIA',np.int8,('lon','lat'))

    outf.variables['longitude'][:]=lons
    outf.varibales['longitude'][:]=lat
    im=src.read()
    SHIA[:,:]=im

    outf.description="IA for"
    longitude.units="degrees east"
    latitude.units='degrees north'

    print ("created empty array")

    outf.close()




outf.close()

ошибка в том, что ошибка индекса: размер массива данных не соответствует срезу. Может кто-нибудь взглянуть и помочь мне, где я сделал не так. Очень признателен!

...