極立体図の地理変換?
現在、CANGRID気候データ(Surfer Grid ascii、「。grd」ファイルとして提供)をArcGISにインポートしています。グリッドのサイズは、95行125列です。メタデータは、緯度/経度の原点(左下隅)、セルサイズ(50km)に加えて、中央子午線(110度W)と原点の緯度(60度N)を備えた極立体図としてのノートの投影を提供します。 最初に.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, input_osr) …