gdal
gis_rc
这个作者很懒,什么都没留下…
展开
-
MODIS MCD19A2气溶胶数据(AOD)利用python rasterio转换成tif,并将投影转换成经纬度
1.先使用命令行将数据转成tif,并转换经纬度语句:gdalwarp -t_srs "EPSG:4326" HDF4_EOS:EOS_GRID:"MCD19A2.A2021215.h10v04.006.2021217045137.hdf"::grid1km:Optical_Depth_055 test.tif2.然后还没有结束,你就会得到一个有6个波段的tif数据。经过查看发现,这个六个波段合成的结果就是这整个的数据3.使用python rasterio将数据进行一个融合输出代码:impor原创 2021-12-11 16:55:23 · 2174 阅读 · 10 评论 -
使用python+gdal实现arcgis的核密度分析,针对点密度
核密度分析的工作原理:核密度分析工具用于计算要素在其周围邻域中的密度。此工具既可计算点要素的密度,也可计算线要素的密度可能的用途包括针对社区规划分析房屋密度或犯罪行为,或探索道路或公共设施管线如何影响野生动物栖息地。可使用 population 字段赋予某些要素比其他要素更大的权重,该字段还允许使用一个点表示多个观察对象。例如,一个地址可以表示一栋六单元的公寓,或者在确定总体犯罪率时可赋予某些罪行比其他罪行更大的权重。对于线要素,分车道高速公路可能比狭窄的土路具有更大的影响。核密度分析用于计算每个输出栅原创 2020-12-23 10:26:18 · 3247 阅读 · 13 评论 -
使用Python的gdal库读取tif格式遥感图像并将其切割为多个小图(按坐标系正确输出)
# -*- coding: utf-8 -*-import osimport numpyfrom osgeo import gdalclass GRID: # 读图像文件 def read_img(self, filename): dataset = gdal.Open(filename) # 打开文件 im_width = dataset.RasterXSize # 栅格矩阵的列数 im_height = dataset.R原创 2020-11-21 10:45:07 · 3157 阅读 · 25 评论 -
使用python gdal将影像重采样之取和值(sum)
我们都知道gdal为我们提供了多种重采样的方法,这里我们不再叙述,有兴趣可以自己去看重采样方法下面我们来介绍怎样进行取和值重采样,针对少数特别的需求。思路把整个图像看成一张网格,一个个的像元,我们要改变格子的大小,将小的格子合并成大的或者将大的分割成小的,我们这里讲合并成大的,因为我们要取和值;使用gdal将图像的矩阵数组读出,取得这些格子的数据,一个矩阵,然后进行计算;使用numpy,将矩阵进行分割,分割成我们想要的想要大小,对应矩阵的大小。均分:np.hsplit(b1, 3) #原创 2020-08-16 11:56:22 · 575 阅读 · 0 评论