首页 文章
  • 1 votes
     answers
     views

    Python中Lat / Lon Meshgrids的GDAL仿射系数

    在创建新的GeoTIFF文件时,我遇到了一些Affine转换系数问题 . 我正在做的是在一个科学数据集上的ETL,它产生一个2D Ndarray以及一组包含Lat和Lon的meshgrid Ndarray . meshgrids和数据集数组都具有645 x 980的相同尺寸 . 据我所知,GeoTIFF需要一个Affine系数列表,当通过 SetGeoTransform() 方法从Python ...
  • 2 votes
     answers
     views

    使用@ rpath / libpoppler.71.dylib找不到图像

    目前我已经安装了anaconda . 我试图在我的脚本中调用osgeo.ogr模块来打开shapefile,但是我有以下错误: ImportError:dlopen(/Applications/anaconda3/lib/python3.6/site-packages/osgeo/_gdal.cpython-36m-darwin.so,2):未加载库:@ rpath / libpoppler.7...
  • 1 votes
     answers
     views

    GDAL2Tiles和向量

    我正在尝试制作一个非常大的仓库的交互式 Map . 我扫描了蓝图,并使用OpenLayers和TMS数据源(使用GDAL2Tiles生成)来显示 Map . 一切看起来都很棒,但我完全不知道如何在它上面创建一个矢量图层! 理想情况下,我的仓库顶部有大约100个左右的向量,每个向量都可以点击以显示其他信息,等等等等等等......但是我用什么来制作向量,我如何参考由于我不使用纬度/经度坐标,因此...
  • 6 votes
     answers
     views

    使用geotiff中的gdal python读取高程

    我正在使用GDAL加载geotiff文件 . 我设法读取坐标X,Y但不是高程 . 有没有人曾经在类似案件上工作过? 问候,
  • 6 votes
     answers
     views

    使用gdal从C#中的现有位图创建地理基准

    我必须在c#项目中使用gdal . 我要做的是将一个简单的位图“转换”为GeoTiff . 我在gdal网站上阅读了一些文档,但是我没有成功完成它 . 事实上,我的位图已成功导出到地理位置,但如果我使用GIS软件(例如QuantumGIS)打开geotiff,则GeoTiff会在y轴上反转: : 而原始位图看起来像这样: 这就是我所做的: 首先我将一个临时文件写入磁盘(即位图),我创建了一个包含...
  • 0 votes
     answers
     views

    Postgres - 将bytea PNG列转换为GeoTiff

    我已经将大量的磁贴加载到postgres数据库中,以便在磁贴服务器中使用 . 这些都被加载到PNG格式的bytea列中 . 我现在发现tile服务器代码需要这些格式为GEOTiff格式 . 命令:- gdal_translate -of GTiff -expand rgb -co COMPRESS = DEFLATE -co ZLEVEL = 6 完美的工作 . 但是,远程服务器上已经加载了大量数...
  • 1 votes
     answers
     views

    显示gdal2tiles输出openlayers 3.15.1

    我用gdal对图像栅格进行了地理配准,并使用EPSG创建了一个新的geotiff:3857 . 现在我想在我的网站上显示它,我有两个问题,可能是因为我不太了解gdal,特别是gdal2tiles的工作原理 . 我使用gdal2tiles.py -p mearcator z 1-8 source.tif Directory_tile,我试图通过openlayers 3.15.1显示这个图块....我...
  • 0 votes
     answers
     views

    为什么gdal无法打开geotiff?

    我是python和地理处理的新手 . 我正在写一些程序来计算ndwi . 为此,我尝试使用gdal打开geotiff数据集,但无法打开数据集 . 我试图打开不同的TIFF文件(Landsat8多个数据,Landsat7复合材料等),但数据集总是无 . 这可能是什么原因?或者我怎么能找到它? 这是代码的一部分: import sys, os, struct import gdal, gdalcons...
  • 1 votes
     answers
     views

    Python-gdal用二进制颜色和NaN编写GeoTiff

    我正在从numpy数组生成一个简单的两类(二进制)geotiff,其中包含3个值:1,2和NaN . 我想用不同的颜色显示geotiff所以我使用了颜色表 . 但是,颜色表仅支持Byte或UInt16数据类型,它将NaN转换为零 . 有没有办法写一个二进制geotiff包括带颜色的NaN?谢谢 . filename = 'test.tif' ds = gdal.Open(filename, gda...
  • 3 votes
     answers
     views

    PySpark 和栅格(GeoTiff)数据

    我正在尝试使用 Spark(PySpark)来对存储在 multi-band GeoTiffs 中的数据进行分析。我还是 Spark 的新手。 设定: Geotiffs 足够小,可以在纯 python 中运行-特别是我正在使用 gdal 读取数据。然后创建数据框并进行分析。 但是分析需要一段时间。而且,在反复的基础上,我将要分析数百个 Geotiff-输入 PySpark。 问题: 我编写了允...
  • 3 votes
     answers
     views

    在Ubuntu 16.04上使用MS SQL的OGR2OGR

    我需要从Ubuntu Server 16.04访问带有OGR2OGR的Microsoft SQL Server . 它正在使用Microsoft,所以基本设置很好 . 我的问题是,当我运行我的OGR2OGR命令时(在使用sudo apt-get install gdal-bin安装GDAL之后): ogr2ogr -overwrite -f MSSQLSpatial "MSSQL:ser...
  • 0 votes
     answers
     views

    在矩阵中查找更高数字的区域

    我正在研究一个项目来检测航拍图像中的某些物体,作为其中的一部分,我正在尝试利用图像的高程数据 . 我正在使用数字高程模型(DEM),基本上是高程值矩阵 . 例如,当我试图检测树木时,我想搜索高于周围地形的树形区域 . 以下是DEM热图中树的示例: https://i.stack.imgur.com/pIvlv.png 我希望能够找到比周围环境更高的小区域 . 我正在使用OpenCV和GDAL进行实...
  • 0 votes
     answers
     views

    JPEG压缩的Geoserver图像金字塔

    我有一个Image Pyramid(使用gdal_retile创建),我想导入Geoserver(2.8.2,使用Image Pyramid插件)并发布为WMS . 我已经能够成功地使用GeoTIFF图像,但是因为我们存储了数TB的图像,所以我们想使用JPEG压缩的TIFF . 如果我使用压缩参数运行gdal_retile,例如: @python%OSGEO4W_ROOT%\ bin \ gda...
  • 2 votes
     answers
     views

    使用gdal压缩GeoTIFF和GeoServer

    我有一个很大的GeoTIFF,我希望通过GeoServer(v.2.11)中的WMS进行流式传输 . 图像的大小约为7GB,包含非常大的高分辨率RGB图像 . 我已经在JVM中允许足够的堆空间来显示图像 . 但是,我想压缩图像,以便在浏览时更具响应性,因此它将分配更少的内存 . 我已经遵循了一些建议here . 我的策略是 compress the GeoTIFF with JPEG compre...
  • 1 votes
     answers
     views

    由gdal2tiles与传单生成的自定义栅格

    我能够使用gdal2tiles.py生成一个tiles文件夹,它允许我在 Map 上显示geotiff文件: 我生成的所有geotiff文件都是灰度级的,我希望能够为它们着色并添加一个图例来查看哪个值对应于哪种颜色 . 这是我的实际代码,它允许我这样做: <div id="map"></div> <button id="popul...
  • 2 votes
     answers
     views

    rgdal有效地读取大型多波段栅格

    我正在使用rgdal包在R中处理图像分类脚本 . 所讨论的光栅是具有28个通道的PCIDSK文件:训练数据通道,验证数据通道和26个光谱数据通道 . 目标是填充包含训练数据信道中不为零的每个像素的值的数据帧,以及26个频带中的相关频谱值 . 在Python / Numpy中,我可以轻松地将整个图像的所有波段导入到多维数组中,但是,由于内存限制,R中的唯一选项似乎是逐块导入此数据,这非常慢: lib...
  • 1 votes
     answers
     views

    来自foreach的神秘错误

    在我第一次涉足foreach时,我收到了这个错误: Error in function (handle) : no function to return from, jumping to top level 我为每个 Worker 得到其中一个 . 我正在使用doMC /多核后端 . 起初我以为这是由我的foreach调用中的函数试图作为副作用发送到控制台的文本进度条引起的,但是在关闭之后...
  • 7 votes
     answers
     views

    在mingw64下使用GetACP链接错误(mingw-builds)

    我试图使用mingw64(来自http://sourceforge.net/projects/mingwbuilds/files/host-windows/ x64-4.8.0-release-posix-seh-rev2.7z)构建gdal-1.10.0(http://trac.osgeo.org/gdal/wiki/DownloadSource) . 我已经在标准的MinGW(32位)版本下编...
  • 1 votes
     answers
     views

    Amazon Linux AMI上的Rgdal软件包安装

    我已经安装了geos,proj-4.8.0和gdal-1.9.2,但是在安装R包rgeos后,rgdal无法安装 . 这是我的主目录: . . . 和我运行的错误 install.packages("rgdal", type="source") : 构建软件包索引安装小插图**测试是否可以加载已安装的软件包dyn.load中的错误(文件,DLLpa...

热门问题