小鱼的博客

一、Windows 环境1、使用condaconda install gdal2、使用whl下载对应版本的gdalhttps://www.lfd.uci.edu/~gohlke/pythonlibs/#gdalpip install GDAL‑3.2.3‑pp37‑pypy37_pp73‑win_amd64.whl注意whl的路径3、使用pip（不推荐）该方法经常出错pip install gdal二、Linux（Ubuntu）1、使用condaconda install -

想直接用的代码，找不着了。。。只能翻出以前写的代码，记录一下# 像素坐标转地理（投影）坐标def pixel2geo(geotrans, px, py): ulX = geotrans[0] ulY = geotrans[3] pixelWidth = geotrans[1] pixelHeight = geotrans[5] rotate1 = geotrans[2] rotate2 = geotrans[4] Xgeo = ulX + px

使用OGR的图层筛选Layer.SetAttributeFilter()和Layer.SetSpatialFilter()，筛选的结果使用for feat in Layer遍历，不能直接使用Layer.GetFeature()获取筛选后的要素

栅格数据赋值矢量数据前段时间给老师写栅格数据赋值矢量的代码，要用到GDAL包，网上搜了一遍没有找到合适的，只在GDAL的一个参考文档中找到一个类似的代码Calculate zonal statisticsimport gdal, ogr, osr, numpyimport sysdef zonal_stats(feat, input_zone_polygon, input_value_raster): # Open data raster = gdal.Open(input_

OGR坐标转换问题之前使用OGR工具做坐标转换，一直没有问题，后来许久不用，再使用就出现无法转换的问题。第一个问题Linux与Windows环境下坐标顺序不同，上图。linux环境下windows环境下两个环境中创建点要素的坐标顺序不同，linux是经度、纬度，而windows是纬度、经度。如果相反，创建后再坐标转换会出现错误。（具体是什么原因，还没有去查，有空闲再看看）第二个问题。（这个问题之前一直没有注意）在Shapefile文件中读取的要素几何类型是带有空间参考的，而自定义的几