RS-GIS
菜鸡的自我拯救
菜,无止境
展开
-
GEE批量下载:控制台代码
function runTaskList() { var tasklist = document.getElementsByClassName('awaiting-user-config'); for (var i = 0; i < tasklist.length; i++) tasklist[i].children[2].click(); } function confirmAll() { var ok = document.getElemen..原创 2022-03-22 17:10:15 · 429 阅读 · 0 评论 -
用R语言实现spatial join
今天实现了一个任务:有一个library(rgdal)library(gstat) # Use gstat's idw routinelibrary(sp) # Used for the spsample functionlibrary(raster) # Used to clip out thiessen polygonsimos=read.csv('D://research//GHM//imos.csv')coordinates(imos) <- c("lon","原创 2021-09-08 23:25:59 · 628 阅读 · 0 评论 -
GEE统计区域气溶胶差值
var shp=tablevar shp1 = shp.filterMetadata('Result_GWR','less_than','30000');var shp2=shp.filterMetadata('Result_GWR','not_less_than','30000');//Map.addLayer(shp1)// var images = imageCollection4.filterBounds(shp1)// .filter(ee.Filter.date('20.原创 2021-01-23 22:14:02 · 475 阅读 · 0 评论 -
使用Arcpy批量进行可达性分析
# Name: MakeODCostMatrixAnalysisLayer_Workflow2.py# Description: Find the travel time to the closest hospital from each census# tract and join the travel time and hospital name to the input# tracts.# Requirements: Network An原创 2020-12-03 01:25:28 · 488 阅读 · 0 评论 -
使用R语言的spgwr包进行地理加权回归(GWR)
今天需要做要给GWR回归,使用了R语言中的spgwr包。做的工作主要有;1)安装R语言直接去R官网就好。2)安装spgwr包> install.packages("spgwr")发现报错Warning in install.packages : 'lib = "C:/Program Files/R/R-4.0.3/library"' is not writableWarning in install.packages : cannot create dir 'D:\onedri原创 2020-11-19 01:04:59 · 6993 阅读 · 3 评论 -
使用arcpy批量计算路网可达性(网路分析方法)
# Name: MakeODCostMatrixAnalysisLayer_Workflow2.py# Description: Find the travel time to the closest hospital from each census# tract and join the travel time and hospital name to the input# tracts.# Requirements: Network An原创 2020-11-18 02:22:07 · 966 阅读 · 0 评论 -
计算各行政区划内的人口加权质心
本研究的目的是使用人口栅格数据计算每个行政区划内的人口加权质心。原创 2020-11-14 05:22:46 · 1787 阅读 · 0 评论 -
问题解决:ArcGis的shp连接表格时,显示不出对应字段名称
如图,想要用CSV匹配shp(含同一字段名称),但是在join时却显示不出字段名称,实际上也匹配不上数据解决方法:将改CSV文件导出为xls文件(老版excel格式),即可匹配原创 2020-08-03 08:56:18 · 8328 阅读 · 0 评论 -
GIS学习(二):插值分析与密度分析的区别
密度分析与插值分析的区别1 定义1 .1插值分析在实际工作中,由于成本的限制、测量工作实施困难大等因素,我们不能对研究区域的每一位置都进行测量(如高程、降雨、化学物质浓度和噪声等级)。这时,我们可以考虑合理选取采样点,然后通过采样点的测量值,使用适当的数学模型,对区域所有位置进行预测,形成测量值表面。插值之所以可称为一种可行的方案,是因为我们假设,空间分布对象都是空间相关的,也就是说,彼此接...原创 2020-04-08 11:28:37 · 5182 阅读 · 1 评论 -
GIS学习(一):密度分析、可达性分析、空间句法
1)密度分析2) 路网可达性分析所需数据:成本栅格、目标点成本栅格一般由路网矢量转换栅格、重分类(nodata为高成本值)权重的设置:人为设置权重、空间句法设置权重3)空间句法人与环境相互影响可以对路网进行空间句法分析,用作可达性分析的权重...原创 2020-04-08 00:10:44 · 8059 阅读 · 0 评论 -
POI高德坐标转WGS84
因为下载的高德POI数据,相比WGS84是有一定的坐标偏移,因此需要进行校正一下。首先,文件夹是这样的:CSV文件格式为:代码如下:#!/usr/bin/python# -*- coding: UTF-8 -*- import pandas as pdimport jsonimport mathimport osimport csv## 转换函数:x_pi = 3...原创 2020-01-07 21:37:01 · 5052 阅读 · 2 评论 -
ArcGIS/ENVI学习(一):近期一些操作总结
1)计算城市比例目的:使用清华10m的全球土地覆盖数据,得到城市像元所占的比例影像Arcgis操作:重分类:将本来属于不透水层的8,重分类为1,其余类别(像元值)重分类为0块统计:对一个33的矩形,计算mean值。即可得到30m30m的城市区块比例影像(注意,此时影像的分辨率仍旧为10m,只是每一个30m*30m的区块都是同样的数值)重采样:重采样至所需的分辨率2)计算某影像的标...原创 2019-08-22 09:10:40 · 1570 阅读 · 0 评论 -
GGE学习(一)统计所有区域的影像均值并输出为shp
任务:在不下载影像的前提下,统计全国每一个县的SAR数据的均值,并且输出这个均值。代码:function exportmean(countyID){ var images = imageCollection2.filterBounds(table2) .filter(ee.Filter.date('2015-12-01', '2015-12-31')) .sel...原创 2019-08-06 16:08:49 · 1049 阅读 · 0 评论