Mapbox studio上传的tiff在反梅里迪亚和不';不要延伸到两极

Mapbox studio上传的tiff在反梅里迪亚和不';不要延伸到两极,mapbox,geotiff,mapbox-studio,Mapbox,Geotiff,Mapbox Studio,我有一个GeoTIFF文件,我正试图在Mapbox Studio中显示它,但它在反日线(日期线)处显示一条空白的垂直线,光栅没有像它应该的那样延伸到极点。GeoTIFF文件是在Python中从RGB numpy阵列标注栏创建的,该标注栏的尺寸范围为-180到180度,纬度为90到-90度。GeoTIFF在QGIS中显示良好 以下是创建文件的代码段: driver = gdal.GetDriverByName("GTiff") ds_out = driver.Create('gfs.t12z.pg

我有一个GeoTIFF文件,我正试图在Mapbox Studio中显示它,但它在反日线(日期线)处显示一条空白的垂直线,光栅没有像它应该的那样延伸到极点。GeoTIFF文件是在Python中从RGB numpy阵列标注栏创建的,该标注栏的尺寸范围为-180到180度,纬度为90到-90度。GeoTIFF在QGIS中显示良好

以下是创建文件的代码段:

driver = gdal.GetDriverByName("GTiff")
ds_out = driver.Create('gfs.t12z.pgrb2.0p25.spd_rgb_10.f012_EPSG_GRIB.tif', 1440, 721, 3, 
gdal.GDT_Byte)
ds_out.SetGeoTransform([-180.125, 0.25, 0, 90.125, 0, -0.25])
ds_out.SetProjection('GEOGCS["Coordinate System imported from GRIB file",DATUM["unknown",SPHEROID["Sphere",6371229,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]')
ds_out.GetRasterBand(1).WriteArray(red)
ds_out.GetRasterBand(2).WriteArray(green)
ds_out.GetRasterBand(3).WriteArray(blue)
ds_out.FlushCache()
ds_out = None
以下是Mac OSX预览中显示的GeoTIFF图像:

以下是它在Mapbox Studio中的显示方式的两个视图:


在上传到Mapbox Studio之前,我还尝试使用gdalwarp将投影更改为EPSG:3857,但它仍然显示相同的内容。

我发现了问题所在。我将数组插值到EPSG:3857坐标系中,并将输出文件的投影设置为EPSG:3857:

driver = gdal.GetDriverByName("GTiff")
ds_out = driver.Create('gfs.t12z.pgrb2.0p25.spd_rgb_10.f012_EPSG_3857.tif', 
x_size, x_size, 3, gdal.GDT_Byte)
cell_dim = 20037508.34 / 720.0
ds_out.SetGeoTransform([-20037508.34, cell_dim, 0, 20037508.34, 0, -cell_dim])
spatialReference = osr.SpatialReference()
spatialReference.ImportFromEPSG(3857)
ds_out.SetProjection(spatialReference.ExportToWkt())
ds_out.GetRasterBand(1).WriteArray(red)
ds_out.GetRasterBand(2).WriteArray(green)
ds_out.GetRasterBand(3).WriteArray(blue)
ds_out.FlushCache()
ds_out = None
ds_in = None