Questions tagged «netcdf»

7
各种栅格数据格式的速度
我在查找不同栅格文件格式的任何讨论或比较基准测试时遇到麻烦(例如,用于R中的数据分析)。是否有人对为什么特定格式可能更快或更慢有任何见解?还是应该使差异最小? 具体来说,我感兴趣的是将栅格(例如GEOTIFF文件)转换为其他格式(例如netCDF)是否值得为了加快读取/写入和其他操作的速度。

2
将GTiff栅格的时间序列堆栈转换为单个NetCDF
从gdal-dev邮件列表移出: 在2013年9月2日星期一晚上7:09,David Shean写道: 嗨,大家好,我正在尝试将具有相同投影/范围/分辨率的GTiff栅格时间序列打包为单个NetCDF文件进行分发。我花了一个小时来咨询在线文档,并且使用gdal_translate,gdalbuildvrt和gdalwarp都没有成功。 有使用现有的gdal命令行实用程序执行此操作的简单方法吗?在使用NetCDF Python API求助于自定义解决方案之前,我想过要问一下。 谢谢。-大卫 在2013年9月3日,星期二,上午10:15,Etienne Tourigny写道: 您想要的可能超出了gdal的范围。它将需要一些聪明的元数据管理,以便gdal_translate将它们放在一个文件中... 我建议您使用gdal_translate将它们全部转换为netcdf,然后使用python-netcdf4(而不是numpy / scipy中的那个)将它们堆叠在时间维度上。 2013年9月3日,星期二,上午7:55,“ Signell,Richard”写道: 大卫,如果您将问题发布到GIS stackexchange组 /gis//上, 我将提供示例代码,该代码应该会有所帮助。 -丰富 ==================== PDT更新9/3/13 17:04 这是我的输入数据集之一的gdalinfo输出: gdalinfo 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Driver: GTiff/GeoTIFF Files: 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Size is 10666, 13387 Coordinate System is: PROJCS["unnamed", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Polar_Stereographic"], …
12 python  gdal  netcdf 

3
使用Python中的shapefile多边形从NetCDF中提取数据[重复]
这个问题已经在这里有了答案: 如何在Python中从shapefile屏蔽NetCDF时间序列数据? (1个答案) 15天前关闭。 我需要使用特定的shapefile从NetCDF子集数据。数据是1/4度分辨率下的海面温度和海洋颜色。我有4个描述美国的多边形。东北大陆架大型海洋生态系统及其提取数据所需的子组件。我正在处理1982-2014年的每月复合文件,因此此数据提取例程需要自动化。文件已被子集为[35、45,-80,-60]的近似工作区域网格。 以前,我们将HDF5数据文件转换为R中的栅格并以这种方式进行处理,但是这种方法的效率很低,我相信使用当前NetCDF文件在Python中有更好的解决方案。 到目前为止,我一直在使用GDAL和Fiona读取shapefile,并使用NetCDF4加载数据文件。我不确定如何对数据进行子集。我找到了这个: GDAL for Python:从NetCDF文件提取子域吗? 但是我不知道如何使用除简单边界框以外的任何东西来对NetCDF文件进行子集化,这些多边形最肯定不是。 多边形例程中的点可能需要花很长时间才能工作,但是也许我可以使用较小的边界框将数据子集化,使其旋转以适合如下所示的这些形状作为初始起点,然后进行多边形点搜索: 使用lon / lat边界框对曲线netCDF文件(ROMS模型输出)进行子设置。 有任何想法吗? 编辑1: 我刚刚遇到了OpenClimateGIS软件包,它似乎很适合该法案...我将继续研究它是否可以正常工作:http: //ncpp.github.io/ocgis/examples。 html#advanced-subsetting

1
我如何解决Gdalwarp错误“太多点无法转换”以将对地静止映射为Lambert保形?
我正在尝试使用gdalwarp将对地静止映射为Lambert保形。我的输入数据在netcdf中,并且在地理坐标(度)中,我想将重新映射的数据输出到netcdf。我为输入的netcdf数据创建了一个相应的vrt文件。Gdalwarp将输出netcdf文件,但输出数据全为零,并且我收到以下错误: Creating output file that is 5120P x 5120L. Processing input file netcdf.vrt. ERROR 1: Too many points (441 out of 441) failed to transform, unable to compute output bounds. Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping. 0...10...20...30...40...50...60...70...80...90...100 - done. 我尝试了以下命令: /usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 …

3
使用GDAL将NetCDF文件转换为GeoTIFF?
我有一个NetCDF文件,其中包含英国下雨天数的数据。我想将其转换为GeoTIFF,但无法正确获取GDAL地理参考。请参阅下面来自gdalinfo的转储。 我想要一个显示“ raindays_tmean_abs”变量的TIFF乐队。“ lat”和“ lon”变量似乎具有合理的最小/最大值,可以将数据放到整个英国,但是当我使用gdal_translate进行转换时,它最终出现在非洲西海岸。 关于第一个尺寸不是“时间或垂直”尺寸的警告。我读到gdal期望数据按一定顺序排列,所以也许这就是问题所在,但是我不知道如何修改它。 如何转换/重新投影此数据? gdalinfo NETCDF:RainDays1_1961-1990_LTA_25km.nc Warning 1: dimension #0 (meaning_period) is not a Time or Vertical dimension. Driver: netCDF/Network Common Data Format Files: RainDays1_1961-1990_LTA_25km.nc RainDays1_1961-1990_LTA_25km.nc.aux.xml Size is 39, 52 Coordinate System is `' Origin = (-17.650008352179277,11.329999312758446) Pixel Size = (0.219999614514803,-0.219999998807907) Metadata: axis_0#units= meaning_period#comment=Meaning periods are all …

3
此netCDF栅格是什么坐标系?
我获得了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 …
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.