QGIS：纬度坐标反转；地图颠倒

Question

我一直在使用 QGIS 来显示荷兰长期降水平均值的地图。然而，当QGIS打开数据时，地图显示倒置

我注意到坐标显示为 0 - 266 (lon) 和 -315 - 0 (lat)。我认为纬度是颠倒的

而不是 -315 - 0 它应该是 0 - 315 并且地图应该看起来不错。但我不知道如何反转这个值。该文件是一个 NetCdf 文件。我用 EmEditor 打开 QGIS 为我制作的 XML 元数据，但它确实显示了正确的坐标（在 lat/lon 中），所以我认为这与 QGIS 设置地图的方式或它将 lat/lon 转换为米的方式。有谁遇到过和我一样的问题吗？提前致谢！

Answer 1

我很确定您可以使用 GDAL 配置选项 GDAL_NETCDF_BOTTOMUP=[YES/NO] 将 NetCDF 转换为 geotiff，并使生成的栅格正确定向为北上。尝试使用 gdal_translate with the above option. See here 了解更多详情。

Answer 2

感谢 Micha（见评论）：

有人告诉我使用 GDAL（地理空间数据抽象库）解决问题，这是一种查看 translate/process 元数据的方法。这很难理解，而我在编程和使用像 GDAL 这样强大的 'languages' 方面相对较新。为了输入 GDAL 代码，我使用了 QGIS 附带的 OSGeo4W Shell。我用来翻转地图的命令是：

gdal_translate -of netCDF -co WRITE_BOTTOMUP=NO 我的 netcdf.nc output.nc

（另见 this 简短 GDAL/netCDF 介绍）。

Answer 3

在R中你可以使用rotate函数



library(raster)
library(gdalUtils)


workdir <- "Your workind dir"
setwd(workdir)

ncfname <- "adaptor.mars.internal-1563580591.3629887-31353-13-1b665d79-17ad-44a4-90ec-12c7e371994d.nc"

# get the variables you want
dname <- c("v10","u10") 


# open using raster 
datasetName <-dname[1]
r <- raster(ncfname,  varname = datasetName)
r2 <- rotate(r)

writeRaster(r2,"wind.tif",driver = "TIFF")

QGIS：纬度坐标反转；地图颠倒

QGIS: latitude coordinates are reversed; map upside down

latitude-longitude

gdal

netcdf

qgis