python - 散点数据插值（纬度、经度和温度）

Question

在此处输入图像描述我设法用下面的代码进行插值，我的问题是投影，当我在GIS software（Erdas 想象）上打开最终输出时，它在加利福尼亚的位置不正确我的输出不合适，有没有办法不affline从我的坐标计算，我不想手动放置它们

lon = data2['Longitude'].values

lon2=data2['Longitude2'].values

lat = data2['Lattitude2'].values

t = data2['Temparature ℃'].values

xp, yp,_= to_proj.transform_points(ccrs.Geodetic(), lon, lat).T

x_m,y_m,t_m=remove_nan_observations(xp, yp, t)


x,y,z=interpolate_to_grid(x_m, y_m, t_m, interp_type='linear', hres=1000, 
                minimum_neighbors=4, gamma=0.25, kappa_star=5.052, search_radius=1000, rbf_func='linear', rbf_smooth=1)


new_transform1 = Affine(1883.50, 0.0, 406063.73814945, 0.0, 2300.95, 4674584.0)

image = rasterio.open(' minus_Central_meridian_ - lat.tiff', 'w', driver='GTiff',
                        height = z.shape[0], width = z.shape[1],
                        count=1, dtype=str(z.dtype),
                  crs='+proj=utm +zone=10 +south+ellps=WGS84+EPSG:857+ +datum=WGS84 +units=m +no_defs',transform=new_transform1
                  )


image.write(z, 1)
image.close()

先感谢您！！！！！

看我的输出它应该在加利福尼亚很好地对齐，但它不合适

标签： pythoninterpolationprojectionrasterio