Python 散乱数据插值(纬度、经度和温度)
我设法用下面的代码来插值,我的问题是投影,当我在Python 散乱数据插值(纬度、经度和温度),python,interpolation,projection,rasterio,Python,Interpolation,Projection,Rasterio,我设法用下面的代码来插值,我的问题是投影,当我在GIS软件(Erdas imagine)上打开最终输出时,它在加利福尼亚的位置不正确,我的输出不合适,有没有办法不从我的坐标计算affline,我不想手动放置它们 lon = data2['Longitude'].values lon2=data2['Longitude2'].values lat = data2['Lattitude2'].values t = data2['Temparature ℃'].values xp, yp,_=
GIS软件
(Erdas imagine)上打开最终输出时,它在加利福尼亚的位置不正确,我的输出不合适,有没有办法不从我的坐标计算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()
提前谢谢你
看我的输出,它本应该在加利福尼亚很好地对齐,但它不合适