SunStrongInChina的博客

关于Matplotlib中文乱码的问题一直困扰着我这个朴朴素素的制图人员，我也是试来试去各种方法无果之后，发现了一个最佳的解决方案。首先我列举一些市面上并不一定有效的解决方案:import matplotlib.pyplot as pltimport seaborn as snsplt.rcParams['font.sans-serif'] = ['SimHei'] # Matplotlib中设置字体-黑体，解决Matplotlib中文乱码问题plt.rcParams['axes.unic

1.安装Qgis或者python的gdal库，确保电脑上有两个工具，gdalbuidlvrt以及gdaltranslateWindows+R,输入cmd，然后输入 cd /d +Qgis的bin目录 导航到Qgis的bin目录下cd /d D:\Qgis\bin2.将所有需要拼接的影像放到一个文件夹下，大致如下图所示：3.构建vrt，前一个参数为vrt路径，后一个为tif文件路径gdalbuildvrt D:\\merged_img\\merged.vrt D:\\merged_im

sen+mk python实现代码免费共享-----赶紧收藏吧python开源社区公布了进行sen+mk趋势性检验的官方包，有关该官方包的主要内容详见：https://github.com/Coder2cdb/pyMannKendall，开源不易，有些人还收费查看，我觉得没必要，本次代码免费分享给大家，有需要的收藏点赞！本次sen+mk python代码输出参数说明'' trend: tells the trend (increasing, decreasing or no trend)

rasterio库中索引像元坐标1.根据像元所在行列号查询坐标以及反查询rasterio库提供了专门用于行列号及像元中心坐标转换的函数，其中dataset.xy可将像元的行列号转化为像元中心坐标，而dataset.index可将像元中心坐标转化为行列号import rasterio as rspath=r"C:\Users\Lenovo\Desktop\test_images\images_doys2010_127_37.tif"dataset=rs.open(path)#获取波段数>&

Python使用gdal实现影像镶嵌如果要对某个文件夹下的多景影像进行镶嵌，我们需要知道镶嵌后影像的行列数目，以及该影像的六个地理变换参数，（值得特别注意的是，无论是影像的重采样还是镶嵌，都需要特别关注影像的6个地理变换参数），关于这六个地理变换参数，请参考我的另一篇博文：https://blog.csdn.net/SunStrongInChina/article/details/1042629...

gdal实现影像的重采样对于tif格式的遥感影像而言，通常使用左上角那个像元真实的x坐标值，像素宽度，x像素旋转角度，原点y坐标，y像素旋转角度，像素高度（负值）共六个参数来表示整幅影像像素点的坐标。当影像分辨率改变后，x坐标值和y坐标值度都不会改变，但像素高度和宽度会改变，而像素高度和宽度改变的话，影像的行列数也会改变。1.gdal实现影像重采样至高分辨率这里，将一幅影像的像元高度和宽...

osr和pyproj库的使用坐标系统的构成无论是矢量数据处理还是栅格数据的处理都需要用到坐标系统，坐标系统的构成主要是三部分，一、椭球体，可以用两个非常简单的参数来表示，a长半轴长度，b短半轴长度（也可以表述为长半轴a，短半轴b，扁率e三者知道任意两个也可表述该坐标系统），二、坐标系，为了描述地球上某一点的位置，人们想出了不同的坐标系统，如非常有名的空间直角坐标系，三、投影，为了制图以及量测的...

Python实现csv文件转换为shapefile文件（最简便方法）python使用ogr实现csv转shapefile文件使用python实现csv转化为shapefile文件的方法有很多，我大致看了下csdn的博主们的分享 ，大部分使用了csv库和shapefile库进行csv文件的读取和shapefile文件的写入，我觉得还是略过于复杂了，这里我提供一个最为简便的方法实现csv文件转换为...