![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
GIS
文章平均质量分 79
撼沧
From zero to hero.
展开
-
使用ArcGIS提取HWSD中的土壤属性数据
基于世界土壤数据库的土壤数据集Siol Map Based Harmonized World Soil Database (HWSD)写在前面:网上关于该数据库的基本介绍以及如何加载到ArcGIS里方法教程较多,此处附上几个链接,读者自行参考。官方链接寒区旱区科学数据中心HWSD数据库的使用孟现勇-如何使用世界土壤数据库(HWSD)建立SWAT模型土壤数据库本文意在将HWSD数据库提取方法进行统一描述,方便自己复用和他人借鉴。HWSD拥有较多土壤属性数据,使用时往往需要将某条土壤属性原创 2022-01-17 16:18:59 · 9976 阅读 · 11 评论 -
基于R语言进行栅格数据统计及Raster包简介
基于R语言进行栅格数据统计及Raster包简介写在前面:日常数据处理中常需要对单个栅格图层进行统计,包括计算其均值、最值和总和等等。使用R语言可方面快速地实现该需求,主要借助的package为Raster包。library(raster)x = raster('data/data.tif')cellStats(x,stat='mean')x表示待计算的栅格;stat表示需要统计的内容,包括sum, mean, min, max, sd(标准差), skew(偏度) 和 rms(均方根);原创 2021-05-25 15:41:58 · 11284 阅读 · 0 评论 -
ArcGIS许可管理器License server administrator无法启动的解决方法
ArcGIS许可管理器License server administrator无法启动的解决方法在安装ArcGIS软件,甚至日常使用启动软件时,ArcGIS许可管理器License Manager经常会无法启动,点多少次都没反应,更别谈往后的破解了。这里谈一下自己遇到问题后用过的能用的方法。原创 2021-01-07 22:50:16 · 117152 阅读 · 25 评论 -
基于R语言读取NC数据并提取单个波段(netcdf转为tiff或其他格式)
写在前面netcdf是气象、生态、地理领域非常常用的一种数据格式,读写的方法很多,arcgis、arcpy、python、matlab等等均有相应的方法。本文基于R语言实现netcdf的数据读取,提取数据中的单个波段(某一时间的数据),展示和写出。注:本文仅展示基础的代码,多个波段或多个数据的读写操作加上相应的循环即可。具体代码加载包library(raster)library(rasterVis)library(ncdf4)library(lattice)打开nc数据,并展示其中原创 2020-12-26 00:33:43 · 10456 阅读 · 2 评论 -
基于R语言或Matlab直接读取压缩包中的数据
写在前面以gosif网站上的tiff数据为例。该数据是地学上常用的一种数据,希望读取该网站上的gosif数据:http://data.globalecology.unh.edu/data/GOSIF_v2/8day/但数据全是压缩包格式的,解压后内存占用又会变大。所以需要直接从压缩包中读取数据。基于R语言library(R.utils)library(stringr)library(raster)library(rasterVis)# download and read the dat原创 2020-12-23 12:47:01 · 3170 阅读 · 0 评论 -
【python】使用python绘制地图时手动添加比例尺
写在前面使用basemap、cartopy等包在python中绘制地图添加比例尺有投影限制,不太方便。本博文使用最基础的matplotlib画线实现比例尺的添加,为大家提供一种可选思路。思路:提供了三种可选风格,根据自己需要进行选择。原理为根据1°=111千米的原理计算需要绘制的长度,因此本代码主要适用于低纬度地区的wgs84投影的绘图。高纬度地区投影变形太大,使用这种方法会不准确。大家也可以根据这种思路自己修改代码。详细代码:#-----------函数:添加比例尺------------原创 2020-12-08 23:15:14 · 5781 阅读 · 0 评论 -
【R语言】使用r语言进行单个或多个栅格数据间的计算
本文介绍如何使用r语言进行单个或多个栅格数据间的计算,类似于ArcGIS中的栅格计算器。单个栅格数据:library(raster)calc(x, fun, filename='', na.rm, forcefun=FALSE, forceapply=FALSE, ...)元素:x:Raster* objectfun:functionfilename:character. Output filename (optional)Examples# NOT RUN {r <- rast原创 2020-12-07 10:08:31 · 7615 阅读 · 0 评论 -
自然资源部标准地图底图转矢量Shapefile并配准
当前许多期刊投稿要求凡涉国界图件(国内部分地区、全国、世界部分地区、全球)必须使用自然资源部标准地图底图(下载网址:http://bzdt.ch.mnr.gov.cn),所用底图边界要完全无修改。但该网站并未提供可编辑的.shp矢量文件,除p图外,无法直接应用于自己的数据,在arcgis或其他软件进行作图。因此,本文介绍了如何将标准地图底图转为Shapefile矢量图层,及如何将自己的数据与标准底图进行配准,使shp边界对齐,达到在arcgis里出图的目的。原创 2020-12-04 10:18:02 · 31691 阅读 · 24 评论 -
【R语言】根据经纬度从tif图像中提取兴趣点
library(stringr)library(raster)library(sp)#set work spacegetwd()setwd("G:/data")#初始设置,提取点lai <- matrix(nr=50,nc=2)point<-cbind(128.5,45.9) #lon,lati <- 1n <- 1repeat { #读取栅...原创 2020-04-03 10:45:28 · 3761 阅读 · 2 评论