Bu netCDF raster hangi koordinat sistemi?


10

Bir netCDF tarama dosyası elde ettim, ancak tarama sisteminin oluşturduğu koordinat sisteminin adını almak için herhangi bir meta veri alamadım. Raster cihazda gömülü koordinat sistemi yoktur. Sadece bir WGS84 olacağını düşündüm ve ilk bakışta böyle görünüyordu, ancak ArcMap'ta daha fazla araştırma yaparak, bunun oldukça yaygın olmayan bir sistem olduğunu gördüm. Nasıl görüntülendiği aşağıda açıklanmıştır: resim açıklamasını buraya girin

Turuncu raster, karşılaştırma amacıyla buraya eklediğim WGS84'te normal bir raster. Mor olan, bilinmeyen koordinat sistemine sahip rasterdir. Bunun ne olabileceğine dair bir fikrin var mı?

Bazı güncellemeler: İşte netCDF raster: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 Görüntü sağlayıcıdan bazı meta veriler de alabildim:

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" ;

Sınırı en geç uzun süre gösterirler, ancak görünüşe göre koordinat sistemi ile ilgili herhangi bir bilgi yoktur.


1
Antarktika genellikle Polar Stereografik Projeksiyondadır, bu nedenle her şey güney kutbundan kuzeye gider. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz

1
Tuhaf, 0-360'a benziyor, ama neden kuzey yarımkürede çiziliyor?
mkennedy

1
NetCDF üstbilgisinde bildirilen uzantılar lat / lon gibi görünüyor, ancak lat: 0 ila?, Lon: 0 ila 360 ise. Bunları lat: -90 ila (? - 90) ve lon: -180 ila + 180?
mkennedy

@Mapperz & mkennedy, bir Polar Stereografik Projeksiyon atamaya çalıştım, ama işe yaramadı. Muhtemelen dikey olarak aşağıya doğru yansıtabilirim, ancak sağ yarısı hala 180E boylamlarının ötesinde olacak. NetCDF'yi not defteri ile açtım, ancak lat, long ile ilgili herhangi bir bilgi bulamadım. Yanınızda açmayı merak ediyorsanız, soruya raster ekliyorum.
multigoodverse

Bu bağlantı çalışmıyor - dropbox "sahibin erişim yetkisi yok" diyor.
BradHards

Yanıtlar:


6

Aşağıdaki GDAL komutlarını kullanarak verileri doğru konuma kaydırabilirsiniz:

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'ın küçük bir yardımı ile: Raster 0 360 ile -180 180 arasında 180 meridyen kesilerek nasıl yeniden üretilir )

stereografik bir projeksiyonda sonuç şöyle görünür:

resim açıklamasını buraya girin


1

(Bu bir cevap değil, yorum yapmak için çok uzun bir süreç önerisi)

Her iki rasterdeki ortak noktaları tanımlayabildiğiniz için GeoReferencer aracını kullanarak QGIS'de coğrafi referans verebilirsiniz. Bu, bir dizi GCP noktası ve dönüşüm parametreleri içeren bir GDAL komut dosyası oluşturur. İki senaryo arasındaki farkın ne olduğunu açıklayabileceğinden, bu komut dosyasını bu noktaların bir örneğiyle birlikte ekleyebilirsiniz.

En kötü senaryo, hedefinize yardımcı olabilecek yeniden konumlandırılmış raster tabanlı değerleri kabaca tahmin edersiniz.


+ 1 ve evet, coğrafi referans son mermi olabilir, ancak kontrol noktalarıyla eşleşmesi zor olduğundan işlem bazı hatalar verebilir. Ek olarak, netCDF rasterinin bölünmesi gerekecektir çünkü aynı ana meridyene sahip değildir. Örneğin, Antarktika kuyruğu sağda, turuncu raster solda gösterir.
multigoodverse

0

Sol alt köşe (0,0) koordinatlarındadır, ancak benim için (0,0) sizin durumunuzdaki sağ üst köşenin koordinatıdır ve kökeni -90 olan herhangi bir coğrafi koordinat sistemi bilmiyorum.

ArcGIS'te çalışırken , verilerinizi doğru konuma taşımak için coğrafi referans araç çubuğunu kullanabilirsiniz . Sağ üst (0,0) ve sol alt (-360, -90) 'a getirin, bu kesin bir dönüşüm sağlamalıdır, ancak yine de datelin ile ilgili bir sorununuz olacaktır.

resim açıklamasını buraya girin

Böylece , özel bir coğrafi koordinat sistemi oluşturmanıza yardımcı olacak bu bağlantıyı deneyebilirsiniz . Özetlemek gerekirse, referans sisteminizde, birincil meridyen -180 ile özel bir WGS 84, ardından boylam dönüşü kullanarak özel bir dönüşüm oluşturun.

resim açıklamasını buraya girin

Bu yeni sistemdeki coğrafi referans, başlangıç ​​noktasının (-180, -90) olarak ayarlanmasını içerecektir. Bir kontrol noktası yeterlidir çünkü sadece çeviriye ihtiyacınız vardır. Sol alt pikseli dikkatlice seçin, sağ tıklayın ve tam XY değerlerini ekleyin). O zaman böyle görünmeli

resim açıklamasını buraya girin

Son bir açıklama olarak, koordinat sisteminin WGS 84'e dayandığı kesin değildir. Ancak kesin olarak bir coğrafi koordinat sisteminin derece cinsinden ve rasterinizin çözünürlüğünün yanlış veri seçiminden kaynaklanabilecek potansiyel hatadan muhtemelen daha kaba olduğu kesin değildir.

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.