我想使用python、gdal和numpy创建一些高程/高度场栅格。我被困在numpy(可能还有python和gdal。)在numpy中，我一直在尝试以下操作:>>>a=numpy.linspace(4,1,4,endpoint=True)>>>b=numpy.vstack(a)>>>c=numpy.repeat(b,4,axis=1)>>>carray([[4.,4.,4.,4.],[3.,3.,3.,3.],[2.,2.,2.,2.],[1.,1.,1.,1.]])#ThisistheelevationdataIwant从osgeo导入gdal从gdalconst导入*>>>for
为了扩展OpenOffice的功能，我在userscriptdirectory的文件中定义了一些Python宏。(~/Library/ApplicationSupport/OpenOffice.org/3/user/Scripts/python/，在我的例子中)。宏在Python宏管理器中可见。但是，在单元格公式中使用这些函数会导致“#NAME?”(OO.orgerror525)。假设我定义了以下函数:defpytype(val):returnstr(type(val))如何在单元格公式中调用pytype(例如=PYTYPE("string"))？背景我正在将一些数据从Authoriz
我正在寻找更改GDAL栅格数据集分辨率的最佳方法。例如，我有一个像素大小为(30,-30)的栅格，我想将像素大小更改为(5,-5)，将给定像素的所有值插值到输出中光栅。因此对于输入栅格的每个像素，我希望输出栅格中有36个像素都共享相同的值。如果我运行gdalwarp-tr5-5inputRaster.tifoutputRaster.tif，我会得到我正在寻找的结果，所以我假设我应该能够复制此功能具有一些GDAL功能。如果可能的话，我宁愿避免调用python的Subprocess类。 最佳答案 您需要重新投影光栅。例如，从交互式Pyt
我一直在尝试设置我的Windows计算机，以便我可以拥有一个带有PostGIS扩展名的本地postgreSQL。安装了这个之后，我希望能够在将其放入云之前在本地使用geodjango创建一个项目。我已经在我的本地机器上使用Django工作了一段时间，现在使用SQLiteDB，但是由于下一个项目将部分基于基于坐标的数据，所以我想设置正确的环境。导入说明:我已经安装了mini-conda以在单独的环境中运行。不过，我在工作时会激活这个“开发”环境我已尝试在线关注大部分geodjango信息/教程，但无法正常工作。我所做的(主要是遵循:https://docs.djangoproject.c
我编写了一个调用gdal_calc.py的R函数来计算RasterStack(一系列输入栅格文件)中的像素最小值。我这样做是因为它比大型栅格的raster::min快得多。该函数适用于最多23个文件，但在传递24个或更多文件时会引发警告，并返回一个用零填充的输出栅格。由于R只是准备对pythongdal_calc.py进行系统调用，所以这个问题并不特定于R，我鼓励python/numpy爱好者继续阅读。这是函数。最终gdal_calc调用的结构显示在本文底部的有问题的用法引发的警告消息中。gdal_min26||nbands>26)stop('Maximumnumberofinputs
如何从python获取已安装的GDAL/OGR版本？我知道gdal-config程序并且目前正在使用以下程序:In[3]:importcommandsIn[4]:commands.getoutput('gdal-config--version')Out[4]:'1.7.2'但是，我怀疑有一种方法可以使用pythonAPI本身来执行此操作。有骰子吗？ 最佳答案 gdal.VersionInfo()做我想做的事:>>>osgeo.gdal.VersionInfo()'1604'这适用于我的Windows机器和Ubuntu安装。gdal.
我尝试在virtualenvwrapper环境中安装gdal(pipinstallgdal)，但出现此错误:error:command'gcc'failedwithexitstatus1----------------------------------------FailedbuildingwheelforgdalFailedtobuildgdal我也试过“pipinstall--no-installGDAL”但是没有选项--no-install我该怎么办!？ 最佳答案 是的，在venv中安装GDAL很麻烦。方便的是，我刚刚为我的
关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题？更新问题，以便editingthispost可以用事实和引用来回答它.关闭8年前。Improvethisquestion我已经为Excel中的各种内容编写了一堆VBA代码。我正在考虑迁移到libreOffice。在Tool->Macros->OrganizeMacros下:两个选择是LibreOfficeBasic和Python。我应该学习其中之一，两者还是其他。我完全是在浪费时间吗？任何建议表示赞赏。
澄清:我以某种方式遗漏了关键方面:不使用os.system或子进程-仅使用pythonAPI。我正在尝试转换NOAAGTX偏移网格的一部分以进行垂直基准转换，而不是完全遵循如何在GDAL中使用python执行此操作。我想要一个网格(在本例中是BathymetryAttributedGrid，但它可能是一个geotif)并将其用作我想要做的模板。如果我能做到这一点，我觉得这将极大地帮助人们利用这种类型的数据。这是我所拥有的，但绝对不起作用。当我在生成的目标数据集(dst_ds)上运行gdalinfo时，它与源网格BAG不匹配。fromosgeoimportgdal,osrbag=gdal
我想使用matplotlibbasemap绘制光栅tiff(download-723Kb)。我的光栅投影坐标以米为单位:In[2]:path=r'albers_5km.tif'raster=gdal.Open(path,gdal.GA_ReadOnly)array=raster.GetRasterBand(20).ReadAsArray()print('RasterProjection:\n',raster.GetProjection())print('RasterGeoTransform:\n',raster.GetGeoTransform())Out[2]:RasterProjec
