摘要: import sklearn sorted(sklearn.metrics.SCORERS.keys()) Out[111]: ['accuracy', 'adjusted_mutual_info_score', 'adjusted_rand_score', 'average_precision', 阅读全文
posted @ 2021-06-29 09:25 vivid_autumn 阅读(368) 评论(0) 推荐(0) 编辑
摘要: z=zscore(active); M=cov(z) ;%协方差 [V,D]=eig(M); %求出协方差矩阵的特征向量、特征根 d=diag(D); %取出特征根矩阵列向量(提取出每一主成分的贡献率) eig1=sort(d,'descend') %将贡献率按从大到小元素排列 v=fliplr(V 阅读全文
posted @ 2020-06-21 11:09 vivid_autumn 阅读(184) 评论(0) 推荐(0) 编辑
摘要: library(installr)copy.packages.between.libraries(from = "F:\\R\\R-3.3.2\\library" , to = "F:\\R\\R-3.5.0\\library" , ask =T, keep_old = TRUE , do_NOT_ 阅读全文
posted @ 2020-06-21 10:46 vivid_autumn 阅读(967) 评论(0) 推荐(0) 编辑
摘要: https://sklearn.apachecn.org/1.离群值(点)识别离群值(outlier)是指在一组数据中出现的与大部分数值相比差异较大的个别值。关键问题:差异有多大、如何判断? ( 1)利用直方图或盒状图直接判断置信水平, 显著偏离直方图主体频数区域的数值;或者距离盒状图的箱体(第 1 阅读全文
posted @ 2019-11-20 14:44 vivid_autumn 阅读(487) 评论(0) 推荐(0) 编辑
摘要: 1、本篇内容包含两个部分,一是使用PROJ4包对点进行投影转换,二是栅格数据投影转换的示例 2、 #3\另外一个投影包PROJ4from pyproj import Proj,Geod,transform#projection1:UTM zone15, grs80 ellipse, NAD83 da 阅读全文
posted @ 2019-10-03 20:12 vivid_autumn 阅读(719) 评论(0) 推荐(0) 编辑
摘要: 1、 (1)重点学习如何应用 gdal.RasterizeLayer 函数 gdal.RasterizeLayer( target_ds, [3, 2, 1],source_lyr, burn_values = [10,10, 55], options = ["ALL_TOUCHED=TRUE"]) 阅读全文
posted @ 2019-10-03 19:36 vivid_autumn 阅读(525) 评论(0) 推荐(0) 编辑
摘要: 1、除了一些转为shapefile文件必要的设置外,python读取csv文件也有一些需要注意的地方。 (1)首先读取了第一行fds作为字段名 (2)之后以此读取文件内容存在data里,这部分代码比较冗长 (3)代码容易出错的地方在于ds = driver.CreateDataSource(file 阅读全文
posted @ 2019-09-30 16:01 vivid_autumn 阅读(729) 评论(0) 推荐(0) 编辑
摘要: 1、这部分的代码需要在前面“读ArcGIS shape文件”的代码后执行,先读后写,需要借用之前fieldlist等变量。 2、代码示例如下 from osgeo import ogr,gdal,osrimport osgeoimport os# 解决中文路径osgeo.gdal.SetConfig 阅读全文
posted @ 2019-09-30 09:15 vivid_autumn 阅读(300) 评论(0) 推荐(0) 编辑
摘要: 1、读shapefile文件主要读以下内容,包括spatialRef投影信息,layerDefn图层定义信息,geomType几何对象类型,fieldDefn字段定义信息。geomlist是得到了每一个feature的geometryRef后转为Wkt形式表示的坐标点位,reclist通过name属 阅读全文
posted @ 2019-09-30 08:03 vivid_autumn 阅读(2070) 评论(0) 推荐(0) 编辑
摘要: #创建 SF 对象from osgeo import ogr#SF 对象的 WKT 表示格式pnt_wkt = "POINT(1 1)"lin_wkt = "LINESTRING(0 0,15 15) "pol_wkt = "POLYGON((1 1,1 10,10 10,10 1,1 1),( 2 阅读全文
posted @ 2019-09-20 15:09 vivid_autumn 阅读(233) 评论(0) 推荐(0) 编辑