极地立体影像的地理转换?
我目前正在将CANGRID气候数据(以Surfer Grid ascii,“。grd”文件提供)导入ArcGIS。网格大小为95行乘125列。元数据提供原始纬度/经度(左下角),像元大小(50公里)以及以中央子午线(纬度110度)和纬度纬度(北纬60度)的极地立体图形式进行的注释投影。 首次尝试将.grd转换为.ascii和.flt栅格失败后,我设法使用GDAL设置范围和投影,但是数据集未正确与目标区域的边界对齐。见下图。 是否存在公认的极地立体地理变换,可以解释这种缺乏对齐的情况? 例如,是否应该使用特定的转换因子或轮换? 数据集中的示例文件在此处:“ t201113.grd” 这是我目前在GDAL中使用的代码 ds = gdal.Open("t201113.grd") array = ds.ReadAsArray() x_rotation = 0 y_rotation = 0 xres = 1 yres = -1 llx = -129.8530 lly = 40.0451 ulx = -175.144 uly = 71.385 input_osr = osr.SpatialReference() input_osr.ImportFromWkt(ds.GetProjection()) wgs84_osr = osr.SpatialReference() wgs84_osr.ImportFromEPSG(4326) wgs_to_nps_trans = osr.CoordinateTransformation(wgs84_osr, …