qq_37935516的博客

1.Assignment 5a：Create an NDVI image Read in data from aster.img Create an NDVI image Write out NDVI to new file Can do entire image at once or block by block  Don't forget to calculate statist...

      利用GDAL读入DEM与Landsat影像，由于DEM是WG84坐标系，Landsat是WGS84坐标系UTM投影，因此处理在实际应用中需要将DEM进行投影转换。大概分为以下几个步骤：读取DEM，读取Landsat影像 获取Landsat影像的投影信息，将其赋给DEM，并对DEM进行重采样 计算dx和dy 计算坡度和坡向 输出坡度和坡向的影像 from osge...

      之前写过一篇【GDAL】python读取txt影像名称文件生成shp，本篇博客是对之前那篇的扩展。之前那篇是把所有.gz文件名写在了一个txt文件里，通过读取txt文件将里面的文件名保存到一个文件名列表中即可。但是这篇博客是读取一个文件夹下所有文件夹中的所有.gz压缩包，将这些压缩包的文件路径名存放在一个列表中解析出文件名路径、XY、传感器类型、日期和产品ID，创建一个点状shape文...

1.在https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal下载GDAL安装包，拷贝至相应的环境中，打开shell 输入以下代码安装GDALpip install GDAL‑2.3.2‑cp35‑cp35m‑win_amd64.whl2.在Anaconda Prompt中输入：activate tensorflowconda instal...

1.属性过滤器>>>import ogr,os>>>os.chdir('E:/data/GDAL/ospy_data3')>>>driver=ogr.GetDriverByName('ESRI Shapefile')>>>ds=driver.Open('sites.shp',0)>>>l

1.创建点状要素：import ogrimport osos.chdir('E:/data/GDAL/ospy_data2')driver = ogr.GetDriverByName('ESRI Shapefile')if os.path.exists('out.shp'): driver.DeleteDataSource('out.shp')ds = driver.Cr...

1.根据坐标读取遥感影像的单个像素值# week 4: get pixel values at a set of coordinates by reading in one pixel at a timeimport os, sys, timefrom osgeo import gdalfrom gdalconst import *# start timingstartTime ...

例如设计一个3 x 3的滑动窗口，写算法执行就有两种方式：1.pixel by piexl每个进行逐像素运算，效率太低，速度慢2.使用 slice切片形式循环，效率高，速度快 两个作业就是分别用pixel和slice方式完成高通滤波操作进行对比1.Assignment 6aUse a 3x3 high pass filter to detect edges in ban...

学习资料：犹他州立大学：https://www.gis.usu.edu/~chrisg/python/2009/lectures/ospy_slides1.pdf 开放地理空间实验室 http://www.osgeo.cn/python_gdal_utah_tutorial/ch02.html GDAL官方网站：https://www.gdal.org/ # Reading and...

读取一个由影像名称组成的txt文件，获得文件名称中的经纬度，并生成shptxt文件结果import ogrimport os# read txt and return a filename listdef readFile(filename): f = open(infile, 'r') lines = f.readlines() point...