SMOS数据处理，投影变换，‘EPSG:6933‘转为‘EPSG:4326‘

最新推荐文章于 2024-11-28 20:26:32 发布

原创

最新推荐文章于 2024-11-28 20:26:32 发布 · 696 阅读

4 ·

CC 4.0 BY-SA版权

文章标签：

#python

文章讲述了在处理SMOS数据时，如何从EPSG:6933投影转换为EPSG:4326，通过Python代码实现读取nc文件、数据处理、地理变换和GDALwarp操作，确保全球数据能重叠显示。

在处理SMOS数据时，遇到了读取nc数据并存为tif后，影像投影无法改变，因此全球数据无法重叠。源数据的投影为'EPSG:6933'，希望转为'EPSG:4326'。

解决代码。

```python
import os
import netCDF4 as nc
import numpy as np
from osgeo import gdal, gdalconst
# *************************导入netCDF4
file_path = r'I:\\'
input_name = file_path + 'SM_SCIE_MIR_L4AGBB_20110101T000000_20211231T235959_100_001_8.nc'
out_path = r'I:\out\\'
res = 0.25 #目标数据的分辨率
top = 57 #目标数据的top坐标
geo_transform = [-179, res, 0.00, top, 0.00, -res] #投影转换

os.makedirs(out_path, exist_ok=True)

with nc.Dataset(input_name) as DS:
    data = DS.variables['AGB'] #nc文件数据
    data_arr1 = np.asarray(data) #转为数组
    src_proj = DS.proj4text #nc原始文件投影信息

for j in range(0, 11):
    year = str(2011 + j)
    data_arr3 = np.where(data_arr1[j] < 0, 0, data_arr1[j])
    data_arr4=np.flipud(data_arr3) #影像上下变换

    # 获取数组的行数和列数