yy_nju的专栏

原创 从USGS批量下载L2级Landsat-4/5/7/8影像

网址：https://espa.cr.usgs.gov/index/第一步：Bulk Ordering① Add Input Products ：输入文件列表.txt文件②Select Product Contents：别点③其他该设置的设置一下④ 提交订单第二步：下载ESPA Downloader网址：https://github.com/USGS-EROS...

原创 《TensorFlow实战Google深度学习架构》——Tensorflow实现迁移学习

参考：https://blog.csdn.net/nnnnnnnnnnnny/article/details/70244232迁移学习就是把一个已经训练好的神经网络隐层的输出当做特征，然后用自己的标记数据训练一个新的全连接层，用于分类。 # %load transfer_flower.py#!/usr/bin/env python3# glob，文件路径查找模块import ...

转载 《TensorFlow实战Google深度学习架构》——LeNet-5勘误

参考：https://zhuanlan.zhihu.com/p/31534286 # 《TensorFlow实战Google深度学习框架》06 图像识别与卷积神经网络# win10 Tensorflow1.0.1 python3.5.3# CUDA v8.0 cudnn-8.0-windows10-x64-v5.1# filename:LeNet5_infernece.py # L...

原创 《TensorFlow实战Google深度学习架构》——MNIST数字识别问题勘误

最近在学习《TensorFlow实战Google深度学习架构》这本书，上面的代码大部分都跑不通，真是一言难尽啊...搜了好多资料，好容易把MNIST数字识别问题这段代码改好了，自己写个笔记纪念一下吧。-----------------------------------------# 导入必要的库import tensorflow as tfimport osfrom tens...

原创 逻辑回归Logistic Regression

原创 梯度下降法

原创 将分类图像转换为ERDAS可识别的专题图像（Continuous to Thematic）

so easy~~~ :D

原创 基于卷积神经网络的图像分类

《Convolutional Neural Networks for Large-Scale Remote-Sensing Image Classification》本文提出了一种全卷积神经网络来产生密集分类（逐像素分类），特点是：位置无关、输出是一系列卷积层的结果。首先将负责分类的全连接层转换为一个卷积层，选定卷积核使它的维度和前一层相同，这样它等价于全连接层，区别是当输入图像的尺寸增加时，输出...

原创 通过BNT Toolbox建立GMM分类模型

参考文档：http://www.media.mit.edu/wearables/mithril/BNT/mixtureBNT.txt代码下载地址：http://www.media.mit.edu/wearables/mithril/BNT/mixtureBNT.zip案例：通过建立一个GMM分类器对人类活动（跑、走）进行分类，实验数据为可穿戴设备获取的31维特征向量。说明数据：

原创 基于Bayes Net ToolBox对层次隐马尔可夫模型进行建模、训练和推理

参考论文“Linear Time Inference in Hierarchical HMMs”，示例文件夹：..\bnt\BNT\examples\dynamic\HHMMSample 1. 建立已知结构的HHMM模型对应的DBN模型如下图所示Qsize = [2 3 2];    % 第1、2、3层的节点个数Qnodes = 1:

转载 MatLab2012b/MatLab2013b 分类器大全(svm,knn,随机森林等)

train_data是训练特征数据, train_label是分类标签。Predict_label是预测的标签。MatLab训练数据, 得到语义标签向量 Scores(概率输出)。1.逻辑回归(多项式MultiNomial logistic Regression)Factor = mnrfit(train_data, train_label);Scores = mnrv

原创 Bayes Net toolbox 的使用实例

Example 1、离散贝叶斯网络模型 Cloudy, Sprinklet, Rain, WetGrass模型，每个节点的条件（先验）概率分布如图所示%% ------------------ 模型设置 ------------------%%% model structureN = 4;dag = zeros(N,N);C = 1; S =

原创 贝叶斯公式及推广

回顾下经典贝叶斯公式：P(A,B)=P(A)*P(B|A)=P(B)*P(A|B)多元情况贝叶斯公式推广： P(A|B,C) = P(A,B,C)/P(B,C)  = P(C|A,B)*P(A,B)/P(B,C)=P(C|A,B)*P(B|A)*P(A)/P(C|B)*P(B)

原创 MOD09A1数据与Landsat 4/5 TM波段的对应关系

MOD09A1数据与Landsat 4/5 TM波段的对应关系Band 1            620-670nm                    TM 3                    620-690nmBand 2            841-876nm                     TM 4                    760-9

原创 利用BP神经网络分类iris数据集

