天外来客的花田

基于流域边界批量裁剪栅格图像

用MATLAB的websave下数据时遇到了上述报错，在外网找到了答案：加一句options.CertificateFilename=('');即可！后记：写博客的初衷是分享经验，同时是算是自己对思路和代码的整理，方便日后处理数据，应该可以帮到很多人。我已免费分享我的心得，如果看官还有其他问题的，那么：知识付费，我的时间和经验正好可以解决你的问题。咨询问题请添加QQ：8193693542022年4月20日...

# -*- coding: utf-8 -*-'''批量实现shp转raster2021/12/24'''import arcpyimport osimport re#检查Arcgis的权限arcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput = 1#设置了输出范围及输出坐标系，这三行也可注释掉arcpy.env.overwriteOutput = Truearcpy.env.extent = arcpy.Ext

%% 修改日期 2021/12/9clcclearclose all%% DEM空间分布图figure(1)worldmap([35 75],[-150 -50]); %研究区范围[DEM,R] = geotiffread('F:\DEM_SRTM\North America\Canada_USA_mianland.tif'); %加载带有坐标信息的图DEM = double(DEM); %需要将

任务描述：已有散点图，想画密度图，效果如下网上很多方法都太麻烦了，废话不多说，直接贴代码：clcclearclose allfigure(1)X = [DI,Er]; %DI Er 是两个列向量N = hist3(X,'CdataMode','auto','Ctrs',{0:0.005:3 0:0.005:1.4}); % N的值正是密度% N是一个二维矩阵，Size大小由Ctrs决定，此处是280*600，因为0.005*600=3 , 0.005*280=1.4N

因为课题处理30年的降雨和蒸发的遥感资料（.NC格式），而想要在Arcgis中处理要求的是raster格式的，所以需要批量转化为tif文件，所以在此分享自己改编之后的代码，可以简洁明了的实现这个过程：版本：MATLAB_2018bclc;clear;%% 批读取NC文件的准备工作datadir = 'G:\Global_P_ET\MSWEP_V2.2\'; %指定批量数据所在的文件夹...

在做地理信息相关的处理时，我们希望的是图片格式（或raster格式）便于在GIS中处理以及可视化，但拿到的数据往往又是TXT文本格式，需要转化，代码如下： data = load('F:\budyko_ds\Zr\Mean_Zr.txt','%s'); data = flipud(data); lat = 89.75:-0.5:-89.75 lon = -179.75:0.5:179.75...

clc; %清除命令行窗口clear; %清除工作区filename='F:\Global_P_ET\ET\GLEAM\1982\E_1982_GLEAM_v3.2a.nc';nc_dump (filename); %展示文件信息%% 提取变量lon = ncread(filename,'lon')';

帮网友处理了个土壤水的hdf（h5）数据，记录一下心得：注：以下只是对单个hdf文件的处理，关于批处理我在另一个博客中写得很清楚，思路大同小异，请看：MATLAB 批量读取NC文件并转为TIF文件代码：clcclear%第一步 展示h5文件的结构（名字 几个group 几个dataset 以及其他属性）h5info('C:\Users\DELL\Desktop\数据处理\SMAP_L3_SM_P_20150331_R16510_001.h5')%第二步 展示所有group的信息，查到我们

任务描述：需要下载NASA发布的GPM（全球降雨观测计划）数据，已得到Data File Links如下：从第3行开始就是所需要的数据的URL链接NASA下载数据是要先登录的，所以我们也需要进行身份识别，否则就会报错（error401，unauthorized），详见：代码如下：（加了一些比较有意思的东西。如waitbar，调教一下MATLAB哈哈哈）%% 读取NASA的Data File Links,批量下载GPM数据url_total = textread('C:\Users\DELL\Do

数字高程模型（Digital Elevation Model)，简称DEM网上有很多中国DEM的下载链接，要么收费，要么地图不对（缺藏南、台湾等等），要么版本太老所以自己做了一个，流程如下：第一步：NASA官网下载STRM 90m文件包（http://srtm.csi.cgiar.org）(i) 第二张图的Tile size就是在选区域下载，所以也可以选5x5degree，但就要选很多区域，而且下载的时候要一个个点，太麻烦了，30x30degree的话，只需要选六块就能包住中国；(ii) 下载速度

有时网速不太行或者下载的文件比较大时就会可能会报错，只需要修改options.Timeout即可options.Timeout = 20;如上，我将时间改为20s，问题解决

任务描述：如题，以2010年月尺度1km的MODIS的植被覆盖度（NDVI）数据为例第一步 ：获得MODIS数据下载链接Earthdata Search(下载地址)可以选择自己想要的时间空间范围，NASA官网会自动生成下载链接，点击保存得到txt文件第二步：下载数据（使用MATLAB）那么如何根据Txt文件中的url链接下载数据呢？%% 代码用于下载MODIS NDVI月尺度数据，全部存在默认文件夹中，需要将系统默认浏览器设置为Chrome，存储路径在浏览器中修改clcclearclos

第一步：将网页保存（html格式），然后用记事本打开（或者直接将文件后缀名改为txt）如图，红框内就是我们想要扒取的文件名手动删掉无用的前十行第二步：提取链接并用url函数开启下载%% 下载Globsnow数据 修改时间2021/5/12clcclearclose all%% 下载链接，将网页保存为txt后需要将前十行及最后一行删去（无用行）fid = fopen('Index of _swe_archive_v3.0_L3A_daily_SWE_NetCDF4.txt');tli.