此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,我尝试为它指定一个Polar Stereographic Projection,但是没有用。我可能可以垂直向下镜像它,但是它的右侧一半仍将超过180E经度。我用记事本打开了netCDF,但是找不到与经纬度有关的任何信息。如果您想打开侧面的栅格,请附加这个栅格。
multigoodverse

该链接无效-投递箱显示“所有者未授权访问”。
BradHards 2014年

Answers:


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脚本。您可以将这个脚本以及这些要点的样本添加到您的问题中,因为它可以解释这两个预测之间的区别。

最坏的情况是,您将根据重新定位的栅格对值进行粗略的估计,这可能会帮助您实现目标。


+ 1,是的,地理配准可能是最后一个项目符号,但是该过程会产生一些错误,因为很难匹配控制点。此外,由于netCDF栅格没有相同的本初子午线,因此必须对其进行分割。例如,南极洲的尾巴在右侧,而橙色栅格则在左侧。
multigoodverse

0

左下角位于(0,0)坐标处,但对我来说(0,0)是您情况下右上角的坐标,我不知道任何原点为-90的地理坐标系。

在ArcGIS中工作时,可以使用地理配准工具栏将数据移动到正确的位置。将右上角移至(0,0),将左下角移至(-360,-90),这将产生精确的变换,但日期线仍然存在问题。

在此处输入图片说明

因此,您可以尝试使用此链接该链接应有助于您创建自定义地理坐标系。综上所述,在参考系统中,创建一个自定义子午线为-180的自定义WGS 84,然后使用经度旋转进行自定义变换。

在此处输入图片说明

在这个新系统中,地理配准包括将原点设置为(-180,-90)。一个控制点就足够了,因为您只需要翻译。仔细选择左下像素,右键单击并添加确切的XY值。那么它应该看起来像这样

在此处输入图片说明

最后,不能确定坐标系是否基于WGS84。但是可以确定,以度为单位的地理坐标系和栅格分辨率可能比由于错误的基准选择引起的潜在误差还要大。

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.