如何更改新栅格的投影

时间:2021-06-05 16:33:16

标签: python pandas gdal coordinate-systems rasterio

我尝试在 Python 中制作一个算法,其中我输入了地理参考栅格(EPSG 32722 投影系统),将其所有 NaN 值转换为零,然后使用 EPSG 4326 投影系统保存了新图像。< /p>

import skimage.io
import pandas as pd
import numpy as np
from sklearn.impute import SimpleImputer

pathhr = 'C:\\Users\\dataset\\S30W051.tif'
HR = skimage.io.imread(pathhr)
df1 = pd.DataFrame(HR)
imputer = SimpleImputer(fill_value=np.nan, strategy='mean')
X = imputer.fit_transform(df1)
X = pd.DataFrame(X, columns=df1.columns)
X.isna().sum()

#save function
savedata = df1.to_numpy()
skimage.io.imsave('C:\\Users\\dataset\\S30W051_TEST.tif', savedata)

但是当我在这个脚本的末尾保存我的光栅时,我得到一个非地理参考的 TIFF 光栅。

如何使用 EPSG 4326 投影系统保存这个新光栅?

我寻求帮助以在 Python 中解决这个问题。非常感谢。

0 个答案:

没有答案