clc;clear;load('irisdata.mat');%读取训练数据%[f1,f2,f3,f4,class] = textread('trainData.txt' , '%f%f%f%f%f',150);f1=[setosa(:,1);versicolour(:,1);virginica(:,1)];f2=[setosa(:,2);versicolour(:,2);

原创 图像分类精度评价

新年第一篇，介绍一个经常混淆的问题——图像分类精度评价，多种精度评价指标的含义1、混淆矩阵2、 用户精度 User's Accuracy (UA)用户（使用者）精度表示在该次分类中，在分类图上，落在该类别上的检验点，被正确分类为该类别的比率。2.1 误判误差 （Commission Error）误判误差指该类别被错误分类的概率。

原创 用ENVI评价分类精度

（1）首先需要将外部程序生成的分类文件转化为ENVI可以识别的分类图，方法如下：Envi打开图像->Tools->Color Mapping->Density Slice->Set default numer of ranges(设置为聚类数目)->Applydefault ranges->output ranges to class image（2）对分类结果进行合并与命名，

原创 ERDAS精度评价时如何避免随机生成的样本点落在背景区域

ERDAS中进行精度评价时生成的随即点如何能都落在影像图的范围内，而不落在影像图周围的黑色区域（在对影像图进行不规则边界裁剪时形成的周围像元值是0的区域），可以在生成随机点时选择随机点在哪些分类中生成，这样就可以避免随即点生成在影像周围的黑色边界中。具体的步骤如下：（1）打开精度评价对话框，然后打开待评价的分类图。（2）在精度评价对话框中选择Edit→Create/Add

转载 MATLAB绘制带图案的柱状图

转自：http://blog.sciencenet.cn/home.php?mod=space&uid=40165&do=blog&id=308439function applyhatch(h,patterns,colorlist)%APPLYHATCH Apply hatched patterns to a figure% APPLYHATCH(H,PATTERNS) creates

转载 Matlab 画图修饰-随机线条和随机颜色

转载自：http://www.zhaoyanpeng.cn/archives/237当需要对同一曲线不同参数下进行模拟时需要不同的颜色来加以区分:1234567x=-2:0.0001:2;fori=0:0.2:8    y=i*x.^2;    plot(x,

转载 MATLAB plot多条曲线

Matlab 提供了几种线型？4种plot 多条曲线 （2）" title="Matlab plot 多条曲线 （2）" style="margin:0px; padding:0px; border:0px; list-style:none">Matlab提供了几种Marker？10几种plot 多条曲线 （2）" title="Matlab plot 多条曲线 （2）" styl

转载 MATLAB多维数组操作

1.一个三维数组由行、列和页三维组成，其中每一页包含一个由行和列构成的二维数组。2.利用标准数组函数创建多维数组A=zeros(4,3,2) 生成一个4行3列2页的三维全0数组，ones，rand和randn等函数有相似的用法。3.利用直接索引方式生成多维数组  A=zeros(2,3)  A(:,:,2)=ones(2,3)  A(:,:,3)=4

转载 Matlab运行时出现提示 too many input arguments问题

装载自：http://blog.csdn.net/begtostudy/article/details/4150185matlab 有时候经常出现“too many input arguments”错误提示。 还有就是程序在别人电脑能运行，在另一台电脑就不行。 经本人研究发现，这些问题大多是因为使用者或者外挂的函数库有与Matlab同名的函数导致。 

原创 用R语言读取Modis时间序列图像

由于直接用R语言处理Modis影像比较麻烦，我的方法是用R语言的MODIS和Raster Packages读取Modis序列影像，并存储为.txt文件，然后用熟悉的MATLAB处理，这样比较方便。1、打开安装了MODIS和Raster Packages的RStudio，启动这两个Packages。2、输入如下代码，读取Modis NDVI图像（我使用的数据集是16天合成的M

原创 安装和配置R语言环境

最近在学习使用R语言处理MODIS时间序列数据处理，有一些

原创 《Sparse and Redundant Representations》第六章 迭代收缩算法

原创 《Sparse and Redundant Representations》第五章 从精确解到近似解

原创 《Sparse and Redundant Representations》第四章 追赶法的保证

原创 《Sparse and Redundant Representations》第三章：追赶法-实验

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

原创 《Sparse and Redundant Representations》第二章：唯一性与不确定性

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

原创 《Sparse and Redundant Representations：From Theory to Applications in Signal and Image Processing》序论

本人是纯小菜鸟一只，最近在精读Elad的《Sparse and Redundant Representations From Theory to Applications in Signal and Image Proces》，写下每一章的读书报告，激励自己不断学习！+U第一章，作者主要探讨了线性系统方程的解的问题。1.1 理解线性系统矩阵，且A为满秩矩阵。有如下方程：