x2434417239的博客

Tips：采用梯度下降法实现线性回归，同时可批量绘制散点图和拟合曲线以及损失图。样例数据：拟合曲线图：损失曲线：Python代码：#!/usr/bin/env python# -*- coding: utf-8 -*-"""Created on Fri Jan 22 15:05:16 2021E-mail: 1321378489@qq.com@author: xiao_gf"""import mathimport numpy as npimport pandas a

Tips：基于GDAL库的gdal.Warp工具实现栅格数据的批量裁剪功能。#!/usr/bin/env python# -*- coding: utf-8 -*-"""Created on Fri Jan 8 16:31:59 2021E-mail: 1321378489@qq.com@author: xiao_gf"""import time,osfrom osgeo import gdalimport numba @numba.jit # 程序加速def clip

Tips: 通过ArcGIS可实现提取，本代码主要是批量处理。仅供参考~代码适应样例说明：数据：2000-2018年逐年的全国县级人口数据，例如：2000_population.shp、2001_population.shp等要求：分布提取每年每个省份的县级人口数据数据属性表：字段包括：FID、population、cityName等最后结果：输出19个文件夹，每个文件夹中包括各省份的人口数据#!/usr/bin/env python# -*- coding: utf-8 -*-"""C

Tips ：本代码实现ArcGIS中的区域统计功能（Zonal），将统计的属性值直接赋予矢量文件属性表，可实现批量统计。import ogr, osfrom rasterstats import zonal_statsdef zonal(input_vector,in_raster,Stats_type,raster_name): shp = ogr.Open(input_vector,1) lyr = shp.GetLayer()

Tips ：本代码实现ArcGIS中的区域统计功能（Zonal），有两点特色：1）能够将统计得到的属性值直接赋予矢量文件中；2）可实现多波段影像区域统计（方法比较笨，先拆分波段，然后进行统计）。区域统计的类型参考ArcGIS：MEAN— Calculates the average of all cells in the value raster that belong to the same zone as the output cell.MAJORITY — Determines the va

#nc数据批量转为tiff数据详细代码（Python代码）# -*- coding: utf-8 -*-"""Created on Wed Dec 2 11:12:50 2020@author: xiao_gf"""import osimport numpy as npimport netCDF4 as ncfrom osgeo import gdal,osrdef strstr(str1, str2): #查找某一字符是否存在于另一字符串中，str2为需要查找的字符，str1

中国气象数据网的.nc数据批量下载（Python代码）# -*- coding: utf-8 -*-"""Created on Thu Dec 10 13:55:39 2020@author: xiao下载链接列表格式样例：http://101.201.177.119/dlorder/orderData/xxx/Z_NAFP_C_BABJ_20191029012028_P_CLDAS_RT_CHN_0P0625_HOR-RSM000010-2019102900.nc?Expires=xxx&