このnetCDFラスターはどの座標系ですか?


10

netCDFラスターファイルを取得しましたが、ラスターが構築されている座標系の名前を取得するためのメタデータを取得できませんでした。ラスター自体には座標系が埋め込まれていません。WGS84だけだと思っていましたが、最初はそのように見えましたが、ArcMapをさらに調査したところ、それほど一般的ではないシステムであることがわかりました。表示方法は次のとおりです。 ここに画像の説明を入力してください

オレンジ色のラスターは、比較のためにここに挿入したWGS84の通常のラスターです。紫色のものは、座標系が不明なラスターです。これが何かの手掛かりはありますか?

一部の更新:ここでのnetCDFラスタがある:https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 私も画像提供者からいくつかのメタデータを取得することができました:

netcdf sic_average_nclimate {
dimensions:
    nlon = 361 ;
    nlat = 90 ;
    nseas = 4 ;
variables:
    float SIC_Change(nlat, nlon) ;
        SIC_Change:Title = "Gridded Multi-Model Ensemble Mean Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change(nseas, nlat, nlon) ;
        SIC_Season_Change:Title = "Gridded Multi-Model Ensemble Mean Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Change_STD(nlat, nlon) ;
        SIC_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change_STD(nseas, nlat, nlon) ;
        SIC_Season_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float LAT(nlat) ;
        LAT:Title = "Latitude" ;
    float LON(nlon) ;
        LON:Title = "Longitude" ;

// global attributes:
        :Title = "Ice Concentration metrics for Model subset as in Figure 1 of NCLIMATE paper" ;

これらは緯度と経度の境界を示していますが、座標系に関する情報は明らかにありません。


1
南極大陸は通常極地ステレオ投影であり、すべて南極から北に向かっています。nsidc.org/data/polar_stereo/ps_grids.html
Mapperz

1
奇妙な、0〜360のように見えますが、なぜ北半球でプロットしているのですか?
mkennedy 2014

1
NetCDFヘッダーで報告されたエクステントがlat / lonのように見えるが、lat:0から?、lon:0から360である場合、それらをlat:-90から(?-90)およびlon:-180から+にハックできますか。 180?
mkennedy 2014

@Mapperzとmkennedy、ポーラーステレオグラフィックプロジェクションを割り当てようとしましたが、機能しませんでした。垂直に下向きにミラーリングすることはできますが、右側の半分は180Eの経度を超えています。netCDFをメモ帳で開きましたが、lat、longに関する情報は見つかりませんでした。あなたがあなたの側でそれを開くことに興味があるなら、私は問題のラスターを添付しています。
multigoodverse

そのリンクは機能しません-Dropboxには「所有者はアクセスを承認していません」と表示されます。
BradHards 2014

回答:


6

次のGDALコマンドを使用して、データを正しい位置にシフトできます。

gdal_translate -a_srs EPSG:4326 -a_ullr 0 0 360 -90 NETCDF:"sic_average_nclimate.nc":SIC_Change change360.tif
gdalwarp -t_srs WGS84 change360.tif change180.tif  -wo SOURCE_EXTRA=1000 --config CENTER_LONG 0

(Frank Warmerdamによる少しの助けを借りて:180子午線をカットして0 360から-180 180にラスターを再投影する方法

結果は、立体投影では次のようになります。

ここに画像の説明を入力してください


1

(これは答えではありませんが、コメントにするには長すぎるプロセスの提案です)

両方のラスタで共通のポイントを識別できるように見えるので、GeoReferencerツールを使用してQGISで地理参照できます。これにより、一連のGCPポイントと変換パラメータを含むGDALスクリプトが生成されます。2つの予測の違いを説明している可能性があるため、このスクリプトにこれらのポイントのサンプルを追加して質問に追加できます。

最悪の場合のシナリオでは、再配置されたラスターに基づく値の概算があり、目標に役立つ可能性があります。


+ 1とはい、ジオリファレンスは最後の弾丸かもしれませんが、コントロールポイントを一致させることが難しいため、プロセスによっていくつかのエラーが発生します。さらに、netCDFラスターは同じ子午線を持たないため、分割する必要があります。たとえば、南極の尾は右側にあり、オレンジ色のラスターは左側に示しています。
multigoodverse 2014

0

左下隅は(0,0)座標にありますが、私にとって(0,0)はあなたのケースの右上隅の座標であり、原点が-90である地理座標系はわかりません。

ArcGISでの作業中に、ジオリファレンスツールバーを使用してデータを適切な場所に移動できます。右上を(0,0)に移動し、左下を(-360、-90)に移動します。これにより、正確な変換が得られますが、日付線に問題があります。

ここに画像の説明を入力してください

したがって、カスタム地理座標系を作成するのに役立つこのリンクを試すことができます。要約すると、参照システムで、-180に本初子午線を持つカスタムWGS 84を作成し、次に経度回転を使用してカスタム変換を行います。

ここに画像の説明を入力してください

この新しいシステムでのジオリファレンスは、原点を(-180、-90)に設定することで構成されます。翻訳が必要なだけなので、1つの制御点で十分です。左下のピクセルを注意深く選択し、右クリックして正確なXY値を追加します)。その後、このようになります

ここに画像の説明を入力してください

最後の注意として、座標系がWGS 84に基づいていることは確実ではありません。ただし、地理座標系が度であり、ラスターの解像度は、誤ったデータムの選択による潜在的なエラーよりもおそらく粗くなります。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.